poj 2478 Farey Sequence

求n以内中所有互质对数的个数

需要用类似于埃筛的方法来找n个数每个数比自己小且互质数的个数  

#include<iostream>
#include<cstring>
#include<cstdio>
using namespace std;
long long int ans[1000002];
int main(){
	memset(ans,0,sizeof(ans));
	for(int j=2;j<=1000000;j++)
   { 
     if(ans[j]==0)
	 {
	    for(int i=j;i<=1000000;i+=j)
		{
			if(ans[i]==0)ans[i]=i;
			ans[i]=ans[i]*(j-1)/j;
			 } 	
	}   
   }
   int n;
   while(scanf("%d",&n))
   {
   	if(n==0)break;
   	   long long int an=0;
   	   for(int i=2;i<=n;i++)
   	   an+=ans[i];
   	   printf("%lld\n",an);
   }
	return 0;
}
//需要用类似于埃筛的方法来找n个数每个数比自己小且互质数的个数