>>925
よく考えたらもっと直接
gcd(i,n)=Σ[j=1,n] (1 - 2{ij/n})
と計算できることがわかりました！