--------莫比乌斯反演
扩展的灰
扩展的灰(Extended Ash)
展开
-
Jzoj5425 数论
聪明的0v0正在学习莫比乌斯反演。她看到了这样的一道题:有n*m个人站成了一个n*m的方阵……剩下的题面,聪明的0v0不记得了。但是,她通过自己高超的数论技巧,给出了一个转化后的模型:给出n和m,求ΣΣmin(n/i,m/j)*[gcd(i,j)=1]{1聪明的0v0当然知道怎么做了,但是她想考考你。是个反演题也,可以用杜教筛也可以分块不会繁衍反演怎么办?我原创 2017-11-07 19:39:25 · 288 阅读 · 0 评论 -
UOJ62 怎样跑的更快
bi=∑j=1n(i,j)c−d∗idjdzjb_i=\sum_{j=1}^n(i,j)^{c-d}*i^dj^dz_jbi=j=1∑n(i,j)c−d∗idjdzjbiid=∑j=1nzjjd(i,j)c−d\frac{b_i}{i^d}=\sum_{j=1}^nz_jj^d(i,j)^{c-d}idbi=j=1∑nzjjd(i,j)c−d记(i,j)=x,Bi=biid,Zi=...原创 2019-01-18 21:29:40 · 453 阅读 · 0 评论 -
51Nod2026 Gcd and Lcm
题目看这里 一个非常好的题! 好的,看到题目就很懵逼 首先这个f不就是ϕϕ\phi吗,认真一看才发现不对 让后问题?f(lcm)*f(gcd)? 肯定有问题,推了一会没有结论,去看看题解: 有这么一个神奇结论f(gcd(x,y))∗f(lcm(x,y))=f(x)∗f(y)f(gcd(x,y))∗f(lcm(x,y))=f(x)∗f(y)f(gcd(x,y))*f(l...原创 2018-08-02 09:32:13 · 235 阅读 · 3 评论 -
51Nod1192 Gcd表中的质数
题目看这里 又到了推式子的时候了,莫比乌斯反演入门题 ∑i=1n∑j=1m[prime(i,j)]∑i=1n∑j=1m[prime(i,j)]\sum_{i=1}^n\sum_{j=1}^m[prime(i,j)] 首先搞一个f(d)f(d)f(d)表示有多少对(i,j)的gcd就是d 那么原式=∑nd=1[prime(d)]∗f(d)=∑d=1n[prime(d)]∗f(d)=...原创 2018-08-02 09:12:06 · 312 阅读 · 0 评论 -
51Nod1239 欧拉函数之和
题目看这里 最近才学的杜教筛,顺便就写一写模板题吧 杜教筛简易教程 直接上代码了 #include<map>#include<stdio.h>#include<string.h>#include<algorithm>#define N 1000原创 2018-08-02 08:49:19 · 291 阅读 · 0 评论 -
51Nod1594 Gcd and Phi
题目看这里 一个简单的反演题目: ∑i=1n∑j=1nϕ(gcd(ϕ(i),ϕ(j)))∑i=1n∑j=1nϕ(gcd(ϕ(i),ϕ(j)))\sum_{i=1}^n\sum_{j=1}^n\phi\bigg(gcd(\phi(i),\phi(j))\bigg) 首先做一下变换 ∑ni=1∑nj=1ϕ(gcd(ϕ(i),ϕ(j)))∑i=1n∑j=1nϕ(gcd(ϕ(i),ϕ(j)))...原创 2018-07-27 20:18:35 · 341 阅读 · 0 评论 -
[笔记]杜教筛核心原理
∑ni=1∑mj=1gcd(i,j)∑i=1n∑j=1mgcd(i,j)\sum_{i=1}^n\sum_{j=1}^mgcd(i,j) =∑nd=1d∗f(d)=∑d=1nd∗f(d)=\sum_{d=1}^nd*f(d) =∑nd=1d∗∑d|jμ(j/d)∗F(d)=∑d=1nd∗∑d|jμ(j/d)∗F(d)=\sum_{d=1}^nd*\sum_{d|j}\mu(j/d)*F(d) ...原创 2018-07-29 18:43:34 · 404 阅读 · 0 评论 -
51Nod1190 最小公倍数之和 V2
题目看这里繁衍反演真好玩来看看这个题的式子求Σlcm(i,b) (a<=i<=b)首先不难发现以下这样的变换让后做一下差分就得到答案,注意这题比较卡时间,需要比较好的优化#pragma GCC opitmize("O3")#pragma G++ opitmize("O3")#include<stdio.h>#include<string.h>#inclu...原创 2018-05-15 22:08:58 · 388 阅读 · 0 评论 -
一些基础算法的模板(持续更新)
更新中//Templates From Extended_Ash/Cooevjnz/JacaJava/Tubbcrafft//To be continued...//Suffix Automationchar str[N]; int s[N][26],mx[N],f[N],sz[N]; int last=1,cnt=1,n,v[N],r[N],ans=0; inlin原创 2017-10-25 22:05:58 · 3734 阅读 · 0 评论 -
一些反演求和过程?
.问题∑i,j&lt;=niaja(i,j)b\sum_{i,j&lt;=n}i^aj^a(i,j)^bi,j<=n∑iaja(i,j)b记f(d)=∑(i,j)=diajaF(d)=∑d∣i,d∣jiajaf(d)=\sum_{(i,j)=d}i^aj^a \\ F(d)=\sum_{d|i,d|j}i^aj^af(d)=(i,j)=d∑iajaF(d)=d∣i,d∣...原创 2019-01-18 22:38:30 · 421 阅读 · 1 评论