Farey Sequence
题意:给定一个数n,求在[1,n]这个范围内两两互质的数的个数。(转化为给定一个数n,比n小且与n互质的数的个数)
知识点:
欧拉函数:
普通求法:
int Euler(int n)
{
int ans=n;
for(int i=0;i<cnt&&prime[i]<=n;i++)
{
if(n%prime[i]==0)
{
ans=ans-ans/prime[i];
while(n%prime[i]==0)
n/=prime[i];
}
}
if(n==1)
return ans;
if(n>1)
return ans-ans/n; }
筛选法:(基于素数筛,跟着代码模拟一遍就懂)
void Init()
{
for(int i=1;i<N;i++)
euler[i]=i;
for(int i=2;i<N;i++) //i=1时,euler[1] 不变
if(euler[i]==i)
for(int j=i;j<N;j+=i)
euler[j]=euler[j]/i*(i-1);
}
题解:筛选法求欧拉函数
#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
const int N=1e6+5; int euler[N]; void Init()
{
euler[1]=1;
for(int i=2;i<N;i++)
euler[i]=i;
for(int i=2;i<N;i++)
if(euler[i]==i)
for(int j=i;j<N;j+=i)
euler[j]=euler[j]/i*(i-1);
} int main()
{
Init();
int n;
while(cin>>n&&n)
{
long long ans=0;
for(int i=2;i<=n;i++)
ans+=euler[i];
cout<<ans<<endl;
}
}