快速求解歐拉函數

原始碼：

void get(){  for(int i=2;i<1001000;i++){      num[i]=1;      euler[i]=i;  }  for(int i=2;i<1001000;i++)   {       if(num[i])       for(int j=i;j<1001000;j+=i)       {           num[j]=0;           euler[j]=euler[j]/i*(i-1);       }   }}

剛開始的時候，第一個迴圈的判斷條件是i*i<1001000，後來發現這樣的話，不能將j的所有質因子都找到，他這樣只枚舉了sqrt(1001000)的素數。

搞得我調了一段時間。