用到了好多技巧的状压DP
我们先统计总数然后除以m的阶乘就可以了
设f[i]表示状态为i的集合造成的贡献数（也就是状态为i的集合 不与集合外的点联通 且 这个集合联通块数是1 的情况数）
不与集合外的点联通的话只用考虑结合i之间连边，集合外那些点之间两边就可以啦
这个集合联通块数是1 就比较难处理了，有个技巧就是枚举包含1号的子集的f造成的贡献即可
#include<bits/stdc++.h>
#define int long long
using namespace std;
int n,m,ans;
const int mod=998244353,N=100010;
int jc[N],inv[N],cnt1[N],cnt2[N],tong1[1<<18],tong2[1<<18],f[1<<18];
int lowbit(int x){return x&(-x);}
int ksm(int a,int b,int mod)
{
int res=1;
for(;b;b>>=1,a=a*a%mod)
if(b&1)res=res*a%mod;
return res;
}
void YYCH()
{
jc[0]=jc[1]=inv[0]=inv[1]=1;
for(int i=2;i<=m;++i)jc[i]=jc[i-1]*i%mod;
inv[m]=ksm(jc[m],mod-2,mod);
for(int i=m-1;i>=1;--i)inv[i]=inv[i+1]*(i+1)%mod;
}
int C(int n,int m)
{
if(m>n)return 0;
return jc[n]*inv[m]%mod*inv[n-m]%mod;
}
signed main()
{
cin>>n>>m;
YYCH();
for(int i=1,x;i<=m;++i)scanf("%lld",&x),++cnt1[x];
for(int i=1,x;i<=m;++i)scanf("%lld",&x),++cnt2[x];
for(int i=0;i<=(1<<n)-1;++i)
for(int j=1;j<=n;++j)
if(!(i&(1<<(j-1))))
{
tong1[i|(1<<(j-1))]=tong1[i]+cnt1[j];
tong2[i|(1<<(j-1))]=tong2[i]+cnt2[j];
}
for(int i=1;i<=(1<<n)-1;++i)
{
if(tong1[i]==tong2[i])
{
f[i]=jc[tong1[i]]*jc[m-tong1[i]]%mod;
for(int j=(i-1)&i;j;j=(j-1)&i)
if(lowbit(i)==lowbit(j))(f[i]-=f[j]*inv[m-tong1[j]]%mod*jc[tong1[i]-tong1[j]]%mod*jc[m-tong1[i]]%mod)%=mod;
//cout<<"-> "<<i<<" "<<f[i]<<endl;
(ans+=f[i])%=mod;
}
else f[i]=0;
}
((ans%=mod)+=mod)%=mod;
cout<<ans*inv[m]%mod;
return 0;
}