由UVA - 11426学到的小知识
由于数据量为4×le6,我们暴力跑答案不算gcd的时间的话还要1e13。大致需要10000s.所以我们需要预处理一些值。
直接就往分块想,我们知道gcd(1,n)、 gcd(2,n)……………gcd(n-1,n)里有很多对的gcd是相同的。而gcd(x,n)=i等价于gcd(x/i,n/i)=1,因此gcd为i的个数为phi(n/i)。因此我们可以通过分块处理它们gcd的和。
#include<bits/stdc++.h>
typedef long long ll;
const int inf = (1<<31)-1;
const ll INF = (1ll<<63)-1;
const int maxn = 4*1e6+1;
bool check[maxn];
ll phi[maxn],Prime[maxn],f[maxn],ans[maxn];
void init(){
int N = maxn-1;
memset(check,false,sizeof(check));
phi[1] = 1;
int tot = 0;
for(int i = 2;i <= N; ++i){
if(!check[i]){
Prime[tot++] = i;
phi[i] = i-1;
}
for(int j = 0;j < tot; ++j){
if(i*Prime[j] > N) break;
check[i*Prime[j]] = true;
if(i%Prime[j] == 0){
phi[i*Prime[j]] = phi[i]*Prime[j];
break;
}
else{
phi[i*Prime[j]] = phi[i]*(Prime[j]-1);
}
}
}
for(int i=1;i<maxn;i++)
for(int j=i+i;j<maxn;j+=i)
f[j]+=i*phi[j/i];
ans[1]=f[1];
for(int i=2;i<maxn;i++)
ans[i]=ans[i-1]+f[i];
}
int main()
{
int n;
init();
while(scanf("%d",&n)!=EOF&&n!=0)
{
printf("%lld\n",ans[n]);
}
return 0;
}
.