莫比乌斯反演
文章平均质量分 68
alan_cty
蒟蒻一只
别打脸(⊙o⊙)哦
展开
-
关于欧拉函数的一个性质
昨天,有一个人过来问我,如何证明∑d|nϕ(d)=n\sum_{d|n}\phi(d)=n 毫无头绪,一群人乱搞乱搞,一个晚上都没证出来。于是决定组队去请(bei)教(diao)PhilipsWeng大犇。大犇愣了三秒(注意,3s),然后给出了一种高级的证法:已知ϕ(n)=∑ni=1e((i,n))\phi(n)=\sum_{i=1}^{n}{e((i,n))}然后e(n)=∑d|nμ(d)e(n原创 2016-02-19 11:29:04 · 1772 阅读 · 7 评论 -
【GDOI2018模拟8.12】求和
Description求∑i=1n∑j=1n∑d=1kfd(gcd(i,j))\sum_{i=1}^{n}\sum_{j=1}^{n}\sum_{d=1}^{k}f_d(\gcd(i,j)) 其中当n=∏ipaiin=\prod_{i}p_i^{a_i}时fd(n)=∏i(−1)ai[ai<=d]f_d(n)=\prod_{i}(-1)^{a_i}[a_i<=d] 答案对2^30取模 n<=1原创 2017-08-21 12:01:31 · 365 阅读 · 0 评论 -
[51nod1190]最小公倍数之和V2
Description给出a,b,求∑i=abgcd(i,b)\sum_{i=a}^{b}\gcd(i,b) a,b<=1e9,数据组数<=1e5,答案对1e9+7取模Solution看到gcd想反演(然而这个是lcm) 这个反演不是正常套路 坑了我好久才跳出来 首先ans=b∑d|b1d∑i=abi[gcd(i,b)=d]ans=b\sum_{d|b}{1\over d}\sum_{i=a原创 2017-06-29 21:12:39 · 894 阅读 · 1 评论 -
[51nod1375]再选数
Description给出n个数,{a},求在其中选k个数使其gcd为1的方案数。 如果k=-1表示任意取,不过至少得取一个。 答案对998244353取模。 n<=10^5,ai<=10^6Solution记得某位dalao说过,看到gcd想到反演=w= 但这道题都能用反演?!涨姿势了~ 开始一直在想容斥,想不出来,然后去Orz了栋爷的题解(数论蒟蒻不要鄙视我QwQ 然后看到第一句话就原创 2016-12-28 20:45:36 · 840 阅读 · 0 评论 -
[51nod1223]分数等式的数量
Description给出一个数L,求所有的x < y <=l且满足1/x+1/y=1/n(n为整数)的(x,y)二元组的数量。 L<=10^11Solution近似暴力算法碾过去了=w= 安利一份更劲的题解 首先我们就相当于求x+y|xyx+y|xy的二元组的数量 提取一个d=gcd(x,y)d=gcd(x,y), x′=x/dx'=x/d, y′=y/dy'=y/d,那么 d(x′+y原创 2016-12-22 22:06:22 · 549 阅读 · 0 评论 -
[51nod1192]gcd表中的质数
Description求∑i=1n∑j=1me(gcd(i,j)是质数)\sum_{i=1}^{n}\sum_{j=1}^{m}e(gcd(i,j)是质数) T<=1000,n,m<=5*10^6Solution经典反演套路题,貌似没有其他做法(其他做法的大爷不要鄙视蒟蒻=w=) 复习一下莫比乌斯反演,所以就来打了这道题 首先约定n<mn<m 设FdF_d表示∑ni=1∑mj=1e(gcd(原创 2016-11-30 20:11:31 · 798 阅读 · 2 评论 -
[51nod1188]最大公约数之和 V2
Description求∑i=1n−1∑j=i+1ngcd(i,j)\sum_{i=1}^{n-1}\sum_{j=i+1}^{n}gcd(i,j) n<=5*10^6,数据组数<=50000Solution一眼看到莫比乌斯反演(感觉中毒了) 但是会超时!!! 想想其他方法。 我们知道,如果要在区间[1..i]中询问gcd(i,x)=t的个数。那么也就相当于询问gcd(i/t,x/t)=1的原创 2016-07-06 15:16:40 · 1293 阅读 · 0 评论 -
[51nod1239]欧拉函数之和
Description求∑i=1nφ(i)\sum_{i=1}^{n}\varphi(i) n<=10^10Solution这道题和莫比乌斯函数一行,都可以通过神奇的推导的出结论。 我们设ϕ(n)=∑i=1nφ(i)\phi(n)=\sum_{i=1}^{n}\varphi(i) 众所周知,∑d|nφ(d)=n\sum_{d|n}\varphi(d)=n 那么,φ(n)=n−∑d|n,d<n原创 2016-07-06 10:37:41 · 2136 阅读 · 0 评论 -
【2011集训队出题】Crash的数字表格
Description求∑i=1n∑j=1mlcm(i,j)\sum_{i=1}^{n}\sum_{j=1}^{m}lcm(i,j) n,m<=10^7Solution(注意,以下内容默认n<=mn<=m) 看着lcm不爽,把它变一变: ∑i=1n∑j=1mijgcd(i,j)\sum_{i=1}^{n}\sum_{j=1}^{m}{ij\over \gcd(i,j)} 莫比乌斯反演常用,枚原创 2016-03-17 21:53:42 · 927 阅读 · 0 评论 -
[bzoj2005][Noi2010]能量采集
Description在平面直角坐标系中,点(x,y)的代价定义为它和原点的连线中经过多少的其他整点个数*2+1。求横坐标在1~n且纵坐标在1~m的所有点的代价和。 n,m<=10^5Solution首先,一个点的代价就是gcd(x,y)*2-1(蒟蒻不会证,各路大犇们带带我呗) 然后,我们就变成了要求∑i=1n∑j=1mgcd(i,j)\sum_{i=1}^{n}\sum_{j=1}^{m}g原创 2016-06-16 19:55:08 · 1567 阅读 · 0 评论 -
[bzoj2818]gcd
Description求∑i=1N∑j=1Ngcd(i,j)为质数的个数\sum_{i=1}^{N}\sum_{j=1}^{N}gcd(i,j)为质数的个数 N<=10^7Solution很显然的莫比乌斯反演~(≧▽≦)/~啦啦啦 然而本蒟蒻只会这种傻逼方法,跑了 WerKeyTom_FTD爷用了机智的phi法,跑的飞起。 好吧,回归正题。 首先,我们知道,根据普通的莫比乌斯反演,原创 2016-05-30 20:53:38 · 991 阅读 · 0 评论 -
GDSOI 2016 T1 互补约数
Description求∑i=1n∑d|igcd(d,id)\sum_{i=1}^{n}\sum_{d|i}gcd(d,{i\over d}) n<=10^11Solution首先,我们发现gcd中的两个东西是所有乘积不超过n的数对,即 Ans=∑i∑j,i∗j<=ngcd(i,j)Ans=\sum_{i}\sum_{j,i*j<=n}gcd(i,j) 然后Ans=∑i=1n∑j=1⌊ni⌋g原创 2016-05-25 20:49:11 · 1006 阅读 · 2 评论 -
线筛模板
忘了的时候来抄模板(我记性差~)scanf("%d%d",&n,&m);if (n>m) swap(n,m);mu[1]=1;fo(i,2,n) { if (!f[i]) p[++tot]=i,mu[i]=-1; fo(j,1,tot) { int k=p[j]*i;if (k>n) break;f[k]=1; if (!(i%p[j])) break;m原创 2016-03-16 20:46:48 · 747 阅读 · 0 评论 -
「MtOI2019」幽灵乐团
Description求∏i=1A∏j=1B∏k=1C([i,j](i,k))f(type)\prod_{i=1}^{A}\prod_{j=1}^{B}\prod_{k=1}^{C}({[i,j]\over (i,k)})^{f(type)}i=1∏Aj=1∏Bk=1∏C((i,k)[i,j])f(type)f(0)=1,f(1)=ijk,f(2)=(i,j,k)f(0)=1,f(1)...原创 2019-09-02 21:23:46 · 534 阅读 · 0 评论