<div class="solution-view">
    <div class="row">
        <div class="col-md-6">
            <p>提交时间：2025-06-16 19:18:57</p>
        </div>
        <div class="col-md-6">
            <p>运行 ID: <a href="/solution/detail?id=38090">38090</a></p>
        </div>
    </div>
    <div class="pre"><p>#include&lt;bits/stdc++.h&gt;
using namespace std;
const int N=20,MS=35000,mod=998244353;
void Add(int &amp;a,int b){a+=b;if(a&gt;=mod) a-=mod;}
void Sub(int &amp;a,int b){a-=b;if(a&lt;0) a+=mod;}
void Mul(int &amp;a,int b){a=1ll*a*b%mod;}
int f[N][MS],tf[N][MS];//count-tmp
int g[N][MS],tg[N][MS];//num-tmp
int h[N][MS],dc[N][N];//decval
int a[N][N];
int n;
bool next_subset(int &amp;a,int b){
    if(!a) return 0;
    a=(a-1)&amp;b;return 1;
}
/*
int solve(int x=0,int y=0){
    if(x&amp;&amp;y) if(!a[x][y]) return 0;
    memset(f,0,sizeof(f));
    memset(t,0,sizeof(t));
    a[x][y]=0;
    for(int i=1;i&lt;=n;i++) t[i][0]=1;
    int U=(1&lt;&lt;n)-1;
    for(int i=1;i&lt;=U;i++){
        for(int j=1;j&lt;=n;j++){
            if(x&amp;&amp;y)
                if(j==y&amp;&amp;((i&gt;&gt;(x-1))&amp;1)) continue;
            if(!((i&gt;&gt;(j-1))&amp;1)) continue;
            //son=1
            int S=i^(1&lt;&lt;(j-1));
            //Add(f[j][i],2ll*t[j][S]%mod);
            //son=2
            int T=S;
            do{
                int ST=S^T;
                //if(!T||!ST) continue;
                if(x&amp;&amp;y) if(((T&gt;&gt;(x-1))&amp;1)&amp;&amp;((ST&gt;&gt;(y-1))&amp;1)) continue;
                if(x&amp;&amp;y) if(((T&gt;&gt;(y-1))&amp;1)&amp;&amp;((ST&gt;&gt;(x-1))&amp;1)) continue;
                Add(f[j][i],1ll*t[j][T]*t[j][ST]%mod);
            }while(next_subset(T,S));
            for(int k=1;k&lt;=n;k++){
                if((i&gt;&gt;(k-1))&amp;1) continue;
                if(a[k][j]) Add(t[k][i],f[j][i]);
            }
            //cout&lt;&lt;i&lt;&lt;&#039; &#039;&lt;&lt;j&lt;&lt;&#039; &#039;&lt;&lt;f[j][i]&lt;&lt;endl;
        }
    }
    a[x][y]=1;
    return f[1][U];
}
*/

int main(){
    cin&gt;&gt;n;
    for(int i=1;i&lt;=n;i++){
        for(int j=1;j&lt;=n;j++){
            cin&gt;&gt;a[i][j];
        }
    }
    int U=(1&lt;&lt;n)-1;
    memset(f,0,sizeof(f));
    memset(tf,0,sizeof(tf));
    for(int i=1;i&lt;=n;i++) tf[i][0]=1;
    for(int i=0;i&lt;=U;i++){
        for(int j=1;j&lt;=n;j++){
            if(!((i&gt;&gt;(j-1))&amp;1)) continue;
            //son=1
            int S=i^(1&lt;&lt;(j-1));
            //Add(f[j][i],2ll*t[j][S]%mod);
            //son=2
            int T=S;
            do{
                int O=S^T;
                Add(f[j][i],1ll*tf[j][T]*tf[j][O]%mod);
            }while(next_subset(T,S));
            for(int k=1;k&lt;=n;k++){
                if((i&gt;&gt;(k-1))&amp;1) continue;
                if(a[k][j]) Add(tf[k][i],f[j][i]);
            }
        }
    }
    memset(g,0,sizeof(g));
    memset(tg,0,sizeof(tg));
    g[1][U]=1;
    for(int i=U;i&gt;=1;i--){
        for(int j=1;j&lt;=n;j++){
            if((i&gt;&gt;(j-1))&amp;1) continue;
            int S=U^i^(1&lt;&lt;(j-1)),T=S;
            do{
                int A=i|T|(1&lt;&lt;(j-1));
                Add(tg[j][i],2ll*tf[j][T]*g[j][A]%mod);
            }while(next_subset(T,S));
            //cout&lt;&lt;&quot;tg:&quot;&lt;&lt;j&lt;&lt;&#039; &#039;&lt;&lt;i&lt;&lt;&#039; &#039;&lt;&lt;tg[j][i]&lt;&lt;endl;
            for(int k=1;k&lt;=n;k++){
                if(!((i&gt;&gt;(k-1))&amp;1)) continue;
                if(a[j][k]) Add(g[k][i],tg[j][i]);
                //cout&lt;&lt;&quot;to:&quot;&lt;&lt;k&lt;&lt;&#039; &#039;&lt;&lt;tg[j][i]&lt;&lt;endl;
            }
        }
    }
    memset(h,0,sizeof(h));
    for(int i=1;i&lt;=U;i++){
        for(int j=1;j&lt;=n;j++){
            if((i&gt;&gt;(j-1))&amp;1) continue;
            int S=U^i^(1&lt;&lt;(j-1)),T=S;
            do{
                int A=i|T|(1&lt;&lt;(j-1));
                Add(h[j][i],2ll*tf[j][T]*g[j][A]%mod);
            }while(next_subset(T,S));
            for(int k=1;k&lt;=n;k++){
                if(!((i&gt;&gt;(k-1))&amp;1)) continue;
                if(a[j][k]) Sub(dc[j][k],1ll*h[j][i]*f[k][i]%mod);
            }
        }
    }
    cout&lt;&lt;f[1][U]&lt;&lt;endl;
    for(int i=1;i&lt;=n;i++){
        for(int j=1;j&lt;=n;j++){
            int ans=dc[i][j];
            for(int k=1;k&lt;=U;k++){
                if((k&gt;&gt;(j-1))&amp;1) Add(ans,1ll*f[i][k]*g[i][k]%mod);
            }
            if(!a[i][j]) ans=0; 
            cout&lt;&lt;ans&lt;&lt;&#039; &#039;;
        }cout&lt;&lt;endl;
    }
    return 0;
}</p></div>
</div>
<script type="text/javascript">
    (function ($) {
        $(document).ready(function () {
            $('.pre p').each(function(i, block) {  // use <pre><p>
                hljs.highlightBlock(block);
            });
        })
    })(jQuery);
</script>