HDU 2824 The Euler function

大概:求a~b段欧拉值和,要用到筛法欧拉函数。

筛选:1)如果n为素数,φ(n)=n-1;

           2)如果m,n互质,φ(n*m)=φ(n)*φ(m);

   

#include<iostream>
#include<cstdio>
#include<ctime>
using namespace std;
typedef  __int64 INT;
#define N 3000001
int prime[216818],phi[N];
bool a[N];
void is_prime()
{
	int i,j,k;
	k=0;
	for(i=2;i<N;i++)
	{
	   if(!a[i])
	   {
	      prime[k++]=i;
		  phi[i]=i-1;
	   }
	   for(j=0;j<k&&i*prime[j]<N;j++)
	   {
	     a[prime[j]*i]=1;
		 if(i%prime[j])
			 phi[i*prime[j]]=phi[i]*(prime[j]-1);
		 else
		 {
			 phi[i*prime[j]]=phi[i]*prime[j];
			 break;
		 }
	   
	   }
	
	}
}
int main()
{
	//freopen("c_in.txt","r",stdin);
	//freopen("myout5.txt","w",stdout);
    int a,b,i;
	INT sum;
	is_prime();
	while(cin>>a>>b)
	{
		sum=0;
		for(i=a;i<=b;i++)
			sum+=phi[i];
		printf("%I64d\n",sum);
	
	}
	//printf("%lfMS\n",(double)clock()/CLOCKS_PER_SEC);
	return 0;
}

       

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值