数论
文章平均质量分 83
LH_Frank
这个作者很懒,什么都没留下…
展开
-
CF1716 F-Bags With Balls
第二类斯特林数原创 2022-08-06 00:05:35 · 165 阅读 · 1 评论 -
SEERC2018 Broken Clock
SEERC2018 Broken ClockProblem Statement时针, 分针, 秒针的程度分别为A,B,CA,B,CA,B,C, 并且假设时钟共有NNN个时刻. 并且由于时钟被破坏时针, 分针和秒针可以指向任意一时刻, 问时针, 分针和秒针能够形成多少个包含时针中心的三角形.Solution我们考虑三角形的生成函数. 假设时针分针和秒针组成三角形的三条边分别为a,b,ca,b,ca,b,c, 那么他们所对应的圆心角应该有θ1+θ2+θ3=2π\theta_1+\theta_2+\the原创 2022-04-13 23:50:01 · 416 阅读 · 0 评论 -
NowCoder LCMs(莫比乌斯反演)
NowCoder LCMsStatement给你一个长度为nnn的序列, A0,A1,…,An−1A_0,A_1,\dots,A_{n-1}A0,A1,…,An−1. ∑i=0n−2∑j=i+1n−1lcm(Ai,Aj) mod 998244353\sum_{i=0}^{n-2}\sum_{j=i+1}^{n-1}\text{lcm}(A_i,A_j)\text{ mod }998244353∑i=0n−2∑j=i+1n−1lcm(Ai,Aj) mod&nbs原创 2022-03-15 18:53:46 · 128 阅读 · 0 评论 -
NowCoder 14648 序列(莫比乌斯反演)
NowCoder 14648 序列Statement两个长度为nnn的序列a,ba,ba,b, 求满足gcd(x,y)=1\gcd(x,y)=1gcd(x,y)=1且abx=baya_{b_x}=b_{a_y}abx=bay的的对(x,y)(x,y)(x,y)数量SolutionAns=∑1≤x≤n且1≤y≤n[abx=bay]⋅[gcd(x,y)=1]Ans=\sum\limits_{1\leq x\leq n且1\leq y\leq n}[a_{b_x}=b_{a_y}]\cdot[原创 2022-03-15 18:49:46 · 147 阅读 · 0 评论 -
[HAOI2011] Problem B(莫比乌斯反演)
[HAOI2011] Problem BStatementAns=∑i=ab∑j=cd[gcd(i,j)=k]Ans=\sum_{i=a}^b\sum_{j=c}^d[\gcd(i,j)=k]Ans=∑i=ab∑j=cd[gcd(i,j)=k]Solution设: s(a,b,c,d,k)=∑i=ab∑j=bc[gcd(i,j)=k]s(a,b,c,d,k)=\sum_{i=a}^b\sum_{j=b}^c[\gcd(i,j)=k]s(a,b,c,d,k)=∑i=ab∑j=bc[gcd(原创 2022-03-15 18:45:27 · 93 阅读 · 0 评论 -
SDOI2015约数个数和(莫比乌斯反演)
SDOI2015约数个数和Statement∑i=1N∑j=1Md(i⋅j)\sum_{i=1}^N\sum_{j=1}^Md(i\cdot j)∑i=1N∑j=1Md(i⋅j).Solutiond(i⋅j)=∑x∣i∑y∣j[gcd(x,y)=1]d(i\cdot j)=\sum_{x|i}\sum_{y|j}[\gcd(x,y)=1]d(i⋅j)=∑x∣i∑y∣j[gcd(x,y)=1]Ans=∑i=1N∑j=1M∑x∣i∑y∣j[gcd(x,y)=1]=∑x=1N∑y=1M⌊原创 2022-03-15 18:33:07 · 80 阅读 · 0 评论 -
YY的GCD(莫比乌斯反演)
YY的GCDStatement求∑i=1N∑j=1M[gcd(i,j)=Prime]\sum_{i=1}^N\sum_{j=1}^M[\gcd(i,j)=Prime]∑i=1N∑j=1M[gcd(i,j)=Prime].Solution设: f(d)=∑i=1N∑j=1M[gcd(i,j)=d]f(d)=\sum_{i=1}^N\sum_{j=1}^M[\gcd(i,j)=d]f(d)=∑i=1N∑j=1M[gcd(i,j)=d].设: g(d)=∑d∣d′f(d′)=∑i=1N∑j=原创 2022-03-15 16:29:46 · 158 阅读 · 0 评论 -
POI2007 ZAP-Queries(莫比乌斯反演)
POI2007 ZAP-QueriesStatement求∑i=1a∑j=1b[gcd(a,b)=d]\sum_{i=1}^a\sum_{j=1}^b[\gcd(a,b)=d]∑i=1a∑j=1b[gcd(a,b)=d]Solution设f(d)=∑i=1N∑j=1M[gcd(i,j)=d]f(d)=\sum_{i=1}^N\sum_{j=1}^M[\gcd(i,j)=d]f(d)=∑i=1N∑j=1M[gcd(i,j)=d]g(d)=∑d∣d′f(d′)=∑i=1N∑j=1M[gcd原创 2022-03-15 16:28:34 · 90 阅读 · 0 评论 -
ABC162E Sum of gcd of Tuples(莫比乌斯反演)
ABC162E Sum of gcd of TuplesStatement∑A1=1K∑A2=1K⋯∑AN=1Kgcd(A1,A2,…,An)\sum_{A_1=1}^K\sum_{A_2=1}^K\dots\sum_{A_N=1}^K\gcd(A_1,A_2,\dots,A_n)∑A1=1K∑A2=1K⋯∑AN=1Kgcd(A1,A2,…,An).SolutionAns=∑d=1Kd∑1≤Ai≤K[gcd(A1,A2,…,AN)=d]Ans=\sum_{d=1}^Kd\sum原创 2022-03-14 16:05:21 · 370 阅读 · 0 评论 -
环计数问题
环计数问题Statement假设nnn个元组成一个环, 每个元都是1,…,r1,\dots,r1,…,r中的一个数, 两个环是不同的当且仅当它们不能通过旋转使得两个环对应的每个元素都相等.Solution如果我们现在先考虑链的情况, 那么一个长度为nnn的链, 每个元素是1,…,r1,\dots,r1,…,r中的一个数的方案数为rnr^nrn.我们设fxf_xfx表示最小循环节为xxx的环的方案数. 则我们可以得到如下性质.性质1: rn=∑d∣ndf(d)r^n=\sum\limits_{原创 2022-03-14 13:59:50 · 306 阅读 · 1 评论 -
CF757E Bash Plays with Functions(积性函数)
CF757E Bash Plays with FunctionsStatementf0(n)=∑p⋅q=n[gcd(p,q)=1]f_0(n)=\sum\limits_{p\cdot q=n}[\gcd(p,q)=1]f0(n)=p⋅q=n∑[gcd(p,q)=1]fr+1(n)=∑u⋅v=nfr(u)+fr(v)2f_{r+1}(n)=\sum\limits_{u\cdot v=n}\frac{f_r(u)+f_r(v)}{2}fr+1(n)=u⋅v=n∑2fr(u)+fr(v)原创 2022-03-14 13:44:32 · 151 阅读 · 0 评论