51Nod 1220 约数之和&&PE439

$∂(i)$表示i的约数和
$∑_{i=1}^{n}∑_{j=1}^{n}∂(ij)$
$=∑_{i=1}^{n}∑_{j=1}^{n}∑_{w|i}∑_{v/j}wv*∑_{d|w,d|(i/v)}u(d)$
我们尝试把d提前 会发现有
$=∑_{d=1}^{n}d*u(d)(∑_{i=1}^{(n/d)}∂(i))^2$
后面$(∑_{i=1}^{(n/d)}∂(i))^2$可以看作$f(n/d)$
对于$∑_{i=1}^{(n/d)}∂(i)$我们可以小数据预处理&&大数据$O(n^{0.5})$暴力
那么我们有
$=∑_{d=1}^{n}d*u(d)*f(n/d)$
现在我们只需要求出$∑_{d=1}^{n}d*u(d)$即可
设$∑_{d=1}^{n}d*u(d)=S(n)$我们给前式卷上一个id
可得
$S(n)=1-∑_{d=2}^{n}d*S(n/d)$
好了 完美解决
如果你要交PE的话。。
祝你好运
由于PE模数原因。。我和另一个同学花了一个晚上找哪里爆了

#include<cstdio>
#include<iostream>
#include<cstring>
#include<cstdlib>
#include<algorithm>
#include<map>
using namespace std;
#define ll long long
char c;
inline void read(ll &a)
{
    a=0;do c=getchar();while(c<'0'||c>'9');
    while(c<='9'&&c>='0')a=(a<<3)+(a<<1)+c-'0',c=getchar();
}

const
    ll Mod=(int)1e9,MAXN=(int)3e7+1;
bool check[MAXN];
int mu[MAXN],prime[MAXN],tot,Min[MAXN],siz[MAXN],s[MAXN];

inline int add(int x,int y)
{
    int res=x+y;
    if(res>=Mod)
        res-=Mod;
    if(res<0)res+=Mod;
    return res;
}
inline int mul(int x,int y){return x*1ll*y%Mod;}
inline 
int A(ll x)
{
    if(x&1)
        return (x%Mod)*((x+1>>1)%Mod)%Mod;
    return ((x>>1)%Mod)*((x+1)%Mod)%Mod;
}
inline 
int cn(ll x,ll y){return A(y)-A(x-1);}
map<ll ,int >M;
int Mu(ll x)
{
    if(x<MAXN)return mu[x];
    if(M[x])return M[x];
    ll i=2,j;
    int res=0;
    while(i<=x)
    {
        j=x/(x/i);
        res=add(res,mul(cn(i,j),Mu(x/i))%Mod);
        i=j+1;
    }
    return M[x]=add(1,Mod-res);
}

map<ll,int >C;
int Calc(ll x)
{
    if(x<MAXN)return s[x];
    if(C[x])return C[x];
    int res=0;
    ll i=1,j;
    while(i<=x)
    {
        j=x/(x/i);
        res=add(res,mul(x/i%Mod,cn(i,j)));
        i=j+1;
    }
    return C[x]=mul(res,res);
}
int S(ll x)
{
    ll i=1,j;
    int res=0;
    while(i<=x)
    {
        j=x/(x/i);
        int tp=add(Mu(j),-Mu(i-1));
        res=add(res,mul(tp,Calc(x/i)));
        i=j+1;
    }
    return res;
}

void out(int x)
{
    if(x<0)putchar('-'),x=-x;
    if(x>9)out(x/10);
    putchar('0'+x%10);
}

int main()
{
    s[1]=mu[1]=siz[1]=Min[1]=1;
    for(ll i=2;i<MAXN;i++)
    {
        if(!check[i])prime[++tot]=i,s[i]=1+i,Min[i]=i,siz[i]=i,mu[i]=-1;
        for(ll j=1,k;j<=tot&&(k=prime[j]*i)<MAXN;j++)
        {
            check[k]=true;
            if(i%prime[j]==0)
            {
                Min[k]=Min[i];
                siz[k]=siz[i]*prime[j];
                mu[k]=0;
                s[k]=s[i/siz[i]]*(s[i]/s[i/siz[i]]+siz[k]);
                break;
            }
            Min[k]=prime[j];
            siz[k]=prime[j];
            mu[k]=-mu[i];
            s[k]=s[i]*(prime[j]+1);
        }
        mu[i]=mul(mu[i],i);
        mu[i]=add(mu[i],mu[i-1]);
    }
    for(ll i=1;i<MAXN;i++)
        s[i]=add(s[i],s[i-1]);
    for(ll i=1;i<MAXN;i++)
        s[i]=mul(s[i],s[i]);
    out(S((ll)1e11));
    return 0;
}