「XXI Opencup GP of Tokyo」 Count Min Ratio

link。


非 GF 做法:https://xyix.gitee.io/posts/?&postname=cf-gym-102978-c


简单转化问题:求所有 \(x \in [1, \lfloor R / B\rfloor]\) 对应的合法方案数之和。

对于某一 \(x\),枚举左边的红蓝球数量,可得:

\[\begin{aligned} ans_x =& \sum_{b = 0}^{B}\sum_{r = bx}^{R - (B-b)x}\binom{r + b}{b}\binom{R - r + B - b}{B - b} \\ =& \sum_{b = 0}^{B}\sum_{d = 0}^{R - Bx}\binom{(bx + d) + b}{b}\binom{R - (bx + d) + B - b}{B - b} \\ =& \sum_{d = 0}^{R - Bx}\sum_{b = 0}^{B}\binom{(x + 1)b + d}{b}\binom{(x + 1)(B - b) + (R - Bx - d)}{B - b} \\ \end{aligned} \]

以下记 \(t = x + 1\)


系数 \(\binom{tb + d}{b}\)广义二项级数 \(\mathcal B_t(z)\) 有关。

\(G(z) = \frac{z}{(1 + z)^t}\),则它的复合逆 \(F(z) = \mathcal B_t(z) - 1\)。于是:

\[\sum_{b=0}\binom{tb + d}{b}z^b = (1 + F)^d\times\frac{1 + F}{1 - (t - 1)F} \]

关于广义二项级数的详细内容可见文末。


代入原式,可得:

\[\begin{aligned} & \sum_{d = 0}^{R - Bx}\sum_{b = 0}^{B}\binom{tb + d}{b}\binom{t(B - b) + (R - Bx - d)}{B - b} \\ =& \sum_{d = 0}^{R - Bx}\sum_{b = 0}^{B}\left([z^b](1 + F)^d\times\frac{1 + F}{1 - (t - 1)F}\right)\left([z^{B-b}](1 + F)^{R - Bx - d}\times\frac{1 + F}{1 - (t - 1)F}\right) \\ =& [z^B]\sum_{d = 0}^{R - Bx}\left((1 + F)^d\times\frac{1 + F}{1 - (t - 1)F}\right)\left((1 + F)^{R - Bx - d}\times\frac{1 + F}{1 - (t - 1)F}\right) \\ =& (R - Bx + 1)\left([z^B](1 + F)^{R - Bx}\times\left(\frac{1 + F}{1 - (t - 1)F}\right)^2\right) \\ \end{aligned} \]

考虑对 \(H(F) = (1 + F)^{R - Bx}\times\left(\frac{1 + F}{1 - (t - 1)F}\right)^2\) 进行拉格朗日反演:

\[\begin{aligned}[] [z^B]H(F) =&[z^0]zH\times G'G^{-B-1} \\ =&[z^0]z(1 + z)^{R - Bx}\times\left(\frac{1 + z}{1 - (t - 1)z}\right)^2\times\frac{1 - (t - 1)z}{(1 + z)^{t + 1}}\times\left(\frac{(1 + z)^{(B + 1)t}}{z^{B + 1}}\right) \\ =&[z^B]\frac{(1 + z)^{(R - Bx) + 2 - (t + 1) + (B + 1)t}}{1 - (t - 1)z} \\ =&[z^B]\frac{(1 + z)^{R + B + 1}}{1 - xz} \\ =&\sum_{i = 0}^{B}\binom{R + B + 1}{i}x^{B - i} \end{aligned} \]

那么答案为 \((R + 1)\sum_{i = 0}^{B}\binom{R + B + 1}{i}\left(\sum_{x = 1}^{\lfloor R / B\rfloor}x^{B - i}\right) - B\sum_{i = 0}^{B}\binom{R + B + 1}{i}\left(\sum_{x = 1}^{\lfloor R / B\rfloor}x^{B - i + 1}\right)\)

多项式求逆算自然数幂和即可。


一些关于 \(\mathcal B_t(x)\) 的补充。

\(\mathcal B_t(x)\) 满足方程 \(\mathcal B_t(x) = x\mathcal B^t_t(x) + 1\)。当 \(t = 2\) 时,\(\mathcal B_2(x)\) 即卡特兰数的 OGF,那么 \(\mathcal B_t\) 其实是卡特兰数的推广。

参考卡特兰数,这个方程至少有如下几种组合意义:

(1)儿子有序的 \(t\) 叉树。

按照定义,显然有方程 \(\mathcal B_t(x) = x\mathcal B^t_t(x) + 1\)

(2)凸多边形的 \(t + 1\) 边形划分(好像去年 FJOI 考了这道题)。

选取一条凸多边形上的边,它所在的 \(t + 1\) 边形将原凸多边形划分成 \(t\) 份,因此有方程 \(\mathcal B_t(x) = x\mathcal B^t_t(x) + 1\)

(3)\((0, 0)\to(n, (t - 1)n)\) 且始终在直线 \(y = (t - 1)x\) 下方的路径数。

仿照卡特兰数,枚举第一次碰到 \(y = (t - 1)x\) 的点,那么有 \(\mathcal B_t(x) = x\mathcal B_t(x)T(x) + 1\),其中 \(T(x)\) 表示完全不碰到 \(y = (t - 1)x\) 的方案数。

考虑这条路径第一次碰到直线 \(y = (t - 1)x - i,i\in (0, t - 1)\) 的时刻(该路径一定会碰到),并按这 \(t - 1\) 个时刻分段,每一段都贡献一个 \(\mathcal B_t(x)\)。因此有 \(T(x) = \mathcal B_t^{t-1}(x)\)

感觉不是很平凡,如果有更好的构造请告诉我(用 Raney 引理直接得到封闭形式并惊讶发现它就是广义二项级数之类的就算了)

这个组合意义可以用于解决本题。

p.s:

尽管上面探讨了 \(\mathcal B_t\) 的组合意义,但是其实 \(t\) 在取任意值时下文中的公式都是成立的(此时定义为广义二项系数)。

不过 \(\mathcal B_t\) 有常数项,不方便拉格朗日反演。记 \(F = \mathcal B_t - 1\),它有复合逆 \(G = x / (1 + x)^t\)。那么:

\[\begin{aligned}[] [x^n]\mathcal B_t^r =& [x^n](1 + F)^r \\ =& [x^0]x(1 + x)^r\frac{1-(t-1)x}{(1 + x)^{t + 1}}\frac{(1+x)^{tn+t}}{x^{n+1}} \\ =& [x^n](1 + x)^{tn+r-1}(1-(t-1)x) \\ =& \binom{tn + r - 1}{n} - \binom{tn + r - 1}{n - 1}(t - 1)\\ =& \binom{tn + r}{n}\times\frac{r}{tn + r} \end{aligned} \]

最后来推导 \(\sum_{n=0}\binom{tn + r}{n}x^n = (1 + F)^r\times\frac{1 + F}{1 - (t - 1)F}\)(当然,右推左是 trivial 的):

\[\begin{aligned} \binom{tn + r}{n} &= [x^n](1 + x)^{tn + r} \\ &= [x^0]\left(\frac{(1 + x)^t}{x}\right)^n(1 + x)^r \\ &= [x^0]G^{-n}(1 + F(G))^r \\ &= [x^0]x\times x^{-n}(1 + F)^r\times \frac{F'}{F} \\ &= [x^n](1 + F)^r \frac{xF'}{F} \end{aligned} \]

之后只需再证明 \(\frac{xF'}{F} = \frac{1 + F}{1 - (t - 1)F}\)

注意到 \(G' = \frac{1 - (t - 1)x}{(1 + x)^{t + 1}}\),那么有 \(\frac{xG'}{G} = \frac{1 - (t - 1)x}{1 + x}\)(这也表明 \(G\) 是微分有限的 废话 它能不微分有限吗),由链式法则求导又有 \(F' = \frac{1}{G'(F)}\)

于是 \(\frac{xF'}{F} = \frac{G(F)}{G'(F)F} = \frac{1 + F}{1 - (t - 1)F}\)

posted @ 2021-04-30 13:13  Tiw_Air_OAO  阅读(533)  评论(2编辑  收藏  举报