UVA11426 GCD - Extreme (II) 小于n的数对的gcd之和

http://www.bnuoj.com/v3/problem_show.php?pid=19844

 

输入正整数n，求gcd(1,2)+gcd(1,3)+...+gcd(1,n) + gcd(2,3)+...+g(2,n) + ...... + gcd(n-1,n)

即所有满足1<=i<j<=n的数对(i,j)对应的gcd(i,j)之和，

n=10答案为67 ， n=100答案为13015。

分析：设f(n)=gcd(1,n)+gcd(2,n)+...+gcd(n-1,n)，则所求答案S(n)=f(2)+f(3)+...+f(n)，

可以看出S(n)=S(n-1)+f(n)

解析f(n)：因为gcd(x,n)=i一定是n的约数，而且很多gcd(x1,n)=gcd(x2,n)=...=i，有多少这样的gcd呢?    

          gcd(x,n)=i  <==>  gcd(x/i,n/i)=1  ,所以满足条件的x/i有phi[n/i]个，即满足gcd(x,n)=i的x有phi[n/i]个。    

        所以f(n)=∑(i*phi[n/i])  | i是n的约数   

        由此我们可以想到类似素数筛法的方式预处理出f数组。     

        p之后S(n)=S(n-1)+f(n)预处理出S数组即可

 

#include<cstring>
#include<cstdio>
#define MT(x,i) memset(x,i,sizeof(x))
using namespace std;
typedef long long LL;
const int MAX = 4000000;

int phi[MAX+5];
void phi_table(int n)///欧拉函数表
{
    for(int i=2;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]) phi[j]=j;
            phi[j]=phi[j]/i*(i-1);
    }
}

LL S[MAX+5],f[MAX+5];
void init()
{
    phi_table(MAX);

    ///初始化f
    MT(f,0);
    for(int i=1;i<=MAX;i++)
        for(int n=i*2;n<=MAX;n+=i) f[n]+=(LL)i*phi[n/i];

    ///初始化S
    S[2]=f[2];
    for(int n=3;n<=MAX;n++) S[n]=S[n-1]+f[n];
}

int main()
{
    init();
    int n;
    while(scanf("%d",&n)&&n){
        printf("%lld\n",S[n]);
    }
    return 0;
}