欧拉函数及欧拉函数表

欧拉函数
Phi（n）=n(1-1/p1) (1-1/p2)…… (1-1/pk)
其中p1， p2 ，pk是n的所有素数因子
Phi（n）：所有小于等于n的且与n互素的数的个数

(1)直接实现

  int oula(int n)
  {
      int rea=n;
      for(int i=2; i<=n; i++)
          if(n%i==0)//第一次找到的必为素因子
          {
              rea=rea-rea/i;
              while(n % i == 0) n/=i;//把该素因子全部约掉
         }
     return rea;
 }

(2)欧拉函数表

void lin()
{
    for(int i=2; i<1000001; i++)//欧拉公式，计算比j小且和j互斥的数的个数
    {
        if(d[i]==0)//如果是素数
        {
            for(int j=i; j<1000001; j+=i)//所有能整除i的j，在下面公式中先运算，最后可以运算完满足j的素数
            {
                if(d[j]==0)
                    d[j]=j;
                d[j]=d[j]/i*(i-1);
            }
        }
    }
    mp[1] = 1;
    for(int i=2; i<1000001; i++)
        mp[i]=mp[i-1]+d[i];
}

例题Farey Sequence Length

欧拉函数表的裸题
但是要注意mp[1]的初值是2,因为规定了0/1也是一个

#include<stdio.h>
#include<string.h>
#include<algorithm>
#include <iostream>
using namespace std;
int d[1000010];
long long mp[1000010];
void lin()
{
    for(int i=2; i<1000001; i++)//欧拉公式，计算比j小且和j互斥的数的个数
    {
        if(d[i]==0)//如果是素数
        {
            for(int j=i; j<1000001; j+=i)//所有能整除i的j，在下面公式中先运算，最后可以运算完满足j的素数
            {
                if(d[j]==0)
                    d[j]=j;
                d[j]=d[j]/i*(i-1);
            }
        }
    }
    mp[1] = 2;
    for(int i=2; i<1000001; i++)
        mp[i]=mp[i-1]+d[i];
}
int main()
{
    lin();
    int n;
    int t;
    cin>>t;
    while(t--)
    {
        int k ;
        scanf("%d%d",&k,&n);
        printf("%d %lld\n",k,mp[n]);
    }
    return 0;
}