Description

給定整數N，求1<=x,y<=N且Gcd(x,y)為素數的
數對(x,y)有多少對.

Input

一個整數N

Output

如題

Sample Input

4

Sample Output

4

Hint

hint

對於樣例(2,2),(2,4),(3,3),(4,2)

1<=N<=10^7

GCD（px,py）==p;
有GCD（x,y）==1,在1-n/p的情況數。即1-n/p的尤拉和。p為質數。
所以直接模板套上篩。。就可以了。。
注意（2,4）（4,2）是兩種情況要*2.。但是（2,2）這種有不是，又要剪掉。。

#include<bits/stdc++.h>
 

using namespace std;
#define LL long long

const int mod=1000000007;
const int maxn=10000005;
bool check[maxn];          //用於打表記錄的中間量
LL sumPhi[maxn];           //前i個的尤拉函式和
int cnt,phi[maxn],prime[maxn];  //素數個數，尤拉表，素數表
//素數表是第幾個素數是什麼，尤拉表是i的尤拉是phi[i];
void init(int n){               //素數+尤拉表
    phi[1]=1;
    cnt=0;
    for
 
(int i=2;i<=n;i++){
        if(!check[i]){
            phi[i]=i-1;
            prime[cnt++]=i;
        }
        for(int j=0;j<cnt;j++){
            if(i*prime[j]>n)break;
            check[i*prime[j]]=true;
            if(i%prime[j]==0){
                phi[i*prime[j]]=phi[i]*prime[j];
                break
 
;
            }
            else{
                phi[i*prime[j]]=phi[i]*(prime[j]-1);
            }
        }
    }
    sumPhi[0]=0;
    for(int i=1;i<=n;i++) sumPhi[i]=(sumPhi[i-1]+phi[i]);
}

int main(){
    int n;
    while(cin>>n){
        LL ans=0;
        init(n);
        for(int i=0;i<cnt;i++){
            ans+=sumPhi[n/prime[i]];
        }
        cout<<2*ans-cnt<<endl;
    }
    return 0;
}