欧拉函数 - GCD(最大公约数) - 洛谷 P2568

欧拉函数 - GCD(最大公约数) - 洛谷 P2568

$给定整数N，求1\le x,y\le N且GCD(x,y)为素数的数对(x,y)有多少对。$

$GCD(x,y)即求x，y的最大公约数。$

输入格式

输入一个整数N

输出格式

输出一个整数，表示满足条件的数对数量。

数据范围

$1\le N\le 10^7$

输入样例：

4

输出样例：

4

---

分析：

问题转化：

$求gcd(x,y)=P<=>gcd(\frac{x}{P},\frac{y}{P})=1，其中P是质数。$

$问题转化为：求1\le x,y\le \frac{N}{P}且GCD(x,y)=1的数对(x,y)有多少对。$

$求某段区间内，互质的点对的数量，参照：$欧拉函数 - Visible Lattice Points - POJ 3090

$所以，我们先预处理出不超过N的所有质数，对于每一个质数P，都需要计算一遍{\sum }_{i=2}^{\frac{N}{P}}\varphi (i)\times 2+1$

$注意：本题需要单独加上y=x上的点(1,1)。$

$由于需要多次查询区间和，故我们另预处理前缀和数组S[i]={\sum }_{j=1}^i\varphi (j)，$

$注意\varphi (1)=0，因为点(1,1)是单独处理累加的。$

$最终答案：{\sum }_{i=1}^{cnt}2\times S[\frac{N}{P_i}]+1，其中cnt是N以内质数的总个数，P_i表示第i个质数。$

代码：

#include<iostream>

#define ll long long

using namespace std;

const int N=1e7+10;

int primes[N],cnt;
bool st[N];
int phi[N];
ll s[N];

void get_prime(int n)
{
    for(int i=2;i<=n;i++)
    {
        if(!st[i])
        {
            primes[cnt++]=i;
            phi[i]=i-1;
        }
        for(int j=0;primes[j]*i<=n;j++)
        {
            st[primes[j]*i]=true;
            if(i%primes[j]==0)
            {
                phi[i*primes[j]]=phi[i]*primes[j];
                break;
            }
            phi[i*primes[j]]=phi[i]*(primes[j]-1);
        }
    }
    
    for(int i=1;i<=n;i++) s[i]=s[i-1]+phi[i];
}

int main()
{
    int n;
    cin>>n;
    get_prime(n);
    
    ll res=0;
    for(int i=0;i<cnt;i++)
    {
        int p=primes[i];
        res+=(s[n/p]*2)+1;
    }
    
    cout<<res<<endl;
    
    return 0;
}