数论计数
Freopen
这个作者很懒,什么都没留下…
展开
-
LOJ #2476. 「2018 集训队互测 Day 3」蒜头的奖杯(三元环计数)
题目 泰勒应天下大雨! 类似于SDOI2018SDOI2018SDOI2018旧试题。 给DDD,EEE,FFF卷上μ\muμ。 得到D′,E′,F′D',E',F'D′,E′,F′。 原式变为∑i,j,kAiBjCk∑a∣i,a∣jDa′∑b∣i,b∣kEb′∑c∣j,c∣kFc′\sum_{i,j,k} A_iB_jC_k \sum_{a|i,a|j} D'_a \sum_{b|i,b|k} E'_b \sum_{c|j,c|k} F'_ci,j,k∑AiBjCka∣i,a∣j∑Da′b∣i原创 2020-07-03 21:12:58 · 503 阅读 · 0 评论 -
【WC2014】时空穿梭(莫比乌斯反演)
题目 枚举每一维上第一个点和最后一个点的坐标差(x1,..xn)(x_1,..x_n)(x1,..xn) ans=∑x1,x2...xn((gcdi=1nxi)−1c−2)∏i=1n(mi−xi−1)=∑d=1min1nmi(d−1c−2)∑x1,x2...xn≤mid[gcdi=1nxi=1]∏i=1n(mi−dxi−1)=∑d=1min1nmi(d−1c−2)∑x1,x2...xn≤mid∑p∣x1,p∣x2...μ(p)∏i=1n(mi−dxi−1)=∑d=1min1nmi(d−1c−2原创 2020-05-29 18:56:12 · 317 阅读 · 0 评论 -
Powerful Number学习笔记
比较好的资源只是可能会转很久 就是对于积性函数F(x)F(x)F(x)求前缀和。 构造出一个满足对于任意质数p,F(p)=G(p)p,F(p) = G(p)p,F(p)=G(p)的积性函数G(x)G(x)G(x)。 那么设H(x)=F(x)G(x)H(x) = \frac {F(x)}{G(x)}H(x)=G(x)F(x),注意这里是狄利克雷卷积意义下的除法。 有F(x)=H(x)G(x)F(x) = H(x) G(x)F(x)=H(x)G(x) 考虑F(p)=H(1)G(p)+H(p)G(1)=G(p)原创 2020-05-20 18:51:55 · 519 阅读 · 0 评论 -
LOJ #138. 类欧几里得算法
题目 推导方式比较牛逼。 消去公式F(k1,k2,n,a,b,c)F(k1,k2,n,a,b,c)F(k1,k2,n,a,b,c)原创 2020-01-22 13:36:58 · 474 阅读 · 1 评论 -
SPOJ DIVCNTK(MIN25筛)
模板题,没啥好讲的。 AC Code\mathrm{AC \ Code}AC Code #include<bits/stdc++.h> #define maxn 2000005 #define LL unsigned long long using namespace std; LL n,K=3,s[maxn],a[maxn]; int id,sn,cnt_pr,...原创 2020-01-14 16:29:42 · 154 阅读 · 0 评论 -
SPOJ 26073 DIVCNT1 - Counting Divisors(SB树上二分折线斜率拟合曲线)
题目 可以看LG的题解 还可以看zzt的论文(但是他好像也没说什么) 其实还是和之前一样,树上二分出当前最优的斜率,然后就按照这个斜率走到不能走, 然后继续二分,用单调栈符合了反比例函数拟合曲线斜率单调的特点,所以可以证明其折线(和二分次数)数量不会太大,为O(n13)O(n^{\frac 13})O(n31),然后就可以算了。 AC Code\mathrm {AC \ Code}AC...原创 2020-01-14 15:50:22 · 403 阅读 · 0 评论 -
#6543. 大葱的神谕(杜教筛)
It’s Obvious that the Chinese Green Onion is the STRONGest MAN in the universe. It’s Such A Universal truth that Even woman and children all know it !! 求∑i=1n(<=1e9)φ(n)(−1)n+1\sum_{i=1}^{n(<=...原创 2019-12-15 11:15:09 · 132 阅读 · 0 评论 -
莫比乌斯与猫与下取整
这式子太丑了我推不来 数论函数求和化简第一步尝试交换和号。 这个题发现没有啥好交换的。 第二步就不知道了,我尝试换为枚举T=dgT = dgT=dg,然而也不好做。 第三步:没有好交换的和号我们就创造和号。 ∑d∣n,d∣S∑g∣ndμ(g)⌊mdg⌋=∑d∣(n,S)∑g∣ndμ(g)∑i=1⌊mdg⌋1=∑d∣(n,S)∑i=1⌊md⌋∑g∣(nd,i)μ(g)=∑d∣(n,S)∑i=1⌊m...原创 2019-11-12 13:45:09 · 165 阅读 · 0 评论 -
BZOJ 3601 一个人的数论
第一次被搞。 这个套路深啊。。。。。。 ∑i=1n−1[gcd(i,n)==1]id=∑i=1n∑p∣i,p∣nμ(p)id=∑p∣nnμ(p)∑i=1np(ip)d=∑p∣nnμ(p)pd∑i=1npid\sum_{i=1}^{n-1} [\gcd(i,n) == 1] i^d \\ =\sum_{i=1}^n \sum_{p|i,p|n} \mu(p) i^d\\ =\sum_{p|n}^...原创 2019-10-16 17:00:24 · 184 阅读 · 0 评论