51nod1188 最大公约数之和V2（筛法 + 欧拉函数）

题目

给出一个数N，输出小于等于N的所有数，两两之间的最大公约数之和。

 

 

相当于计算这段程序（程序中的gcd(i,j)表示i与j的最大公约数）：

 

G=0;

for(i=1;i<N;i++)

for(j=i+1;j<=N;j++)

{

    G+=gcd(i,j);

}

 收起

输入

第1行：1个数T，表示后面用作输入测试的数的数量。(1 <= T <= 50000)
第2 - T + 1行：每行一个数N。(2 <= N <= 5000000)

输出

共T行，输出最大公约数之和。

输入样例

3
10
100
200000

输出样例

67
13015
143295493160

解题思路:做之前先看一下51nod1040,。对于可以用的时间打表出来。然后可以o(1)查询了。

代码:

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int maxn = 5e6 + 10;
ll phi[maxn],sum[maxn];
int prime[maxn],mark[maxn];

void getphi(){
	phi[1] = mark[1] = 1;
	for(int i = 2;i < maxn;++i){
		if(!mark[i])	prime[++prime[0]] = i,phi[i] = i - 1;
		for(int j = 1;j <= prime[0] && prime[j] * i < maxn;++j){
			mark[i * prime[j]] = 1;
			if(i % prime[j] == 0){
				phi[i * prime[j]] = phi[i] * prime[j];
				break;
			}
			phi[i * prime[j]] = phi[i] * (prime[j] - 1);
		}
	}
	
	for(int i = 1;i < maxn;++i)
		for(int j = i;j < maxn;j += i)
			sum[j] += (ll)i * phi[j / i];
	for(int i = 1;i < maxn;++i)	sum[i] += sum[i - 1];//前缀和
}

int main(){
	int T,N;
	getphi();

	cin>>T;
	while(T--){
		cin>>N;
		cout<<sum[N] - (ll)N * (N + 1) / 2 <<endl;
		
	}
	
	return 0;
}