正题

题目链接:https://www.luogu.com.cn/problem/P1829

题目大意

给出n,mn,mn,m求∑i=1n∑j=1mlcm(i,j)\sum_{i=1}^n\sum_{j=1}^mlcm(i,j)i=1∑nj=1∑mlcm(i,j)

解题思路

∑i=1n∑j=1mijgcd(i,j)\sum_{i=1}^n\sum_{j=1}^m\frac{ij}{gcd(i,j)}i=1∑nj=1∑mgcd(i,j)ij
∑x=1nx×(∑i=1⌊nx⌋∑j=1⌊mx⌋ij×[gcd(i,j)==1])\sum_{x=1}^nx\times (\sum_{i=1}^{\lfloor\frac{n}{x}\rfloor}\sum_{j=1}^{\lfloor\frac{m}{x}\rfloor}ij\times [gcd(i,j)==1])x=1∑nx×(i=1∑⌊xn⌋j=1∑⌊xm⌋ij×[gcd(i,j)==1])
∑x=1nx×∑x∣dμ(dx)d2∗∑i=1⌊nx⌋∑j=1⌊mx⌋ij\sum_{x=1}^nx\times \sum_{x|d}\mu(\frac{d}{x})d^2*\sum_{i=1}^{\lfloor\frac{n}{x}\rfloor}\sum_{j=1}^{\lfloor\frac{m}{x}\rfloor}ijx=1∑nx×x∣d∑μ(xd)d2∗i=1∑⌊xn⌋j=1∑⌊xm⌋ij
单独拎出后面的来看
∑x∣dμ(dx)d2∗∑i=1⌊nx⌋∑j=1⌊mx⌋ij\sum_{x|d}\mu(\frac{d}{x})d^2*\sum_{i=1}^{\lfloor\frac{n}{x}\rfloor}\sum_{j=1}^{\lfloor\frac{m}{x}\rfloor}ijx∣d∑μ(xd)d2∗i=1∑⌊xn⌋j=1∑⌊xm⌋ij
∑x∣dμ(dx)d2∗∑i=1⌊nd⌋i∗∑j=1⌊md⌋j\sum_{x|d}\mu(\frac{d}{x})d^2*\sum_{i=1}^{\lfloor\frac{n}{d}\rfloor}i*\sum_{j=1}^{\lfloor\frac{m}{d}\rfloor}jx∣d∑μ(xd)d2∗i=1∑⌊dn⌋i∗j=1∑⌊dm⌋j
然后整除分块处理这一个
之后再在外面套一个整除分块即可。

时间复杂度O(n+m)O(n+m)O(n+m)

codecodecode

#include<cstdio>
#include<cstring>
#include<algorithm>
#define ll long long
using namespace std;
const ll N=1e7+10,XJQ=20101009;
ll n,m,mu[N],pri[N],g[N],cnt,ans;
bool v[N];
void ycl(ll n){mu[1]=1;for(ll i=2;i<=n;i++){if(!v[i])pri[++cnt]=i,mu[i]=-1;for(ll j=1;j<=cnt&&i*pri[j]<=n;j++){v[i*pri[j]]=1;if(i%pri[j]==0)break;mu[i*pri[j]]=-mu[i];}}for(ll i=1;i<=n;i++)mu[i]=(mu[i-1]+mu[i]*i*i%XJQ)%XJQ;return;
}
ll solve(ll n,ll m){ll ans=0;for(ll l=1,r;l<=n;l=r+1){r=min(n/(n/l),m/(m/l));ll q=n/l,p=m/l;ans=(ans+(q*(q+1)/2%XJQ)*(p*(p+1)/2%XJQ)%XJQ*(mu[r]-mu[l-1]+XJQ)%XJQ)%XJQ;}return ans;
}
int main()
{scanf("%lld%lld",&n,&m);if(n>m)swap(n,m);ycl(m);for(ll l=1,r;l<=n;l=r+1){r=min(n/(n/l),m/(m/l));(ans+=(r+l)*(r-l+1)/2%XJQ*solve(n/l,m/l)%XJQ)%=XJQ;}printf("%lld\n",ans);
}

P1829-[国家集训队]Crash的数字表格/JZPTAB【莫比乌斯反演】相关推荐

P1829 [国家集训队]Crash的数字表格 / JZPTAB
P1829 [国家集训队]Crash的数字表格 / JZPTAB 题意: 求∑i=1n∑j=1mlcm(i,j)\sum_{i=1}^{n}\sum_{j=1}^{m}lcm(i,j)∑i=1n∑j ...
P1829 [国家集训队]Crash的数字表格 / JZPTAB（莫比乌斯反演）
[国家集训队]Crash的数字表格 / JZPTAB 题目描述今天的数学课上,Crash 小朋友学习了最小公倍数(Least Common Multiple).对于两个正整数 a a a 和 b b ...
BZOJ 2154 [国家集训队]Crash的数字表格 / JZPTAB（莫比乌斯反演，经典好题）（Luogu P1829）
整理的算法模板合集: ACM模板点我看算法全家桶系列!!! 实际上是一个全新的精炼模板整合计划 P1829 [国家集训队]Crash的数字表格 / JZPTAB(反演,经典好题) Problem S ...
P1829 [国家集训队]Crash的数字表格（推了好久的mobius反演）
P1829 [国家集训队]Crash的数字表格 / JZPTAB 推导过程 ∑i=1n∑j=1mlcm(i,j)\sum_{i = 1} ^{n} \sum_{j = 1} ^{m} lcm(i, j ...
【数论】Crash的数字表格 / JZPTAB（P1829）
正题 P1829 题目大意给出n,m,求∑i=1n∑j=1mlcm(i,j)\sum_{i=1}^n\sum_{j=1}^mlcm(i,j)i=1∑nj=1∑mlcm(i,j) 解题思路 ∑i= ...
P3704-[SDOI2017]数字表格【莫比乌斯反演】
正题题目链接:https://www.luogu.com.cn/problem/P3704 题目大意 TTT组询问,给出n,mn,mn,m求∏i=1n∏j=1mFbigcd(i,j)\prod_{i ...
BZOJ 4816[SDOI2017]数字表格（莫比乌斯反演）
题目链接 \(Description\) 用\(f_i\)表示\(fibonacci\)数列第\(i\)项,求\(\prod_{i=1}^{n}\prod_{j=1}^{m}f[gcd(i,j)]\) ...
bzoj2154 Crash的数字表格
2154: Crash的数字表格 Time Limit: 20 Sec Memory Limit: 259 MB Submit: 4549 Solved: 1643 [Submit][Status ...
【BZOJ2117】 [2010国家集训队]Crash的旅游计划
[BZOJ2117] [2010国家集训队]Crash的旅游计划 Description 眼看着假期就要到了,Crash由于长期切题而感到无聊了,因此他决定利用这个假期和好友陶陶一起出去旅游. Cra ...

P1829-[国家集训队]Crash的数字表格/JZPTAB【莫比乌斯反演】

正题

题目大意

解题思路

codecodecode

P1829-[国家集训队]Crash的数字表格/JZPTAB【莫比乌斯反演】相关推荐

最新文章

热门文章