文章标题 51nod 最大公约数之和

51nod 最大公约数之和

先是简单版本：
51nod 1040 最大公约数之和
题目链接：最大公约数之和
Description
给出一个n，求1-n这n个数，同n的最大公约数的和。(n<=10^9)
Solution
最大公约数为n的约数。对于每个约数i，找[1,n]区间GCD(x,n)==i
的x个数即可，转化为 [1,n/i]区间内与n互质的数个数，即欧拉函数的
值，累加即可。

Code

#include <iostream>
#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <algorithm>
#define ll long long
using namespace std;
int phi(int v);
int main()
{
    int n;
    cin>>n;
    ll ans=0;
    for(int i=1;i*i<=n;++i)
        if(n%i==0){
            if(i*i==n){
                ans+=i*phi(i);
                break;
            }
            else{
                ans+=i*phi(n/i);
                ans+=n/i*phi(i);
            }
        }
    cout<<ans<<endl;
    return 0;
}
int phi(int v)
{
    int ans=v;
    for(int i=2;i*i<=v;++i)
        if(v%i==0){
            ans=ans-ans/i;
            v/=i;
            while(v%i==0)
                v/=i;
        }
    if(v!=1)
        ans=ans-ans/v;
    return ans;
}

第二题：
51nod 1188 最大公约数之和 V2

题目链接：
1188 最大公约数之和 V2

Solution

题目要求i < j，对于每一个j，即求[1,j-1]中每一个数与j的最大公约数之和，这正是上题的版本，
只是不取上界j，于是：

$ans=$$\sum_{j=2}^n \sum_{i=1}^{j-1} \gcd (i,j)$

$=\sum_{j=2}^n \sum_{i|j,i<j} i*phi(\frac{j}{i})$

$=\sum_{i=1}^{\lfloor n/2\rfloor} \sum_ {k=2}^{\lfloor n/i\rfloor}i*phi(k)$

即预处理N=5000000以内的欧拉函数值，枚举最大公约数i，显然只需到N/2，内层从i的2倍开始枚举，
这样，用一个数组保存答案，询问的时候直接输出。

Code

#include <iostream>
#include <cstdio>
#define ll long long
using namespace std;
const int N=5000000;
const int M=348513;
int phi[N+2],prime[M+2];
bool use[N+2];
ll ans[N+2];
void phi_table();
int cnt;
int main()
{
    int n,t;
    phi_table();
    int v=N>>1;
    for(int i=1;i<=v;++i)
        for(int j=2;j<=N/i;++j)
            ans[i*j]+=phi[j]*i;
    for(int i=1;i<=N;++i)
        ans[i]+=ans[i-1];
    // cout<<cnt<<endl;
    cin>>t;
    while(t--){
        cin>>n;
        cout<<ans[n]<<endl;
    }
    return 0;
}
void phi_table()
{
    phi[1]=1;
    for(int i=2;i<=N;++i){
        if(!use[i]){
            prime[cnt++]=i;
            phi[i]=i-1;
        }
        for(int j=0;j<cnt&&i*prime[j]<=N;++j){
            use[i*prime[j]]=1;
            if(i%prime[j]==0){
                phi[i*prime[j]]=phi[i]*prime[j];
                break;
            }
            else
                phi[i*prime[j]]=phi[i]*(prime[j]-1);
        }
    }
}

第三题再补，第一次写博客啊，请见谅。