【NOI2018模拟4.10】随机定向

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;
}