莫比乌斯反演
文章平均质量分 51
Icefox_zhx
这个作者很懒,什么都没留下…
展开
-
CIA3 Sometimes Naive(莫比乌斯反演+积性函数前缀和)
题意:求∑x=1n∑y=1mφ(gcd(x,y))\sum\limits_{x=1}^n\sum\limits_{y=1}^m\varphi(gcd(x,y)) 根据各种老套路化简得 ∑i=1n⌊n/i⌋⌊m/i⌋∑d|iφ(i/d)∗μ(d)\sum\limits_{i=1}^n\lfloor n/i\rfloor \lfloor m/i\rfloor\sum\limits_{d|i}\varp原创 2017-10-18 21:23:54 · 870 阅读 · 0 评论 -
bzoj5332/loj2565 [Sdoi2018]旧试题(Mobius反演+图论)
求∑x=1A∑y=1B∑z=1Cd(xyz)∑x=1A∑y=1B∑z=1Cd(xyz)\sum\limits_{x=1}^A\sum\limits_{y=1}^B\sum\limits_{z=1}^Cd(xyz) 类似这道题的结论,我们有 d(xyz)=∑i|x∑j|y∑k|z[(i,j)==1][(j,k)==1][(i,k)==1]d(xyz)=∑i|x∑j|y∑k|z[(i,j)==1]...原创 2018-05-24 15:35:30 · 448 阅读 · 0 评论 -
bzoj3561 DZY Loves Math VI(Mobius反演)
枚举gcd反演一下,化简得到答案: Ans=∑d=1ndd∑t=1⌊nd⌋μ(t)t2d∑i=1⌊ndt⌋∑j=1⌊mdt⌋idjdAns=∑d=1ndd∑t=1⌊nd⌋μ(t)t2d∑i=1⌊ndt⌋∑j=1⌊mdt⌋idjdAns=\sum\limits_{d=1}^nd^d\sum\limits_{t=1}^{\lfloor\frac{n}{d}\rfloor}\mu(t)t^{2d}\s...原创 2018-05-22 22:06:05 · 120 阅读 · 0 评论 -
bzoj3994 [SDOI2015]约数个数和(Mobius反演)
求∑ni=1∑mj=1σ0(ij)∑i=1n∑j=1mσ0(ij)\sum_{i=1}^n\sum_{j=1}^m\sigma_0(ij) 我们有σ0(xy)=∑d1|x∑d2|y[gcd(d1,d2)=1]σ0(xy)=∑d1|x∑d2|y[gcd(d1,d2)=1]\sigma_0(xy)=\sum_{d_1|x}\sum_{d_2|y}[gcd(d_1,d_2)=1] 然后化式子就好了。...原创 2018-05-22 21:17:29 · 378 阅读 · 0 评论 -
bzoj4652 [Noi2016]循环之美(Mobius反演+杜教筛+Hash表)
挖了很久的坑qaq 还记得第一次xtx神犇给我们讲这道题,已经过去快整整一年了呢owo首先怎么算是纯循环小数呢? 因为要求数值不同的数对个数,因此我们只考虑gcd(x,y)==1的数对(x,y)。 什么时候会出现循环呢?当出现相同余数时。那纯循环呢?就要求第一次出现的相同余数为x%y。即只需要满足xkl≡xmody,l>0xkl≡xmody,l>0xk^l \equiv x ...原创 2018-05-15 11:32:21 · 236 阅读 · 0 评论 -
bzoj3944 Sum(杜教筛)
杜教筛裸题。。。 mu函数的前缀和怎么搞呢,我们有μ∗1=ϵ\mu*1=\epsilon 令g=1就好了。原创 2018-05-07 21:09:12 · 238 阅读 · 0 评论 -
bzoj2005 [Noi2010]能量采集(Mobius反演,欧拉函数)
我们发现点(x,y)到原点的连线上有gcd(x,y)-1个点(不含始终点),所以我们要求的答案就是∑x=1n∑y=1m2∗(gcd(x,y)−1)+1\sum\limits_{x=1}^n\sum\limits_{y=1}^m2*(gcd(x,y)-1)+1也就是2∑x=1n∑y=1mgcd(x,y)−n∗m2\sum\limits_{x=1}^n\sum\limits_{y=1}^mgcd(x,y原创 2018-03-30 08:57:12 · 257 阅读 · 0 评论 -
bzoj3434 [Wc2014]时空穿梭(莫比乌斯函数+多项式乘法)
蒟蒻我第一步共线就没想出来x暴力dp T飞。直接爆零qaq 我们记C个点分别为CiC_i,V为{v1,v2...vn}\{v_1,v_2...v_n\},令Bi=Ci−Ci−1B_i=C_i-C_{i-1}。则我们发现Bi=(ti−ti−1)VB_i=(t_i-t_{i-1})V,也就是说每一个B的各维坐标都是成比例的。我们设第i维的坐标极差为Δxi\Delta x_i(Cn−C1C_n-C_1)原创 2018-01-29 17:53:02 · 481 阅读 · 0 评论 -
bzoj3529 [Sdoi2014]数表(莫比乌斯反演+树状数组)
%%%popoqqq大爷的题解:传送门 又掌握了线性素数筛处理sigma函数的技能。。。复杂度O(T*logT+T*sqrt(n)*logn)原创 2017-10-13 22:14:36 · 367 阅读 · 0 评论 -
bzoj2440 [中山市选2011]完全平方数(二分答案+莫比乌斯反演)
类似vijos1889。原创 2017-10-13 11:43:26 · 403 阅读 · 0 评论 -
bzoj2301 [HAOI2011]Problem b(求gcd==k的个数)(莫比乌斯反演+容斥原理)
首先我们搞掉下界,怎么搞呢,用容斥原理即可。(看做矩形区间),然后我们需要求∑x=1n∑y=1ngcd(x,y)==k\sum\limits_{x=1}^n\sum\limits_{y=1}^ngcd(x,y)==k。 ∑x=1⌊n/k⌋∑y=1⌊m/k⌋gcd(x,y)==1\sum\limits_{x=1}^{\lfloor{n/k}\rfloor}\sum\limits_{y=1}^{\lf原创 2017-10-10 14:39:41 · 1371 阅读 · 1 评论 -
vijos1889 天真的因数分解(莫比乌斯反演+二分答案)
求第k个mu为0的数,k<=1e10,所以O(n)肯定也挂。考虑不合格的数,就是含有一个平方因子的数。设F(x)表示1~x中不合格的数,则根据容斥原理,利用莫比乌斯函数,我们可以得到F[x]=−∑i=2x√⌊x/i2]∗μ(i)F[x]=-\sum\limits_{i=2}^{\sqrt x}\lfloor x/{i^2}]*\mu(i),然后就是喜闻乐见的分块啦~(直接暴力算也是可以的。。。)复杂原创 2017-10-11 16:02:34 · 441 阅读 · 0 评论 -
bzoj2820 [bzoj2820]YY的GCD(线性素数筛+莫比乌斯反演)
求gcd(i,j)为质数的个数。即 ∑p∑i=1n∑j=1mgcd(i,j)==p\sum\limits_p\sum\limits_{i=1}^n\sum\limits_{j=1}^mgcd(i,j)==p 首先我们可以像bzoj2301一样把后一部分化简,得到(假定n< m) ∑p∑d=1⌊n/p⌋μ(d)⌊n/pd⌋⌊m/pd⌋\sum\limits_p\sum\limits_{d=1}^原创 2017-10-10 21:58:56 · 567 阅读 · 0 评论 -
bzoj2154 Crash的数字表格(莫比乌斯反演)
求lcm(i,j)的和. 这个题解写的挺详细的呢:传送门。 算答案时分块套路复杂度O(sqrt(n)),算f[x][y]时利用了莫比乌斯反演,分块计算,复杂度也是O(sqrt(n)),所以总的复杂度是O(n) tips:注意因为取模带来的负数问题。。。原创 2017-10-10 22:49:52 · 326 阅读 · 0 评论 -
bzoj1101 [POI2007]ZAP-Queries(莫比乌斯反演)
简单化的bzoj2301原创 2017-10-10 14:46:37 · 386 阅读 · 0 评论 -
bzoj2986 Non-Squarefree Numbers(Mobius函数+容斥+二分答案)
其实和bzoj2440一模一样。 就是利用Mobius函数来容斥就好了#include <bits/stdc++.h>using namespace std;#define ll long long#define inf 0x3f3f3f3f#define N 500010inline char gc(){ static char buf[1<<16...原创 2018-06-22 16:57:52 · 212 阅读 · 0 评论