Description
给出一个n个点的竞赛图,其中每条边(ui,vi)有pi的概率是(ui,vi),有1-pi的概率是(vi,ui)。
给出m条关键边,关键边的pi给出,其余边的pi均为0.5,求这个竞赛图的期望强联通分量个数。
n<=38,m<=19
Solution
n这么大经典的套路显然是做不了了,我们不妨换个思路
众所周知,竞赛图缩点之后的拓扑序是唯一的,而且强联通分量的个数就等于”间隙”的个数+1
那么我们不妨统计间隙的个数。
这个很好统计,只需要枚举被分成哪两个集合,然后把中间边的概率乘起来,最后累加。
这样我们得到了一个O(2^n)的做法,发现m很小,说明原图中特殊边所形成的连通块的大小最大为20
于是我们可以对每个连通块分别考虑,最后用一个卷积卷起来就好了。
Code
#include <cstdio>
#include <cstring>
#include <algorithm>
#define fo(i,a,b) for(int i=a;i<=b;i++)
#define fd(i,a,b) for(int i=a;i>=b;i--)
#define rep(i,a) for(int i=lst[a];i;i=nxt[i])
using namespace std;
typedef long long ll;
const int N=40,Mo=998244353;
int pwr(int x,int y) {
int z=1;
for(;y;y>>=1,x=(ll)x*x%Mo)
if (y&1) z=(ll)z*x%Mo;
return z;
}
int t[N<<1],nxt[N<<1],lst[N],l;
void add(int x,int y) {t[++l]=y;nxt[l]=lst[x];lst[x]=l;}
int n,m,tot,x,y,z,P,a[N],b[N],c[N],f[N],g[N],h[N],p[N][N];
bool vis[N],bz[N];
void dfs(int x) {
vis[x]=1;a[++tot]=x;
rep(i,x) if (!vis[t[i]]) dfs(t[i]);
}
void calc(int x,int y) {
if (x>tot) {
b[0]=c[0]=0;
fo(i,1,tot)
if (bz[i]) b[++b[0]]=a[i];
else c[++c[0]]=a[i];
int res=1;
fo(i,1,b[0])
fo(j,1,c[0])
res=(ll)res*p[b[i]][c[j]]%Mo;
(g[y]+=res)%=Mo;
return;
}
bz[x]=1;calc(x+1,y+1);
bz[x]=0;calc(x+1,y);
}
int main() {
freopen("random.in","r",stdin);
freopen("random.out","w",stdout);
scanf("%d%d",&n,&m);
int inv=pwr(10000,Mo-2);P=(ll)5000*inv%Mo;
fo(i,1,n) fo(j,1,n) if (i!=j) p[i][j]=1;
fo(i,1,m) {
scanf("%d%d%d",&x,&y,&z);
z=(ll)z*inv%Mo;
add(x,y);add(y,x);
p[x][y]=z*2%Mo;p[y][x]=(1+Mo-z)*2%Mo;
}
f[0]=1;
fo(i,1,n)
if (!vis[i]) {
tot=0;
dfs(i);
fo(j,0,n) g[j]=0;
calc(1,0);
fo(j,0,n) {
h[j]=0;
fo(k,0,j) (h[j]+=(ll)f[k]*g[j-k]%Mo)%=Mo;
}
fo(j,0,n) f[j]=h[j];
}
int ans=0;
fo(i,1,n) (ans+=(ll)f[i]*pwr(P,i*(n-i))%Mo)%=Mo;
ans=(ll)ans*pwr(10000,n*(n-1))%Mo;
printf("%d\n",ans);
return 0;
}