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