51nod1227 平均最小公倍数

$
Ans(l,r)=ans(r)-ans(l-1)
\\
ans(n)=\sum\limits_{i=1}^n
\sum\limits_{j=1}^i
\frac{j}{gcd(i,j)}
\\=
\sum\limits_{g=1}^n
\sum\limits_{i=1}^{n/g}
\sum\limits_{j=1}^n
j\cdot [gcd(i,j)=1]
\\=
\sum\limits_{g=1}^n
\sum\limits_{i=1}^{n/g}
\sum\limits_{j=1}^n
j\sum\limits_{d|i\wedge d|j}\mu(d)
\\=
\sum\limits_{g=1}^n
\sum\limits_{d=1}^{n/g}
\mu(d)d\sum\limits_{i=1}^{n/gd}
\sum\limits_{j=1}^n
j
\\=
\sum\limits_{t=1}^n
\sum\limits_{d|t}\mu(d)d
\sum\limits_{i=1}^{n/t}
\sum\limits_{j=1}^n
j
\\=
\sum\limits_{t=1}^n
a(t)S(n/t)
,\;
a(n)=\sum\limits_{d|n}\mu(d)d
,\;
S(n)=n(n+1)(n+2)/6=n(n+1)(2n+1)/12+n(n+1)/4
\\
A(n)=\sum\limits_{i=1}^na(i)
\\
A(n)=n-\sum\limits_{i=2}^nA(n/i)i
$

时间复杂度$O(n^{2/3})$

#include<bits/stdc++.h>
typedef unsigned long long i64;
const int P=1e9+,I2=(P+)/,I6=(P+)/;
int B;
const int M=3e6+;
int ps[M/],pp=;
bool np[M];
int va[M],vA[],n0;
int S3(int n){
    return i64(n)*(n+)%P*(n+)%P;
}
int A(int n){
    if(n<=B)return va[n];
    int&w=vA[n0/n];
    if(w)return w;
    i64 s=;
    for(int l=,r,c;l<n;l=r){
        r=n/(c=n/(l+));
        if((s+=i64(r+l+)*(r-l)%P*A(c))>i64(1.5e19))s%=P;
    }
    s=s%P*I2%P;
    return w=(n-s+P)%P;
}
int F(int n){
    i64 s=;
    for(int l=,r,c,s0=,s1;l<n;l=r){
        r=n/(c=n/(l+));
        s1=A(r);
        if((s+=i64(s1-s0+P)*S3(c))>i64(1.5e19))s%=P;
        s0=s1;
    }
    return s%P*I6%P;
}
void pre(){
    va[]=;
    for(int i=;i<=B;++i){
        if(!np[i]){
            ps[pp++]=i;
            va[i]=P+-i;
        }
        for(int j=,k;j<pp&&(k=i*ps[j])<=B;++j){
            np[k]=;
            if(i%ps[j]){
                va[k]=va[i]*i64(P+-ps[j])%P;
            }else{
                va[k]=va[i];
                break;
            }
        }
    }
    for(int i=;i<=B;++i)if((va[i]+=va[i-])>=P)va[i]-=P;
}
int cal(int n){
    n0=n;
    memset(vA,,sizeof(vA));
    return F(n);
}
int main(){
    int l,r;
    scanf("%d%d",&l,&r);
    B=pow(r,./.)*2.5;
    pre();
    printf("%d\n",(cal(r)-cal(l-)+P)%P);
    return ;
}