Farey Sequence POJ - 2478

点击打开链接

这个题是数论题 一开始 感觉 emmmm 挺简单的  后来才发现 这道题 贼 恐怖 超时超的 怀疑人生  最后我看了别人的 思路  别人的代码 总是那么厉害。。。 自己理解了一下 发现这道题 还是比较好理解的 主要 还是欧拉函数  主要思路是 先用一个 初始化 所有ans[i]=i  ans[i] 代表 i之前的所有的与它最大公约数为一的所有数字 然后  最后只需要 将n前所有的数组加起来就行 怎么处理ans 这个是重要的点  其实 只要将减去所有的素数的即 ans[j]-ans[j]/素数  就是ans[j],这个 就是ans[j]的所有数值 。。

下面就是代码

#include<stdio.h>
#include<iostream>
#include<string.h>
#include<algorithm>
#include<string>
#include<math.h>
using namespace std;
const int n=1e6+10;
bool pem[n+100];
int ans[n+100];
void init()
{
    for(int i=0;i<n;i++)
        ans[i]=i;
        pem[0]=pem[1]=0;
    long long int i,j,k=sqrt(n);
    for(i=2; i<=k; i++)
    {
         if(!pem[i])
         {
              for(j=i*i;j<n;j+=i)
              {
                    pem[j]=1;
              }
         }
    }
     for(i=2;i<n;i++)
     {
           if(!pem[i])
           {
                for(j=i;j<n;j+=i)
                {
                      ans[j]=ans[j]/i*(i-1);
                }
                //printf("%d\n",1);
           }
     }
    // printf("1\n");
}
int main()
{
    int t,i;
    long long int sum;
    init();
    while(~scanf("%d",&t)&&t)
    {
       sum=0;
       for(i=2;i<=t;i++)
          sum+=ans[i];
       printf("%lld\n",sum);
    }
    return 0;
}