hdu 2824欧拉函数(晒法模板)

点击打开链接

#include <iostream>
#include <algorithm>
#include <cmath>
using namespace std;
typedef long long ll;
const int N=3e6+20;
ll phi[N+20];
void Euler()
{
	//Euler(n)=n*(1-1/p1)*(1-1/p2)...(1-1/pn)
	//通分=n*((p1-1)/p1)*((p2-1)/p2)...((pn-1)/pn)
	for(int i=2;i<=N;i++)
	phi[i]=i;
	
	phi[1]=1;
	for(int i=2;i<=N;i++)
	{
		if(phi[i]==i)//i是素因子
		{
			phi[i]=i-1;//
			for(int j=i+i;j<=N;j+=i)//j含有素因子i 
			{
				phi[j]=phi[j]/i*(i-1); 
			}	
		} 
 	}

}
int main()
{
	ll a,b;
	Euler();
	while(cin>>a>>b)
	{
		ll ans=0;
		for(int i=a;i<=b;i++)
		ans+=phi[i];

		cout<<ans<<endl;
	}
	return 0;
}