HDU - A Simple Stone Game (数论，质因子)

题目传送
题意:
给你n堆石头，现在你每一次能将一堆中的一个石头转移到其他任意一堆石头中。如果存在一个x > 1 使得现在的每一堆的石头总数都是x的倍数那么是完美的，问现在最少移动多少次，使得这n堆石头完美?

思路:
1.很显然每一堆石头都是x的倍数，那么也就是说，这n堆石头有共同的一个因子(约数)。

2.那么如果n堆石头有一个共同大于1的因子x，那么很显然，这些石头加起来的总数(一共的石头数)的因子也有x。

3.那么我们现在就分解总数的因子，又因为唯一分解定理，所以我们只用分解出总数的质因子，然后来枚举每一堆石头所改变的次数，再更新最小值。

具体代码操作

#include <bits/stdc++.h>
inline long long read(){char c = getchar();long long x = 0,s = 1;
while(c < '0' || c > '9') {if(c == '-') s = -1;c = getchar();}
while(c >= '0' && c <= '9') {x = x*10 + c -'0';c = getchar();}
return x*s;}
using namespace std;
#define NewNode (TreeNode *)malloc(sizeof(TreeNode))
#define Mem(a,b) memset(a,b,sizeof(a))
#define lowbit(x) (x)&(-x)
const int N = 1e5 + 10;
const long long INFINF = 0x7f7f7f7f7f7f7f;
const int INF = 0x3f3f3f3f;
const double EPS = 1e-7;
const unsigned long long mod = 1e9+7;
const double II = acos(-1);
const double PP = (II*1.0)/(180.00);
typedef long long ll;
typedef unsigned long long ull;
typedef pair<int,int> pii;
typedef pair<ll,ll> piil;
int main()
{
    std::ios::sync_with_stdio(false);
    cin.tie(0),cout.tie(0);
    ll t;
    t = read();//由于数据较大，所以我选择了快读
    while(t--)
    {
        vector<ll> k,ans;
        ll n,sum = 0,Min = 1e18;
        n = read();
        ll arr[n+5] = {0};
        for(ll i = 0;i < n;i++) {arr[i] = read();sum += arr[i];}//sum为总的石头数
        for(ll i = 2;i * i <= sum;i++)//找出所有的质因子，这个是标准的模板代码
        {
            if(sum % i == 0)
            {
                k.push_back(i);//i为质因子
                while(sum % i == 0)
                    sum /= i;
            }
        }
        if(sum > 1) k.push_back(sum);//如果剩下的不是1，那么sum也是一个质因子
        for(ll i = 0;i < k.size();i++)
        {
            ans.clear();
            ll cnt = 0,sumsum = 0;
            for(ll j = 0;j < n;j++)
            {
                ll a = arr[j] % k[i];//a为差多少为k[i]这个因子的倍数
                if(a != 0) ans.push_back(a),sumsum += a;//sumsum为总共差多少，ans记录每一个差多少
            }
            sort(ans.begin(),ans.end());//排序，
            for(ll j = ans.size()-1;j >= 0;j--)//我们一定要从大到小枚举，因为大的，一定是获得其他堆的石子才最小
            {
                cnt += (k[i]-ans[j]);//获得其他堆的石子所需要的操作次数
                sumsum -= k[i];//一定是减少了k[i]差的多少，可以想想
                if(sumsum <= 0) break;//如果不差了，可以结束了
            }
            Min = min(Min,cnt);//更新
        }
        cout << Min << endl;
    }
}