莫比乌斯反演
文章平均质量分 69
_beginend
这个作者很懒,什么都没留下…
展开
-
bzoj 2440: [中山市选2011]完全平方数 莫比乌斯函数
题目传送门分析莫比乌斯函数μ(i)=(−1)k(-1)^k(i为k个互不相等的素数的成绩),否则μ(i)=0 那么对于这题,可以想到用容斥原理,那么答案就是n-奇数个质数的平方的倍数的个数+偶数个质数的平方的倍数的个数 那么ans=Σμ[i]*(n/i^2) (i<=sqrt(n))记得开long long代码#include<iostream>#include<cstdio>#includ原创 2016-11-08 12:20:15 · 334 阅读 · 0 评论 -
51nod 1222 最小公倍数计数 莫比乌斯反演+杜教筛
题意定义F(n)表示最小公倍数为n的二元组的数量。 即:如果存在两个数(二元组)X,Y(X <= Y),它们的最小公倍数为N,则F(n)的计数加1。 例如:F(6) = 5,因为[2,3] [1,6] [2,6] [3,6] [6,6]的最小公倍数等于6。 给出一个区间[a,b],求最小公倍数在这个区间的不同二元组的数量。 例如:a = 4,b = 6。符合条件的二元组包括: [1,4]原创 2017-10-30 20:03:23 · 364 阅读 · 0 评论 -
51nod 1355 斐波那契的最小公倍数 莫比乌斯反演+数学
题意F(0) = 0 F(1) = 1 F(n) = F(n-1) + F(n-2) 给出n个正整数a1, a2,…… an,求对应的斐波那契数的最小公倍数,由于数字很大,输出Mod 1000000007的结果即可。 2 <= N <= 50000,1 <= ai <= 1000000。分析跪zyz大爷写的题解。首先有个结论就是 gcd(fn,fm)=fgcd(n,m)gcd(f_n,f_m原创 2017-12-24 15:44:17 · 550 阅读 · 0 评论 -
51nod 1675 序列变换 莫比乌斯反演
题意lyk有两序列a和b。 lyk想知道存在多少对x,y,满足以下两个条件。 1:gcd(x,y)=1。 2: abx=baya_{bx} = b_{ay} 。 例如若a={1,1,1},b={1,1,1}。那么存在7对,因为除了x=2,y=2或x=3,y=3外都满足条件。 n<=100000,1<=a[i],b[i]<=n分析一开始没什么思路,后来看到评论有人说用反演可做,于是就会了。原创 2017-08-12 20:17:08 · 497 阅读 · 0 评论 -
51nod 1227 平均最小公倍数 莫比乌斯反演+杜教筛
题意Lcm(a,b)表示a和b的最小公倍数,A(n)表示Lcm(n,i)的平均数(1 例如:A(4) = (Lcm(1,4) + Lcm(2,4) + Lcm(3,4) + Lcm(4,4)) / 4 = (4 + 4 + 12 + 4) / 4 = 6。 F(a, b) = A(a) + A(a + 1) + …… A(b)。(F(a,b) = ∑A(k), a 例如:F(2, 4)原创 2018-01-05 07:53:51 · 291 阅读 · 0 评论 -
51nod 1584 加权约数和 莫比乌斯反演
题意去年的 tangjz 非常喜欢做数论题,但是一年以后的 tangjz 却不那么会做了。 在整理以前的试题时,他发现了这样一道题目:“求 ∑σ(i) ,其中 1≤i≤N , σ(i) 表示 i 的约数之和。” 现在他长大了,题目也变难了,所以麻烦你来帮他解决一道数论题吧。 他需要你求如下表达式的值: ∑ni=1∑nj=1max(i,j)∗σ(i∗j)\sum_{i=1}^n\sum_原创 2018-01-08 21:53:41 · 366 阅读 · 0 评论 -
51nod 1594 Gcd and Phi 莫比乌斯反演
题意F(n)=∑ni=1∑nj=1φ(φ(i),φ(j))F(n)=\sum_{i=1}^n\sum_{j=1}^n\varphi_{(\varphi(i),\varphi(j))} 其中φ\varphi表示欧拉函数。欧拉函数φn\varphi_n是不超过n的数中与n互质的数的数目。 φ(φ(i),φ(j))\varphi_{(\varphi(i),\varphi(j))}表示i,j欧拉函原创 2018-01-09 10:09:13 · 289 阅读 · 0 评论 -
bzoj 4176: Lucas的数论 莫比乌斯反演+杜教筛
题意设f(i)f(i)表示i的约数个数。给出n,求∑ni=1∑nj=1f(ij)\sum_{i=1}^n\sum_{j=1}^nf(ij)。答案模1e9+7. n分析按套路来,有f(nm)=∑i|n∑j|m[(i,j)=1]f(nm)=\sum_{i|n}\sum_{j|m}[(i,j)=1] 把这条式子代进原式,反演一波,就可以得到ans=∑d=1nμ(d)(∑i=1⌊nd⌋原创 2018-01-09 20:18:30 · 323 阅读 · 0 评论 -
Codeforces 235E Number Challenge 莫比乌斯反演+数论
题意设d(s)表示s的约数个数,给出a,b,c,求∑i=1a∑j=1b∑k=1cd(ijk)∑i=1a∑j=1b∑k=1cd(ijk)\sum_{i=1}^a\sum_{j=1}^b\sum_{k=1}^cd(ijk) a,b,c<=2000分析题解貌似是一个很鬼畜的做法。。。 约数个数函数σ0(d)σ0(d)\sigma_0(d)有一个小结论,就是σ0(ij)=∑p|i...原创 2018-02-28 20:08:40 · 469 阅读 · 0 评论 -
bzoj 4174: tty的求助 数学+莫比乌斯反演
题意 N,M<=500000,0<x<=100000N,M<=500000,0<x<=100000N,M∑k=0m−1⌊nk+xm⌋∑k=0m−1⌊nk+xm⌋\sum_{k=0}^{m-1}\lfloor\frac{nk+x}{m}\rfloor =∑k=0m−1⌊nkmodm+xm⌋+nk−nkmodmm=∑k=0m−1⌊nkmodm+xm⌋+nk−...原创 2018-03-18 20:50:07 · 279 阅读 · 0 评论 -
luogu 3768 简单的数学题 莫比乌斯反演+杜教筛
题意给出n,pn,pn,p,求(∑i=1n∑j=1nijgcd(i,j))modp(∑i=1n∑j=1nijgcd(i,j))modp(\sum_{i=1}^n\sum_{j=1}^nijgcd(i,j))\bmod p n≤1010n≤1010n\le10^{10}分析不多说,直接上式子∑i=1n∑j=1nijgcd(i,j)∑i=1n∑j=1nijgcd(i,j)\sum_{...原创 2018-04-26 14:43:06 · 248 阅读 · 0 评论 -
LibreOJ #2563.「SDOI2018」反回文串 pollard_rho+莫比乌斯反演
题意求所有长度为n字符集大小为k的字符串中有多少个字符串可以由某个回文串移位若干次后得到。 n,k≤1018n,k≤1018n,k\le10^{18}分析不难发现一个回文串移位若干位后一定可以得到一个合法的串,关键在于如何不重不漏的计数。 若一个回文串移位x次后得到的字符串也是一个回文串,则继续往后移位则必然会算重。 那移位多少次后会得到一个新的回文串呢?这显然跟原串的周期有...原创 2018-06-20 19:54:04 · 401 阅读 · 0 评论 -
Codeforces 809E Surprise me! 莫比乌斯反演+点分治
题意给一棵n个节点的数,点权为1到n的整数且互不相同。求∑i≠jφ(vali∗valj)∗dis(i,j)n∗(n−1)\frac{\sum_{i=\not j}\varphi(val_i*val_j)*dis(i,j)}{n*(n-1)}n∗(n−1)∑i≠jφ(vali∗valj)∗dis(i,j)n≤2∗105n\le2*10^5n≤2∗105分析注意到φ(ab)=φ(...原创 2019-07-02 11:11:01 · 202 阅读 · 0 评论 -
51nod 1192 Gcd表中的质数 莫比乌斯反演
题意有一个M * N的表格,行与列分别是1 - M和1 - N,格子中间写着行与列的最大公约数Gcd(i, j)(1 <= i <= M, 1 <= j <= N)。 给出M和N,求这张表中有多少个质数。 T<=1000,n,m<=5000000分析比较好的思路。[p][p]表示若p为素数则为1否则为0. ∑p∑i=1⌊np⌋∑j=1⌊mp⌋[gcd(i,j)=1][p]\sum_p\sum_原创 2017-08-16 22:32:58 · 383 阅读 · 0 评论 -
bzoj 4815: [Cqoi2017]小Q的表格 分块+莫比乌斯反演
题意小Q是个程序员。 作为一个年轻的程序员,小Q总是被老C欺负,老C经常把一些麻烦的任务交给小Q来处理。每当小Q不知道如何解决时,就只好向你求助。为了完成任务,小Q需要列一个表格,表格有无穷多行,无穷多列,行和列都从1开始标号。为了完成任务,表格里面每个格子都填了一个整数,为了方便描述,小Q把第a行第b列的整数记为f(a,b),为了完成任务,这个表格要满足一些条件:(1)对任意的正整数a,b,都要原创 2017-04-14 17:00:00 · 1395 阅读 · 0 评论 -
bzoj 2820: YY的GCD 莫比乌斯反演
题意给出n,m,求有多少数对(x,y)满足gcd(x,y)为质数 n,m<=10000000 多组数据T<=10000分析设g(d)表示1<=x<=n,1<=y<=m中gcd(x,y)=d的数对数量设f(d)表示1<=x<=n,1<=y<=m中d|gcd(x,y)的数对数量显然f(d)=∑d|ng(n)=⌊nd⌋∗⌊md⌋f(d)=\sum_{d|n}g(n)=\lfloor \frac{n}{原创 2017-01-15 20:54:05 · 438 阅读 · 1 评论 -
bzoj 3529: [Sdoi2014]数表 莫比乌斯反演+树状数组
题意有一张N×m的数表,其第i行第j列(1 < =i < =礼,1 < =j < =m)的数值为能同时整除i和j的所有自然数之和。给定a,计算数表中不大于a的数之和。 1 < =N.m < =10^5 , 1 < =Q < =2×10^4分析多做了几道反演题终于找到了一点所谓的套路了,也就是把式子化简之后求前缀和然后用分块处理。直接上PoPoQQQ大佬%%%的题解代码#include<iostr原创 2017-01-16 11:10:24 · 391 阅读 · 0 评论 -
bzoj 2154: Crash的数字表格 莫比乌斯反演
题意求∑ni=1∑mj=1lcm(i,j)\sum_{i=1}^n\sum_{j=1}^mlcm(i,j) n,m<=10000000分析题解见Po神的论文 或者看这篇比较详细的博客 恩全OJ倒数第二……代码#include<iostream>#include<cstdio>#include<cstdlib>#include<cstring>#include<algorithm>#d原创 2017-01-17 20:49:11 · 535 阅读 · 0 评论 -
bzoj 2693: jzptab 莫比乌斯反演
题意 n,m<=10000000,T<=10000分析bzoj 2154的加强版 再次膜拜Po神 题解代码#include<iostream>#include<cstdio>#include<cstdlib>#include<cstring>#include<algorithm>#define ll long longusing namespace std;const int MA原创 2017-01-17 22:30:40 · 309 阅读 · 0 评论 -
bzoj 4407: 于神之怒加强版 莫比乌斯反演
题意给出n,m,k求 多组数据 n,m<=5000000,T<=2000分析设d=gcd(i,j)设d=gcd(i,j)那么ans=∑nd=1dk∗f(d)那么ans=\sum_{d=1}^{n}d^k*f(d)f(d)表示1<=i<=n,1<=j<=m中gcd(i,j)=d的(i,j)对数f(d)表示1<=i<=n,1<=j<=m中gcd(i,j)=d的(i,j)对数反演后可得ans=∑n原创 2017-01-20 11:54:12 · 437 阅读 · 0 评论 -
bzoj 3309: DZY Loves Math 莫比乌斯反演
题意对于正整数n,定义f(n)为n所含质因子的最大幂指数。例如f(1960)=f(2^3 * 5^1 * 7^2)=3, f(10007)=1, f(1)=0。 给定正整数a,b,求sigma(sigma(f(gcd(i,j)))) (i=1..a, j=1..b)。 T<=10000,1<=a,b<=10^7分析%%%PoPoQQQ大佬题解本蒟蒻自己一直推到了最后一步,然后就是不会求那个前缀和原创 2017-03-20 21:22:54 · 404 阅读 · 0 评论 -
bzoj 3994: [SDOI2015]约数个数和 莫比乌斯反演
题意设d(x)为x的约数个数,给定N、M,求∑ni=1∑mj=1d(ij)\sum_{i=1}^n\sum_{j=1}^m d(ij) n,m,T<=50000分析一开始推了半天还是停留在原来的式子,于是无奈地点开了题解。。。首先模一下PoPoQQQ大佬%%% 要把这个讨厌的d(ij)给去掉的话,怎么换元或交换主题都没用,得先知道一个神奇的结论: d(nm)=∑i|n∑j|m[gcd(i,j原创 2017-03-28 14:15:44 · 439 阅读 · 0 评论 -
bzoj 4816: [Sdoi2017]数字表格 莫比乌斯反演
题意定义fibonacci数列。用f[i]f[i]表示数列的第ii项,那么 f[0]=0f[0]=0 f[1]=1f[1]=1 f[n]=f[n−1]+f[n−2],n≥2f[n]=f[n-1]+f[n-2],n\geq 2 有一个n×mn×m的表格,第ii行第jj列的格子中的数是f[gcd(i,j)]f[\gcd(i,j)],其中gcd(i,j)\gcd(i,j)表示i,ji,j的最大公约原创 2017-04-11 14:05:46 · 436 阅读 · 0 评论 -
bzoj 3930: [CQOI2015]选数 莫比乌斯反演+杜教筛
题意我们知道,从区间[L,H](L和H为整数)中选取N个整数,总共有(H-L+1)^N种方案。小z很好奇这样选出的数的最大公约数的规律,他决定对每种方案选出的N个整数都求一次最大公约数,以便进一步研究。然而他很快发现工作量太大了,于是向你寻求帮助。你的任务很简单,小z会告诉你一个整数K,你需要回答他最大公约数刚好为K的选取方案有多少个。由于方案数较大,你只需要输出其除以1000000007的余数即可原创 2017-04-17 21:35:11 · 617 阅读 · 0 评论 -
【hdu 6588 Function】【莫比乌斯反演】
题意求∑i=1ngcd(⌊i3⌋,i)\sum_{i=1}^ngcd(\lfloor\sqrt[3]i\rfloor,i)i=1∑ngcd(⌊3i⌋,i)T≤11,n≤1021T\le 11,n\le10^{21}T≤11,n≤1021分析直接推式子∑i=1ngcd(⌊i3⌋,i)\sum_{i=1}^ngcd(\lfloor\sqrt[3]i\rfloor,i)i=1∑ngcd(...原创 2019-07-28 10:30:36 · 300 阅读 · 0 评论