uvaoj 10820 Send a Table 欧拉函数打表

uvaoj 10820 Send a Table 欧拉函数打表

这个题目就是让求小于等于n的数中，互质二元组(x,y)的个数。因为(x,y)和(y,x)不同，我们假定只考虑x<y的情况，一种特殊的情况就是x和y都是1。那么问题就变成了小于等于n，满足x<y的二元组的个数，这个数目其实就是1到n的欧拉函数的和。

下面先介绍一下欧拉函数：

在数论，对正整数n，欧拉函数是少于或等于n的数中与n互质的数的数目。

此函数以其首名研究者欧拉命名，它又称为Euler's totient function、φ函数、欧拉商数等。

欧拉函数的公式为：

phi(n) = n*(1-1/p1)*(1-1/p2)*...*(1-1/pk)。

其中，p1到pk是n的素因子，即n=p1^a1*p2^a2*...*pk^ak。

欧拉函数的性质：

对于质数p，φ(p) = p - 1。注意φ(1)=1.
欧拉定理：对于互质的正整数a和n，有aφ(n) ≡ 1 mod n。
欧拉函数是积性函数——若m,n互质，φ(mn)=φ(m)φ(n)。
                                 若n是质数p的k次幂，φ(n)=p^k-p^(k-1)=(p-1)p^(k-1)，因为除了p的倍数外，其他数都跟n互质。
特殊性质：当n为奇数时，φ(2n)=φ(n)。

求欧拉函数的方法是，每次找到一个素因子，把它除干净，即可保证找到的因子都是素数，不可能有两个质因子都大于根号n。

方法如下：

int euler_phi(int n)
{
	int m = (int)sqrt(n + 0.5);
	int ans = n;
	for (int i = 2; i <= m; ++i) {
		if (n % i == 0) {
			ans = ans / i * (i - 1);
			while (n % i == 0) {
				n /= i;
			}
		}
	}
	if (n > 1) {
		ans = ans / n * (n - 1);
	}
}

如果要求多个1到n的所有数的欧拉函数，有更快的方法，类似筛选素数，就是先筛选素数，则这个素数就是一个素因子。复杂度为O(nloglogn)。

题目代码如下：

/*************************************************************************
	> File Name: 10820.cpp
	> Author: gwq
	> Mail: gwq5210@qq.com 
	> Created Time: 2015年01月09日 星期五 12时28分23秒
 ************************************************************************/

#include <cmath>
#include <ctime>
#include <cctype>
#include <climits>
#include <cstdio>
#include <cstdlib>
#include <cstring>

#include <map>
#include <set>
#include <queue>
#include <stack>
#include <string>
#include <vector>
#include <sstream>
#include <iostream>
#include <algorithm>

#define INF (INT_MAX / 10)
#define clr(arr, val) memset(arr, val, sizeof(arr))
#define pb push_back
#define sz(a) ((int)(a).size())

using namespace std;
typedef set<int> si;
typedef vector<int> vi;
typedef map<int, int> mii;
typedef long long ll;

const double esp = 1e-5;

#define N 50010

int phi[N], sum[N];

void phi_table(void)
{
	clr(phi, 0);
	phi[1] = 1;
	for (int i = 2; i < N; ++i) {
		if (!phi[i]) {
			for (int j = i; j < N; j += i) {
				if (!phi[j]) {
					phi[j] = j;
				}
				phi[j] = phi[j] / i * (i - 1);
			}
		}
	}
}

int main(int argc, char *argv[])
{
	phi_table();
	clr(sum, 0);
	for (int i = 1; i < N; ++i) {
		sum[i] += sum[i - 1] + phi[i];
	}
	int n;
	while (scanf("%d", &n) != EOF) {
		if (n == 0) {
			break;
		}
		printf("%d\n", 2 * sum[n] - 1);
	}

	return 0;
}

参考：

1）http://baike.baidu.com/link?url=kD_VMAkBJ0sLzngBoiUK8geaSdDFNeQ2GCMRC8f4C5Gq0Wyml4jRJtrGq34DfcAo3Xs9GHqboNyje13jDahgNK

2）http://blog.csdn.net/leolin_/article/details/6642096