HYSBZ - 2818

problem

给定整数N，求1<=x,y<=N且Gcd(x,y)为素数的
数对(x,y)有多少对.
N1e7

思路

枚举素数（大概几e5个） 对于每个素数暴力算（先预处理筛出欧拉值）
时间复杂度 case∗O(素数个数) c a s e ∗ O ( 素 数 个 数 )

代码示例

#include
using namespace std;
typedef long long ll;

const int maxn=1e7+10;
bool valid[maxn];
int phi[maxn];
int ans[maxn/10];
ll sum[maxn];

void getprime(int n,int &tot,int ans[])
{
    tot=0;
    phi[1]=1;
    //memset(valid,true,sizeof(valid));
    for(int i=1;itrue;
    for(int i=2;i<=n;++i){
        if(valid[i]){
            tot++;
            ans[tot]=i;
            phi[i]=i-1;//i为素数
        }
        //下面的主角是小于等于i的每个质数
        for(int j=1;(j<=tot) && (i*ans[j]<=n);++j){
            valid[i*ans[j]]=false;
            if(i%ans[j]==0){//ans[j]是i的素因子
                phi[i*ans[j]]=phi[i]*ans[j];
                break;//如果整除就break;
            }
            else phi[i*ans[j]]=phi[i]*(ans[j]-1);
        }
    }
}

int main()
{
    int tot;
    getprime(1e7,tot,ans);
    //cout<
    //cout<
    for(int i=1;i<=int(1e7);++i){
        sum[i]=sum[i-1]+phi[i];
    }
    ios::sync_with_stdio(false);
//    int t;
//    cin>>t;
//    while(t--)
//    {
        ll n;
        ll res=0;
        cin>>n;
        int num=0;
        for(int i=1;ans[i]<=n;++i){
            res+=2*sum[n/ans[i]]-1;
        }
        cout<//    }
    return 0;
}