专栏文章

目光呆滞

个人记录参与者 1已保存评论 0

文章操作

快速查看文章及其快照的属性,并进行相关操作。

当前评论
0 条
当前快照
1 份
快照标识符
@mipgmrmd
此快照首次捕获于
2025/12/03 11:41
3 个月前
此快照最后确认于
2025/12/03 11:41
3 个月前
查看原文
来点 Adjacent Binomial Coefficients 的无组合意义做法。
直接起手一个二元 GF fn(x,y)f_n(x,y),其中 xiyjx^iy^j 项的含义是长为 nn 的数列,和为 ii,最后一项为 jj 的答案。边界条件是 f0(x,y)=11xyf_0(x,y)=\frac{1}{1-xy},我们最后需要求的就是 [xmyk]fn(x,y)[x^my^k]f_n(x,y)
考虑转移 xiyjxi+kykx^iy^j\rightarrow x^{i+k}y^k 的系数是 (j+kj)\binom{j+k}{j},这不禁让我们想到 [yk]1(1y)j+1[y^k]\frac{1}{(1-y)^{j+1}}。凑一下即可得到下面的式子:
fn(x,y)=fn1(x,11xy)1xyf_{n}(x,y)=\frac{f_{n-1}(x,\frac{1}{1-xy})}{1-xy}
我们发现 fnf_n 其实就是,每次令 yy 变成 11xy\frac{1}{1-xy},把迭代 1n+11\sim n+1 次的结果乘起来得到的分式。
设迭代 ii 次的结果是 AiBi\frac{A_i}{B_i},那么有:
AnBn=Bn1Bn1xAn1\frac{A_n}{B_n}=\frac{B_{n-1}}{B_{n-1}-xA_{n-1}}
为了方便直接让 An=Bn1A_n=B_{n-1}Bn=Bn1xAn1B_{n}=B_{n-1}-xA_{n-1}
然后我们可以发现,因为 An=Bn1A_n=B_{n-1},所以连乘的上一项分母和下一项分子总是能抵消的,因此最后有 fn=1Bn+1f_n=\frac{1}{B_{n+1}}
An=Bn1A_n=B_{n-1} 带入 BB 的递推式,有 Bn=Bn1xBn2B_n=B_{n-1}-xB_{n-2},边界为 B0=1,B1=1xyB_0=1,B_1=1-xy
如果 B0=B1=1B_0=B_1=1,那么 [xi]B[x^i]B 相当于是递推时选了 iin2n-2,剩下的 n2in-2i 次选了 n1n-1 的方案数乘以 (1)i(-1)^i,因此就是 (1)i(nii)(-1)^i\binom{n-i}{i}
B1=1xyB_1=1-xy 只需要再考虑 xy-xy 项的贡献即可,可以发现最后有 Bn+1=11PyQB_{n+1}=\frac{1}{1-P-yQ},其中 P,QP,Q 都是关于 xx 的多项式。
那么最后问题就是求出 [xmyk]11PyQ[x^my^k]\frac{1}{1-P-yQ},直接展开可得:
[yk]11PyQ=[yk]i0(P+yQ)i=i0(ik)PikQk=Qk(1P)k+1\begin{aligned} [y^k]\frac{1}{1-P-yQ}&=[y^k]\sum_{i\geq 0}(P+yQ)^i\\ &=\sum_{i\geq 0}\binom{i}{k}P^{i-k}Q^k\\ &=\frac{Q^k}{(1-P)^{k+1}} \end{aligned}
因此使用多项式 ln,exp\ln,\exp 即可做到 O(mlogm)O(m\log m)

评论

0 条评论,欢迎与作者交流。

正在加载评论...