HDOJ-2824The Euler function(筛选法求欧拉函数)

最新推荐文章于 2021-04-23 20:33:30 发布

原创最新推荐文章于 2021-04-23 20:33:30 发布 · 722 阅读

0 ·

本内容遵循CC 4.0 BY-SA版权协议

标签

#acm #数论

收录于

ACM 专栏收录该内容

76 篇文章

订阅专栏

本文介绍了一种利用筛选法求解欧拉函数的方法，详细解释了其背后的数学原理，并通过C++代码实现了一个高效的计算过程。该方法不仅适用于初学者理解欧拉函数的概念，也为实际应用提供了优化解决方案。

欧拉函数的定义 f(n)=[1,n]中与n互质的数的个数

筛选法求欧拉函数原理，利用了欧拉的一下性质：

(1) if (b%a==0 && (b/a)%a==0) then f(b)=f(b/a)*a;
(2) if (b%a==0 && (b/a)%a!=0) then f(b)=f(b/a)*(a-1);

主要思想和筛选法求素数是一样的，枚举约数，然后对倍数进行计算。

#include <cstdio>
#include <cstdlib>
#include <cstring>
using namespace std;
const int N=3E6;

long long euler[N+5];

void get_euler()
{
	int i,j;
	memset(euler,0,sizeof(euler));
	euler[1]=1;
	for(i=2;i<=N;i++)
		if(!euler[i])
			for(j=i;j<=N;j+=i)
			{
				if(!euler[j]) euler[j]=j;
				euler[j]=euler[j]*(i-1)/i;
			}
}

int main()
{
	int i,a,b;
	get_euler();
	for(i=1;i<=N;i++)
		euler[i]+=euler[i-1];
	while(scanf("%d%d",&a,&b)==2)
		printf("%I64d\n",euler[b]-euler[a-1]);
	return 0;
}