UVa 11426--算法竞赛入门经典训练指南

题目大意：输入正整数n，求所有满足1<=i<j<=n的数对(i,j)所对应的gcd(i,j)之和。

题解：设f(n)=gcd(1,n)+gcd(2,n)+……+gcd(n-1,n),则所求答案S(n)=f(2)+f(3)+……+f(n)

注意到gcd(x,n)的值都是n的约数，可以按照这个约数进行分类，用g(n,i)表示满足gcd(n,x)=i且x<n的正整数x的个数，则f(n)=sum(i*g(n,i)|i是n的约数)。注意到gcd(n,x)=i的充要条件是gcd(n/i,x/i)=1,因此满足条件的x/i有phi(n/i)个没说明g(n,i)=phi(n/i);

求f(n)时，对于每个i枚举它的倍数n(并更新n的值)，类似于素数筛法求解。

#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <cmath>
#include <algorithm>

using namespace std;

typedef long long LL;

const int N=4000001;

int phi[N],n;
LL f[N],s[N];
void phi_table(int n)
{
    for (int i=1;i<=n;i++)phi[i]=0;
    phi[1]=1;
    for (int i=2;i<=n;i++)if (!phi[i])
    for (int j=i;j<=n;j+=i){
        if (phi[j]==0)phi[j]=j;
        phi[j]=phi[j]/i*(i-1);
    }
    //for (int i=1;i<=100;i++)printf("%d %d\n",i,phi[i]);
}
int main()
{
    phi_table(N-1);
    for (int i=1;i<N;i++)for (int j=i+i;j<N;j+=i)f[j]+=i*phi[j/i];
    for (int i=2;i<N;i++)s[i]=s[i-1]+f[i];

    while (scanf("%d",&n)!=EOF && n){
        cout<<s[n]<<endl;
    }
    return 0;
}