hdu6237-2017哈尔滨CCPC-质因子&贪心&暴力-A Simple Stone Game

http://acm.hdu.edu.cn/showproblem.php?pid=6237
当我们枚举素因子的时候，
我们统计mod的和。sumb
并且对mod的结果搞一个b数组。
然后队友的思路是直接分为两半，然后左边那一半小的往右边那一半靠。
以下疑问：
①可能转移的时候中介线并不是一半。。
因为我们怎么确定b数组在范围内是线性的分布？
但是当时写的也有问题。没有考虑到贪心，
然后 最优的方案是让最大的变成pi（那个素因子）。那种花费最小。
可以直接倒着算，计算以下各个mod的和，后面的大数（升序）每搞一个，sumb就减去 因子，ans加上 因子和大数之差。

就是让小的填充到绿色部分。（红线是mod）。 这样花费最低（如果让大的取填小的，看着就多。）
有一个简短的数学证明（思维）。这样处理是肯定可以得到结果的。
因为sumb 是pi的倍数，所以b数组中的每一个数mod pi，和也一定是倍数。因此这种处理方法是有效的，对所有情况都适用的。
有一个优化，直接从大到小，sumb-=pi。ans+=vec[i]-p[j]。这样少一个循环。并且很多代码都是这样qwq
g++比c++快。
但是g++得用i64d。。。。。！！！

  #include <iostream>
#include <cstring>
#include <cstdio>
#include <vector>
#include <algorithm>
using namespace std;
/* 为什么我这么菜。
   cf上做过一道类似的。和这一道题意思差不多。
   都有贪心的一些意思。我怎么就没想到。
   不过那个是修改（本质也是让其花费尽可能的小）

   当我们枚举素因子的时候，
   我们统计mod的和。sumb
   并且对mod的结果搞一个b数组。
   然后队友的思路是直接分为两半，然后左边那一半小的往右边那一半靠。
   以下疑问：
   ①可能转移的时候中介线并不是一半。。
   因为我们怎么确定b数组在范围内是线性的分布？
   但是当时写的也有问题。没有考虑到贪心，
   然后 最优的方案是让最大的变成pi（那个素因子）
*/
typedef long long ll;
ll sum;
const int maxn=1e5+200;
ll a[maxn];
vector<ll>vec;
vector<ll>q;
int main()
{   int t,m;
     scanf("%d",&t);
     while(t--){
           sum=0;
           scanf("%d",&m);
           for(int i=0;i<m;i++){
               scanf("%I64d",&a[i]);
               sum+=a[i];
           }
         vec.clear();
         q.clear();
         //素因子这个也没有想到。。
          for(ll i=2;i*i<=sum;i++)           //求出总和的质因子
        {
            if(sum%i==0)
            {
                vec.push_back(i);
                while(sum%i==0)
                {
                    sum/=i;
                }
            }
        }
        if(sum>1)
            vec.push_back(sum);
        ll res=1e17;
        for(int i=0;i<vec.size();i++){

            q.clear();
            ll ans=0;
             sum=0;
            for(int j=0;j<m;j++){
                 ll b=a[j]%vec[i];

                 q.push_back(b);
                 sum+=b;
            }
            int siz=sum/vec[i];
            int  siz2=0;
            ll sum2=0;
            sort(q.begin(),q.end());
            for(int j=q.size()-1;j>=0;j--){
                sum2+=(vec[i]-q[j]);
                siz2++;
                if(siz2==siz) break;
            }
            for(int j=0;j<q.size();j++){
                ans+=q[j];
                sum2-=q[j];
                if(sum2<=0) break;
            }
            res=min(res,ans);
        }
        printf("%I64d\n",res);//i64d。！！！
     }
    return 0;
}