GCD - Extreme (II) UVA - 11426 最大公约数之和

题目链接

输入正整数n ， 求gcd(1 ,2)+gcd(1 ,3)+gcd(2,3)+. . . +gcd(n-1，n)时，即所有满足1 <= i <j<=n 的数对(i,j)所对应的gcd(i,j)之和。比如n= l0 时答案为67 ， n=100 时答案为13 015 , n=200 000时答案为143295493 160 。
【输入格式】
输入包含不超过100 组数据，每组数据占一行，包含正整数n (2 <= n <=4000000) 。输入结束标志为n=0 。
【输出格式】
对于每组数据，输出一行，即所求和。答案保证在64 位带符号整数范围内。
【分析】
设f(n)=gcd(1 ,n)+gcd(2 ,n)+gcd(3 ,n)+...+gcd(n-1, ,n)，则所求答案为S(n)=f(2)+f(3)+...f(n) 。只需求出f(n) ， 就可以递推出所有答案: S(n)=S(n-1)+f(n) 。因此，本题的重点在于如何计算f(n) 。
注意到所有gcd(x,n)的值都是n 的约数，可以按照这个约数进行分类，用g(n， i)表示满足gcd(x，n)=i 且x<n 的正整数x 的个数，则f(n)=sum{ixg(n, i) | i 是n 的约数}。注意到gcd(x，n)=i的充要条件是gcd(x/i，n/i)=1(此时x/i和n/i互质）， 因此满足条件的x/i 有phi(n/i)个，说明g(n ， i)=phi(n/i) 。问题到这里还没有结束。如果依次计算j例，需要对每个n 枚举它的约数i ， 速度较慢，但如果把思路逆转过来，对于每个i 枚举它的倍数n (并且更新j(n) 的值) ，时间复杂度将降为与素数筛法同阶。至此，问题得到了完整的解决。代码如下。

#include <cstdio>
const int N = 4000000;
typedef long long LL;
int phi[N+5];
LL f[N+5], s[N+5];

void phi_table(){
	phi[1] = 1;
	for(int i = 2; i <= N; i++){
		if(!phi[i]){
			for(int j = i; j <= N; j += i){
				if(!phi[j]) phi[j] = j;
				phi[j] = phi[j] / i * (i-1);
			}
		}
	}
} 

int main(int argc, char** argv) {
	phi_table();
	for(int i = 1; i <= N; i++) // 预处理f
		for(int n = i*2; n <= N; n += i)
			f[n] += i * phi[n/i];
	s[2] = f[2]; // 预处理s
	for(int n = 3; n <= N; n++) 
		s[n] = s[n-1] + f[n];
	int n;
	while(~scanf("%d",&n) && n)
		printf("%lld\n",s[n]);
	return 0;
}