hdu6237 素因子

传送门

题意：给定n堆石子，每次可以将一个石子从一堆移到另一堆，问最少移动几次，使得最后每一堆石子的个数都为一个大于等于2的数x的倍数

思路：首先x一定是n堆石子数和s的因子，所以可以枚举s的因子来对每堆石子进行调整，再找最小的调整数。

但是其实不需要找到s的所有因子，只需要找s的素因子即可

我们可以枚举每个素因子p[i]，记录下每堆石子数距离p[i]的倍数的求模后的值（也就是余数），并对这些余数求和，然后对余数再排序，我们从余数最大值开始补（这样需要补的数一定最少），每补一次，余数和都减去这个素因子，直到余数和减到0

#include <cstdio>
#include <iostream>
#include <algorithm>
#define INF 0x3f3f3f3f3f3f3f3f
using namespace std;
typedef long long ll;
const int maxn=1e6+5;
int a[100005],b[100005];
long long prim[maxn/10];
bool vis[maxn+10];
long long p[maxn];
int cnt;
int ct;
void init()
{
    cnt=0;
    for(int i=2;i*i<=maxn;i++)
        if(!vis[i])
          for(int j=i*i;j<=maxn;j+=i)  vis[j]=1;
    for(int i=2;i<=maxn;i++)
        if(!vis[i])  prim[cnt++]=i;
}
void fun(ll a)
{
    ct=0;
    for(int i=0;i<cnt&&prim[i]*prim[i]<=a;i++)
    {
        if(a%prim[i]==0)
        {
            p[ct++]=prim[i];
            while(a%prim[i]==0)
                a/=prim[i];
        }
    }
    if(a>1) p[ct++]=a;
}
int main()
{
    int w;
    scanf("%d",&w);
    init();
    while (w--)
    {
        int n,big=-1;
        long long s=0,ans=INF;
        scanf("%d",&n);
        for (int i=0;i<n;i++)
        {
            scanf("%d",&a[i]);
            s+=a[i];
            big=max(big,a[i]);
        }
        fun(s);
        if (n==1)
        {
            printf("0\n");
            continue;
        }
        for (int i=0;i<ct;i++)
        {
            long long t=0;
            long long ys=0;
            for (int j=0;j<n;j++)
            {
                b[j]=a[j]%p[i];
                ys+=b[j];
            }
            sort(b,b+n);
            for (int j=n-1;j>=0;j--)
            {
                t+=p[i]-b[j];
                ys-=p[i];
                if (ys==0)
                {
                    ans=min(ans,t);
                    break;
                }
            }
        }
        printf("%lld\n",ans);
    }
    return 0;
}