Luogu P4321 随机漫游

期望DP要倒着推

Luogu P4321

---

题意

LOJ #2542

不一定是树，询问点不一定均为1

---

$Solution$

设计一个巧妙的DP状态

设$ F(S,x)$表示当前在点$ x$已经走遍了$ S$，走完剩下所有点的期望步数

这样推转移$ DP$的时候一定是从$ F(S|y,y)$转移过来

容易发现$ S|y$->$S$是不可能会变大的，即这维不可能成环

因此从大到小枚举$ S$，对当前$ S$，显然比$ S$大的状态已经被计算，暴力$ n^3$高斯消元消出这维就好了

时间复杂度$ O(2·n·n^3)$

但

为什么不直接$ Min-Max$容斥呢

随机游走

每次枚举$ S$消出当前$ Min(S,x)$中的x这维

复杂度不变

而且很好写啊...

---

$ my \ code$

#include<ctime>
#include<cmath>
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<queue>
#include<vector>
#define p 998244353
#define rt register int
#define ll long long
using namespace std;
inline ll read(){
    ll x=;char zf=;char ch=getchar();
    while(ch!='-'&&!isdigit(ch))ch=getchar();
    if(ch=='-')zf=-,ch=getchar();
    while(isdigit(ch))x=x*+ch-'',ch=getchar();return x*zf;
}
void write(ll y){if(y<)putchar('-'),y=-y;if(y>)write(y/);putchar(y%+);}
void writeln(const ll y){write(y);putchar('\n');}
int k,m,n,x,y,z,cnt,S;
bool link[][];int a[][],d[];
int inv[],ans[][<<];
int ksm(int x,int y=p-){
    int ans=;
    for(rt i=y;i;i>>=,x=1ll*x*x%p)if(i&)ans=1ll*ans*x%p;
    return ans;
}
void gauss(){
    for(rt i=;i<=n;i++)if(!(S>>i-&)){
        for(rt j=i;j<=n;j++)if(a[j][i]){
            if(j!=i)for(rt k=;k<=n+;k++)swap(a[i][k],a[j][k]);
            break;
        }
        if(!a[i][i])continue;
        inv[i]=ksm(a[i][i]);
        for(rt j=i+;j<=n;j++)if(a[j][i]){
            const int x=1ll*a[j][i]*inv[i]%p;
            for(rt k=i;k<=n+;k++)if(a[i][k])(a[j][k]-=1ll*a[i][k]*x%p)%=p;
        }
    }
    for(rt i=n;i>=;i--)if(!(S>>i-&)){
        ans[i][S]=1ll*a[i][n+]*inv[i]%p;
        for(rt j=i-;j>=;j--)if(a[j][i])(a[j][n+]-=1ll*ans[i][S]*a[j][i]%p)%=p;
    }
}
int ret[][<<];
int main(){
    n=read();m=read();
    for(rt i=;i<=m;i++){
        x=read();y=read();d[x]++;d[y]++;
        link[x][y]=link[y][x]=;
    }
    //Max(S)走遍集合S的时间
    //Min(S)第一次走到S的时间
    for(S=;S<(<<n);S++){
        memset(a,,sizeof(a));
        for(rt i=;i<=n;i++){
            if(S>>i-&)a[i][i]=;
            else {
                for(rt j=;j<=n;j++)if(link[i][j]&&!(S>>j-&))a[i][j]=-;
                a[i][i]=d[i];a[i][n+]+=d[i];
            }
        }
        gauss();
    }
    //Max(S)=sigma Min(T) (-1)|T|+1
    for(rt i=;i<=n;i++){
        for(rt j=;j<(<<n);j++)if(!(__builtin_popcount(j)&))ans[i][j]*=-;
        for(rt j=;j<n;j++)
        for(rt k=;k<(<<n);k++)if(k>>j&)(ans[i][k]+=ans[i][k^(<<j)])%=p;
    }
    for(rt T=read();T;T--){
        n=read();S=;
        for(rt i=;i<=n;i++)S|=(<<read()-);
        x=read();writeln((ans[x][S]+p)%p);
    }
    return ;
}