挖掘机技术哪家强

$$挖掘机技术哪家强$$

题目链接：jzoj 3858

题目

有人问现实中为什么总是男生追求女生，反过来很少。实际上女生也是想主动追求男生的，但是世俗中对于主动追求男生的女生有种歧视，这样就使得女生不大敢主动追求男生。但是面对喜欢的男生，难道就不出手么？女生只能步步为营，挖坑来引诱男生往里跳。这时候问题就来了，挖掘机技术到底哪家强？
被热血沸腾的广告洗脑了若干天后，$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）=n\ast \phi （n）/2$$
而在这其中，$\phi （n）$是欧拉函数：

（$p_i$为$x$的所有质因数）

我们只要先处理出素数，在枚举$N$能被那些数整除，加上它的难挖指数，就是答案了。

代码

#include<cmath>
#include<cstdio>
#define ll long long
#define max(x, y) (x) > (y) ? (x) : (y)

using namespace std;

ll T, n, ans, su_num[100001];
bool susu[100001];

ll phi(ll nn) {
	if (nn < 3) return 1;
	ll an = nn, temp_n = nn;
	for (ll i = 1; i <= su_num[0]; i++)
		if (nn % su_num[i] == 0) {
			an = an / su_num[i] * (su_num[i] - 1);//求出phi(nn)的值
			while (nn % su_num[i] == 0)
				nn /= su_num[i];
			if (nn == 1) break;
		}
	if (nn > 1) an = an / nn * (nn - 1);
	
	return an * temp_n / 2;//求出nn的难挖指数
}

void get_susu() {//素数筛选法
	susu[1] = 1;
	for (ll i = 2; i <= 100000; i++)
		if (!susu[i]) {
			su_num[++su_num[0]] = i;
			for (ll j = i + i; j <= 100000; j += i)
				susu[j] = 1;
		}
}

int main() {
	get_susu();//事前处理出素数表
	
	scanf("%lld", &T);//读入
	for (ll i = 1; i <= T; i++) {
		scanf("%lld", &n);//读入
		
		ans = 0;//初始化
		for (ll j = 1; j * j <= n; j++)//找到能整除n的数
			if (n % j == 0) {
				ans += phi(j);//求出这个数的难挖指数
				if (j * j != n) ans += phi(n / j);//这个数对应的另一个因子也求出难挖指数
			}
		
		printf("%lld\n", ans);//输出
	}
	
	return 0;
}