【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);
    }
}
  • 1
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
1、资源项目源码均已通过严格测试验证,保证能够正常运行; 2、项目问题、技术讨论,可以给博主私信或留言,博主看到后会第一时间与您进行沟通; 3、本项目比较适合计算机领域相关的毕业设计课题、课程作业等使用,尤其对于人工智能、计算机科学与技术等相关专业,更为适合; 4、下载使用后,可先查看README.md文件(如有),本项目仅用作交流学习参考,请切勿用于商业用途。1、资源项目源码均已通过严格测试验证,保证能够正常运行; 2、项目问题、技术讨论,可以给博主私信或留言,博主看到后会第一时间与您进行沟通; 3、本项目比较适合计算机领域相关的毕业设计课题、课程作业等使用,尤其对于人工智能、计算机科学与技术等相关专业,更为适合; 4、下载使用后,可先查看README.md文件(如有),本项目仅用作交流学习参考,请切勿用于商业用途。1、资源项目源码均已通过严格测试验证,保证能够正常运行; 2、项目问题、技术讨论,可以给博主私信或留言,博主看到后会第一时间与您进行沟通; 3、本项目比较适合计算机领域相关的毕业设计课题、课程作业等使用,尤其对于人工智能、计算机科学与技术等相关专业,更为适合; 4、下载使用后,可先查看README.md文件(如有),本项目仅用作交流学习参考,请切勿用于商业用途。1、资源项目源码均已通过严格测试验证,保证能够正常运行; 2、项目问题、技术讨论,可以给博主私信或留言,博主看到后会第一时间与您进行沟通; 3、本项目比较适合计算机领域相关的毕业设计课题、课程作业等使用,尤其对于人工智能、计算机科学与技术等相关专业,更为适合; 4、下载使用后,可先查看README.md文件(如有),本项目仅用作交流学习参考,请切勿用于商业用途。1、资源项目源码均已通过严格测试验证,保证能够正常运行; 2、项目问题、技术讨论,可以给博主私信或留言,博主看到后会第一时间与您进行沟通; 3、本项目比较适合计算机领域相关的毕业设计课题、课程作业等使用,尤其对于人工智能、计算机科学与技术等相关专业,更为适合; 4、下载使用后,可先查看README.md文件(如有),本项目仅用作交流学习参考,请切勿用于商业用途。1、资源项目源码均已通过严格测试验证,保证能够正常运行; 2、项目问题、技术讨论,可以给博主私信或留言,博主看到后会第一时间与您进行沟通; 3、本项目比较适合计算机领域相关的毕业设计课题、课程作业等使用,尤其对于人工智能、计算机科学与技术等相关专业,更为适合; 4、下载使用后,可先查看README.md文件(如有),本项目仅用作交流学习参考,请切勿用于商业用途。1、资源项目源码均已通过严格测试验证,保证能够正常运行; 2、项目问题、技术讨论,可以给博主私信或留言,博主看到后会第一时间与您进行沟通; 3、本项目比较适合计算机领域相关的毕业设计课题、课程作业等使用,尤其对于人工智能、计算机科学与技术等相关专业,更为适合; 4、下载使用后,可先查看README.md文件(如有),本项目仅用作交流学习参考,请切勿用于商业用途。
1、资源项目源码均已通过严格测试验证,保证能够正常运行; 2、项目问题、技术讨论,可以给博主私信或留言,博主看到后会第一时间与您进行沟通; 3、本项目比较适合计算机领域相关的毕业设计课题、课程作业等使用,尤其对于人工智能、计算机科学与技术等相关专业,更为适合; 4、下载使用后,可先查看README.md文件(如有),本项目仅用作交流学习参考,请切勿用于商业用途。1、资源项目源码均已通过严格测试验证,保证能够正常运行; 2、项目问题、技术讨论,可以给博主私信或留言,博主看到后会第一时间与您进行沟通; 3、本项目比较适合计算机领域相关的毕业设计课题、课程作业等使用,尤其对于人工智能、计算机科学与技术等相关专业,更为适合; 4、下载使用后,可先查看README.md文件(如有),本项目仅用作交流学习参考,请切勿用于商业用途。1、资源项目源码均已通过严格测试验证,保证能够正常运行; 2、项目问题、技术讨论,可以给博主私信或留言,博主看到后会第一时间与您进行沟通; 3、本项目比较适合计算机领域相关的毕业设计课题、课程作业等使用,尤其对于人工智能、计算机科学与技术等相关专业,更为适合; 4、下载使用后,可先查看README.md文件(如有),本项目仅用作交流学习参考,请切勿用于商业用途。1、资源项目源码均已通过严格测试验证,保证能够正常运行; 2、项目问题、技术讨论,可以给博主私信或留言,博主看到后会第一时间与您进行沟通; 3、本项目比较适合计算机领域相关的毕业设计课题、课程作业等使用,尤其对于人工智能、计算机科学与技术等相关专业,更为适合; 4、下载使用后,可先查看README.md文件(如有),本项目仅用作交流学习参考,请切勿用于商业用途。1、资源项目源码均已通过严格测试验证,保证能够正常运行; 2、项目问题、技术讨论,可以给博主私信或留言,博主看到后会第一时间与您进行沟通; 3、本项目比较适合计算机领域相关的毕业设计课题、课程作业等使用,尤其对于人工智能、计算机科学与技术等相关专业,更为适合; 4、下载使用后,可先查看README.md文件(如有),本项目仅用作交流学习参考,请切勿用于商业用途。1、资源项目源码均已通过严格测试验证,保证能够正常运行; 2、项目问题、技术讨论,可以给博主私信或留言,博主看到后会第一时间与您进行沟通; 3、本项目比较适合计算机领域相关的毕业设计课题、课程作业等使用,尤其对于人工智能、计算机科学与技术等相关专业,更为适合; 4、下载使用后,可先查看README.md文件(如有),本项目仅用作交流学习参考,请切勿用于商业用途。1、资源项目源码均已通过严格测试验证,保证能够正常运行; 2、项目问题、技术讨论,可以给博主私信或留言,博主看到后会第一时间与您进行沟通; 3、本项目比较适合计算机领域相关的毕业设计课题、课程作业等使用,尤其对于人工智能、计算机科学与技术等相关专业,更为适合; 4、下载使用后,可先查看README.md文件(如有),本项目仅用作交流学习参考,请切勿用于商业用途。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值