———反演第11页

bzoj 2820 mobius反演

学了一晚上mobius，终于A了一道了。。。。假设枚举到i，质数枚举到p(程序里的prime[j])，要更新A=i*p的信息。1.p|i这时A的素数分解式中，p这一项的次数>=2。考虑g(A)的求和式：如果枚举的质数p'不等于p，A/p'就也会有p这一项，次数>=2，这时候miu(A/p')=0如果枚举的质数p'=p，A/p=i，这一项就是miu(i)因此g(A)=miu(i)2.p!|i(即i%

weixin_34292402·2020-09-16 17:24

[SDOI2018]旧试题

[SDOI2018]旧试题Bzoj5332:[Sdoi2018]旧试题反演+三元环计数+卡常反演部分：第一步还是“约数个数和”那个结论的扩展：（考虑质数的次数的分配即可证明）最后一步就是枚举u,v,w考虑满足能凑出

weixin_34072458·2020-09-16 17:20

于神之怒加强版

然后大力反演一波：其实很套路。见e就miu，见两项du就D=du（这样的好处是把枚举两个变成枚举一个，另一个的枚举用枚举约数来代替。这样就有可能凑出卷积形式）后面的F函数是积性函数。直接sie

weixin_33909059·2020-09-16 17:19

数学(论)里的一些定理（莫比乌斯反演，傅立叶变换，数论变换...）

莫比乌斯反演莫比乌斯反演在数论中占有重要的地位，许多情况下能大大简化运算。那么我们先来认识莫比乌斯反演公式。

巷中人·2020-09-16 17:18

BZOJ2154: Crash的数字表格 & BZOJ2693: jzptab

【传送门：BZOJ2154&BZOJ2693】简要题意：给出n,m，求$\sum_{i=1}^{n}\sum_{j=1}^{m}LCM(i,j)$题解：莫比乌斯反演（因为BZOJ2693是多组数据，数据强一点

weixin_30954265·2020-09-16 17:47

BZOJ4407 于神之怒加强版 - 莫比乌斯反演

题解非常裸的莫比乌斯反演。

weixin_30337251·2020-09-16 17:45

【莫比乌斯反演】[HYSBZ/BZOJ2693]jzptab

题目大意就是求∑ni=1∑mj=1lcm(i,j)，但是有多组输入数据.如果之前做过【莫比乌斯反演】[HYSBZ\BZOJ2154]Crash的数字表格,就会发现，对于每一个询问，有O(n)的做法，但显然不够快

bo o ya ka·2020-09-16 17:14

bzoj 2301 Problem b 莫比乌斯反演+容斥

题意：对于给出的n个询问，每次求有多少个数对(x,y)，满足a≤x≤b，c≤y≤d，且gcd(x,y)=k，gcd(x,y)函数为x和y的最大公约数思路：在hdu1695的基础上加上容斥，即:ans=solve(b/k,d/k)-solve((a-1)/k,d/k)-solve((c-1)/k,b/k)+solve((a-1)/k,(c-1)/k),详见代码：/******************

Excelsior_kereo·2020-09-16 17:37

BZOJ2005:能量采集（Mobius反演）

传送门题意：求∑i=1n∑j=1mgcd(i,j)题解：Mobius反演。首先,n=∑d|nϕ(d)然后phi(d)可线性筛出。

DZYO·2020-09-16 17:18

【莫比乌斯反演】BZOJ2154Crash的数字表格&BZOJ2693jzptab

首先是BZOJ2154题目大意给出n,m求∑i=1n∑j=1mlcm(i,j)满足（n,m#include#include#include#include#defineSFscanf#definePFprintf#defineMAXN10000010intMAXPRIME;#defineMOD20101009usingnamespacestd;vectorprimes;intisprime[MAX

616156·2020-09-16 17:17

Bzoj2820:YY的GCD:莫比乌斯反演

题目链接：YY的GCD公式和上一题一样地推最后得到答案为令T=pd'，化简公式得到后面那个sigma可以线性筛筛出来，枚举素数贡献答案即可#include#include#include#include#include#include#include#include#include#include#defineLLlonglongusingnamespacestd;constintmaxn=100

TheWolfWhistlingSong·2020-09-16 17:16

【CJOJ2512】gcd之和（莫比乌斯反演）

题面给定n,m(n,m#include#include#include#include#include#include#include#include#includeusingnamespacestd;#defineMOD998244353#defineMAX10000000inlineintread(){intx=0,t=1;charch=getchar();while((ch'9')&&ch!

小蒟蒻yyb·2020-09-16 17:44

BZOJ_P2820 YY的GCD(数论+莫比乌斯反演)

BZOJ传送门TimeLimit:10SecMemoryLimit:512MBSubmit:1227Solved:636[Submit][Status][Discuss]Description神犇YY虐完数论后给傻×kAc出了一题给定N,M,求1#includeusingnamespacestd;#defineN10000005typedeflonglongll;intn,m,cnt,t;intp

BeiYu-oi·2020-09-16 17:11

【莫比乌斯反演】[HYSBZ\BZOJ2154]Crash的数字表格

题目题目大意就是求∑ni=1∑mj=1lcm(i,j).分析：ans=∑i=1n∑j=1mlcm(i,j)=∑i=1n∑j=1mi∗jgcd(i,j)枚举d=gcd(i,j),令f(n,m,k)=∑1#includeusingnamespacestd;#defineMAXN10000000#defineMOD20101009#defineSum(x,y)(1ll*x*(x+1)/2%MOD*(1l

outer_form·2020-09-16 17:39

bzoj 3601: 一个人的数论高斯消元&莫比乌斯反演

AC代码如下：#include#include#include#definelllonglong#definemod1000000007#defineinv(x)ksm(x,1000000005)usingnamespacestd;intm,n,a[105][105],p[1005],q[1005];intksm(intx,inty){intt=1;if(y>=1,x=(ll)x*x%mod)if

lych_cys·2020-09-16 17:33

bzoj3994 [SDOI2015]约数个数和（Mobius反演）

求∑ni=1∑mj=1σ0(ij)∑i=1n∑j=1mσ0(ij)我们有σ0(xy)=∑d1|x∑d2|y[gcd(d1,d2)=1]σ0(xy)=∑d1|x∑d2|y[gcd(d1,d2)=1]然后化式子就好了。最后答案就是∑d=1nμ(d)∑t1=1⌊nd⌋⌊⌊nd⌋t1⌋∑t2=1⌊md⌋⌊⌊md⌋t2⌋∑d=1nμ(d)∑t1=1⌊nd⌋⌊⌊nd⌋t1⌋∑t2=1⌊md⌋⌊⌊md⌋t2⌋O

Icefox_zhx·2020-09-16 17:21

洛谷4449 BZOJ4407 于神之怒加强版莫比乌斯反演线性筛

题目链接题意：给定n,m,kn,m,kn,m,k，求∑i=1n∑j=1mgcd(i,j)k\sum_{i=1}^n\sum_{j=1}^mgcd(i,j)^k∑i=1n∑j=1mgcd(i,j)k，有TTT组数据，每组数据的kkk是相同的。T<=2000,n,m,k<=5000000T<=2000,n,m,k<=5000000Tusingnamespacestd;const

forever_shi·2020-09-16 17:49

bzoj4407 于神之怒加强版

当你莫比乌斯反演发现复杂度不对的话怎么办？继续瞎**化简！令x=pd，那么有哎这个式子看起来很和善后面那个好像是积性函数（逃那么线性筛就行啦！！！问题是怎么筛呢

aklm45097·2020-09-16 17:35

「SDOI2018」旧试题（莫比乌斯反演）...

题目链接LOJ2476：https://loj.ac/problem/2476LOJ2565：https://loj.ac/problem/2565题解参考照搬了wxh的博客。为了方便，下文用$(x,y)$表示${\rmgcd}(x,y)$。先分析LOJ2476。注意到对于任意一个数组$a$，第$x$项的值$a_x$可以展开写成\(\sum_\limits{i=1}^{x}a_

anzi3457·2020-09-16 17:04

【bzoj2693】jzptab【反演】

反演是不是就是拿莫比乌斯函数乱搞……如果我说错了请回复--倒数第三行d’变成了倒数第四行的dd’……真是精妙。然后观察到d∑d′|dd′μ(d′)是积性函数，线性筛出来即可。

ZMOIYNLP·2020-09-16 17:29

BZOJ 4176 [莫比乌斯反演][杜教筛]

Description求∑i=1n∑j=1nd(ij)Solution通过陈老师r老师等式可以的得到该式子就等于∑i=1n∑j=1n⌊ni⌋⌊nj⌋[(i,j)=1]一波反演以后就得到∑d=1nμ(d)

Vectorxj·2020-09-16 17:58

BZOJ 2693 jzptab 莫比乌斯反演

题目大意：同2154多组数据后面那坨东西由于积性函数的约数和仍是积性函数因此只需要线性筛一下就行i%prime[j]==0那部分由于多出来的因数都不是无平方因子数因此μ值都为0增加的只有原先的D/i#include#include#include#include#defineM10001000#defineMOD100000009usingnamespacestd;typedeflonglongl

PoPoQQQ·2020-09-16 17:57

莫比乌斯反演详解

前言由于莫比乌斯反演的应用非常广泛，内容很多但是结论却并不复杂。然而如果没有接触过的话是很难得到超过暴力分的分数，最近的省选也经常考到，所以开单篇记一下公式。

Ripped·2020-09-16 17:57

BZOJ 4174 tty的求助莫比乌斯反演

题目大意：求∑Nn=1∑Mm=1∑m−1k=0⌊nk+xm⌋mod998244353假设n和m都已经确定了，现在要求这坨玩应：∑m−1k=0⌊nk+xm⌋=∑m−1k=0(⌊nk%m+xm⌋+nk−nk%mm)=∑m−1k=0(⌊nk%m+xm⌋+nkm−nk%mm)我们一项一项考虑令d=gcd(n,m)，那么有∑m−1k=0⌊nk%m+xm⌋=d∗∑md−1k=0⌊kd+xm⌋=d∗(md∗x−

PoPoQQQ·2020-09-16 17:26

YY的GCD，洛谷P2257，莫比乌斯反演+狄利克雷卷积

正题题目要求这个东西其实就是求根据换出来变成枚举d变成设，那么就变成发现做不了。看一下后面的东西，我们把它设成那么很明显发现prime函数指的是第i位是否为质数，是则为1，否则为0.然后整除分块一下，预处理狄利克雷卷积。然后n才1e7，按定义枚举#include#include#include#includeusingnamespacestd;intT;constintmaxn=1e7;intmu

Deep_Kevin·2020-09-16 17:50

【莫比乌斯反演】BZOJ4018 小Q的幻想之乡

【题目】原题地址题目可以转化为给定N,MN,M,求∑Ni=1∑Mj=1|i−j|gcd(i,j)∑i=1N∑j=1M|i−j|gcd(i,j)【题目分析】这种题就是反演辣，不过我不是很会，直接搬过来了大佬的

Dream_Lolita·2020-09-16 17:19

[BZOJ5276] Skyfall & [CF235E] Number Challenge [莫比乌斯反演]

LinkBZOJ-https://www.lydsy.com/JudgeOnline/problem.php?id=5276Codeforces-http://codeforces.com/contest/235/problem/EDescription∑i=1A∑j=1B∑k=1Cd(ijk)mod 230\sum\limits_{i=1}^A\sum\

*éphia·2020-09-16 17:19

【数论】Mobius反演学习笔记

于是我不得不开始治疗于是在写了n天某个半平面交题未果之后我决定先来搞搞mobius反演换换口味。。。说正文啦魂淡！好吧我不扯淡了开始正文时间。——————

CreationAugust·2020-09-16 17:18

莫比乌斯反演公式（数论）

转自：http://blog.csdn.net/f_zyj/article/details/51852672莫比乌斯反演公式则莫比乌斯函数µ另一种更常用的形式：在某一个范围内：则线性筛法求解/**莫比乌斯反演公式

BBHHTT·2020-09-16 17:17

BZOJ2154+BZOJ2693 Crash的数字表格&jzptab-莫比乌斯反演

1n∑j=1nlcm(i,j)=∑i=1n∑j=1nijgcd(i,j)∑i=1n∑j=1nlcm(i,j)=∑i=1n∑j=1nijgcd(i,j)然后开始拼命化简，未果听完dalao的提示：莫比乌斯反演题一定要想办法搞掉分母

Fizzmy·2020-09-16 17:46

BZOJ-2301-莫比乌斯

zyz_3_14159·2020-09-16 17:12

bzoj4659 Lcm （莫比乌斯反演）

题意求∑i=1n∑j=1mlcm(i,j)(gcd(i,j)无平方质因子)∑i=1n∑j=1mlcm(i,j)(gcd(i,j)无平方质因子)题解无平方质因子，即μ2=1μ2=1因此式子可以转化为∑i=1n∑j=1mlcm(i,j)(gcd(i,j)无平方质因子)=∑i=1n∑j=1mlcm(i,j)μ2(gcd(i,j))=∑i=1n∑j=1mijgcd(i,j)μ2(gcd(i,j))=∑d=

sunshiness_s·2020-09-16 17:05

BZOJ-2820-YY的GCD-（Mobius反演）

神犇YY虐完数论后给傻×kAc出了一题给定N,M,求1#defineINF0x3f3f3f3f#defineLLlonglong#definebugcoutMAXN)break;check[i*prime[j]]=1;if(i%prime[j]==0){mu[i*prime[j]]=0;break;}elsemu[i*prime[j]]-=mu[i];}}for(inti=0;im)swap(n,

EIKY·2020-09-16 17:03

BZOJ2154,BZOJ2693 (Mobius反演)

BZOJ2154:Crash的数字表格BZOJ2693:jzptabBZOJ2154给n,m,求∑i=1n∑j=1n[lcm(i,j)]将原式适当变形得：∑i=1n∑j=1n[lcm(i,j)]====∑i=1n∑j=1ni⋅jgcd(i,j)∑d∑i=1⌊nd⌋∑i=1⌊md⌋d2⋅i⋅jd[gcd(i,j)==1]∑dd∑i=1⌊nd⌋∑i=1⌊md⌋i⋅j⋅∑x|i∧x|jμ(x)∑dd∑x

DZYO·2020-09-16 17:32

BZOJ-2818-Gcd（欧拉函数/Mobius反演）

Description给定整数N，求1#defineINF0x3f3f3f3f#defineLLlonglong#definebugcout<<"bug\n"usingnamespacestd;constintMAXN=1e7+7;constintMAXM=1e9+7;longlongphi[MAXN];longlongprime[MAXN];intnum_prime;voidget_euler(

EIKY·2020-09-16 17:31

bzoj5332/loj2565 [Sdoi2018]旧试题（Mobius反演+图论）

求∑x=1A∑y=1B∑z=1Cd(xyz)∑x=1A∑y=1B∑z=1Cd(xyz)类似这道题的结论，我们有d(xyz)=∑i|x∑j|y∑k|z[(i,j)==1][(j,k)==1][(i,k)==1]d(xyz)=∑i|x∑j|y∑k|z[(i,j)==1][(j,k)==1][(i,k)==1]于是带进去，就是求∑x=1A∑y=1B∑z=1C∑i|x∑j|y∑k|z[(i,j)==1][

Icefox_zhx·2020-09-16 17:53

影像出现条带斜纹式的乱码——reshape操作中按行与按列取数据的差异

Case2：做反演的时候，矩阵变向量输入模型，输出向量再变回矩阵时，reshape使用错误时，出现上图结果。如，对于输入的矩阵，用以下代码生成输入到模型里面的列向量：f

-Y-Q-Q-·2020-09-16 11:51

莫比乌斯反演

莫比乌斯反演在数论中占有重要的地位，许多情况下能大大简化运算。那么我们先来认识莫比乌斯反演公式。

XTU_wcb·2020-09-15 23:15

20200802 T3 我永远喜欢【生成函数容斥，拉格朗日反演】

题目描述有nnn种颜色的石子，每种cic_ici个，记一个石子序列首尾相接后极长连续段的长度为lil_ili，求所有石子序列的1∏li!\frac1{\prodl_i!}∏li!1的和。n≤105,∑ci≤2∗105n\le10^5,\sumc_i\le2*10^5n≤105,∑ci≤2∗105题目分析先考虑去掉首尾相接的情况怎么做。因为限制了极长，容易想到把每种颜色分成几段，然后合并，但是不好保

Master.Yi·2020-09-15 13:44

51nod(1220 约数之和推式子反演+线性筛+杜教筛)

学习链接2.慢慢推式子需要用到一个小结论(挺好yy证明的)+反演+积累技巧把(比如dd前缀和可以转化成别的式子进而分块求)学习链接第三个式子是反演得到的。

Helium_wild·2020-09-15 03:37

[莫比乌斯反演] bzoj2301: [HAOI2011]Problem b

id=2301一看题目模板题模板题：caioj1280:[视频]【莫比乌斯反演模板题】GCDhttp://caioj.cn/problem.php?

kkkGIGi_qtt·2020-09-15 02:03

【容斥+莫比乌斯反演】BZOJ2301 [HAOI2011]Problem b

题面在这里首先容斥，把问题转化为求∑i=1n∑j=1m[gcd(i,j)=k]⇒∑i=1⌊nk⌋∑j=1⌊mk⌋[gcd(i,j)=1]由u∗I=e，得到：∑i=1⌊nk⌋∑j=1⌊mk⌋∑d|(i,j)μ(d)⇒∑dμ(d)∑d|i∑d|j1⇒∑dμ(d)⌊nkd⌋⌊mkd⌋由于⌊nkd⌋⌊mkd⌋是可以分块处理的于是每个询问都可以O(n√)得到答案示例程序：#include#includeus

linkfqy·2020-09-15 02:56

BZOJ 2301 HAOI2011 Problem b 容斥原理+莫比乌斯反演

题目大意：多次询问有多少个数对(x,y)满足a#include#include#include#defineM100100usingnamespacestd;typedeflonglongll;llmu[M],prime[M],tot;boolnot_prime[M];voidLinear_Shaker(){inti,j;mu[1]=1;for(i=2;i>T;T;T--){scanf("%d%d

PoPoQQQ·2020-09-15 02:55

[51nod][积性函数][杜教筛]最大公约数之和 V3

简单反演∑ndd∑ndi∑ndj(i,j)=1∑dnd∑ind∑jnd(i,j)=1∑ndd(2(∑ndjϕ(i)

Gzb1128·2020-09-15 02:47

基于飞桨实现高光谱反演：通过遥感数据获取土壤某物质含量

高光谱反演是什么？高光谱反演是使用遥感卫星拍摄的高光谱数据以及实地采样化验的某物质含量数据来建立一个反演模型。简单来说就是：有模型以后卫星一拍，就能得知土壤中某物质的含量，不用实地采样化验了。

百度大脑·2020-09-15 00:26

F的ACM暑期集训计划

暑假的知识计划(补充中...)1.数论相关(7days)待完成多项式同余方程/高次同余方程/欧拉函数/克莱姆法则/高斯消元/莫比乌斯反演/伪素数判定/baby-step-gaint-step2.组合数学相关

weixin_30512785·2020-09-14 23:34

CodeForces 900D-Unusual Sequences(快速幂，莫比乌斯反演)

题解练习题里的题，我原来根本看不出是莫比乌斯反演

Terry: ）·2020-09-14 21:49

Unusual Sequences(容斥莫比乌斯反演)

题意题目链接Sol首先若y%x不为0则答案为0否则，问题可以转化为，有多少个数列满足和为y/x，且整个序列的gcd=1考虑容斥，设$g[i]$表示满足和为$i$的序列的方案数，显然$g[i]=2^{i-1}$(插板后每空位放不放)同时还可以枚举一下gcd，设$f[i]$表示满足和为$i$且所有数的gcd为1的方案，\(g[i]=\sum_{d|i}f[\frac{n}{d}]\

weixin_30716141·2020-09-14 20:22

牛客练习赛69 F.解方程（莫比乌斯反演 + 迪利克雷卷积性质 + 欧拉筛）

n,p,qn,p,qn,p,q均为10710^7107，这个规模下可以考虑求出所有的f(i)f(i)f(i)，题目给出的式子很明显要用莫比乌斯反演，因为莫比乌斯反演的形式为：形式一：设F(n)=∑d∣nf

猝死在学ACM的路上·2020-09-14 20:29

HDU 6868：Absolute Math（莫比乌斯反演 + 离线暴力）

f(n)=∑d∣nμ(d)2f(n)=\displaystyle\sum_{d|n}\mu(d)^2f(n)=d∣n∑μ(d)2，根据μ(d)\mu(d)μ(d)的定义，当ddd的质因子的幂次不高于一次时有贡献，设ddd的质因子共有w(d)w(d)w(d)种，f(n)=∑d∣nμ(d)2=∑i=0w(d)C(w(d),i)=2w(d)f(n)=\displaystyle\sum_{d|n}\mu(