3233: [Ahoi2013]找硬币 线性筛+DP

人太弱了什么题也不会。。。

---

用$f_i$表示当前最大面值的硬币为$i$，零头所需的硬币数量的最小值。
那么有$f_j=min\{f_i+\sum_{k=1}^n \lfloor \frac {a_k \ mod \ j}{i} \rfloor\}(i\mid j)$
则有$ans=min\{f_i+\sum_{k=1}^n\lfloor\frac{a_k}{i}\rfloor\}$
优化：要保证$\frac{j}{i}$是质数，否则存在更优解，所以筛一下质数枚举就行了。

#include<iostream>
#include<cstdio>
#include<cstring>
#define N 100005
#define inf 1000000007
using namespace std;
int n,a[55];
int f[N],prime[N],ans=inf;
bool flag[N];
inline int read()
{
    int a=0,f=1; char c=getchar();
    while (c<'0'||c>'9') {if (c=='-') f=-1; c=getchar();}
    while (c>='0'&&c<='9') {a=a*10+c-'0'; c=getchar();}
    return a*f;
}
inline void get_prime()
{
    for (int i=2;i<=100000;i++)
    {
        if (!flag[i]) prime[++prime[0]]=i;
        for (int j=1;j<=prime[0]&&i*prime[j]<=100000;j++)
        {
            flag[i*prime[j]]=1;
            if (i%prime[j]==0) break;
        }
    }
}
int main()
{
    get_prime();
    n=read();
    for (int i=1;i<=n;i++) a[i]=read();
    memset(f,0x3f,sizeof(f)); f[1]=0;
    for (int i=1;i<=100000;i++)
    {
        for (int j=1;j<=prime[0]&&prime[j]*i<=100000;j++)
        {
            int val=prime[j]*i;
            int tmp=0;
            for (int k=1;k<=n;k++)
                tmp+=a[k]%val/i;
            f[val]=min(f[val],f[i]+tmp);
        }
        int tmp=f[i];
        for (int k=1;k<=n;k++)
            tmp+=a[k]/i;
        ans=min(ans,tmp);
    }
    cout << ans << endl;
    return 0;
}