BZOJ 2818 Gcd

传送门
题目大意就是求在给定范围内
$gcd(x,y)=prime$
的所有的有序$(x,y)$的对数

分析
$\text{我们令}gcd(x,y) =k,k\text{是质数，那么就有}\\gcd(x/k,y/k)=1\\$
$\text{不妨先让}(x,y)\text{中的}y>x\\\text{问题就变成了求}y/k\text{中与其互质的数的个数}$
$\text{显然是}\phi(y/k)\\\text{所以我们只需要对}\phi\text{求一个前缀和，然后扫一遍}N\text{之前的全部质数统计下就行了}$
下面是代码

#include <stdio.h>
#include <string.h>
#include <iostream>
#include <algorithm>
#define N 10000000+5
using namespace std;

int prime[N],cnt;
bool f[N];
int phi[N];
long long sum[N];
int n;
void init()
{
    phi[1]=1;
    register int i;
    for(i=2;i<=n;++i)
    {
        if(!f[i])
        {
            prime[++cnt] = i;
            phi[i] = i-1;
        }
        for(int j=1;(long long)prime[j]*i <= n; ++j)
        {
            f[prime[j]*i] = 1;
            if(i%prime[j]==0)
            {
                phi[prime[j]*i] = phi[i] * prime[j];
                break;
            }
            else
                phi[prime[j]*i] = phi[i] * (prime[j]-1);
        }
    }
}

int main()
{
    cin>>n;
    init();
//  for(int i=1;i<=10;++i)
//s     cout<<phi[i]<<" ";
    for(int i=1;i<=n;++i)
        sum[i] = sum[i-1] + phi[i];
    long long ans = 0;
    for(int i=1;i<=cnt;++i)
        ans += sum[n/prime[i]]*2-1;//1 1 ->-1
    cout<<ans;
}