3233: [Ahoi2013]找硬币

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

题解：

f[i]表示以i为最大值的最少硬币数，

f[i]=min(f[j]-sigma(a[k]/i)*(i/j-1))

a[k]/i就表示这个兔子需要几个i值的硬币，这些硬币值(i/j)个值为j的硬币，所以可以少(i/j-1)

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
const int N=100001;
int n;
int a[60];
int prime[N],top=0;
bool v[N];
void ge_pr()
{
	for(int i=2;i<N;i++)
	{
		if(!v[i]) prime[++top]=i;
		for(int j=1;j<=top&&prime[j]*i<N;j++)
		{
			v[i*prime[j]]=1;
			if(i%prime[j]==0)break;
		}
	}
}
int tot=0,ans,mx=0;
int f[N];
int main()
{
	ge_pr();
	scanf("%d",&n);
	for(int i=1;i<=n;i++)
	scanf("%d",&a[i]),tot+=a[i],mx=max(mx,a[i]);
	memset(f,63,sizeof(f));
	ans=f[1]=tot;
	for(int i=2;i<=mx;i++)
	{
		f[i]=tot;
		int p=i;
		for(int j=1;j<=top;j++)
		{
			if(p%prime[j]==0)
			{
				int x=i/prime[j];
				int sum=0;
				for(int u=1;u<=n;u++) sum+=(a[u]/i)*(prime[j]-1);
				f[i]=min(f[i],f[x]-sum);
				while(p%prime[j]==0) p/=prime[j];
				if(p==1) break;
			}
		}
		ans=min(ans,f[i]);
	}
	printf("%d\n",ans);
}