![](https://img-blog.csdnimg.cn/20201014180756928.png?x-oss-process=image/resize,m_fixed,h_64,w_64)
莫比乌斯反演
romiqi_new
这个作者很懒,什么都没留下…
展开
-
191106CSP(NOI?)模拟及NOI(CSP?)模拟题解
CSP模拟:T1:求∑i=1n∑j=1mCgcd(i,j)B%mod\sum_{i=1}^n\sum_{j=1}^mC_{gcd(i,j)}^B \%mod∑i=1n∑j=1mCgcd(i,j)B%modn,m≤1e10,B≤mod=9990017n,m\le1e10,B\le mod=9990017n,m≤1e10,B≤mod=9990017莫比乌斯反演,设f(i)f(i)f(i)表...原创 2019-11-06 19:04:24 · 503 阅读 · 2 评论 -
[BZOJ3994][莫比乌斯反演][SDOI2015]约数个数和
BZOJ3994------------------------------------这里是留给证明的地方-------------------------------------------------------------------------这里是留给证明的地方-------------------------------------预处理莫比乌斯函数前缀和 和 约数个数和的前缀...原创 2019-02-13 10:21:46 · 164 阅读 · 0 评论 -
[莫比乌斯反演]YY的GCD[证明已补]
BZOJ2820莫比乌斯反演模板Code:#include<bits/stdc++.h>#define int long longusing namespace std;inline int read(){ int res=0,f=1;char ch=getchar(); while(!isdigit(ch)) {if(ch=='-') f=-f;ch=getchar(...原创 2019-01-19 16:38:31 · 194 阅读 · 0 评论 -
[莫比乌斯反演][BZOJ2301]problem b
BZOJ2301先把询问拆成四个部分然后分别询问(别告诉我你不会容斥)然后问题就是求∑i=1n∑j=1m[gcd(i,j)=d]\sum_{i=1}^n\sum_{j=1}^m[gcd(i,j)=d]i=1∑nj=1∑m[gcd(i,j)=d]=∑i=1⌊nd⌋∑j=1⌊md⌋[gcd(i,j)=1]=\sum_{i=1}^{\lfloor \frac{n}{d} \rfloor} \s...原创 2019-02-19 09:20:05 · 110 阅读 · 0 评论 -
[BZOJ2440][二分][数学][中山市选2011]完全平方数
BZOJ2440一句话:求第k个不含完全平方因子的数(k<=1e9)首先可以二分,然后问题就转化为求区间[1,mid]中有多少个不含完全平方因子的数考虑容斥原理,1-mid中不含完全平方因子的数的个数=mid-有 由一个质数平方构成的完全平方数(如4,9) 这个因子的数 的个数+有 由两个不同质数相乘的数 平方 构成的完全平方数(如36) 这个因子的数 的个数…所以可以发现ans=...原创 2019-02-13 11:40:57 · 147 阅读 · 0 评论 -
[BZOJ4816][莫比乌斯反演]SDOI2017:数字表格
BZOJ4816∏i=1n∏j=1mfib[gcd(i,j)]\prod_{i=1}^n\prod_{j=1}^m{fib[gcd(i,j)]}i=1∏nj=1∏mfib[gcd(i,j)]∏k=1min(n,m)∏i=1n∏j=1mfib[gcd(i,j)]\prod_{k=1}^{min(n,m)}\prod_{i=1}^n\prod_{j=1}^m{fib[gcd(i,j)]}k=1...原创 2019-08-09 19:09:31 · 191 阅读 · 0 评论 -
[BZOJ3529][莫比乌斯反演][树状数组]SDOI2014:数表
BZOJ3529要求的是约数和如果没有a的限制,直接反演就完了有a的限制的话,把询问按a排序,用树状数组维护约数和乘上μ\muμ的前缀和即可Code:#include<bits/stdc++.h>#define mod 2147483648#define ll long longusing namespace std;inline int read(){ int r...原创 2019-08-09 19:21:31 · 159 阅读 · 0 评论