bzoj2818 Gcd

题目链接：http://www.lydsy.com/JudgeOnline/problem.php?id=2818

题目大意：给你一个n，让你求使得gcd(x,y)为素数的有序数对(x,y）的个数，其中1<=x,y<=n。n<=1e7

题目分析：我们枚举1到n内的素数(gcd(x,y)的值)，对于每一个素数p，我们需要找到数对(a,b)的个数，其中a和b互质且p*a和p*b均小于n。那么我们先用筛法求出phi和素数。并求出phi的前缀和，这时，我们枚举素数p，给答案ans加上sum[n/phi[p]]。这样计算出来的ans不是最终答案，因为我们要求的是有序数对x,y的个数，而之前的计算方法a,b 和b,a之计算了一次。我们首先排除掉a=b的情况，这一共有素数个数个，然后乘2，再加上a=b的个数即可。

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<cmath>
#include<iostream>
#include<algorithm>
#define maxn 11000000
using namespace std;
typedef long long LL;
int read()
{
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
    return x*f;
}
int minDiv[maxn],prim[1000000];
LL sum[maxn],phi[maxn];
bool mark[maxn];
int n,primm=0; 
void getphi()
{
	phi[1]=1;
	for(int i=2;i<=n;i++)
	{
		if(!mark[i])
		{
			phi[i]=i-1;
			prim[++primm]=i;
		}
		for(int j=1;j<=primm;j++)
		{
			int x=prim[j];
			if(i*x>n)break;
			mark[i*x]=1;
			if(i%x==0)
			{
				phi[i*x]=phi[i]*x;
				break;
			}
			else phi[i*x]=phi[i]*phi[x];
		}
	}
}
int main()
{
 
    scanf("%d",&n); 
    getphi();
    sum[0]=0;
    for (int i=1;i<=n;i++)
    { 
        sum[i]=sum[i-1]+phi[i];
	}
	LL ans=0;
	for (int i=1;i<=primm;i++)
	{
		ans+=sum[n/prim[i]];
	}
	cout<<(ans-primm)*2+primm<<endl;
	return 0;
}