![](https://img-blog.csdnimg.cn/20190918140053667.png?x-oss-process=image/resize,m_fixed,h_224,w_224)
数论题
文章平均质量分 89
数论
凌乱之风
ACM蒟蒻
展开
-
[2021 CCPC 广州 A] Math Ball (生成函数 多项式)
[2021 CCPC 广州 A] Math Ball (生成函数 多项式)原创 2022-10-16 00:33:06 · 700 阅读 · 0 评论 -
[2021CCPC 威海G] Shinyruo and KFC (下降幂多项式乘法+下降幂转普通幂+多项式多点求值)
[2021CCPC 威海G] Shinyruo and KFC (下降幂多项式乘法+下降幂转普通幂+多项式多点求值)原创 2022-10-04 15:22:49 · 814 阅读 · 0 评论 -
[TJOI2019] 唱、跳、rap和篮球 (二项式反演 NTT)
[TJOI2019] 唱、跳、rap和篮球 (二项式反演 NTT)原创 2022-09-22 20:07:24 · 305 阅读 · 0 评论 -
[NOI2020统一省选 A] 组合数问题 (推式子)
[NOI2020统一省选 A] 组合数问题 (推式子)原创 2022-09-01 16:09:34 · 303 阅读 · 0 评论 -
[2021ICPC济南 L] Strange Series (Bell 数 多项式exp)
[2021ICPC济南 L] Strange Series (Bell 数 多项式exp)原创 2022-08-30 13:27:21 · 384 阅读 · 0 评论 -
[2022CCPC华为云1005] 带权子集和 (NTT)
[2022CCPC华为云1005] 带权子集和 (NTT)原创 2022-08-22 14:18:01 · 348 阅读 · 0 评论 -
[Educational Codeforces Round 133 F] Bags with Balls (组合计数 推式子)
[Educational Codeforces Round 133 F] Bags with Balls (组合计数 推式子)原创 2022-08-05 15:38:39 · 199 阅读 · 0 评论 -
[2022 杭电多校5] Count Set (生成函数 分治NTT)
[2022 杭电多校5] Count Set (生成函数 分治NTT)原创 2022-08-03 11:22:33 · 421 阅读 · 4 评论 -
[2022 牛客多校2 E] Falfa with Substring (二项式反演 NTT)
[2022 牛客多校2 E] Falfa with Substring (二项式反演 NTT)原创 2022-08-01 10:49:35 · 262 阅读 · 0 评论 -
[2022 牛客多校4 C] Easy Counting Problem (生成函数 NTT)
[2022 牛客多校4 C] Easy Counting Problem (生成函数 NTT)原创 2022-08-01 10:48:51 · 266 阅读 · 0 评论 -
[2016 CCPC 杭州J] Just a Math Problem (莫比乌斯反演)
[2016 CCPC 杭州J] Just a Math Problem (莫比乌斯反演)原创 2022-07-15 09:02:02 · 146 阅读 · 0 评论 -
[2022 广东省赛M] 拉格朗日插值 (多元函数极值 分治NTT)
gdcpc m 题解原创 2022-07-01 21:59:25 · 337 阅读 · 0 评论 -
[AtCoder arc140 D] One to One (图论 计数 分治NTT)
题意初始有 nnn 个点,给定一个长度为 nnn 的数组 aia_iai,若 ai≠−1a_i \ne -1ai=−1,则有无向边 (i,ai)(i, a_i)(i,ai),若 ai=−1a_i = -1ai=−1,则点 iii 可以连向 1∼n1 \sim n1∼n 任意点,求所有图的联通块个数之和1≤n≤2×103,ai∈[1,n]∪{−1}1 \le n \le 2 \times 10 ^ 3, a_i \in [1, n] \cup \{-1\}1≤n≤2×103,ai∈[1,n]原创 2022-05-18 12:49:24 · 380 阅读 · 0 评论 -
[Wannafly 28] msc的背包 (生成函数 组合数)
题意有 nnn 种体积为 111 的物品和 mmm 种体积为 222 的物品,求选择物品的体积为 kkk 的方案数对 998244353998244353998244353 取模(1≤n,m≤106,1≤k≤9×108)(1 \le n, m \le 10 ^ 6,1 \le k \le 9 \times 10 ^ 8)(1≤n,m≤106,1≤k≤9×108)分析:所有体积为 111 的生成函数为F(x)=(∑i=0∞xi)nF(x) = \left ( \sum_{i = 0} ^ {\i原创 2022-05-04 21:22:27 · 483 阅读 · 1 评论 -
[2020 BSUIRPC] Function analysis (EGF NTT)
题意给定三个正整数 n,d,kn, d, kn,d,k ,现有排列 p=(1,2,3,⋯ ,n)p = (1,2,3,\cdots,n)p=(1,2,3,⋯,n) ,有 n−d+1n - d + 1n−d+1 个询问,对于每个询问有正整数 m(d≤m≤n)m (d \le m \le n)m(d≤m≤n) ,现从 ppp 中随机可重复地选取 mmm 个数构成序列 qqq ,求 qqq 中第 ddd 小数大于 kkk 的概率,对 998244353998244353998244353 取模。分析:对于原创 2022-04-13 12:58:30 · 197 阅读 · 0 评论 -
[2021ICPC上海 B] Strange Permutations (容斥 分治NTT)
题意给定一个长度为 nnn 的 1∼n1 \sim n1∼n 排列 PPP,找到有多少个 1∼n1 \sim n1∼n 的排列 QQQ 使得 ∀i∈[1,n−1],Qi+1≠PQi\forall i \in[1, n - 1], Q_{i + 1} \ne P_{Q_i}∀i∈[1,n−1],Qi+1=PQi对 998244353998244353998244353 取模(1≤n≤105,1≤Pi≤n)(1 \le n \le 10 ^ 5, 1 \le P_i \le n)(1≤n≤105原创 2022-03-28 11:37:18 · 856 阅读 · 3 评论 -
[AtCoder abc234 F] Reordering (EGF NTT)
题意给定一个长度为 nnn 的字符串 SSS,计算有多少种非空子序列的排列,对 998244353998244353998244353 取模。1≤n≤5×1031 \le n \le 5\times 10^31≤n≤5×103分析:设 262626 个英文字母每个字母 uuu 的生成函数为(因为要计算排列,所以是 EGF\textbf{EGF}EGF)1+11!x+12!x2+⋯+1cntu!xcntu1 + \frac{1}{1!}x + \frac{1}{2!}x^2 + \cdots +原创 2022-01-09 16:59:50 · 889 阅读 · 0 评论 -
[2021CCPC 威海M] 810975 (组合计数 容斥原理 多项式快速幂)
题意夜吹下了 nnn 盘酒馆战旗,吃了 mmm 次鸡,并且 kkk 连鸡,求方案数对 998244353998244353998244353 取模。分析:已知下了 nnn 盘棋,吃了 mmm 次鸡,所以有 n−mn - mn−m 局是未吃鸡的,所以考虑把 mmm 次鸡插空,n−mn - mn−m 局未吃鸡局有 n−m+1n - m + 1n−m+1 个空位,那么每个空的生成函数为1+x+x2+x3+⋯+xk1+x+x^2+x^3+\cdots+x^k1+x+x2+x3+⋯+xk所以所有的方案数原创 2021-12-08 14:59:43 · 1897 阅读 · 6 评论 -
[2021CCPC威海热身赛] Number Theory (打表 推公式)
题意求∑k=1n∑i∣k∑j∣iλ(i)λ(j)\sum_{k = 1}^{n}\sum_{i \mid k} \sum_{j \mid i} \lambda(i) \lambda(j)k=1∑ni∣k∑j∣i∑λ(i)λ(j)对 998244353998244353998244353 取模其中 λ(x)=(−1)∑iei,x=∏ipiei\lambda(x) = (-1)^{\sum\limits_{i}e_i},x=\prod\limits_{i}p_i^{e_i}λ(x)=(−1)i∑原创 2021-11-20 18:52:14 · 952 阅读 · 2 评论 -
[CQOI2015] 选数 (莫比乌斯反演 杜教筛)
题意求从区间 [L,R][L,R][L,R] 选出 nnn 个数使得最大公约数为 kkk 的方案数,对 109+710^9 + 7109+7 取模1≤n,k≤1091 \le n,k \le 10^91≤n,k≤1091≤L≤R≤1091 \le L \le R \le10^91≤L≤R≤109分析:根据题意∑a1=LR∑a2=LR⋯∑an=LR[gcd(a1,a2,⋯ ,an)=k]\sum_{a_1=L}^{R}\sum_{a_2=L}^{R}\cdots\sum_{a_n=L}^{R}[原创 2021-10-19 17:24:55 · 100 阅读 · 0 评论 -
[洛谷 P5221] Product (莫比乌斯反演)
题意求∏i=1n∏j=1nlcm(i,j)gcd(i,j)\prod_{i=1}^{n}\prod^{n}_{j=1}\frac{\text{lcm}(i,j)}{\gcd(i,j)}i=1∏nj=1∏ngcd(i,j)lcm(i,j)对 104857601104857601104857601 取模1≤n≤1061 \le n \le 10^61≤n≤106分析:将 lcm(i,j)=i⋅jgcd(i,j)\text{lcm}(i,j)=\dfrac{i \cdot j}{\gcd(原创 2021-10-19 10:14:31 · 113 阅读 · 1 评论 -
[2021 CCCC天梯赛] 可怜的简单题 (概率期望 莫比乌斯反演 杜教筛)
题意每次从 [1,n][1,n][1,n] 中选择一个数加到一个序列末尾,当 gcd(a1,⋯ ,an)=1\gcd(a_1,\cdots,a_n)=1gcd(a1,⋯,an)=1 时停止,求期望长度,对 ppp 取模1≤n≤1011,n<p≤10121\le n \le 10^{11},n< p \le 10 ^{12}1≤n≤1011,n<p≤1012分析:设 E(len)E(len)E(len) 为期望长度,那么根据期望定义E(len)=∑i=1∞P(len=i)⋅i原创 2021-09-23 12:36:41 · 1707 阅读 · 3 评论 -
[NC 200008] Lady Layton with Math (杜教筛)
题意求∑i=1n∑j=1nφ(gcd(i,j))\sum_{i=1}^{n}\sum_{j=1}^{n} \varphi(\gcd(i,j))i=1∑nj=1∑nφ(gcd(i,j))1≤n≤1091 \le n \le 10^91≤n≤109,对 109+710^9+7109+7 取模分析:枚举 gcd(i,j)\gcd(i,j)gcd(i,j)∑d=1nφ(d)∑i=1n∑j=1n[gcd(i,j)=d]\sum_{d=1}^{n}\varphi(d)\sum_{i=1}^{n}\原创 2021-09-23 08:59:04 · 103 阅读 · 0 评论 -
[2018 icpc徐州网络赛] Easy Math (杜教筛)
题意求∑i=1mμ(in)\sum_{i=1}^{m} \mu(in)i=1∑mμ(in)m≤2×109,n≤1012m \le 2×10^9,n\le 10^{12}m≤2×109,n≤1012分析:首先分析 nnn 有没有平方因子,如果有那么答案就是 000原创 2021-09-18 09:14:09 · 118 阅读 · 0 评论 -
[洛谷 P6055] [RC-02] GCD (莫比乌斯反演 杜教筛)
题意求∑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][gcd(p,q)=1]对 99824原创 2021-09-17 20:34:22 · 125 阅读 · 0 评论 -
[2019 icpc西安邀请赛] Product (莫比乌斯反演 杜教筛)
题意求∏i=1n∏j=1n∏k=1nmgcd(i,j)[k∣gcd(i,j)] mod p\prod_{i=1} ^{n} \prod_{j=1}^{n}\prod_{k=1}^{n}m^{\gcd(i,j)[k \mid \gcd(i,j)]} \bmod pi=1∏nj=1∏nk=1∏nmgcd(i,j)[k∣gcd(i,j)]modpn≤109,m≤2×109,p≤2×109n \le10^9,m \le 2 ×10^9,p\le 2×10^9n≤109,m≤2×109,p≤2×109原创 2021-09-17 20:13:13 · 150 阅读 · 0 评论 -
[洛谷 P4318] 完全平方数 (杜教筛)
题意TTT 组询问,回答第 KiK_iKi 个不是完全平方数的正整数倍的数。1≤Ki≤109,T≤501\le K_i \le 10^9,T \le 501≤Ki≤109,T≤50分析:法一:如果一个数 nnn 不是完全平方数,那么 n=p1α1p2α2⋯pkαkn=p_1^{\alpha_1}p_2^{\alpha_2} \cdots p_k^{\alpha_k}n=p1α1p2α2⋯pkαk 中 0≤αi≤10 \le \alpha_i \le 10≤αi≤1,所以就想到了莫原创 2021-09-15 20:26:31 · 168 阅读 · 0 评论 -
[2019 CCPC网络赛] huntian oy (杜教筛)
题意TTT 组输入,给定 n,a,bn,a,bn,a,b 求f(n,a,b)=∑i=1n∑j=1igcd(ia−ja,ib−jb)[gcd(i,j)=1]f(n,a,b)=\sum_{i=1} ^{n}\sum_{j=1}^{i} \gcd(i^a-j^a, i ^ b - j ^ b)[\gcd(i,j)=1]f(n,a,b)=i=1∑nj=1∑igcd(ia−ja,ib−jb)[gcd(i,j)=1]对 109+710^9+7109+7 取模,1≤n,a,b≤1091 \le n,a,b原创 2021-09-14 10:50:41 · 128 阅读 · 0 评论 -
[洛谷 P3768] 简单的数学题 (莫比乌斯反演 狄利克雷卷积 杜教筛)
题意:求∑i=1n∑j=1nijgcd(i,j)\sum_{i=1}^{n}\sum_{j=1}^{n}ij\gcd(i,j)i=1∑nj=1∑nijgcd(i,j)对 ppp 取模,n≤1010,5×108≤p≤1.1×109n \le10^{10}, 5 ×10^8\le p \le1.1 ×10^{9}n≤1010,5×108≤p≤1.1×109分析:∑i=1n∑j=1nijgcd(i,j)\sum_{i=1}^{n}\sum_{j=1}^{n}ij\gcd(i,j)i=1∑nj=原创 2021-09-13 21:21:46 · 192 阅读 · 2 评论 -
[洛谷 P1829] Crash的数字表格 (莫比乌斯反演)
题意:求∑i=1n∑j=1mlcm(i,j)\sum_{i = 1} ^{n} \sum_{j = 1} ^{m} \text{lcm}(i, j)i=1∑nj=1∑mlcm(i,j)对 201010092010100920101009 取模分析:首先 lcm(i,j)=i⋅jgcd(i,j)\text{lcm}(i, j) = \dfrac{i \cdot j}{\gcd(i,j)}lcm(i,j)=gcd(i,j)i⋅j 代入:∑i=1n∑j=1mi⋅jgcd(i,j)\sum_{原创 2021-09-09 11:48:44 · 99 阅读 · 0 评论 -
[洛谷 P1891] 疯狂 LCM (欧拉函数 莫比乌斯反演)
题意:求∑i=1nlcm(i,n)\sum_{i=1} ^{n} \text{lcm}(i,n)i=1∑nlcm(i,n)分析:法一:拆一下 lcm(i,n)=i⋅ngcd(i,n)\text{lcm}(i,n) = \dfrac{i \cdot n}{\gcd{(i,n)}}lcm(i,n)=gcd(i,n)i⋅n 变为:∑i=1ni⋅ngcd(i,n)\sum_{i=1} ^{n} \frac{i \cdot n}{\gcd{(i,n)}}i=1∑ngcd(i,n)i⋅n枚举原创 2021-09-08 20:02:50 · 268 阅读 · 0 评论 -
[NOI2010] 能量采集 (莫比乌斯反演 欧拉反演)
题意:求2∗∑i=1n∑j=1mgcd(i,j)−n∗m2*\sum_{i=1}^{n}\sum_{j=1}^{m}\gcd(i,j)-n*m2∗i=1∑nj=1∑mgcd(i,j)−n∗m分析:莫比乌斯反演:那么只需要求∑i=1n∑j=1mgcd(i,j)\sum_{i=1}^{n}\sum_{j=1}^{m}\gcd(i,j)i=1∑nj=1∑mgcd(i,j)枚举每个 gcd(i,j)\gcd(i,j)gcd(i,j) 的值∑d=1nd∑i=1n∑j=1m[gcd(i,原创 2021-08-24 08:27:51 · 117 阅读 · 0 评论 -
[SDOI2015] 约数个数和 (莫比乌斯反演)
题意:设 d(x)d(x)d(x) 为 xxx 的约数个数,求∑i=1N∑j=1Md(ij)\sum_{i=1}^{N}\sum_{j=1}^{M} d(ij)i=1∑Nj=1∑Md(ij)分析:首先 i=p1α1⋯pkαk,j=p1β1⋯pkβk,i×j=p1α1+β1⋯pkαk+βki=p_1 ^{\alpha_1}\cdots p_k^{\alpha_k},j=p_1 ^{\beta_1}\cdots p_k^{\beta_k},i×j=p_1 ^{\alpha_1+\beta_1}\cd原创 2021-07-23 16:08:06 · 171 阅读 · 0 评论 -
[HAOI2011] problem b (莫比乌斯反演)
题意:求∑i=ab∑j=cd[gcd(i,j)=k]\sum_{i=a}^{b}\sum_{j=c}^{d}[\gcd(i,j)=k]i=a∑bj=c∑d[gcd(i,j)=k]分析:用二维前缀和的思想,把答案分为444个部分令S(x,y)=∑i=1x∑j=1y[gcd(i,j)=k]令S(x,y)=\sum_{i=1}^{x}\sum_{j=1}^{y}[\gcd(i,j)=k]令S(x,y)=i=1∑xj=1∑y[gcd(i,j)=k]则∑i=ab∑j=bd[gcd(i,j)原创 2021-07-22 17:33:32 · 76 阅读 · 0 评论