题目大意:
用 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;
}