莫比乌斯反演
Facico
生命是闪耀的此刻,不是过程,就像芳香不需要道路一样。
展开
-
【2011集训队出题】Crash的数字表格
Description求 \sum_{i=1}^n\sum_{j=1}^mlcm(i,j)\mod\sum_{i=1}^n\sum_{j=1}^mlcm(i,j)\mod原创 2016-03-18 20:38:34 · 974 阅读 · 0 评论 -
线性筛μ的模板
时效O(n)O(n) 因为质数的μ都是-1,然后一个数t里面如果有多个质数=zhi[j]*i,i是有很多个与zhi[j]互异的质数相乘的,那么μ[t]=-μ[i]。miu[1]=1;qian[1]=1; fo(i,2,n){ if(!bz[i]){ zhi[++zhi[0]]=i; miu[i]=-1; }原创 2016-03-18 21:46:34 · 610 阅读 · 0 评论 -
【51Nod1244】莫比乌斯函数之和
Description莫比乌斯函数,由德国数学家和天文学家莫比乌斯提出。梅滕斯(Mertens)首先使用μ(n)(miu(n))作为莫比乌斯函数的记号。具体定义如下: 如果一个数包含平方因子,那么miu(n) = 0。例如:miu(4), miu(12), miu(18) = 0。 如果一个数不包含平方因子,并且有k个不同的质因子,那么miu(n) = (-1)^k。例如:miu(2), miu原创 2016-07-06 21:58:41 · 2494 阅读 · 0 评论 -
【GDSOI 2016】互补约数
DescriptionSolution看到对于gcd的求和就知道是关于莫比乌斯繁衍的题目。 先来看一看要求什么:所有数分解成两个相乘的数的gcd,那么我们可以直接枚举那些相乘的数。 所以 ans=∑i=1n∑j=1⌊ni⌋gcd(i,j)ans=\sum_{i=1}^n\sum_{j=1}^{\lfloor{n\over i}\rfloor}gcd(i,j) 和一般的莫比乌斯繁衍的长得很像,原创 2016-08-03 09:14:45 · 809 阅读 · 0 评论 -
【51Nod 1180】方格射击游戏
DescriptionM*N的方格矩阵,一个人在左下角格子的中心,除他所站位置外,其他格子的中心都有一个敌人,他一次可发射一枚子弹干掉一条直线上的所有敌人,问至少要发射多少子弹才能干掉所有敌人。Solution在一条直线上如果经过了两个点及以上,那么第一个经过的点一定是互质的,因为经过的第一个点的斜率就可以确定,如果不是互质,那么肯定是在前面确定的。 所以答案等于2+∑n−1i=1∑m−1j=1(原创 2016-10-20 19:55:55 · 1121 阅读 · 0 评论 -
【SDOI2017】数字表格
DescriptionSolution这道题明显可以用莫比乌斯反演来做,非常的裸的繁衍。 假设n < m ans=∏ni=1∏mj=1f[gcd(i,j)]ans=\prod_{i=1}^n\prod_{j=1}^mf[gcd(i,j)] 那么直接枚举gcd=d ans=∏nd=1f[d]g[d]ans=\prod_{d=1}^nf[d]^{g[d]} g[d]表示gcd=d的个数g[d]原创 2017-04-18 21:09:11 · 1287 阅读 · 2 评论 -
JZOJ 5224【GDOI2018模拟7.12】C(数论&&杜教筛)
DescriptionSolution这题直接用反演可以过70分。 但是一看到这种题就反演还是不行。 我们表示gcd(i,j)=dgcd(i,j)=d还可以用gcd(id,jd)=1=φ(id),i≥jgcd({i\over d},{j\over d})=1=\varphi({i\over d}),i≥j来搞 ∑i=1n∑j=1ngcd(i,j)k\sum_{i=1}^n\sum_{j=1}^原创 2017-07-14 22:13:25 · 499 阅读 · 0 评论