BZOJ传送门
题目描述
小蛇是金融部部长。最近她决定制造一系列新的货币。假设她要制造的货币的面值为 x 1 , x 2 , x 3 … x_1,x_2,x_3… x1,x2,x3… 那么 x 1 x_1 x1必须为 1 1 1, x b x_b xb必须为 x a x_a xa的正整数倍( b > a b>a b>a)。例如 1 , 5 , 125 , 250 1,5,125,250 1,5,125,250就是一组合法的硬币序列,而 1 , 5 , 100 , 125 1,5,100,125 1,5,100,125就不是。不知从哪一天开始,可爱的蛇爱上了一种萌物——兔纸!从此,小蛇便走上了遇上兔纸娃娃就买的不归路。某天,小蛇看到了 N N N只可爱的兔纸,假设这 N N N 只兔纸的价钱分别是 a 1 , a 2 … a N a_1,a_2…a_N a1,a2…aN。现在小蛇想知道,在哪一组合法的硬币序列下,买这 N N N只兔纸所需要的硬币数最少。买兔纸时不能找零。
输入输出格式
输入格式
第一行,一个整数 N N N,表示兔纸的个数
第二行, N N N个用空格隔开的整数,分别为 N N N只兔纸的价钱
输出格式
一行,一个整数,表示最少付的钱币数。
输入输出样例
输入样例#1:
2
25 102
输出样例#1:
4
提示
样例解释:共有两只兔纸,价钱分别为 25 25 25和 102 102 102。现在小蛇构造 1 1 1, 25 25 25, 100 100 100这样一组硬币序列,那么付第一只兔纸只需要一个面值为 25 25 25的硬币,第二只兔纸需要一个面值为 100 100 100的硬币和两个面值为 1 1 1的硬币,总共两只兔纸需要付 4 4 4个硬币。这也是所有方案中最少所需要付的硬币数。
1 ≤ N ≤ 50 , 1 ≤ a i ≤ 100 , 000 1\le N\le 50, 1\le ai\le 100,000 1≤N≤50,1≤ai≤100,000
解题分析
一开始没啥思路, 然后发现这个倍数有个很好的性质: 无论更大面值的硬币怎么选, 选到面值为 x x x的硬币的时候价格为 a i a_i ai的兔子一定会剩下 a i m o d    x a_i\mod x aimodx的零头无法付清, 所以我们直接从大到小 d p dp dp, 枚举质因数转移即可。
总复杂度? 貌似是 O ( n l o g l o g n ) O(nloglogn) O(nloglogn)的?
代码如下:
#include <cstdio>
#include <cstring>
#include <cmath>
#include <cstdlib>
#include <cctype>
#include <algorithm>
#define R register
#define IN inline
#define W while
#define gc getchar()
#define MX 100500
template <class T>
IN void in(T &x)
{
x = 0; R char c = gc;
for (; !isdigit(c); c = gc);
for (; isdigit(c); c = gc)
x = (x << 1) + (x << 3) + c - 48;
}
template <class T> IN T max(T a, T b) {return a > b ? a : b;}
template <class T> IN T min(T a, T b) {return a < b ? a : b;}
int n, mx, pcnt;
int pri[MX], low[MX], dp[MX], dat[55];
bool npr[MX];
void Sieve()
{
int tar;
for (R int i = 2; i <= mx; ++i)
{
if (!npr[i]) pri[++pcnt] = i, low[i] = i;
for (R int j = 1; j <= pcnt; ++j)
{
tar = i * pri[j];
if (tar > mx) break;
npr[tar] = true; low[tar] = pri[j];
if (!(i % pri[j])) break;
}
}
}
int main(void)
{
in(n); int foo, bar, mn, tar, tot;
for (R int i = 1; i <= n; ++i) in(dat[i]), mx = max(mx, dat[i]);
Sieve(); std::memset(dp, 63, sizeof(dp));
for (R int i = mx; i >= 1; --i)
{
mn = 0;
for (R int j = 1; j <= n; ++j) mn += dat[j] / i;
dp[i] = min(dp[i], mn);
foo = i;
W (foo ^ 1)
{
bar = low[foo];
tar = i / bar; tot = 0;
for (R int j = 1; j <= n; ++j) tot += dat[j] % i / tar;
dp[tar] = min(dp[tar], dp[i] + tot);
W (!(foo % bar)) foo /= bar;
}
}
printf("%d", dp[1]);
}