题目链接:http://codeforces.com/gym/101982/attachments 题目大意:有区间[a,b]和区间[c,d],求gcd(x,y)=1,其中x属于[a,b],y属于[c,d],求这样的x,y有多少对。 解题思路: 第一种反演思路: 把下界变换一下 代码: 第二种反演思 ...
分类:
其他好文 时间:
2019-04-29 12:49:31
阅读次数:
201
P3455 [POI2007]ZAP Queries 最简单的求gcd==d的个数,直接分块瞎搞。 cpp include using namespace std; define ll long long define MAXN 50000 int mu[MAXN+5]; int pri[MAXN+ ...
分类:
其他好文 时间:
2019-04-18 01:02:00
阅读次数:
195
#include #define ll long long using namespace std; const int maxn=1e7+10; int vis[maxn]; int mu[maxn]; int prime[maxn]; int tot=0; int sum1[maxn]; int... ...
分类:
其他好文 时间:
2019-04-14 12:28:06
阅读次数:
122
数论卷积: 对于两个数论函数f(x),g(x) f(n)g(n)=∑ f(d)g(n/d) d|n 莫比乌斯函数: 设一个数n=(p1^k1)*(p2^k2)*(p3^k3)*..........*(pn^kn) / 1,n=1 则μ(n)= 0,max(k1,...kn)=1 \ -1^r ,ma ...
分类:
其他好文 时间:
2019-04-12 20:46:04
阅读次数:
239
莫比乌斯函数 定义 对$d$进行质因数分解:$d=p_1^{r1}p_2^{r2}p_3^{r3}····p_k^{rk}$ $r=max\{r_1,r_2,r_3···r_k\}$ 莫比乌斯函数的定义为 $$\mu(d) = \begin{cases}1\qquad d=1\\ 0\qquad r ...
分类:
其他好文 时间:
2019-04-06 23:14:13
阅读次数:
220
Codeforces 1139D. Steps to One 题目大意: 给出$m$,一个空的数列,每次可以$rand$一个数$x\in[1,m]$放到数列的末尾,若整个数列的$gcd==1$则停止加入数 求数列的期望长度 思路: 考虑当前整个数列的$gcd$为质数$p$的倍数,则若下一个数还是$p ...
分类:
其他好文 时间:
2019-03-30 10:40:25
阅读次数:
174
题目链接: https://jzoj.net/senior/#main/show/6084 https://www.luogu.org/problemnew/show/P4916 题目: 题解: 我们设$f(x)$表示最小循环节长度为x的合法序列数,那么有$ans=\sum_{d|gcd(n,m)} ...
分类:
其他好文 时间:
2019-03-25 17:47:25
阅读次数:
266
反演套 DP 的好题(不用反演貌似也能做 Description Vivek initially has an empty array $a$ and some integer constant $m$. He performs the following algorithm: 1. Select ...
分类:
其他好文 时间:
2019-03-25 12:15:37
阅读次数:
221
几个有用的结论: 记 $(f*g)(n)=\sum\limits_{d|n} f(d) * g(\frac{n}{d})$ 则有$\mu * 1 = [n=1]$与$\phi *1 = id$ 若$F(n)=\sum\limits_{d|n} f(d)$ 则$f(n)=\sum\limits_{d| ...
分类:
其他好文 时间:
2019-03-19 10:22:34
阅读次数:
121
Description 设 $d(x)$ 为 $x$ 的约数个数,给定 $N$ 、$M$,求 $\sum\limits_{i=1}^N \sum\limits_{j=1}^M d(ij)$ Input 输入文件包含多组测试数据。 第一行,一个整数 $T$ ,表示测试数据的组数。 接下来的 $T$ 行 ...
分类:
其他好文 时间:
2019-03-15 23:18:32
阅读次数:
289