![](https://img-blog.csdnimg.cn/20201014180756926.png?x-oss-process=image/resize,m_fixed,h_64,w_64)
莫比乌斯反演
zsyz_lb2003
这个作者很懒,什么都没留下…
展开
-
YY的GCD[BZOJ2820]
欢迎大家访问我的老师的OJ———caioj.cn题面描述因为没钱,所以思路由Zap的结论我们可以知道我们这题应该求的是∑isprime(p)∑i=1min(⌊n/p⌋,⌊m/p⌋)μ(i)∗⌊⌊n/p⌋/i⌋∗⌊⌊m/p⌋/i⌋\large\sum_{isprime(p)}\sum_{i=1}^{\min(\left\lfloor n/p\right\rfloor,\left\lflo...原创 2019-05-03 10:45:25 · 84 阅读 · 0 评论 -
[中山市选2011]完全平方数[BZOJ2440]
题面描述传送门思路很明显地,一个数不属于完全平方数或完全平方数的倍数,它的μ\muμ一定不等于0.因此我们只要二分枚举xxx,计算[1,x][1,x][1,x]之中的非完全平方数根据容斥原理:xxx以内非完全平凡数的个数为:xxx-[1,x][1,x][1,x]中一个质数(2,3,5,...2,3,5,...2,3,5,...)的平方的倍数的数量+[1,x][1,x][1,x]中两个...原创 2019-04-30 12:51:57 · 168 阅读 · 0 评论 -
GCD4[caioj1283]\[Zoj3435]
欢迎大家访问我的老师的OJ———caioj.cn题面描述传送门思路算是仪仗队的进阶版吧。因为这题的左下角为(1,1,1)(1,1,1)(1,1,1),我们不妨将其变为(0,0,0)(0,0,0)(0,0,0),右上角为(a−1,b−1,c−1)(a-1,b-1,c-1)(a−1,b−1,c−1)方便讨论。首先我们可以观察到,三条棱上的(1,0,0),(0,1,0),(0,0,1)(1,...原创 2019-04-30 10:57:06 · 146 阅读 · 0 评论 -
仪仗队[SDOI2008]\[LuoguP2158]
题面描述传送门思路这道题和可见点数没有差别。欧拉函数做法就不写了。主要讲讲莫反做法吧。回顾一下Zap(这真的是一道模板题啊!)设D(a,b,k)D(a,b,k)D(a,b,k)表示满足x≤a,y≤bx\le a,y\le bx≤a,y≤b且k∣gcd(x,y)k \mid \gcd(x,y)k∣gcd(x,y)的二元组有多少对。由于D(a,b,k)=⌊a/k⌋∗⌊b/k⌋D(a,...原创 2019-04-30 10:38:10 · 224 阅读 · 0 评论 -
GCD3[caioj1282]\[BZOJ2301]莫反
题面描述caiojbzoj注意,caioj题面略有不同思路与Zap很像。设D(a,b,k)D(a,b,k)D(a,b,k)表示满足x≤a,y≤bx\le a,y\le bx≤a,y≤b且k∣gcd(x,y)k \mid \gcd(x,y)k∣gcd(x,y)的二元组有多少对。由于D(a,b,k)=⌊a/k⌋∗⌊b/k⌋D(a,b,k)=\left\lfloor a/k\right\...原创 2019-04-30 09:08:49 · 187 阅读 · 0 评论 -
GCD2[caioj1281]\[BZOJ2818]莫反
欢迎大家访问我的老师的OJ———caioj.cn题面描述caiojBZOJ思路这是一道不用莫函的莫反(伪)题啊。我们观察一下,1≤x,y≤n,gcd(x,y)=p,p是质数1\le x,y\le n,\gcd(x,y)=p,p是质数1≤x,y≤n,gcd(x,y)=p,p是质数可以转化为:1≤x,y≤n/p,gcd(x,y)=1,也就是x,y互质1\le x,y\le n/p...原创 2019-04-30 07:55:25 · 129 阅读 · 0 评论 -
GCD[caioj1280 莫反]
欢迎大家访问我的老师的OJ———caioj.cn前情提要建议先去看看Zap,数据加强了,但思路较为简单。题面描述传送门思路与Zap不同的就是多了注意:2,3和3,2是一种情况那么我们就要去重了。设D(a,b,k)D(a,b,k)D(a,b,k)表示满足x≤a,y≤bx\le a,y\le bx≤a,y≤b且k∣gcd(x,y)k \mid \gcd(x,y)k∣gcd(x,y)...原创 2019-04-29 20:07:28 · 557 阅读 · 0 评论 -
Calculating[LuoguP3935]
题面描述传送门思路根据算数基本定理,正约数个数为:(c1+1)∗(c2+1)∗......∗(cm+1)=∏i=1m(ci+1)(c_1+1)*(c_2+1)*......*(c_m+1)={\large\prod\limits_{i=1}^{m}}(c_i+1)(c1+1)∗(c2+1)∗......∗(cm+1)=i=1∏m(ci+1)则f(i)f(i)f(i)实际上就是ii...原创 2019-04-29 18:22:56 · 174 阅读 · 0 评论 -
Crash的数字表格[BZOJ2154]
题面描述传送门思路快去%PoPoQQQAC code#include<cstdio>#include<cmath>#include<cstdlib>#include<cstring>#include<algorithm>#define ll long long#define gc getchar()using na...原创 2019-05-05 07:46:40 · 82 阅读 · 0 评论 -
Zap[BZOJ1101]\[POI2007]
题面描述传送门思考题目要求等价于求有多少二元组(x,y)(x,y)(x,y)满足x≤a/k,y≤b/kx\le a/k,y\le b/kx≤a/k,y≤b/k并且x,yx,yx,y互质(因为gcd(x,y)=k\operatorname{gcd}(x,y)=kgcd(x,y)=k)。设D(a,b,k)D(a,b,k)D(a,b,k)表示满足x≤a,y≤bx\le a,y\le bx≤a,...原创 2019-04-29 11:55:13 · 179 阅读 · 0 评论 -
莫比乌斯函数
设正整数N按照算数基本定理分解质因数为N=p1c1p2c2⋯pmcmN=p_1^{c_1}p_2^{c_2}\cdots p_m^{c_m}N=p1c1p2c2⋯pmcm,定义函数μ(N)={0∃i∈[1,m],ci>11m≡0(mod2),∀i∈[1,m],ci=1−1m≡1(mod2),∀i∈[1,m],ci=1\mu(N)=\begin{cases}0 &am...原创 2019-04-29 12:55:26 · 466 阅读 · 1 评论 -
数表[BZOJ3529]\[SDOI2014]
题面描述传送门思路先%一波dalao——PoPoQQQ话说这道题真的毒瘤啊!一句话题意,设f(i)f(i)f(i)为iii的约数和,求ans=∑i=1n∑j=1m并且f(gcd(i,j))<=af(gcd(i,j))mod231\large ans=\sum_{i=1}^n\sum_{j=1}^{m并且f(\gcd(i,j))<=a}f(\gcd(i,j...原创 2019-05-03 15:54:30 · 133 阅读 · 0 评论 -
[NOI2010]能量采集
题目描述传送门思路暂时咕咕咕了。AC code#include<cstdio>#include<algorithm>#include<cstdlib>#include<cstring>#include<cmath>#define ll long longusing namespace std;const int N...原创 2019-06-27 08:55:18 · 203 阅读 · 0 评论