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;
}

       

你可能感兴趣的:(HDU 2824 The Euler function)