未知：挖掘机技术哪家强——题解

题目描述
有人问现实中为什么总是男生追求女生，反过来很少。实际上女生也是想主动追求男生的，但是世俗中对于主动追求男生的女生有种歧视，这样就使得女生不大敢主动追求男生。但是面对喜欢的男生，难道就不出手么？女生只能步步为营，挖坑来引诱男生往里跳。这时候问题就来了，挖掘机技术到底哪家强？

被热血沸腾的广告洗脑了若干天后，Matt终于下定决心，毅然登上了开往泉城的列车，决心寻找生活的希望。

来到布鲁谢特学院后，Matt逐渐地了解了各种型号的挖掘机。在这里我们可以认为有大挖掘机和小挖掘机两种。

今天Matt的任务很简单：首先他要用大挖掘机挖出恰好N单位体积的砂土。由于小挖掘机比较笨拙，它每次挖的砂土体积是固定的。也就是说，设每次挖x单位体积砂土，那么N需要被x整除。在挖出若干堆体积为x的砂土后，Matt需要计算x的“难挖指数”。体积x的“难挖指数”定义如下：对于某个不超过x的体积y，如果x与y的最大公约数为1，那我们认为体积y是“难挖的”，x的“难挖指数”就要加上y。

由于Matt之后还需要用小挖掘机处理被大挖掘机挖出的砂土，他希望知道所有可能的x的难挖指数的和，这样他好估算今天要干多久，不然作为布鲁谢特的高才生，他出门要被笑话的。

输入
第一行一个整数T，表示数据组数。

接下来T行每行一个整数表示N。

输出
对于每个数据输出一行一个整数表示难挖指数的和。

样例输入
3
2
3
4
样例输出
2
4
6
提示
对于30%的数据有T<=20，N<=10^4。

对于60%的数据有T<=100，N<=10^7。

对于100%的数据有1<=T<=1000，1<=N<=10^9。

通过打表发现规律。
答案=（（n的所有的因数i的欧拉函数值*i）+1）/2（这里的+1是为了防止除不开（当然实际上一定会除不开所以加一没毛病））
证明：若a与n互质，那么n-a与n也互质。我们把a与n-a结成一对，那么n以内就存在phi(n)/2对，每一对的和都为n。
然后欧拉函数优化——欧拉筛法
然后就以飞速AC啦！

#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
using namespace std;
long long sum;
int su[1000000]={0},he[1000000]={0};
int cnt=0;
void eular(int n){
    if(n==1||n==2){
        sum++;
        return;
    }
    long long ans=n,a=n;
    for(int i=1;i<=cnt;i++){
        if(n%su[i]==0){
            ans=ans/su[i]*(su[i]-1);
            while(n%su[i]==0){
                n/=su[i];
            }
        }
        if(n==1)break;
    }
    if(n>1){
        ans=ans/n*(n-1);
    }
    sum+=ans*a/2;
    return;
}
void Euler(int n){
    for(int i=2;i<=n;i++){
        if(he[i]==0){
            cnt++;
            su[cnt]=i;  
        }
        for(int j=1;j<=cnt&&i*su[j]<=n;j++){
            he[su[j]*i]=1;
            if(i%su[j]==0)break;
        }
    }
}
int main(){
    Euler(100000);
    int t;
    scanf("%d",&t);
    while(t--){
        int n;
        scanf("%d",&n);
        sum=0;
        for(int i=1;i*i<=n;i++){
            if(n%i==0){
                eular(i);
                if(i!=n/i){
                    eular(n/i);
                }
            }
        }
        printf("%lld\n",sum);
    }
    return 0;
}