bzoj 2818 Gcd（莫比乌斯反演）

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

裸莫比乌斯反演。

设F(n) 为gcd是n的倍数的个数，f(n)为gcd是n的个数,

直接就求  sigma(f(p[i]))就行了，p[i]是质数。

莫比乌斯反演有两种形式，这里我们使用第二种：

一、

二、

#include
#include
#include
#include
#define LL long long
using namespace std;


const int maxn = 1e7 + 10;
int p[maxn/10];
int flag[maxn];
int mu[maxn];
int cnt = 0;
void init()
{
    int i,j;
    mu[1] = 1;
    for(i=2;i