poj 2407

最新推荐文章于 2020-06-29 17:42:19 发布

傻笨

最新推荐文章于 2020-06-29 17:42:19 发布

阅读量488

点赞数

分类专栏：欧拉函数 poj

本文链接：https://blog.csdn.net/u010195743/article/details/8985551

版权

poj 同时被 2 个专栏收录

22 篇文章 0 订阅

订阅专栏

欧拉函数

5 篇文章 0 订阅

订阅专栏

这是一道裸的欧拉函数题目：

#include <iostream>
#include <cstdio>
#include <cmath>
using namespace std;

int euler(int x)
{

int unsigned i, res=x;
for (i = 2; i < (int)sqrt(x * 1.0) + 1; i++)
if(x%i==0)
{
res = res / i * (i - 1);
while (x % i == 0) x /= i; // 保证i一定是素数
}
if (x > 1) res = res / x * (x - 1);
return res;
}

int main()
{
int n,m,i;
while(scanf("%d",&n,&m)!=EOF)
{
int sum=0;
for(i=n; i<=m; i++)
{
sum+=euler(i);
}
printf("%d\n",sum);
}
return 0;
}

下面是两种求欧拉函数的不同编程方法：

/*==================================================*\
|递推求欧拉函数phi(i)
\*==================================================*/
for (i = 1; i <= maxn; i++) phi[i] = i;
for (i = 2; i <= maxn; i += 2) phi[i] /= 2;
for (i = 3; i <= maxn; i += 2) if(phi[i] == i) {
for (j = i; j <= maxn; j += i)
phi[j] = phi[j] / i * (i - 1);

/*==================================================*\
|单独求欧拉函数phi(x)
\*==================================================*/
unsigned euler(unsigned x)
{// 就是公式
unsigned i, res=x;
for (i = 2; i < (int)sqrt(x * 1.0) + 1; i++)
if(x%i==0) {
res = res / i * (i - 1);
while (x % i == 0) x /= i; // 保证i一定是素数
}
if (x > 1) res = res / x * (x - 1);
return res;
}