[CF848D]Shake It!

题目

传送门 to CF

思路

说出来你可能不信，我在考场上的思路是正解，但是复杂度算错了根本没打……

首先最开始只有一种操作，就是 $s\to w\to t$ 。然后你会发现，这两条边之上建造的图是独立的。所以肯定考虑 $\mathtt{d}\mathtt{p}$ 啊！

设 $f(n,m)$ 为所求的答案，$g(n,m)$ 为 “两截式” 的方案数（即 $s\to w\to t$ 之上建造 $n-1$ 条点），则
$$g(n,m)=\sum _{x_1+x_2=n-1}\sum _{\min (y_1,y_2)=m}f(x_1,y_1)\cdot f(x_2,y_2)$$
如果 $f$ 已知，预处理后缀和，枚举 $n,m,x_1$ 即可求出 $g(n,m)$，复杂度 $\mathcal{O}(n^3)$ 。现在考虑怎么用 $g$ 求出 $f$ 呢？

很容易受到 $e^{F(x)}$ 这种想法的干扰。可是这里并不是那样的。$e^{F(x)}$ 往往是 有标号 的指数型生成函数，因为它是依次分配 $id$，所以交换分配 $id$ 的顺序也可能分配出相同的 $id$，进而导致重复。

这里没有什么好的方法，干脆 组合数硬算。让 $g(x,y)$ 对应的方案中选出 $c_{x,y}$ 个（可重复选择），显然方案数是隔板法 $\left(\genfrac{}{}{0ex}{}{g+c-1}{g-1}\right)$，于是
$$f(n,m)=\sum _{⟨c⟩}\left[\prod \left(\genfrac{}{}{0ex}{}{g(x,y)+c_{x,y}-1}{g(x,y)-1}\right)\right]$$
当然要满足 $\sum c_{x,y}\cdot x=n$ 与 $\sum c_{x,y}\cdot y=m-1$，这里 $y$ 缺的 $1$ 就是 $s\to t$ 直接提供 $1$ 的流量。

考虑生成函数。记 $G(n,m)={\sum }_{t=0}^{+\infty }\left(\genfrac{}{}{0ex}{}{g(n,m)+t-1}{t}\right)\cdot x^{tn}\cdot y^{tm}$，则 $f(n,m)=[x^n{y}^{m-1}]\prod G(a,b)$ 。注意这里 $a,b$ 可以是没有任何范围的。

现在请问，所有 $G(n,m)$ 的项数之和是什么？我粗略计算了一下，是 $(n\ln n)(m\ln m)$，两个都是调和级数。但是实际上它是
$$\sum \min \left(\frac{n}{a},\frac{m}{b}\right)⩽\sum _{a,b}\frac{n}{a}+\sum _{a,b}\frac{m}{b}=nm(\ln n+\ln m)$$
如果进行暴力卷积，每一项乘 $n^2$ 个数，复杂度 $\mathcal{O}(n^4\ln n)$ 。感觉过不了，然后它过了……我也想不出什么方法优化它。

代码

#include <cstdio>
#include <iostream>
#include <cstring>
#include <vector>
#include <algorithm>
using namespace std;
typedef long long int_;
# define rep(i,a,b) for(int i=(a); i<=(b); ++i)
# define drep(i,a,b) for(int i=(a); i>=(b); --i)
inline int readint(){
	int a = 0, c = getchar(), f = 1;
	for(; '0'>c||c>'9'; c=getchar())
		if(c == '-') f = -f;
	for(; '0'<=c&&c<='9'; c=getchar())
		a = (a<<3)+(a<<1)+(c^48);
	return a*f;
}
inline void writeint(int x){
	if(x > 9) writeint(x/10);
	putchar((x-x/10*10)^48);
}

const int MAXN = 205;
const int Mod = 1e9+7;
int f[MAXN][MAXN], suffix_sum[MAXN][MAXN];
int g[MAXN][MAXN]; // product of all G
int tmp[MAXN][MAXN], inv[MAXN];

int main(){
	rep(i,(inv[1]=1)<<1,MAXN-1)
		inv[i] = int(int_(Mod-Mod/i)*inv[Mod%i]%Mod);
	int n = readint(), m = readint();
	f[0][1] = 1; // no operation, 1 flow
	suffix_sum[0][1] = 1; // maybe useful
	int len_x = 0, len_y = 0;
	g[0][0] = 1; // init G = 1
	rep(i,1,n){ // operation times
		rep(j,1,i+1){ // amount of flow
			int now = 0; // current g[i][j] (DP)
			rep(x,0,i-1) // enumerate one of the parts
				now = static_cast<int>((now // accumulate
					+int_(f[x][j])*suffix_sum[i-1-x][j+1]
					+int_(f[i-1-x][j])*suffix_sum[x][j])%Mod);
			int p = 0, q = 0; // exponent
			for(int c=1,t=0; p<=n&&q<=n; p+=i,q+=j,++t){
				rep(x,0,len_x) rep(y,0,len_y) tmp[x+p][y+q] =
					int((tmp[x+p][y+q]+int_(c)*g[x][y])%Mod);
				c = int(int_(now+t)*inv[t+1]%Mod*c%Mod);
			}
			p -= i, q -= j; // the last tranfer
			len_x = min(len_x+p,n), len_y = min(len_y+q,n);
			rep(x,0,len_x) rep(y,0,len_y) // copy back
				g[x][y] = tmp[x][y], tmp[x][y] = 0;
			f[i][j] = g[i][j-1]; // transfer
		}
		suffix_sum[i][i+1] = f[i][i+1];
		drep(j,i,1) suffix_sum[i][j] =
			(suffix_sum[i][j+1]+f[i][j])%Mod;
	}
	printf("%d\n",f[n][m]);
	return 0;
}