BZOJ 3233: [Ahoi2013]找硬币( dp )

时间:2022-08-31 19:30:49

BZOJ 3233: [Ahoi2013]找硬币( dp )

dp(x)表示最大面值为x时需要的最少硬币数.

枚举x的质因数p,  dp(x) = min( dp(x/p) - (p-1) * sigma[a[i]/x] ).

----------------------------------------------------------------------------------

#include<cstdio>
#include<cstring>
#include<algorithm>
 
using namespace std;
 
const int maxn = 59;
const int maxm = 100009;
 
int N, M, w[maxn], dp[maxm];
int p[maxm], minp[maxm], pn;
bool F[maxm];
 
template<class T>
inline void Min(T &x, T t) {
if(t < x) x = t;
}
template<class T>
inline void Max(T &x, T t) {
if(t > x) x = t;
}
 
void Init() {
scanf("%d", &N);
memset(dp, -1, sizeof dp);
M = dp[1] = 0;
for(int i = 0; i < N; i++) {
scanf("%d", w + i);
Max(M, w[i]);
dp[1] += w[i];
}
memset(F, 0, sizeof F);
pn = 0;
for(int i = 2; i <= M; i++) {
if(!F[i])
minp[i] = p[pn++] = i;
for(int j = 0; j < pn && i * p[j] <= M; j++) {
F[i * p[j]] = true;
minp[i * p[j]] = p[j];
if(i % p[j] == 0) break;
}
}
}
 
void Work() {
int ans = dp[1];
for(int i = 0; i < pn; i++) {
dp[p[i]] = 0;
for(int j = 0; j < N; j++)
dp[p[i]] += w[j] / p[i] + w[j] % p[i];
Min(ans, dp[p[i]]);
}
for(int i = 2; i <= M; i++) if(!~dp[i]) {
dp[i] = dp[1];
for(int t = i; t != 1; t /= minp[t]) {
int v = 0;
for(int j = 0; j < N; j++)
v += w[j] / i;
Min(dp[i], dp[i / minp[t]] - v * (minp[t] - 1));
}
Min(ans, dp[i]);
}
printf("%d\n", ans);
}
 
int main() {
Init();
Work();
return 0;
}

----------------------------------------------------------------------------------

3233: [Ahoi2013]找硬币

Time Limit: 10 Sec  Memory Limit: 64 MB
Submit: 617  Solved: 275
[Submit][Status][Discuss]

Description

小蛇是金融部部长。最近她决定制造一系列新的货币。假设她要制造的货币的面值为x1,x2,x3… 那么x1必须为1,xb必须为xa的正整数倍(b>a)。例如 1,5,125,250就是一组合法的硬币序列,而1,5,100,125就不是。不知从哪一天开始,可爱的蛇爱上了一种萌物——兔纸!从此,小蛇便走上了遇上兔纸娃娃就买的不归路。某天,小蛇看到了N只可爱的兔纸,假设这N 只兔纸的价钱分别是a1,a2…aN。现在小蛇想知道,在哪一组合法的硬币序列下,买这N只兔纸所需要的硬币数最少。买兔纸时不能找零。
 
 
 

Input

第一行,一个整数N,表示兔纸的个数
第二行,N个用空格隔开的整数,分别为N只兔纸的价钱
 
 
 

Output

一行,一个整数,表示最少付的钱币数。

 
 
 

Sample Input

2
25 102

Sample Output

4

HINT

样例解释:共有两只兔纸,价钱分别为25和102。现在小蛇构造1,25,100这样一组硬币序列,那么付第一只兔纸只需要一个面值为25的硬币,第二只兔纸需要一个面值为100的硬币和两个面值为1的硬币,总共两只兔纸需要付4个硬币。这也是所有方案中最少所需要付的硬币数。

1<=N<=50, 1<=ai<=100,000

Source

BZOJ 3233: [Ahoi2013]找硬币( dp )的更多相关文章

  1. BZOJ 3233&colon; &lbrack;Ahoi2013&rsqb;找硬币

    BZOJ 3233: [Ahoi2013]找硬币 标签(空格分隔): OI-BZOJ OI-DP Time Limit: 10 Sec Memory Limit: 64 MB Description ...

  2. BZOJ3233&colon;&lbrack;AHOI2013&rsqb;找硬币&lpar;DP&rpar;

    Description 小蛇是金融部部长.最近她决定制造一系列新的货币.假设她要制造的货币的面值为x1,x2,x3… 那么x1必须为1,xb必须为xa的正整数倍(b>a).例如 1,5,125, ...

  3. &lbrack;Bzoj3233&rsqb;&lbrack;Ahoi2013&rsqb;找硬币&lbrack;基础DP&rsqb;

    3233: [Ahoi2013]找硬币 Time Limit: 10 Sec  Memory Limit: 64 MBSubmit: 924  Solved: 482[Submit][Status][ ...

  4. &lbrack;AHOI2013&rsqb;找硬币(搜索)

    [Ahoi2013]找硬币 Time Limit: 10 Sec  Memory Limit: 64 MBSubmit: 348  Solved: 114[Submit][Status] Descri ...

  5. 【bzoj 3233】&lbrack;Ahoi2013&rsqb;找硬币 ——搜索

    Description 小蛇是金融部部长.最近她决定制造一系列新的货币.假设她要制造的货币的面值为x1,x2,x3… 那么x1必须为1,xb必须为xa的正整数倍(b>a).例如 1,5,125, ...

  6. 【BZOJ 3233】 &lbrack;Ahoi2013&rsqb;找硬币

    [题目 描述] 小蛇是金融部部长. 最近她决定制造一系列新的货币. 假设她要制造的货币 的面值为 x1, x2, x3… 那么 x1 必须为 1, xb 必须为 xa 的正整数倍(b>a). 例 ...

  7. &lbrack;bzoj3233&rsqb; &lbrack;Ahoi2013&rsqb;找硬币

    一开始没什么思路...后来想到确定最大硬币面值就知道其他面值能取多少了..而且结果是可以由较小的面值转移过来的. f[i]表示最大面值为i时的最小硬币数.a[i]表示第i个物品的价钱. f[i]=mi ...

  8. &lbrack;BZOJ 3233&rsqb; 找硬币

    Link: BZOJ 3233 传送门 Solution: 在本蒟蒻看来算是一道比较神的$dp$了 一开始转移方程都没看出来…… 首先,如果确定了最大面值,是能推出其他面值的所有可能值的 从而发现最大 ...

  9. BZOJ 3235&colon; &lbrack;Ahoi2013&rsqb;好方的蛇

    BZOJ 3235: [Ahoi2013]好方的蛇 标签(空格分隔): OI-BZOJ OI-DP OI-容斥原理 Time Limit: 10 Sec Memory Limit: 64 MB Des ...

随机推荐

  1. 一行R代码来实现繁琐的可视化

    ggfortify 有着简单易用的统一的界面来用一行代码来对许多受欢迎的R软件包结果进行二维可视化的一个R工具包.这让许多的统计学家以及数据科学家省去了许多繁琐和重复的过程,不用对结果进行任何处理就能 ...

  2. 关于python的最大递归层数详解

    在阅读http://www.cnblogs.com/skabyy/p/3451780.html这篇文章的时候,实验yield的流式迭代素数的时候发现有个问题,故详细记录下来. 首先来看看python默 ...

  3. 《淘宝技术这十年》之LAMP架构的网站

    本文节选自<淘宝技术这十年>一书,子柳(赵超)著,由电子工业出版社出版.作者的系列博文:从P1到P7--我在淘宝这7年 2003年4月7日,马云在杭州成立了一个神秘的组织.他叫来十位员工, ...

  4. IOS第17天&lpar;2&comma;Quartz2D图片剪裁变圆行图,和截屏图片)

    **** #import "HMViewController.h" #import "UIImage+Tool.h" @interface HMViewCont ...

  5. java调用&period;net asmx服务

    有时候,在java下开发会调用一下.net下写的服务,看网上有各种方法,但总是不成功,总结下自己测试过能调用成功的方式: 1. 原始方式http-soap public static String p ...

  6. 【转】【cocos2d-x教程】如何使用WebSocket

    介绍 WebSocket是HTML5开始提供的一种浏览器与服务器间进行全双工通讯的网络技术.在WebSocket API中,浏览器和服务器只需要做一个握手的动作,然后,浏览器和服务器之间就形成了一条快 ...

  7. C&plus;&plus;const使用&lpar;06&rpar;

    可以在类中使用const关键字定义数据成员和成员函数或修饰一个对象.一个const对象只能访问const成员函数,否则将产生编译错误. 常量成员 常量成员包括常量数据成员.静态常数据成员和常引用.静态 ...

  8. CSS制作波浪线

    建议先去了解清楚了径向渐变,线性渐变的用法先 这个作者的css制作波浪线讲解很不错额:https://www.jianshu.com/p/8570433e3669不理解的可以看看这个链接的额 可以去菜 ...

  9. jS弹出新窗口被拦截的解决方法

    使用ajax处理数据,在回调中跳转到或打开新页面,这时就会被浏览器拦截 解决方案:先用window.open打开一个窗口,然后修改该窗口地址 var w = window.open('about:bl ...

  10. AppleID的双重认证

    [链接]AppleID的双重认证https://support.apple.com/zh-cn/HT204915