poj 3090 http://write.blog.csdn.net/postedit

//嘿嘿....做数论的题目现在有点感觉了。
//把学到的内容举一反三 然后 在宽展一部分 你一定会提高
//题意：给你一个n*n的点阵 你站在（0 0）点的位置 问 你能看到几个点
//思路：图对称 所以求出一半乘以2然后加上中间的点即可。
#include<stdio.h>
#include<string.h>
#include<math.h>
#define max 1010
int is_prime[max];
int Euler_num[max];
void if_prime()
{
    memset(is_prime,0,sizeof(is_prime));
    int tem=(int)sqrt(max*1.0);
    for(int i=2; i<=tem; i++)
        if(is_prime[i]==0)
            for(int j=i*i; j<max; j+=i)
                is_prime[j]=1;
}
void Euler()//利用欧几里得的性质求解。
{
    for(int i=1; i<=max; i++)
        Euler_num[i]=i;
    for(int j=2; j<=max; j++)
    {
        if(is_prime[j]==0)
            for(int i=j; i<=max; i+=j)
            {
                Euler_num[i]=Euler_num[i]/j*(j-1);
            }
    }
}
int main()
{
    int n;
    int nn;
    __int64 ans;
    if_prime();
    Euler();
    int step;
    scanf("%d",&n);
    for(int step=1; step<=n; step++)
    {
        scanf("%d",&nn);
        ans=0;
        for(int i=1; i<=nn; i++)
        {
            ans+=Euler_num[i];
        }
        printf("%d %d %I64d\n",step,nn,ans*2+1);
    }
    return 0;
}