POJ 3090 Visible Lattice Points 【欧拉函数】

题目: http://poj.org/problem?id=3090

题意: 如题, 实质是求给出一个n, 求出ψ(1) + ψ(2) + ψ(3) + ... ψ(n).

　　　　ψ(i) 是欧拉函数, 表示在 1..i内是i的质因子的个数.

　　　　要求ψ(i), 由 i = p1^(a1)*p2^(a2)*p3(a3)...*pi^(ai)

　　　　可推出 ψ(i) = i(1 - 1/p1)(1 - 1/p2)...(1 - 1/pi);

　　　　pi表示i的素因子.

代码:

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cstring>

using namespace std;
const int MAX = 1000 + 1;
int D[MAX];
void phi_table(int n)
{
int i, j;
memset(D, 0, sizeof(D));
D[1] = 1;
for(i = 2; i <= n; ++i)
{
if(!D[i])
{
for(j = i; j <= n; j += i)
{
if(!D[j])
D[j] = j;
D[j] =D[j] / i * (i - 1);
}
}
}
D[1] = 3;
//printf("2 = %d\n", D[5]);
for(i = 2; i <= n; ++i)
{
D[i] = D[i - 1] + D[i] * 2;
}
}
int main()
{
phi_table(1000);
int N, n, i;
scanf("%d", &N);
for(i = 1; i <= N; ++i)
{
scanf("%d", &n);
printf("%d %d %d\n", i, n, D
);
}
return 0;
}