題目不多說了，看了一會用了最弱智的暴力求解果然TLE了，最後上網查了一下正確解法是尤拉函式，下面就簡單的總結一下尤拉函式。

尤拉函式：在數論，對正整數n，尤拉函式是少於或等於n的數中與n互質的數的數目。此函式以其首名研究者尤拉命名，它又稱為Euler's totient function、φ函式、尤拉商數等。

它的一般計算式是：phi(n) = n*(1-1/p1)*(1-1/p2)*...*(1-1/pk)，其中p1...pk為n的所有質因子。

關於尤拉函式有如下幾點性質：

1、phi(1) = 1

2、若n是質數，那麼phi(n) = n-1

3、若n是質數x的k次冪，phi(n) = (x-1)*x^(k-1)

4、若m，n互質，那麼phi(m*n) = phi(m)*phi(n)

5、若n是奇數，那麼phi(2*n) = phi(n)

6、若x，y是質數，且n = x*y，那麼phi(n) = (x-1)*(y-1)

7、小於n且與n互質的數的和為：n/2 * phi(n)

那麼給出兩種求解尤拉函式的方法：

//直接求解尤拉函式
int euler(int n){ //返回euler(n) 
     int res=n,a=n;
     for(int i=2;i*i<=a;i++){
         if(a%i==0){
             res=res/i*(i-1);//先進行除法是為了防止中間資料的溢位 
             while(a%i==0) a/=i;
         }
     }
     if(a>1) res=res/a*(a-1);
     return res;
}

//篩選法打尤拉函式表 
#define Max 1000001
int euler[Max];
void Init(){ 
     euler[1]=1;
     for(int i=2;i<Max;i++)
       euler[i]=i;
     for(int i=2;i<Max;i++)
        if(euler[i]==i)
           for(int j=i;j<Max;j+=i)
              euler[j]=euler[j]/i*(i-1);//先進行除法是為了防止中間資料的溢位 
}