SPOJ 5971 LCM Sum 欧拉函数 (或 莫比乌斯反演?)

题目大意：

用 LCM(i, n) 代表 i 和 n 的最小公倍数， 求sigma(LCM(i, n)， 1 <= i <= n)的值

数据范围：

测试数据组数 1 <= T <= 300000 

1 <= n <= 1000000

大致思路：

首先这道题刚开始听学长说的思路是莫比乌斯反演...虽然我做到最后还是没用上...不过因为这些还是取了解了很多和莫比乌斯反演相关的东西， 这样子的话以后碰到必须用莫比乌斯反演的题应该能接受的快一些...不过数论真是高大上 T^T

首先用gcd( a, b)表示a和b的最大公约数有：

算了公式多还是都写在图里吧...

这里欧拉函数用的O( n*loglogn )的复杂度计算,之后存储答案的数组计算用的也是O(n*loglogn)

递归计算即可

代码如下：

Result  :  Accepted     Memory  :  14 M     Time  :  4.54 s (险过= =)

/*
 * Author: Gatevin
 * Created Time:  2014/8/4 18:37:22
 * File Name: test.cpp
 */
#include<iostream>
#include<sstream>
#include<fstream>
#include<vector>
#include<list>
#include<deque>
#include<queue>
#include<stack>
#include<map>
#include<set>
#include<bitset>
#include<algorithm>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cctype>
#include<cmath>
#include<ctime>
#include<iomanip>
using namespace std;
const double eps(1e-8);
typedef long long lint;
typedef unsigned long long ulint;
#define maxn 1000010

int phi[maxn];
lint f[maxn];
int t;
int n;

void initial()
{
    memset(phi, 0, sizeof(phi));
    phi[1] = 1;
    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);
            }
        }
    }
    memset(f, 0, sizeof(f));
    f[1] = 1;
    for(int i = 2; i < maxn; i++)
    {
        for(int j = i; j < maxn; j += i)
        {
            f[j] += i*(lint)phi[i] / 2 * j;
        }
    }
    for(int i = 2; i < maxn; i++)
    {
        f[i] += i;
    }
}

int main()
{
    initial();
    scanf("%d",&t);
    for(int cas = 1; cas <= t; cas++)
    {
        scanf("%d",&n);
        printf("%lld\n", f[n]);
    }
    return 0;
}