![](https://img-blog.csdnimg.cn/20201014180756738.png?x-oss-process=image/resize,m_fixed,h_64,w_64)
数学
文章平均质量分 73
数学
是羊驼鸭
这个作者很懒,什么都没留下…
展开
-
CodeForces - 615D Multipliers (欧拉降幂)
题目链接:点击这里题目大意:给定一个正整数 nnn ,求其所有约数的乘积,即求 ∏d∣ndmod 1e9+7\prod\limits_{d|n}d \mod 1e9+7d∣n∏dmod1e9+7题目分析:将 nnn 用唯一分解定理展开 n=∏i=1xpiain=\prod_{i=1}^xp_i^{a_i}n=∏i=1xpiai ,逐个质数考虑贡献:对于 piaip_i^{a_i}piai 来说,自己有 1+2+...+ai=ai⋅(ai+1)21+2+...+a_i=\frac{a_原创 2022-03-21 22:06:29 · 268 阅读 · 0 评论 -
CodeForces - 1646E Power Board (思维,数学)
题目链接:点击这里题目大意:给定一个 n×mn\times mn×m 的矩阵,其中 a[i][j]=ija[i][j]=i^ja[i][j]=ij ,求这个 n×mn\times mn×m 的矩阵有多少个不同的元素一个 3×33\times 33×3 的矩阵如下图所示:题目分析:不难想到,如果两行有相同的元素,那么两行中必然分别存在一个 aba^bab 和 aca^cac例如第 2,4,82,4,82,4,8 行其元素分别为:21,22,...,2m2^1,2^2,...,2^m21,22,原创 2022-03-05 13:24:52 · 637 阅读 · 1 评论 -
P4931 [MtOI2018]情侣?给我烧了(加强版) (错位排列)
题目链接:点击这里题目大意:一个电影院有 nnn 排座椅,每排座椅有两个座位,现在有 nnn 对情侣(每对情侣是一男一女两人),求有 kkk 对情侣坐在一起的方案数题目分析:有 kkk 对情侣坐在一起,就说明有 n−kn-kn−k 对情侣没有坐在一起,我们可以先将 kkk 对情侣安排入座,然后再对剩下的 n−kn-kn−k 对情侣进行一个错位排列首先考虑将 kkk 对情侣安排入座,其贡献为:Cnk⋅Ank⋅2kC_n^k·A_n^k·2^kCnk⋅Ank⋅2k其意义为先在 nnn 排座椅中原创 2022-02-16 11:01:16 · 622 阅读 · 0 评论 -
P1891 疯狂 LCM(欧拉函数)
题目链接:点击这里题目大意:求 ∑i=1nlcm(i,n)\sum_{i=1}^nlcm(i,n)∑i=1nlcm(i,n)题目分析:思路同 月月给华华出题我们知道 lcm(i,n)=i⋅ngcd(i,n)lcm(i,n)=\frac{i·n}{\gcd(i,n)}lcm(i,n)=gcd(i,n)i⋅n所以题目所求就转化成了∑i=1ni⋅ngcd(i,n)=n⋅∑i=1nigcd(i,n)\sum_{i=1}^n\frac{i·n}{\gcd(i,n)}=n·\sum_{i=1}^n原创 2022-02-14 17:38:35 · 280 阅读 · 0 评论 -
P7884 【模板】Meissel–Lehmer 算法(亚线性求素数个数)
题目链接:点击这里题目大意:给定整数 nnn ,求出 π(n)\pi(n)π(n) 的值。其中 π(n)\pi(n)π(n) 表示 1∼n1 \sim n1∼n 的整数中质数的个数。题目分析:min25min25min25 在求素数个数上被此算法完爆(此题 min25min25min25 求 101310^{13}1013 的数据会 ttt 飞)Meissel–Lehmer\text{Meissel–Lehmer}Meissel–Lehmer 算法时间复杂度 101310^131013 以内为原创 2022-02-13 14:26:27 · 639 阅读 · 0 评论 -
POJ - 3090 Visible Lattice Points (数学,四种做法)
题目链接:点击这里题目大意:题目分析:本题需要一个重要的转化:仔细观察上图,我们会发现:当一个点的横坐标和纵坐标的 gcd\gcdgcd 大于 111 时就会被遮挡,即 gcd(x,y)≠1\gcd(x,y)≠1gcd(x,y)=1 的点会被遮住所以题目所求就是:2+∑i=1n∑j=1n[gcd(i,j)=1]2+\sum_{i=1}^n\sum_{j=1}^n[\gcd(i,j)=1]2+∑i=1n∑j=1n[gcd(i,j)=1](其中2是x=0 || y=0的贡献)所以便有原创 2022-02-11 20:03:26 · 592 阅读 · 3 评论 -
POJ - 2480 Longge‘s problem (欧拉函数)
题目链接:点击这里题目大意:给定正整数 n(1<n<231−1)n(1<n<2^{31}-1)n(1<n<231−1) ,求:∑i=1ngcd(i,n)\sum_{i=1}^n\gcd(i,n)i=1∑ngcd(i,n)题目分析:∑i=1ngcd(i,n)\sum_{i=1}^n\gcd(i,n)i=1∑ngcd(i,n)=∑d∣n∑i=1n[gcd(i,n)=d]⋅d=\sum_{d|n}\sum_{i=1}^n[\gcd(i,n)=d]·d=d∣原创 2022-02-11 15:13:03 · 376 阅读 · 0 评论 -
2022牛客寒假第一场 B站与各唱各的 (期望)
题目链接:点击这里题目大意:有 nnn 个人和一个长度为 mmm 的序列。每个人对序列的每个元素都可以选或不选,如果某个元素被 000 或 nnn 个人选中就算该元素被删除掉,在人们不互相沟通且做出最优选择的前提下,求该序列剩余元素个数的期望题目分析:我们可以发现,对于序列的不同位置,是相互独立的,因此期望就是单个元素的期望乘上 mmm 。因为每个都是在不能互相交流的情况下做出最优选择,所以他们对于同一个元素的选取概率是相同的不妨设第 iii 个人选该元素的概率为 ppp ,因此该元素被删掉的期原创 2022-01-24 20:06:16 · 279 阅读 · 3 评论 -
upc 数字游戏 (质因子分解)
题目大意:Alice的父亲是一个伟大的数学家,他很喜欢和Alice一起玩数学游戏,这次他写下一系列的数,告诉Alice可以进行以下操作:选择序列中的任意两个数A和B,再选择一个能整除A的素数X,然后用A/X替代A,用B*X替代B。上述操作可以进行任意次,最终得分为数列中所有数的最大公约数。请你帮助Alice获得最大得分。题目分析:不难发现这个操作就是把质因子的指数进行了均分,所以可以把这 nnn 个数都质因子分解一下,然后再遍历一遍存指数的桶,把值均分给这 nnn 个数即可具体细节见代码:原创 2022-01-21 23:05:30 · 2278 阅读 · 0 评论 -
SP34096 DIVCNTK - Counting Divisors (general)(min25)
题目链接:点击这里题目大意:题目分析:好久以前的代码,发现一直 waiting\text{waiting}waiting 就 rejudge\text{rejudge}rejudge 了一下,然后 ac\text{ac}ac 了,特此记录一下具体细节见代码://#pragma GCC optimize(2)//#pragma GCC optimize("Ofast","inline","-ffast-math")//#pragma GCC target("avx,sse2,sse3,sse4原创 2021-10-20 16:18:46 · 90 阅读 · 0 评论 -
P6055 [RC-02] GCD(莫比乌斯反演+杜教筛)
题目链接:点击这里题目大意:给定正整数 nnn 求:∑i=1n∑j=1n∑p=1⌊nj⌋∑q=1⌊nj⌋[gcd(i,j)=1][gcd(p,q)=1]\sum_{i=1}^n\sum_{j=1}^n\sum_{p=1}^{\lfloor \frac{n}{j} \rfloor}\sum_{q=1}^{\lfloor \frac{n}{j} \rfloor}[\gcd(i,j)=1][\gcd(p,q)=1]i=1∑nj=1∑np=1∑⌊jn⌋q=1∑⌊jn⌋[gcd(i,j)=1][原创 2021-09-18 23:25:21 · 99 阅读 · 0 评论 -
CodeForces 348D Turtles(LGV引理)
题目链接:点击这里题目大意:一张 nnn 行 mmm 列的网格图,图中的有些格子上面有障碍物,但保证 (1,1)(1,1)(1,1) 和 (n,m)(n,m)(n,m) 上面都没有障碍物。在 (1,1)(1,1)(1,1) 处有两只乌龟,都想要去(n,m)(n,m)(n,m) 。乌龟每次都可以向下或者向右走一格,前提是格子上没有任何障碍物。要求两只乌龟在前往 (n,m)(n,m)(n,m) 的路途中不可以相遇,即除了起点和终点,他们的路径没有其他公共点。求出从起点到终点的不同路径对数。答案对 109+7原创 2021-09-17 18:27:39 · 545 阅读 · 2 评论 -
P6657 【模板】LGV 引理
题目链接:点击这里题目大意:有一个 n×nn\times nn×n 的棋盘,左下角为 (1,1)(1,1)(1,1),右上角为 (n,n)(n,n)(n,n),若一个棋子在点 (x,y)(x,y)(x,y),那么走一步只能走到 (x+1,y)(x+1,y)(x+1,y) 或 (x,y+1)(x,y+1)(x,y+1)。现在有 mmm 个棋子,第 iii 个棋子一开始放在 (ai,1)(a_i,1)(ai,1),最终要走到 (bi,n)(b_i,n)(bi,n) 。问有多少种方案,使得每个棋子都能从原创 2021-09-17 17:15:48 · 254 阅读 · 0 评论 -
P7112 【模板】行列式求值
题目链接:点击这里题目大意:给定一个 nnn 阶行列式 AAA,求 ∣A∣|A|∣A∣ 。结果对 ppp ( ppp 不一定是质数)取模题目分析:高斯消元,类似于求解线性方程组,消成三角矩阵,然后求三角矩阵的行列式即可注意交换行列要对结果乘 −1-1−1具体细节见代码://#pragma GCC optimize(2)//#pragma GCC optimize("Ofast","inline","-ffast-math")//#pragma GCC target("avx,sse2,ss原创 2021-09-16 23:07:34 · 229 阅读 · 1 评论 -
P2455 [SDOI2006]线性方程组
题目链接:点击这里题目大意:给定一个线性方程组,对其求解题目分析:高斯-约旦消元法利用矩阵的初等变换,将矩阵变成对角矩阵即可无解输出 −1-1−1 ,无穷多解输出 000时间复杂度 O(n3)O(n^3)O(n3)注意:洛谷模板题数据极弱,建议用此题测试模板具体细节见代码://#pragma GCC optimize(2)//#pragma GCC optimize("Ofast","inline","-ffast-math")//#pragma GCC target("avx,ss原创 2021-09-16 17:46:39 · 181 阅读 · 0 评论 -
AcWing 3069. 圆的面积并(自适应辛普森积分)
题目链接:点击这里题目大意:给出 nnn 个圆的圆心和半径,求这 nnn 个圆的面积并题目分析:自适应辛普森积分考虑构造 f(x)f(x)f(x) 表示横坐标为 xxx 时的纵坐标长度然后套上辛普森积分的板子即可具体细节见代码://#pragma GCC optimize(2)//#pragma GCC optimize("Ofast","inline","-ffast-math")//#pragma GCC target("avx,sse2,sse3,sse4,mmx")#inclu原创 2021-09-14 23:01:44 · 212 阅读 · 0 评论 -
P3172 [CQOI2015]选数(莫比乌斯反演+杜教筛)
题目链接:点击这里题目大意:给出一个区间 [l,r][l,r][l,r] ,求从区间中选出 kkk 个数(选了的数不会被删除),其 gcd\gcdgcd 等于 kkk 的方案数题目分析:看到 gcd=k\gcd =kgcd=k 不难想到莫比乌斯反演,令 F(k)F(k)F(k) 表示 k∣gcdk|\gcdk∣gcd 的方案数,f(k)f(k)f(k) 表示 gcd=k\gcd =kgcd=k 的方案数则有:F(k)=(⌊rk⌋−⌊l−1k⌋)n=∑n∣df(d)F(k)=(\lfloo原创 2021-09-11 22:22:13 · 100 阅读 · 0 评论 -
月月给华华出题(欧拉函数)
题目链接:点击这里题目大意:给定一个正整数 NNN ,对于所有的 n∈[1,N]n\in[1,N]n∈[1,N] 输出 ∑i=1nigcd(i,n)\sum_{i=1}^n\frac{i}{\gcd(i,n)}∑i=1ngcd(i,n)i题目分析:∑i=1nigcd(i,n)\sum_{i=1}^n\frac{i}{\gcd(i,n)}i=1∑ngcd(i,n)i=∑d=1n1d∑i=1ni[gcd(i,n)=d]=\sum_{d=1}^n\frac{1}{d}\sum_{i=1}^原创 2021-09-08 21:53:40 · 197 阅读 · 0 评论 -
CodeForces 1288C Two Arrays(dp/组合数学+lucas)
题目链接:点击这里题目大意:给出两个数 n,mn,mn,m ,求满足下列条件的 ai,bia_i,b_iai,bi 序列数量:两个序列长度为 mmm两个序列元素取值都在 [1,n][1,n][1,n]对于任意的 iii 均满足 ai≥bia_i\ge b_iai≥biaia_iai 单调不下降bib_ibi 单调不上升题目分析:观察可发现其本质就是构造一个长度为 2m2m2m 的序列且该序列单调不下降dp:容易想到二维的状态 dp[i][j]dp[i][j]dp[i][原创 2021-09-08 21:34:19 · 191 阅读 · 0 评论 -
2021牛客暑期多校训练营6 G Hasse Diagram(min25筛)
题目链接:点击这里题目大意:哈斯图,设 HnH_nHn 为 nnn 的因子组成的哈斯图,对于任意两个因子 u,v(u>v)u,v(u>v)u,v(u>v) ,u,vu,vu,v 之间有边当且仅当 u=p⋅vu=p·vu=p⋅v例如 H12H_{12}H12 如下图所示:设 f(x)f(x)f(x) 为 HxH_xHx 的边数和,现在给定一个正整数 nnn 求 ∑i=1nf(i)\sum_{i=1}^nf(i)∑i=1nf(i)题目分析:...原创 2021-08-23 22:17:33 · 306 阅读 · 2 评论 -
Powerful Number 筛
本文主要讲解 Powerful Number 筛,并为之前的杜教筛进行一定的补充powerful number\text{powerful number}powerful number 定义:如果一个正整数 xxx 所有质因子的幂次均大于等于二,即 x=∏pieix=\prod p_i^{e_i}x=∏piei 且 ei≥2e_i\geq 2ei≥2 ,则称 xxx 是一个 powerful number\text{powerful number}powerful&nb原创 2021-08-23 17:55:27 · 384 阅读 · 0 评论 -
LibreOJ - 6053 简单的函数(min25筛)
题目链接:点击这里题目大意:定义一个正整数域函数 f(x)f(x)f(x) 如下:f(1)=1f(1)=1f(1)=1f(pc)=p⨁cf(p^c)=p\bigoplus cf(pc)=p⨁c (ppp 为质数)f(ab)=f(a)⋅f(b)f(ab)=f(a)·f(b)f(ab)=f(a)⋅f(b) (a,ba,ba,b 互质)现在给定一个 nnn ,求 ∑i=1nf(i)mod 1e9+7\sum_{i=1}^nf(i) \mod 1e9+7∑i=1nf(i)mod1e9+7题目分原创 2021-08-18 17:19:44 · 115 阅读 · 2 评论 -
LibreOJ - 6202 叶氏筛法(min25筛)
题目分析:点击这里题目大意:求区间 [l,r][l,r][l,r] 范围内的质数和题目分析:min25min25min25 筛阉割版令完全积性函数 F(x)=id(x)=xF(x)=id(x)=xF(x)=id(x)=x同样的做 s(n,j)s(n,j)s(n,j) 转移得到 s(n,∣P∣)s(n,|P|)s(n,∣P∣)因为我们要求质数和,所以对 S(n,j)S(n,j)S(n,j) 求质数贡献即可即 ans=g(n,∣p∣)−∑i=1j−1f(pi)ans=g(n,|p|)-\sum_原创 2021-08-17 00:08:41 · 160 阅读 · 0 评论 -
LibreOJ - 6235 区间素数个数(min25筛)
题目链接:点击这里题目大意:给定一个数 nnn ,求 111 到 nnn 的质数个数题目分析:阉割版 min25min25min25 筛不妨令完全积性函数 F(x)=id(x)F(x)=id(x)F(x)=id(x)然后做 s(n,j)s(n,j)s(n,j) 的转移由于 f(pk)=0f(p^k)=0f(pk)=0 ,所以合数贡献就可以不算了,故求得的 s(n,0)s(n,0)s(n,0) 就是所求答案具体细节见代码://#pragma GCC optimize(2)//#pragma原创 2021-08-16 22:56:47 · 153 阅读 · 0 评论 -
min_25筛
笔记写在前面:min25min25min25 筛是用来求积性函数前缀和 ∑i=1nf(i)\sum_{i=1}^nf(i)∑i=1nf(i) 的一个筛法,其时间复杂度为 O(n34logn)O(\frac{n^{\frac 34}}{\log n})O(lognn43) 大约 1s1s1s 能跑 1e101e101e10。该积性函数需要满足:f(p)f(p)f(p) 是一个关于质数 ppp 的项数较小的多项式 或者 能快速求值f(pk)f(p^k)f(pk) 可以快速求值我们把 f(i原创 2021-08-15 20:19:53 · 93 阅读 · 0 评论 -
2021牛客暑期多校训练营2 J Product of GCDs(数学)
题目链接:点击这里题目大意:给定一个整数集合,求所以大小为 kkk 的子集的 gcd\gcdgcd 之积题目分析:我们发现求的是所有 gcd\gcdgcd 的乘积,所以我们可以单独考虑每一个数字 ggg 的贡献,若有 xxx 个数是 ggg 的倍数, yyy 个数是 k⋅g(k≥1)k·g(k\ge 1)k⋅g(k≥1) 的倍数,那么 ggg 产生的贡献就是 gCx−ykg^{C_{x-y}^k}gCx−yk 。在操作容斥时,我们可以先用 CxkC_x^kCxk 减去 ∑g∣iCdpik\原创 2021-07-30 19:40:59 · 176 阅读 · 0 评论 -
hdu-6975 Forgiving Matching(卷积处理模糊串匹配)
题目链接:点击这里题目大意:给出一个长度为 nnn 的字符串 SSS 和一个长度为 mmm 的字符串 TTT 。我们定义 kkk 匹配为两个长度相同的字符串至多有 kkk 个位置不同,求当 k∈[0,m]k \in[0,m]k∈[0,m] 时,TTT 可以和 SSS 中的多少个子串匹配题目分析:类似于 P4173 残缺的字符串 的思路。我们发现字符集很小,所以考虑枚举字符集的元素。如果 Si=TjS_i=T_jSi=Tj ,那么当前字符在 SSS 中以 i+m−ji+m-ji+m−j 结尾的原创 2021-07-28 17:21:36 · 347 阅读 · 0 评论 -
P4173 残缺的字符串(卷积处理模糊串匹配)
题目链接:点击这里题目大意:给你一个长度为 mmm 的字符串 AAA ,和一个长度为 nnn 的字符串 BBB,求对于 BBB 的每一个位置 iii ,从这个位置开始连续 mmm 个字符形成的子串是否可能与 AAA 串完全匹配,输出所有可能匹配的位置题目分析:构造匹配函数 C(x,y)=[A(x)−B(y)]2A(x)B(y)C(x,y)=[A(x)-B(y)]^2A(x)B(y)C(x,y)=[A(x)−B(y)]2A(x)B(y)故完全匹配函数为 P(x)=∑i=0m−1C(x,x−(m−1)原创 2021-07-28 11:29:06 · 144 阅读 · 0 评论 -
2021牛客暑期多校训练营1 H Hash Function(卷积)
题目链接:点击这里题目大意:给定 nnn 个互不相同的书,找一个最小的模域,使得他们在这个模域下互不相同题目分析:我们发现 a≡bmod ma\equiv b\mod ma≡bmodm 当且仅当 ∣a−b∣|a-b|∣a−b∣ 能被 mmm 整除。这样问题就转化成求一个最小的 mmm ,其满足他不是任何一个 ∣ai−aj∣|a_i-a_j|∣ai−aj∣ 的约数。直接枚举时间复杂度是 O(n2)O(n^2)O(n2) 的,显然会超时。考虑如何优化该过程,我们发现 1≤∣ai−aj∣≤5e5原创 2021-07-27 20:14:07 · 104 阅读 · 0 评论 -
hdu-6989 Pass!(广义斐波那契数列求通项+bsgs)
题目链接:点击这里题目大意:有 nnn 个人在传球,从 111 号开始传球,每次只能传给相邻的人,设传了 ttt 次球后球回到 111 的方案数为cntcntcnt ,求最小的 ttt 其满足此时的 cntmod 998244353=xcnt \mod 998244353=xcntmod998244353=x ,如果无解输出 −1-1−1题目分析:设 dp[i]dp[i]dp[i] 为经过 iii 次传球回到 111 号的方案数,容易发现有转移方程 dp[i]=(n−1)∗dp[i−1]+(n−2原创 2021-07-25 13:56:32 · 571 阅读 · 5 评论 -
CodeForces - 1499D The Number of Pairs(数学+埃式筛)
题目链接:点击这里题目大意:给定三个正整数 c,d,xc,d,xc,d,x ,求有多少对 a,ba,ba,b 满足如下式子 c∗lcm(a,b)+d∗gcd(a,b)=xc*lcm(a,b)+d*gcd(a,b)=xc∗lcm(a,b)+d∗gcd(a,b)=x题目分析:原创 2021-07-12 22:11:23 · 203 阅读 · 0 评论 -
CodeForces - 1538D Another Problem About Dividing Numbers(数学+构造)
题目链接:点击这里题目大意:给出两个正整数 a,ba,ba,b ,每轮可以选择二者之一,不妨令为 xxx ,用一个正整数 c(c>1)c(c>1)c(c>1) 来让其替换为 xc\frac x ccx问经过 kkk 轮操作后是否能使 a,ba,ba,b 恰好相等题目分析:显然操作为给 xxx 去除掉他的一个不等于其本身的因子,最终都可以通过操作使其等于 111先考虑一种特殊情况 k=1k=1k=1 ,其特殊在只能对一个数进行操作,此时只有在 a∣ba|ba∣b 或 b∣ab原创 2021-07-09 19:40:05 · 348 阅读 · 0 评论 -
CodeForces - 1475G Strange Beauty(数学+dp)
题目链接:点击这里题目大意:给定一个长度为 nnn 的序列 aia_iai ,求一个该序列的子序列 bib_ibi ,满足对于任意的 i,ji,ji,j 有 bi∣bjb_i|b_jbi∣bj 或 bj∣bib_j|b_ibj∣bi ,求该子序列的最大大小题目分析:...原创 2021-07-09 16:09:49 · 121 阅读 · 0 评论 -
简单组合数学
题目:To monitor cryptocurrency exchange rates trader William invented a wonderful device consisting of n lights arranged in a row. The device functions in the following way:Initially, all lights on William’s device are turned off. At the beginning of a new原创 2021-06-10 20:44:40 · 415 阅读 · 1 评论 -
快速数论变换(NTT)
前置知识:阶:对于互质整数 a,na,na,n ,满足 ar≡1mod na^r \equiv 1\mod nar≡1modn 的最小的 rrr 就是 aaa 模 nnn 的阶。原根:对于正整数 nnn ,若整数 aaa 模 nnn 的阶等于 φ(n)\varphi(n)φ(n) ,则称 aaa 为模 nnn 的一个原根。快速数论变换NTTNTTNTT,实际上就是利用一些特殊的质数(例如 998244353=119×223+1998244353=119\times2^{23}+1998244353原创 2021-05-27 22:04:52 · 577 阅读 · 0 评论 -
快速傅里叶变换(FFT)
前言:老笨蛋羊驼了,发现之前学的多项式技术基本都忘干净了,开始重新学习,顺便在此记录一下快速傅里叶变换可以用来在近似于 O(nlogn)O(nlogn)O(nlogn) 的时间复杂度下计算多项式卷积预备知识复数单位根:定义 nnn 次单位复根为满足 ωn=1\omega^n=1ωn=1 的 ω\omegaω ;我们考虑用欧拉公式(eix=cosx+sinixe^{ix}=\cos x+\sin ixeix=cosx+sinix)取单位根令 x=2πx=2\pix=2π ,此时有 e2πi=c原创 2021-05-27 13:04:20 · 792 阅读 · 0 评论 -
P4213 【模板】杜教筛 (筛欧拉函数和莫比乌斯函数前缀和)
题目链接:点击这里题目大意:给定正整数 nnn ,求:ans1=∑i=1nφ(i)ans_1=\sum_{i=1}^n\varphi(i)ans1=i=1∑nφ(i)ans2=∑i=1nμ(i)ans_2=\sum_{i=1}^n\mu(i)ans2=i=1∑nμ(i)题目分析:杜教筛是一种用来筛积性函数前缀和的神奇筛法构造两个积性函数 h,gh,gh,g ,使之满足 h=f∗gh=f*gh=f∗g ,设 S(n)=∑i=1nf(i)S(n)=\sum_{i=1}^nf(i)S(n)=原创 2021-04-16 22:26:28 · 289 阅读 · 0 评论 -
数学笔记(数论函数、迪利克雷卷积)
数论函数定义:定义域为正整数,值域复数的一类函数积性函数:规定 f(1)=1f(1) = 1f(1)=1 ,当 gcd(a,b)=1\gcd(a,b)=1gcd(a,b)=1 时,有 f(ab)=f(a)f(b)f(ab)=f(a)f(b)f(ab)=f(a)f(b)完全积性函数:对于任意的 a,ba,ba,b ,有 f(ab)=f(a)f(b)f(ab) = f(a)f(b)f(ab)=f(a)f(b)常见的积性函数:ϵ(n)=[n=1] (完全积性函数)\epsilon (n)=[原创 2021-04-16 21:29:03 · 413 阅读 · 0 评论 -
常见数论函数线性筛
筛素数:int cnt,pri[maxn];bool vis[maxn];void get_pri(int n){ for(int i = 2;i <= n;i++) { if(!vis[i]) pri[++cnt] = i; for(int j = 1;j <= cnt && i*pri[j] <= n;j++) { vis[i*pri[j]] = true; if(i%pri[j] == 0) break; } }}筛莫比原创 2021-04-16 20:42:57 · 372 阅读 · 0 评论 -
P2257 YY的GCD (莫比乌斯反演)
题目链接:点击这里题目大意:ttt 组数据,每组数据给出 N,MN , MN,M ,求 1≤x≤N1 \leq x \leq N1≤x≤N ,1≤y≤M1 \leq y \leq M1≤y≤M 且 gcd(x,y)\gcd(x, y)gcd(x,y) 为质数的二元组 (x,y)(x, y)(x,y) 的数目题目分析:题目所求式子为:∑i=1N∑j=1M[gcd(i,j)=prime]\sum_{i=1}^N\sum_{j=1}^M[\gcd(i,j)=prime]i=1∑Nj=1∑M[g原创 2020-11-16 18:22:16 · 186 阅读 · 4 评论