gcd的小特性

由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;
}

.