【题解】
本蒻一直在想二维dp,看了题解才发现竟然一维就可以
设f[i]为最大面值为i时,买下所有兔纸花费的最小硬币数f[i] = min{ f[j] - sigma(a[k]/i*(i/j-1)) } , j|i,其中,j为次大面值,这个方程考虑的是选了i能减小多少j的使用
注意, 如果硬币种类很多,是不影响最优答案的(不用就行了) -----------> 重要的性质
所以,可以在枚举j的这一步下手优化:
规定i/j必须为质数,否则可以在它们之间加一个面值k,使 i/k,k/j都是质数,不影响答案
可以在线性筛的过程中求出数组min[i]:i的最小质因子,用它实现:仅枚举i的 质因子作为i/j的值
n的质因子个数是log(n)级别的,复杂度:O( n*Max*log(Max) )
【代码】
#include<stdio.h>
#include<stdlib.h>
int a[100005],f[100005],pri[100005],min[100005];
int main()
{
int n,i,j,k,Ma=0,cnt=0,t,ans;
scanf("%d",&n);
for(i=1;i<=n;i++)
{
scanf("%d",&a[i]);
f[1]+=a[i];
if(Ma<a[i]) Ma=a[i];
}
for(i=2;i<=Ma;i++)
{
if(min[i]==0)
{
pri[++cnt]=i;
min[i]=i;
}
for(j=1;j<=cnt&&pri[j]*i<=Ma;j++)
{
min[pri[j]*i]=pri[j];
if(i%pri[j]==0) break;
}
}
for(i=2;i<=Ma;i++)
f[i]=f[1];
ans=f[1];
for(i=2;i<=Ma;i++)
{
j=i;
while(j>1)
{
t=f[i/min[j]];
for(k=1;k<=n;k++)
t-=a[k]/i*(min[j]-1);
if(f[i]>t) f[i]=t;
while(min[j]==min[j/min[j]]) j/=min[j];
j/=min[j];
}
if(ans>f[i]) ans=f[i];
}
printf("%d",ans);
return 0;
}