hdu2824(the Euler Function)

用欧拉筛求出欧拉函数。

http://blog.csdn.net/nk_test/article/details/46242311
欧拉函数：在数论中，对正整数n，欧拉函数是少于或等于n的数中与n互质的数的数目。

先给出一个结论：

设P是素数，

若p是x的约数，则E(x*p)=E(x)*p.

若p不是x的约数，则E(x*p)=E(x)E(p)=E(x)(p-1).

证明如下：

E(x)表示比x小的且与x互质的正整数的个数。
*若p是素数，E(p)=p-1。
*E(p^k)=p^k-p^(k-1)=(p-1)*P^(k-1)
证：令n=p^k,小于n的正整数数共有n-1即(p^k-1)个,其中与p不质的数共[p^(k-1)-1]个(分别为1*p,2*p,3*p…p(p^(k-1)-1))。
所以E(p^k)=(p^k-1)-(p^(k-1)-1)=p^k-p^(k-1).得证。
*若ab互质，则E(a*b)=E(a)*E(b)，欧拉函数是积性函数.
对任意数n都可以唯一分解成n=p1^a1*p2^a2*p3^a3…*pn^an(pi为素数).
则E(n)=E(p1^a1)E(p2^a2)*E(p3^a3)…*E(pn^an)
=(p1-1)p1^(a1-1)(p2-1)p2^(a2-1)…*(pn-1)*pn^(an-1)
=(p1^a1*p2^a2*p3^a3*…pn^an)[(p1-1)(p2-1)(p3-1)…(pn-1)]/(p1*p2*p3*…*pn)
=n*(1-1/p1)(1-1/p2)…*(1-1/pn)
* E(p^k) =(p-1)*p^(k-1)=(p-1)*p^(k-2)*p
E(p^(k-1))=(p-1)*p^(k-2)
->当k>1时，E(p^k)=E(p*p^(k-1))=E(p^(k-1))*p.
(当k=1时，E(p)=p-1.)
由上式: 设P是素数，
若p是x的约数，则E(x*p)=E(x)*p.
若p不是x的约数，则E(x*p)=E(x)E(p)=E(x)(p-1). 证明结束。

#include <stdio.h>
#include <string.h>
#include <iostream>
#include <algorithm>
using namespace std;
#define ms(X) memset(X,0,sizeof(X))
typedef long long ll;
bool isnP[4000000];
int phi[4000000],t=0,mprime[4000000];
void euler_ex(void)
{
    ms(isnP);
    phi[1]=1;
    for(int i=2;i<4000000;i++)
    {
        if(!isnP[i]) {mprime[t++]=i;phi[i]=i-1;}
        for(int j=0;i*mprime[j]<4000000;j++)
        {
            isnP[i*mprime[j]]=true;
            if(i%mprime[j]) phi[i*mprime[j]]=(mprime[j]-1)*phi[i];
            else {
                phi[i*mprime[j]]=mprime[j]*phi[i];
                break;
            }
        }
    }
}
int main(int argc, char const *argv[])
{
    int a,b;
    euler_ex();
    while(scanf("%d %d",&a,&b)!=EOF)
    {
        ll sum=0ll;
        for(int i=a;i<=b;i++)
            sum+=(ll)phi[i];
        printf("%lld\n",sum );
    }
    return 0;
}