2016夏季练习——数论

来源：POJ3090

直接就是求欧拉函数

直接上代码：

#include <iostream>#include <cstdio>#include <cstring>using namespace std;const int MAXN = 1010;long long p[MAXN+10];void ini(){memset(p,0,sizeof(p));p[1]=1;for(int i=2;i<MAXN;i++){if(!p[i]){for(int j=i;j<MAXN;j+=i){if(!p[j])p[j]=j;p[j]=p[j]/i*(i-1);}}}}int main(){ini();int T,nc=1;scanf("%d",&T);while(T--){int n;scanf("%d",&n);int sum=3;for(int i=2;i<=n;i++)sum+=2*p[i];cout<<nc++<<" "<<n<<" "<<sum<<endl;}return 0;}