原根

http://zh.wikipedia.org/wiki/%E5%8E%9F%E6%A0%B9

欧拉函数：

void phi_table()

{

    int i,j;

    memset(phi,0,sizeof(phi));

    phi[1]=1;

    for(i=2;i<=n;i++)

    {

         if(!phi[i])

             for(j=i;j<=n;j+=i)

                 {

                     if(!phi[j])  phi[j]=j;

                     phi[j]=phi[j]/i*(i-1);

                 }

    }

}