uva 11426 欧拉函数的运用

题意： 求 n 以内gcd( i , j )  i < j 的和的个数。

n 的欧拉函数值 表示 小于 n 与 n 互质的数的个数。

对于互质的两个数 a,b (a<b<=n),  gcd（a,b） = 1 是显然的。

对 b 而言互质的数有 fai(b) = m个 ，那么b 的贡献就是 m。如果是 2*b呢？ 那贡献就是2*m

因为将与 b 互质的m个数都乘2， 就得到 m 个 gcd(2*a,2*b) = 2的数

但2*b不仅仅有 2*m个，因为它也可以通过别的数相乘得到 。那么只需要求出每个数的欧拉函数，

再对 i * j 下标累加 j * fai(i) 个贡献就行，j从1倍一直到小于maxn。

求出的是 i 以内的gcd()和，前 i 项gcd() 和累加就行。

//uva 11426 欧拉函数的运用 
#include <iostream>
#include <algorithm>
#include <cstring>
#define ll long long
using namespace std;
const int maxn = 4000051;
bool is_prime[maxn];
int prime[maxn],phi[maxn];
ll a[maxn];
void getphi(){
	memset(is_prime,0,sizeof(is_prime));
	phi[1] = 1;
	for(int i=2;i<maxn;i++){
		if(!is_prime[i]){
			prime[++prime[0]] = i;
			phi[i] = i-1;
		}
		for(int j=1;j<=prime[0]&&i*prime[j]<maxn;j++){
			is_prime[prime[j]*i] = 1;
			if(i%prime[j]==0){
				phi[i*prime[j]] = phi[i]*prime[j];
			}else{
				phi[i*prime[j]] = phi[i]*(prime[j]-1);
			}
		}
	}
}
void solve(){
	getphi();
	for(int i=2;i<maxn;i++){
		for(int j=1;j*i<maxn;j++){
			a[j*i]+=j*phi[i];
		}
	}
	for(int i=1;i<maxn;i++)
		a[i]+=a[i-1];
}
int main() {
	solve();
	int n;
	while(scanf("%d",&n)&&n){
		printf("%lld\n",a[n]);
	}
	return 0;
}