【NKOJ3958/BZOJ2693】【莫比乌斯反演】jzptab

题目要求:

i=1nj=1mlcm(i,j)

看到lcm,首先想到转换为做过的gcd:

i=1nj=1mijgcd(i,j)

相同的套路,设:

d=gcd(i,j)

N=min(n,m)

f(n,m,d)=i=1nj=1mij[gcd(i,j)=d]

=d=1Nf(n,m,d)d

f(d)i,jdd=1Ndf(nd,md,1)

F(n,m,d)=i=1nj=1mij[d|gcd(i,j)

F(n,m,d)=d|eNf(n,m,e)

f(n,m,d)=d|eNμ(ed)F(n,m,e)

sum(n,m)=i=1nj=1mij=n(n+1)2m(m+1)2

F(n,m,d)=d2sum(nd,md)

f(n,m,1)=i=1Nμ(i)i2sum(ni,mi)

Ans=d=1Ndi=1Ndμ(i)i2sum(ndi,mdi)

到目前为止,我们可以利用分块优化,在第一个和式利用Nd进行分块,在第二个和式利用ndi进行分块,在O(nn)=O(n)的复杂度下求出答案

然而题目还有T组数据的要求OTZ,卡O(n)

于是回忆起董神的教诲,尝试改变枚举顺序,把比较烦的di直接转成枚举量

Ans=d=1min(n,m)dsum(nd,md)i|ddμ(i)i

前面求sum()肯定是logn复杂度,于是我们需要在后面那一坨上面进行优化

g(d)=di|dμ(i)i

显然μ(i)i为积性函数,根据莫比乌斯反演性质,g(d)也为积性函数

d拆为pk11pk22pk33....pkii (pi为质因子)

对于每个pkii,g(pkii)=1pi(显然)

g(d)=Πdi=1g(pkii)=Πki=1(1pi)

很像phi,放在线性筛中处理

总复杂度O(Tn)

代码:

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N=10000005,Inf=1e9;
const ll LInf=1e15,Mod=100000009;
inline ll _R(){
    ll d=0;bool b=1;char t=getchar();
    while(t<'0'||t>'9'){if(t=='-')b=0;t=getchar();}
    for(;t>='0'&&t<='9';t=getchar())d=(d<<3)+(d<<1)+t-'0';
    return b?d:-d;
}

ll T,n,m,t,Ans;
ll Cnt,Prime[N],Vis[N],Sum[N];
void Euler(){
     Sum[1]=1;
     for(ll i=2;i<=N-5;i++){
        if(!Vis[i]){
            Sum[i]=(i-i*i%Mod+Mod)%Mod;
            Prime[++Cnt]=i;
        }
        for(ll j=1;Prime[j]*i<=N-5&&j<=Cnt;j++){
            ll k=Prime[j]*i;
            Vis[k]=1;
            if(i%Prime[j])Sum[k]=Sum[Prime[j]]*Sum[i]%Mod;
            else {
                Sum[k]=Sum[i]*Prime[j]%Mod;
                break;
            }
        }
     }
     for(int i=1;i<=N-5;i++)Sum[i]=(Sum[i-1]+Sum[i])%Mod;
}
ll Cal(ll x,ll y){
    return (x*(x+1)/2%Mod)*(y*(y+1)/2%Mod)%Mod;
}
int main(){
    T=_R();
    Euler();
    while(T--){
        n=_R(),m=_R();
        Ans=0;
        for(ll i=1;i<=min(n,m);i++){
            t=min(n/(n/i),m/(m/i));
            Ans=(Ans+((Sum[t]-Sum[i-1]+Mod)%Mod)*Cal(n/i,m/i)%Mod)%Mod;
            i=t;
        }
        printf("%lld\n",Ans);
    }
}
发布了5 篇原创文章 · 获赞 3 · 访问量 1499
展开阅读全文

没有更多推荐了,返回首页

©️2019 CSDN 皮肤主题: 大白 设计师: CSDN官方博客

分享到微信朋友圈

×

扫一扫,手机浏览