[UVA 10820]Send a Table[欧拉函数][nloglog(n)]

题目链接： [UVA 10820]Send a Table[欧拉函数][nloglog(n)]

题意分析：

友人A想要打表过题，每个表内数据都是以数对的形式(x,y)出现，但是呢，表太大了，OJ不让上交。机智的他发现了一个规律：数对(x*k, y*k)可以由数对(x,y)推出来，所以他只用打(x,y)上交就行了。比如(1,2)被打入表，那么(2,4) (3,6)....这些数据就不用打入表了。现在给出数字N，求给这些数字两两打表，要打多少个数据？

解题思路：

根据题意，最后留下来的那些数对都有一个特点gcd(x,y) == 1。问题转化为了，以某个数作为数对的第一个数，那么与它互成素数的小于它的整数有多少个？把这些数加起来乘以二（因为两两可以互换），最后加上数对(1, 1)就是最终答案。那么怎么快速求一个整数小于他而且与它互为素数的数的个数呢？这里我们引入欧拉函数求解，关于欧拉函数，详见这篇博文：http://blog.csdn.net/catglory/article/details/47674961

个人感受：

一直在补其它题，没去解决欧拉函数这个遗留的问题，今天个人赛开了这个陈题，然而我并不会（之前压根没理解，也就懒得记忆了）。今晚详细看了下，算了理解了。

具体代码如下：

#include<iostream>
using namespace std;
const int INF = 0x7f7f7f7f, MAXN = 50001;
int phi[MAXN], n;

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

int main()
{
    init();
    while (cin >> n && n)
    {
        cout << phi[n] * 2 + 1 << '\n';
    }
    return 0;
}