[BZOJ 3233] [Ahoi2013]找硬币

BZOJ传送门

题目描述

小蛇是金融部部长。最近她决定制造一系列新的货币。假设她要制造的货币的面值为$x_1,x_2,x_3\dots$ 那么$x_1$必须为$1$，$x_b$必须为$x_a$的正整数倍（$b>a$）。例如 $1,5,125,250$就是一组合法的硬币序列，而$1,5,100,125$就不是。不知从哪一天开始，可爱的蛇爱上了一种萌物——兔纸！从此，小蛇便走上了遇上兔纸娃娃就买的不归路。某天，小蛇看到了$N$只可爱的兔纸，假设这$N$ 只兔纸的价钱分别是$a_1,a_2\dots a_N$。现在小蛇想知道，在哪一组合法的硬币序列下，买这$N$只兔纸所需要的硬币数最少。买兔纸时不能找零。

输入输出格式

输入格式

第一行，一个整数$N$，表示兔纸的个数

第二行，$N$个用空格隔开的整数，分别为$N$只兔纸的价钱

输出格式

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

输入输出样例

输入样例#1：

2 
25 102 

输出样例#1：

4

提示

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

$1\le N\le 50,1\le ai\le 100,000$

解题分析

一开始没啥思路， 然后发现这个倍数有个很好的性质： 无论更大面值的硬币怎么选， 选到面值为$x$的硬币的时候价格为$a_i$的兔子一定会剩下$a_i\mathrm{m}\mathrm{o}\mathrm{d}x$的零头无法付清， 所以我们直接从大到小$dp$， 枚举质因数转移即可。

总复杂度？ 貌似是$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]);
}