<div class="solution-view">
    <div class="row">
        <div class="col-md-6">
            <p>提交时间：2025-06-23 21:04:21</p>
        </div>
        <div class="col-md-6">
            <p>运行 ID: <a href="/solution/detail?id=38149">38149</a></p>
        </div>
    </div>
    <div class="pre"><p>#include&lt;bits/stdc++.h&gt;
using namespace std;
#define ll long long
#define pii pair&lt;int,int&gt;
#define fr first
#define sc second
#define mk make_pair
#define pb push_back
#define popcnt __builtin_popcount
const int N=18,Mod=998244353,M=133333;
inline int read(){int x=0,f=1;char c=getchar();while(c&lt;&#039;0&#039;||c&gt;&#039;9&#039;){if(c==&#039;-&#039;)f=-1;c=getchar();}while(c&gt;=&#039;0&#039;&amp;&amp;c&lt;=&#039;9&#039;)x=x*10+c-&#039;0&#039;,c=getchar();return x*f;}
inline void add(int &amp;x,int y){x+=y;if(x&gt;=Mod)x-=Mod;}
int n,P,ans,a[N][N],as[N][N],f[N][1&lt;&lt;N],w[N][1&lt;&lt;N],p[N][1&lt;&lt;N],g[N][1&lt;&lt;N],pt[1&lt;&lt;N];
int fp[N][N][1&lt;&lt;N],h[N][1&lt;&lt;N],gp[N][N][1&lt;&lt;N],tr[N][M],tr2[N][M];
// string ejz(int x){string res=&quot;&quot;;for(int i=0;i&lt;n;i++)res+=(x&amp;(1&lt;&lt;i))?&quot;1&quot;:&quot;0&quot;;return res;}
void slv(){
    n=read(),P=(1&lt;&lt;n)-1;
    for(int i=0;i&lt;n;i++)for(int j=0;j&lt;n;j++)a[i][j]=read();
    for(int i=0;i&lt;1&lt;&lt;n;i++)pt[i]=popcnt(i);
    for(int i=0;i&lt;n;i++)w[i][0]=fp[i][0][0]=1;
    for(int i=0;i&lt;n;i++){
        int cnt=0,cnt2=0;
        for(int j=0;j&lt;1&lt;&lt;n;j++)if(j&amp;(1&lt;&lt;i))tr[i][cnt++]=j;else tr2[i][cnt2++]=j;
    }
    for(int o=0;o&lt;n;o++){
        for(int i=0;i&lt;n;i++)
        for(int s=0;s&lt;1&lt;&lt;n-1;s++)if(pt[tr2[i][s]]==o)
        for(int j=0;j&lt;n;j++)
            if(tr2[i][s]&amp;(1&lt;&lt;j))
                if(a[i][j])add(fp[i][o][tr2[i][s]],f[j][tr2[i][s]]),add(w[i][tr2[i][s]],f[j][tr2[i][s]]);

        for(int k=0;k&lt;n;k++)for(int i=0;i&lt;n;i++)for(int j=0;j&lt;1&lt;&lt;n-1;j++)if(pt[tr[i][j]]&gt;=o)add(fp[k][o][tr[i][j]],fp[k][o][tr[i][j]^(1&lt;&lt;i)]);
        for(int i=0;i&lt;n;i++){
            for(int s=0;s&lt;1&lt;&lt;n-1;s++)if(pt[tr2[i][s]]&lt;=o){
                ll res=fp[i][o][tr2[i][s]]*(o?2:1);
                for(int k=1;k&lt;o;k++)res+=(ll)fp[i][k][tr2[i][s]]*fp[i][o-k][tr2[i][s]];
                h[i][tr2[i][s]]=res%Mod;
            }
        }
        for(int k=0;k&lt;n;k++)for(int i=0;i&lt;n;i++)for(int j=0;j&lt;1&lt;&lt;n-1;j++)if(pt[j]&lt;=o)add(h[k][tr[i][j]],Mod-h[k][tr[i][j]^(1&lt;&lt;i)]);
        for(int s=0;s&lt;1&lt;&lt;n;s++)if(pt[s]==o)
            for(int j=0;j&lt;n;j++)if(!(s&amp;(1&lt;&lt;j)))f[j][s|(1&lt;&lt;j)]=h[j][s];
    }
    // return;
    g[0][0]=gp[0][0][0]=1;
    for(int o=0;o&lt;n;o++){
        for(int i=0;i&lt;n;i++)
            for(int s=0;s&lt;1&lt;&lt;n-1;s++)if(pt[tr[i][s]]==o){
            for(int j=0;j&lt;n;j++)if(!(tr[i][s]&amp;(1&lt;&lt;j))){
                if(a[i][j])add(gp[j][o][tr[i][s]],p[i][tr[i][s]]),add(g[j][tr[i][s]],p[i][tr[i][s]]);
            }
        }
        for(int k=0;k&lt;n;k++)for(int i=0;i&lt;n;i++)for(int j=0;j&lt;1&lt;&lt;n-1;j++)add(gp[k][o][tr[i][j]],gp[k][o][tr[i][j]^(1&lt;&lt;i)]);
        for(int i=0;i&lt;n;i++){
            for(int s=0;s&lt;1&lt;&lt;n-1;s++)if(pt[tr2[i][s]]&lt;=o){
                ll res=0;
                for(int k=0;k&lt;=o;k++)res+=(ll)gp[i][k][tr2[i][s]]*fp[i][o-k][tr2[i][s]]*2;
                h[i][tr2[i][s]]=res%Mod;
            }
        }
        for(int k=0;k&lt;n;k++)for(int i=0;i&lt;n;i++)for(int j=0;j&lt;1&lt;&lt;n-1;j++)if((tr[i][j]&amp;(1&lt;&lt;i)))add(h[k][tr[i][j]],Mod-h[k][tr[i][j]^(1&lt;&lt;i)]);

        for(int s=0;s&lt;1&lt;&lt;n;s++)if(pt[s]==o)for(int j=0;j&lt;n;j++)if(!(s&amp;(1&lt;&lt;j)))p[j][s|(1&lt;&lt;j)]=h[j][s];
    }
    for(int i=0;i&lt;n;i++)for(int j=0;j&lt;n;j++)if(a[i][j]){
        for(int s=0;s&lt;1&lt;&lt;n;s++)if((s&amp;(1&lt;&lt;i))&amp;&amp;(s&amp;(1&lt;&lt;j)))add(as[i][j],(ll)f[i][s]*g[i][s^P]%Mod);
        for(int s=0;s&lt;1&lt;&lt;n;s++)if((s&amp;(1&lt;&lt;i))&amp;&amp;(s&amp;(1&lt;&lt;j)))add(as[i][j],Mod-(ll)f[j][s^(1&lt;&lt;i)]*p[i][s^P^(1&lt;&lt;i)]%Mod);
    }
    printf(&quot;%d\n&quot;,f[0][(1&lt;&lt;n)-1]);
    for(int i=0;i&lt;n;i++){
        for(int j=0;j&lt;n;j++)printf(&quot;%d &quot;,as[i][j]);
        puts(&quot;&quot;);
    }
}
signed main(){
 //  freopen(&quot;1.in&quot;,&quot;r&quot;,stdin);freopen(&quot;1.out&quot;,&quot;w&quot;,stdout);
    slv();
   //cerr&lt;&lt;(clock()*1.0)/CLOCKS_PER_SEC&lt;&lt;&quot;s&quot;&lt;&lt;endl;
    // system(&quot;grep VmPeak /proc/$PPID/status&quot;);
    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>