bzoj2226[Spoj 5971] LCMSum

题目大意：给定T（10^5）组数据，每次给出n（10^6），求lcm（i,n）{1<=i<=n}的和

一眼看上去就是一道数学题……蒟蒻数学太差思路略微诡异……我想着能不能nlogn直接把所有答案算出来，然后yy了一下……

设F[n]为所求的值，f[n]为小于n的质数的和，我们考虑到lcm(a,b)=ab/gcd(a,b)，那么我们可以去枚举gcd，对于一个n的因数d，在1到n中和n最大公因数为d的数对于F[n]的贡献就显而易见了，即F[n]+=n*d*f[n/d]，这样可以nlogn求出

对于f[n]求法也非常的简单，可以看出f[n]=sigma(i*[gcd(i,n)=1])，再把[gcd(i,n)=1]写成莫比乌斯函数的形式，套路一样的重新组合啊什么的，最终就是可以变成这样的形式：f[n]+=u[k]*(k*((n/k)+1)(n/k)/2){k|n}，也是可以在nlogn时间求出。

听说有根号大法……不过不会……数学太弱QAQ

为了避免常数被卡特意用C++写的

#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
#define ll long long
using namespace std;
const int N=1000086;
ll f[N],F[N],pri[N],vis[N],mu[N],tot;
int read()
{
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
    return x*f;
}
void getmu()
{
    mu[1]=1;tot=0;
    for(int i=2;i<=N;i++)
    {
        if(!vis[i]){pri[++tot]=i;mu[i]=-1;}
        for(int j=1;j<=tot&&pri[j]*i<=N;j++)
        {
            vis[i*pri[j]]=1;
            if(i%pri[j]==0){mu[pri[j]*i]=0;break;}
            else mu[pri[j]*i]=-mu[i];
        }
    }
}
void getf()
{

    for(int i=1;i<=N;i++)
        for(int j=1;j*i<=N;j++)
            f[i*j]+=mu[i]*((ll)i*(j+1)*j/2);
}
void getF()
{
    for(int i=1;i<=N;i++)
        for(int j=1;j*i<=N;j++)
            F[i*j]+=f[j]*i*j;
}
int main()
{
    getmu();getf();getF();
    int T=read();
    while(T--)
    {
        int n=read();printf("%lld
",F[n]);
    }
    return 0;
}