【计数DP】BZOJ5219 [Lydsy2017省队十连测] 最长路径

【题目】
BZOJ
给定$n$，求一幅$n$个点的竞赛图（所有点对之间有一条有向边），对于$i=1\sim n$，从$1$出发的最长简单路径经过点数恰好为$i$的竞赛图个数模$P$。$n\le 2000,P\le 10^9$

【解题思路】
对于竞赛图有一些性质：

• 竞赛图一定有哈密尔顿路径（即沿边访问每个顶点恰好一次
• 缩点之后按拓扑序形成一条“链”
  那么起点$1$一定在这条路径上，且从$1$出发最长路径上的点数等于$1$所在$scc$点数$+$拓扑序在$1$以后的$scc$的点的总数。

令$f_n$表示$n$个节点竞赛图个数，那么$f_n=2^{\frac{n(n-1}{2}}$

令$g_n$表示有$n$个节点的竞赛图$scc$，我们枚举拓扑序最小的$scc$进行容斥，那么有：$g_n=f_n-\sum _{i=1}^{n-1}\left(\genfrac{}{}{0ex}{}{n}{i}\right)g_i{f}_{n-i}$

初始我们有$f_0=1,g_1=1$，然后递推即可。

那么我们枚举$1$所在最后有答按，令$an{s}_k$表示从$1$出发最长路径为$k$的方案数
$$an{s}_k=\sum _{i=1}^n\sum _{j=1}^n\left(\genfrac{}{}{0ex}{}{n-1}{i-1}\right)\left(\genfrac{}{}{0ex}{}{n-i}{j}\right)g_i{f}_j{f}_{n-i-j}$$
具体意义为枚举$1$中$scc$点数以及拓扑序在$1$后面的$scc$总点数。

预处理组合数可以得到一个$O(n^2)$的做法。

【参考代码】

#include<bits/stdc++.h>
using namespace std;

typedef long long ll;
const int N=2005;
int n,mod,c[N][N],f[N],g[N],ans[N];

int upm(int x){return x>=mod?x-mod:(x<0)?x+mod:x;}
void up(int &x,int y){x=upm(x+y);}
int mul(int x,int y){return 1ll*x*y%mod;}
int qpow(int x,int y){int res=1;for(;y;y>>=1,x=mul(x,x))if(y&1)res=mul(res,x);return res;}

int main()
{
#ifndef ONLINE_JUDGE
	freopen("BZOJ5219.in","r",stdin);
	freopen("BZOJ5219.out","w",stdout);
#endif
	scanf("%d%d",&n,&mod);
	for(int i=0;i<=n;++i)
	{
		c[i][0]=c[i][i]=1;
		for(int j=1;j<i;++j) c[i][j]=upm(c[i-1][j]+c[i-1][j-1]);
	}
	for(int i=0;i<=n;++i) f[i]=qpow(2,i*(i-1)/2);
	for(int i=1;i<=n;++i) 
	{
		g[i]=f[i];
		for(int j=1;j<i;++j) up(g[i],-mul(c[i][j],mul(g[j],f[i-j])));
	}
	for(int i=1;i<=n;++i) for(int j=0;j<=n-i;++j)
		up(ans[i+j],1ll*mul(c[n-1][i-1],c[n-i][j])*mul(g[i],f[j])%mod*f[n-i-j]%mod);
	for(int i=1;i<=n;++i) printf("%d\n",ans[i]);
	return 0;
}