NOI模拟 最短路

内网传送门

由于求最短路，边权为$1$，那么考虑在$bfs$树上$dp$，类似于模拟$bfs$树的生成过程。用$f[i][j]$表示当前的$bfs$树上一共有$i$个点，最后一层有$j$个点的状态的概率（$i$中包含了$j$），这里最后一层是指只考虑到对应的层，后面的连边还不确定。$dep[i][j]$代表这样状态的期望深度。我们现在假设选中某一个确定的点，也就是说它的标号确定。因为除了一号点的期望是0，其它点的期望是相同的，那么我们最后给这个点的期望乘上$\frac{n-1}{n}$就是答案了。
考虑从$f[i][j]$转移，枚举下一层有几个点，假设有$k$个点。在维护$bfs$树时，有哪些限制条件？
①这$k$个点的每个点和最后一层的$j$个点至少有一条连边，不然$bfs$不到它。
②除去$bfs$树上的$i$个点以及这$k$个点，剩下的所有点和$bfs$树上的这$i$个点不能有连边，不然就会$bfs$到这些剩下的点，导致$bfs$树最后一层不止$k$个点。
③给这$k$个点分配标号时，从$n-i-1$个点中选$k$个出来，减一是因为那个确定的点。
那么有：$$f[i+k][k]=①\ast ②\ast ③\ast f[i][j]$$
由于这一步期望深度=（上一步的期望深度+转移贡献）×这个转移的概率。
于是有：
$$dep[i+k][k]=(dep[i][j]+f[i][j]\ast 1)\ast ①\ast ②\ast ③$$

①，②，③可以预处理出来。

然后怎么统计答案呢，考虑$f[i][j]$的贡献：
①令最后一层的$j$个点和选定点连通的概率为$P$，那么贡献就是：$P\ast (dep[i][j]+f[i][j]\ast 1)$
②从这个地方开始，后面就不连通了，令最后一层的$j$个点和剩下的$n-i$个点都不联通概率为$Q$，那么贡献就是：$f[i][j]\ast Q\ast inf$。
为什么是和$n-i$个点不连通，而不是和选定点不连通？因为$f$并没有把后面的状态确定下来，所以现在是强制把后面断掉，$bfs$树就确定了。不然会算重的。

#include<bits/stdc++.h>
#define cs const
#define re register
cs int N=4e2+10,mod=998244353,oo=1e9;
int f[N][N],dep[N][N],plink[N][N],pcut[N][N],n,p;
int fac[N],ifac[N],ans=0;
inline int add(int x,int y){return x+y>=mod?x+y-mod:x+y;}
inline int dec(int x,int y){return x-y<0?x-y+mod:x-y;}
inline int mul(int x,int y){return 1ll*x*y%mod;}
inline void Inc(int &x,int y){x=(x+y>=mod?x+y-mod:x+y);}
inline int quickpow(int a,int b,int ret=1){for(;b;b>>=1,a=mul(a,a))if(b&1)ret=mul(ret,a);return ret;}
inline void prework(){
	fac[0]=1;
	for(int i=1;i<N;++i) fac[i]=mul(fac[i-1],i);
	ifac[N-1]=quickpow(fac[N-1],mod-2);
	for(int i=N-2;~i;--i) ifac[i]=mul(ifac[i+1],i+1);
	for(int i=1;i<N;++i)
		for(int j=0;j<N;++j)
			pcut[i][j]=quickpow(dec(1,p),i*j),
			plink[i][j]=quickpow(dec(1,quickpow(dec(1,p),i)),j);
}
inline int C(int n,int m){return n<m?0:mul(fac[n],mul(ifac[m],ifac[n-m]));}
int main(){
	scanf("%d%d",&n,&p),p=mul(p,quickpow(1000000,mod-2));
	prework(),f[1][1]=1,dep[1][1]=0;
	for(int i=1;i<=n-1;++i){
		for(int j=1;j<=i;++j) if(f[i][j]){
			for(int k=1;k+i+1<=n;++k){
				int tmp=mul(plink[j][k],mul(C(n-i-1,k),pcut[j][n-i-k]));
				Inc(f[i+k][k],mul(tmp,f[i][j]));
				Inc(dep[i+k][k],mul(tmp,add(dep[i][j],f[i][j])));
			}
		}
	}
	for(int i=1;i<=n-1;++i){
		for(int j=1;j<=i;++j) if(f[i][j]){
			Inc(ans,mul(plink[j][1],add(dep[i][j],f[i][j])));
			Inc(ans,mul(mul(pcut[j][n-i],f[i][j]),oo%mod));
		}
	}printf("%d\n",mul(mul(ans,n-1),quickpow(10,6*n*n)));
}