欧拉函数
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];
}
欧拉函数表的裸题
但是要注意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;
}