![](https://img-blog.csdnimg.cn/20190927151053287.png?x-oss-process=image/resize,m_fixed,h_224,w_224)
莫比乌斯反演,杜教筛,迪利克雷卷积
文章平均质量分 70
qq_43520313
这个作者很懒,什么都没留下…
展开
-
简单的函数
题目https://loj.ac/p/6053求∑i=1nf(i) mod 109+7\sum_{i=1}^{n}f(i)\ mod\ 10^9+7i=1∑nf(i) mod 109+7f(1)=1f(1)=1f(1)=1f(pc)=p⊕cf(p^c)=p\oplus cf(pc)=p⊕c,ppp为质数f(ab)=f(a)f(b)f(ab)=f(a)f(b)f(ab)=f(a)f(b),(a,b)=1(a,b)=1(a,b)=1思路Min25Mi原创 2021-03-24 19:22:33 · 64 阅读 · 0 评论 -
Sanrd
题目https://uoj.ac/problem/188求∑i=lrf(i)\sum_{i=l}^{r}f(i)i=l∑rf(i)f(i)f(i)f(i)是iii的次大质因子,如果是素数或111,则为000。l+r≤1011l+r\le 10^{11}l+r≤1011思路考虑Min25Min25Min25筛时的SSS函数的dpdpdp,素数部分贡献是000,则只有合数的贡献S(n,j)=∑i=j+1π(n)∑e=1pie≤n(f(pie)[e>1]+∑k=2⌊npie⌋f(piek原创 2021-03-23 20:01:14 · 108 阅读 · 0 评论 -
Min_25筛
题目一般要求∑i=1nF(i)n≤1010\sum_{i=1}^{n}F(i)\quad n\le 10^{10}i=1∑nF(i)n≤1010其中F(x)F(x)F(x)是积性函数。Min25Min25Min25筛能用的前提:质数处的f(p)f(p)f(p)值是关于ppp的低阶多项式,质数次方处的f(pe)f(p^e)f(pe)值可以快速计算。约定pip_ipi表示第iii个素数,下标从111开始P\mathbb{P}P表示素数的集合π(n)\pi(n)π(n)表示1∼n1\sim原创 2021-03-21 21:48:23 · 79 阅读 · 0 评论 -
Min_25筛模板
题目https://www.luogu.com.cn/problem/P5325定义积性函数f(x)f(x)f(x),且f(pk)=pk(pk−1)(p为质数)f(p^k)=p^k(p^k-1)(p为质数)f(pk)=pk(pk−1)(p为质数),求∑i=1nf(i)\sum_{i=1}^{n}f(i)i=1∑nf(i)对109+710^9+7109+7取模。n≤1010n\le 10^{10}n≤1010。思路把f(p)=p2−pf(p)=p^2-pf(p)=p2−p拆成两个完全积性函原创 2021-03-18 20:41:14 · 83 阅读 · 0 评论 -
[SDOI2008]仪仗队
题目:https://ac.nowcoder.com/acm/problem/20313作为体育委员,C君负责这次运动会仪仗队的训练。仪仗队是由学生组成的N * N的方阵,为了保证队伍在行进中整齐划一,C君会跟在仪仗队的左后方,根据其视线所及的学生人数来判断队伍是否整齐(如下图)。思路:对于一个点(x,y)(x≥1,y≥1)(x,y)(x\ge 1,y\ge 1)(x,y)(x≥1,y≥1),如果gcd(x,y)≠1gcd(x,y)\neq 1gcd(x,y)=1,则会被点(xd,yd)(d原创 2020-08-26 00:42:41 · 133 阅读 · 0 评论 -
特殊情况迪利克雷卷积的化简
题目:证明:若f,gf,gf,g为积性函数,n=∏i=1mpiαin=\prod_{i=1}^{m}p_i^{\alpha_i}n=∏i=1mpiαi,有f∗g (n)=∏i=1mf∗g (piαi)f*g\,(n)=\prod\limits_{i=1}^{m}f*g\,(p_i^{\alpha_i})f∗g(n)=i=1∏mf∗g(piαi)思路:f∗g (n)=∑d∣nf(d)⋅g(nd)=∑a1=0α1∑a2=0α2...∑am=0αmf(∏i=1mpiai)⋅g(∏i=1mp原创 2020-08-16 00:34:55 · 173 阅读 · 0 评论 -
互素的数的和
题目:求∑i=1ni[(i,n)=1]\sum_{i=1}^{n}i[(i,n)=1]∑i=1ni[(i,n)=1]思路:∑i=1ni[(i,n)=1]=∑i=1n∑d∣(i,n)μ(d)⋅i=∑d∣n∑i=1ndμ(d)⋅i⋅d=∑d∣nμ(d)⋅d⋅(1+nd)⋅nd2=∑d∣nμ(d)⋅n2+μ(d)⋅n22d=n2[n=1]+n2⋅(μ∗id)=n2[n=1]+n2⋅φ(n)\begin{aligned}&\sum_{i=1}^{n}i[(i,n)=1]\\=&\sum原创 2020-08-15 23:40:51 · 253 阅读 · 0 评论 -
A Very Easy Math Problem
题目:http://acm.hdu.edu.cn/showproblem.php?pid=6833Given you nnn,xxx,kkk , find the value of the following formula:∑a1=1n...∑ax=1n(∏j=1xajk)f(gcd(a1,a2...ax))⋅gcd(a1,a2...ax)\sum_{a_1=1}^{n}...\sum_{a_x=1}^{n}\left(\prod_{j=1}^{x}a_j^k\right)f(gcd(a_1,a_原创 2020-08-08 11:52:57 · 150 阅读 · 0 评论 -
LCM Sum
题目:求∑i=1nlcm(i,n)n≤106\sum_{i=1}^{n}lcm(i,n)\quad n\le 10^6∑i=1nlcm(i,n)n≤106思路:∑i=1nlcm(i,n)=∑i=1ni∗n(i,n)=∑d=1n∑i=1ni∗nd[(i,n)=d]=∑d∣n∑k=1ndk∗n[(k,nd)=1]另i=k∗d=n∑d∣nφ(nd)∗nd2∑i=1ni[(i,n)=1]=φ(n)∗n2=n∑d∣nφ(d)∗d2\begin{aligned}&\sum_{i=1}^{n}lcm(原创 2020-07-24 22:33:44 · 163 阅读 · 0 评论 -
牛牛的最大公约数
题目:https://ac.nowcoder.com/acm/problem/21546牛牛有一个区间[L,R][L,R][L,R],需要选择NNN个数,这NNN个数都在这个区间范围内,那么我们知道一共有(R−L+1)N(R - L + 1) ^ N(R−L+1)N种选法,假如我们想要这NNN个数的最大公约数恰好是KKK.请问一共有多少种选法,输出答案对109+710^9+7109+7取模.N,K,L,R(1≤N,K,L≤109,L≤R≤109)0≤R−L≤105N,K,L,R (1\le N, K,原创 2020-07-24 22:34:56 · 348 阅读 · 0 评论 -
迪利克雷卷积,莫比乌斯反演,杜教筛
常见积性函数μ(n)={1,i=10,n无2次方以上的素因数(−1)k,k为不同素因数的种类\mu(n)=\begin{cases}1,&i=1\\0,&n无2次方以上的素因数\\(-1)^k,&k为不同素因数的种类\end{cases}μ(n)=⎩⎪⎨⎪⎧1,0,(−1)k,i=1n无2次方以上的素因数k为不同素因数的种类φ(n)=∑i=1n[(n,i)=1]\varphi(n)=\sum_{i=1}^{n}[(n,i)=1]φ(n)=∑i=1n[(n,i)=原创 2020-10-15 23:17:54 · 152 阅读 · 0 评论