Visible Lattice Points POJ - 3090(欧拉函数)

题目：

#include<iostream>
#include<string.h>
#include<stdio.h>
#include<algorithm>
#include<queue>
#include<stack>
#include<map>
#include<set>
#include<vector>
#include<cmath>
#define LL long long 
#define mem(a, b) memset(a, b, sizeof(a))
#define N 1010
#define MOD 1000000007
using namespace std;
int phi[N];

void euler(int n) {
	for(int i=2; i<=n; i++)
		phi[i]=i;
	for(int i=2; i<=n; i++) {
		if(phi[i]==i) {
			for(int j=i; j<=n; j+=i) {
				phi[j]=phi[j]/i*(i-1);
			}
		}
	}
}

int main() {
	int n, t, ans;
	int cas=1;
	scanf("%d", &t);
	while(t--) {
		scanf("%d", &n);
		euler(n);
		ans=0;
		for(int i=2; i<=n; i++) {
			ans+=phi[i];
		} 
		cout<<cas++<<" "<<n<<" "<<3+2*ans<<endl;
	} 
}

#include<iostream>
#include<string.h>
#include<stdio.h>
#include<algorithm>
#include<queue>
#include<stack>
#include<map>
#include<set>
#include<vector>
#include<cmath>
#define LL long long 
#define mem(a, b) memset(a, b, sizeof(a))
#define N 1010
#define MOD 1000000007
using namespace std;
int phi[N], v[N], prime[N], m;


void euler(int n) {
	mem(v, 0); //最小质因子
	m=0;//数量 
	for(int i=2; i<=n; i++) {
		if(v[i]==0) {
			v[i]=i, prime[++m]=i;
			phi[i]=i-1;
		}
		for(int j=1; j<=m; j++) {
			if(prime[j]>v[i] || prime[j]>n/i) break;//如果prime[j]不是最小质因子或者超过范围
			v[i*prime[j]]=prime[j];
			phi[i*prime[j]]=phi[i]*(i%prime[j] ? prime[j]-1 : prime[j]); 
		}
	}
}

int main() {
	int n, t, ans;
	int cas=1;
	scanf("%d", &t);
	while(t--) {
		scanf("%d", &n);
		euler(n);
		ans=0;
		for(int i=2; i<=n; i++) {
			ans+=phi[i];
		} 
		cout<<cas++<<" "<<n<<" "<<3+2*ans<<endl;
	} 
}