bzoj 2818 gcd 线性欧拉函数

2818: Gcd

Time Limit: 10 Sec  Memory Limit: 256 MB
[Submit][Status][Discuss]

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（x,y)=p;p为素数；

　　　gcd（x/p,y/p)==1;

　　　一个p的贡献为1-（N/p)，求前缀和意思；

#include<bits/stdc++.h>
using namespace std;
#define ll long long
#define mod 1000000009
#define inf 999999999
#define esp 0.00000000001
//#pragma comment(linker, "/STACK:102400000,102400000")
const int N=1e5,M=1e7+;
ll p[M],ji;
bool vis[M];
ll phi[M];
void get_eular(int n)
{
    ji = ;
    memset(vis, true, sizeof(vis));
    for(int i = ; i <= n; i++)
    {
        if(vis[i])
        {
            p[ji ++] = i;
            phi[i] = i - ;
        }
        for(int j = ; j < ji && i * p[j] <= n; j++)
        {
            vis[i * p[j]] = false;
            if(i % p[j] == )
            {
                phi[i * p[j]] = phi[i] * p[j];
                break;
            }
            else
            phi[i * p[j]] = phi[i] * phi[p[j]];
        }
    }
}
ll sumphi[M];
int main()
{
    ll x,y,z,i,t;
    get_eular();
    sumphi[]=;
    for(i=;i<=;i++)
    sumphi[i]=sumphi[i-]+phi[i];
    while(~scanf("%lld",&x))
    {
        ll ans=;
        for(i=;p[i]<=x;i++)
        ans+=*sumphi[x/p[i]]-;
        printf("%lld\n",ans);
    }
    return ;
}