poj2478 Farey Sequence (欧拉函数)

Farey Sequence

题意:给定一个数n,求在[1,n]这个范围内两两互质的数的个数。(转化为给定一个数n,比n小且与n互质的数的个数)

知识点:

欧拉函数:poj2478 Farey Sequence (欧拉函数)

普通求法:

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;
}
}
上一篇:【Linux】Linux 目录结构


下一篇:转载:Candy? 在线性时间内求出素数与欧拉函数