BZOJ 2749 && Comet OJ C1416 [HAOI2012]外星人

题目链接：BZOJ，Comet OJ

根据唯一分解定理和欧拉函数在自变量互质时的积性性质
我们可以推出$\varphi (n)=\varphi (p_1^{q_1})\ast \varphi (p_2^{q_2})\ast \dots \ast \varphi (p_m^{q_m})$
现在可以单独考虑每个$\varphi (p_i^{q_i})$对答案的贡献
最后把$1-m$每个$p_i^{q_i}$的贡献加起来就可以
根据题目下方的$Hint$
$\varphi (p^q)$就等于$(p-1)\ast \varphi (p^{q-1})$
题目要求多少开次欧拉函数之后结果为$1$
而只有$\varphi (2)$等于$1$
我们至少需要开$q$次才能把后面的那部分开完
对于前面的那部分，如果$(p-1)$是质数，那么$\varphi (p-1)=p-2$
否则就会被分成其他的合数
最后的要开的次数就是能开出多少个$2$，因为$\varphi (2)$才等于$1$
那么我们现在需要统计每个数能这么“开”出多少个$2$
设$f[i]$表示$i$能“开”出$f[i]$个$2$，
$a$与$b$互质时，$f[a\ast b]=f[a]+f[b]$：$i$为质数时，$f[i]=f[i-1]$

#include <bits/stdc++.h>
#define A 100010
#define B 100000

using namespace std;
typedef long long ll;
int T, m, pri[A], cnt;
bool bz[A]; ll f[A], p, q;

int main(int argc, char const *argv[]) {
	f[1] = 1;
	for (int i = 2; i <= B; i++) {
		if (!bz[i]) pri[++cnt] = i, f[i] = f[i - 1];
		for (int j = 1; j <= cnt and i * pri[j] <= B; j++) {
			bz[i * pri[j]] = 1;
			f[i * pri[j]] = f[i] + f[pri[j]];
			if (i % pri[j] == 0) break;
		}
	}
	cin >> T;
	while (T--) {
		scanf("%d", &m); ll ans = 1;
		while (m--) {
			scanf("%lld%lld", &p, &q);
			if (p == 2) ans--;
			ans += f[p] * q;
		}
		printf("%lld\n", ans);
	}
	return 0;
}