![](https://img-blog.csdnimg.cn/20201014180756918.png?x-oss-process=image/resize,m_fixed,h_64,w_64)
学习笔记
文章平均质量分 75
ruizhangj
这个作者很懒,什么都没留下…
展开
-
拉格朗日插值
拉格朗日插值原创 2022-08-11 16:11:08 · 163 阅读 · 1 评论 -
Powerful Number筛
Powerful Number筛原创 2022-08-04 10:27:22 · 174 阅读 · 0 评论 -
q-二项式定理
q-二项式定理原创 2022-08-04 08:29:03 · 240 阅读 · 0 评论 -
自然数幂和
自然数幂和O(k2)O(k^2)O(k2)递推做法首先令S(n,k)=∑i=1nikS(n,k)=\sum_{i=1}^{n}i^kS(n,k)=∑i=1nik,那么我们把S(n,k+1)S(n,k+1)S(n,k+1)给拆开,可以得到如下的方程:$$\begin{aligned}S(n,k+1) & =\sum_{i=1}{n}i{k+1} \S(n,k+1) & =\sum_{i=2}{n}(i-1+1){k+1}+1\S(n,k+1) & =\sum_{j=0}原创 2022-03-29 21:45:44 · 319 阅读 · 0 评论 -
多项式全家桶
多项式全家桶FFT和NTT首先最基本的FFT和NTT什么的应该能够默写了吧。牛顿迭代法假如有几个关于A(x)A(x)A(x)的函数ggg,现在要求解方程g(A(x))≡0mod xng(A(x))\equiv 0 \mod x^ng(A(x))≡0modxn,那么不妨先求出方程g(A0(x))≡0mod x⌈n2⌉g(A_0(x))\equiv 0\mod x^{\lceil \frac{n}{2}\rceil}g(A0(x))≡0modx⌈2n⌉的解A0(x)A_0(x)A0(x),那原创 2022-03-29 10:13:41 · 300 阅读 · 0 评论 -
折线法学习笔记
折线法从卡特兰数入手,卡特兰数是没有上界要求的,所以只需要对触碰到y=−1y=-1y=−1的情况去掉即可。在卡特兰数中,把折线从第一个碰到y=−1y=-1y=−1的位置开始翻折其实是一种一一映射的想法,对于所有不合法的方案都可以映射到一条到(2n,−2)(2n,-2)(2n,−2)的折线,而每一条这样的折线有可以唯一映射回一条不合法的路径。在有上界的情况,也是用翻折+映射的办法,但是需要再用到容斥。一个比较普通的想法是把碰到y=my=my=m的也按同样的方法翻折然后计数,但是这样会把即碰到y=−1y=原创 2022-03-19 19:39:55 · 785 阅读 · 0 评论 -
拓扑序下求原排列和逆排列的字典序最大/小的等价性
在PPT的Page8可以看到“排列在拓扑序的限制下,原排列和逆排列的字典序最大/小的等价性”这里证明一下假设现在有一个DAG,容易知道这个DAG的任意一个拓扑序都是排列,不妨假设要最大化拓扑序。把这个排列记作PPP也就是要maxP\max PmaxP,记P−1P^{-1}P−1表示PPP的逆排列,逆排列定义为∀1≤i≤n,P−1[P[i]]=i\forall 1\le i\le n,P^{-1}[P[i]]=i∀1≤i≤n,P−1[P[i]]=i。不妨用归纳法证明,显然∣V∣=1|V|=1∣V∣=原创 2022-03-19 19:38:04 · 100 阅读 · 0 评论 -
莫队二次离线学习笔记
莫队二次离线一般套路还是正常的莫队,但是如果在移动l,rl,rl,r的同时计算答案复杂度会炸掉,那么考虑莫队二次离线,当然要用的话还要满足答案能够差分。下面是算法的流程。以rrr向右移动为例,当rrr移向r+1r+1r+1时,其实就是把r+1r+1r+1在[l,r][l,r][l,r]上的贡献加入,不妨把这个记作([l,r],r+1)([l,r],r+1)([l,r],r+1),更一般的([l,r],x)([l,r],x)([l,r],x)表示xxx在区间[l,r][l,r][l,r]上的贡献,根据原创 2022-03-19 19:36:03 · 194 阅读 · 0 评论 -
拉格朗日反演学习笔记
拉格朗日反演证明对于两个函数f(x),g(x)f(x),g(x)f(x),g(x),若f(g(x))=xf(g(x))=xf(g(x))=x且f(x),g(x)f(x),g(x)f(x),g(x)的常数项为000,一次项不为000,那么称这两个函数互为复合逆。当然,根据f(g(x))=xf(g(x))=xf(g(x))=x可以推出g(f(x))=xg(f(x))=xg(f(x))=x,证明:f(g(f(x)))=f(x)f(g(f(x)))=f(x)f(g(f(x)))=f(x),令y=f(x)y=f(原创 2022-03-19 19:35:23 · 386 阅读 · 0 评论 -
二项式反演学习笔记
二项式反演引入记AcA^cAc表示AAA的补集。已知简单的容斥原理为:∣⋃i=1nAi∣=∑p1<p2<⋯<pm∣⋂j=1mApj∣(−1)m+1|\bigcup_{i=1}^{n} A_i|=\sum_{p_1<p_2<\dots<p_m}|\bigcap_{j=1}^{m}A_{p_j}|(-1)^{m+1}∣i=1⋃nAi∣=p1<p2<⋯<pm∑∣j=1⋂mApj∣(−1)m+1根据集合交并补的意义可知:∣⋂i=原创 2022-03-19 19:34:41 · 279 阅读 · 0 评论 -
杜教筛学习笔记
杜教筛杜教筛用于计算函数的前缀和。若需要计算前缀和的函数为fff,前缀和函数为SSS,那么有S(n)=∑i=1nf(i)S(n)=\sum_{i=1}^{n}f(i)S(n)=∑i=1nf(i)。现在通过构造一个合适的函数ggg,优化其计算过程。对于任意和两个函数f,gf,gf,g,以下均成立:∑i=1n∑d∣ig(d)f(id)=∑i=1ng(i)∑i=1⌊ni⌋f(i)=∑i=1ng(i)S(⌊ni⌋)g(1)S(n)=∑i=1n(g∗f)(i)−∑i=2ng(i)S(⌊ni⌋)\sum_{原创 2022-03-19 19:33:52 · 96 阅读 · 0 评论 -
单纯形学习笔记
单纯形一般形式单纯形有如下形式:max∑i=1naixi∀1≤i≤m,∑j=1nbi,jxj≤ci∀1≤i≤n,xi≥0\max \sum_{i=1}^{n} a_ix_i\\\forall 1\le i\le m,\sum_{j=1}^{n}b_{i,j}x_j\le c_i\\\forall 1\le i\le n,x_i \ge 0maxi=1∑naixi∀1≤i≤m,j=1∑nbi,jxj≤ci∀1≤i≤n,xi≥0不妨引入mmm个变量xn+1,…,xn+mx_{n原创 2022-03-19 19:33:19 · 104 阅读 · 0 评论 -
从单纯形联想到的一种简单线性规划的费用流做法
从单纯形联想到的一种简单线性规划的费用流做法使用条件观察一下单纯形中各项系数,假如对于所有xix_ixi,其在线性规划中的系数只有+1+1+1或−1-1−1,并且两者至多出现一次,那么就可以用费用流来解决。当然,如果在本身的线性规划中并不成立,但是在对偶问题中成立,也是可以做的。使用方法用费用流中的一个点来描述一个限制,一个限制形如:∑i=1naixi[≥/=/≤]c\sum_{i=1}^{n}a_ix_i[\ge/=/\le]c∑i=1naixi[≥/=/≤]c(中间的符号三选一);用一原创 2022-03-19 19:32:28 · 198 阅读 · 0 评论 -
Min_25筛学习笔记
Min_25筛学习笔记一般套路前提和目的对一个积性函数f(x)f(x)f(x),求其前缀和函数S(n)=∑i=1nf(i)S(n)=\sum_{i=1}^{n}f(i)S(n)=∑i=1nf(i),并且f(p)f(p)f(p)是一个关于ppp的多项式,f(pk)f(p^k)f(pk)能够快速计算。实现令PiP_iPi表示第iii个质数,PPP是质数集。定义如下函数:s(n,j)=∑i=2nf(n)[minp∣ip>Pj]g(n,j)=∑i=2nf(i)[i∈P∨minp∣ip>原创 2022-03-19 19:31:51 · 82 阅读 · 0 评论