hdu-6237

A Simple Stone Game

Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 262144/262144 K (Java/Others)
Total Submission(s): 1247    Accepted Submission(s): 287

Problem Description

After he has learned how to play Nim game, Bob begins to try another stone game which seems much easier.

The game goes like this: one player starts the game with N piles of stones. There is ai stones on the i th pile. On one turn, the player can move exactly one stone from one pile to another pile. After one turn, if there exits a number x(x>1) such that for each pile bi is the multiple of x where bi is the number of stone of the this pile now), the game will stop. Now you need to help Bob to calculate the minimum turns he need to stop this boring game. You can regard that 0 is the multiple of any positive number.

 

Input

The first line is the number of test cases. For each test case, the first line contains one positive number N(1≤N≤100000) , indicating the number of piles of stones.

The second line contains N positive number, the i th number ai(1≤ai≤100000) indicating the number of stones of the i th pile.


The sum of N of all test cases is not exceed 5∗105 .

 

Output

For each test case, output a integer donating the answer as described above. If there exist a satisfied number x initially, you just need to output 0 . It's guaranteed that there exists at least one solution.

 

Sample Input

2 5 1 2 3 4 5 2 5 7

 

Sample Output

2 1

 

Source

2017中国大学生程序设计竞赛-哈尔滨站-重现赛（感谢哈理工）

 

铜牌题坑就是多呀，注意读懂题意，sum是会爆int的，还有sum本身就有可能是个素数，其他的就没啥坑了。

 

解题思路：求出石子总数sum，能整除sum的质数肯定满足x，遍历1-1e5的所有质数，统计满足x的有多少个。然后遍历满足x的质数，求出最小的移动次数，与sum减去最大堆（sum为质数的情况）的值作比较，取小的便是最终答案。具体操作看代码，代码比文字好理解(●'◡'●)

 

#include<iostream>
#include<string.h>
#include<algorithm>
#include <cstdio>
using namespace std;
typedef long long ll;
const int maxn = 1e5+10;
int nu[maxn];
int prim[maxn];
bool isprim[maxn];
int res[maxn];
int cop[maxn];
int len;
void euler(int maxx)
{
    memset(isprim,true,sizeof(isprim));
    isprim[1]=false;
    len=0;
    for(int i=2;i<=maxx;++i)
    {
        if(isprim[i]) prim[++len]=i;
        for(int j=1;j<=len && i*prim[j]<=maxx;++j)
        {
            isprim[i*prim[j]]=false;
            if(i%prim[j]==0) break;
        }
    }
}
int main()
{
    euler(1e5);
    int t;
    scanf("%d",&t);
    while(t--)
    {
        int n;
        scanf("%d",&n);
        ll sum=0;
        for(int i=0;i<n;++i)
        {
            scanf("%d",&nu[i]);
            sum+=nu[i];
        }
        sort(nu,nu+n);
        ll ans = sum-nu[n-1];
        int le=0;
        for(int i=1;i<=len;++i)
        {
            if(sum%prim[i]==0)
            {
                res[++le]=prim[i];
                while(sum%prim[i]==0)
                    sum/=prim[i];
            }
        }

        for(int i=1;i<=le;++i)
        {
      //      printf("%di ",res[i]);
            for(int j=0;j<n;++j)
            {
                cop[j]=nu[j]%res[i];
            //    printf("%d ",cop[j]);
            }
            sort(cop,cop+n);
            ll cnt=0;
            int l=0,r=n-1;
            while(l<r)
            {
       //         printf("@");
                if(cop[l]==0)
                {
                    l++;
                    continue;
                }
                if(cop[l]+cop[r]==res[i])
                {
                    cnt+=cop[l];
                 //   printf("
%d##
",cnt);
                    l++;
                    r--;
                }
                else if(cop[l]+cop[r]<res[i])
                {
                    cop[r]+=cop[l];
                    cnt+=cop[l];
                    l++;
                }
                else if(cop[l]+cop[r]>res[i])
                {
                    cnt+=res[i]-cop[r];
                    cop[l]-=(res[i]-cop[r]);
                    r--;
                }
            }
         //   printf("
%d##
",cnt);
            ans=min(cnt,ans);
         //   printf("
%d###
",ans);
        }
        printf("%lld
",ans);
    }
    return 0;
}