这他妈什么玩意儿……
转化模型，对于第(i)句话：“第(p)句话为真话”，将(i)，(p)连一条白边，“第(p)句话为假话”，将(i)，(p)连一条黑边。
显然我们的图会是一片基环树森林，并且边为无向边，白边连的两点真假相同，黑边相反。
那么要使森林符合要求，每个基环树的环上必定有偶数个黑边。（证明口胡如下）
我们考虑树（没有环），对于一棵树无论如何一定是满足条件的，（这显然）所以重点在环上。
对于黑边，我们可以看做将一个环分为了几段，每一段中点真假性相同，相邻的段真假性不同。
所以要在环上有偶数个黑边。
考虑一颗有(i)条白边，(j)条黑边的基环树方案数。
先枚举环上的黑白边条数：(a)条白边，(b)条黑边，设(n=a+b,m=i+j-a-b)。
每个点都有编号，所以环上的方案数为((n-1)!)，而环上加点的方案数为(n*(n+m)^{m-1})。（实际我也不会证，当结论记着吧……）
我们预处理一个(g[i][j])
[g[i][j]=left{egin{array}{ll}(i-1)!*i*(i+j)^{j-1}&i>0\(i-1)!&i=0 end{array}
ight.
]
所以(f[i][j])就有：
[f[i][j]=sum_{a+b<=i+j,!b&1}{i choose a}{j choose b}g[a+b][i+j-a-b]
]
最后统计答案时需要钦定(1)的位置，不然会计算重复。
#include<bits/stdc++.h>
using namespace std;
#define int long long
inline int read()
{
int f=1,w=0;char x=0;
while(x<'0'||x>'9') {if(x=='-') f=-1; x=getchar();}
while(x!=EOF&&x>='0'&&x<='9') {w=(w<<3)+(w<<1)+(x^48);x=getchar();}
return w*f;
}
const int N=100;
const int p=998244353;
int Blk,Wht,n;
int f[N][N],g[N][N],C[N][N],ans[N][N];
inline int ksm(int b,int k)
{
int a=b;if(!k) return 1;k--;
while(k)
{
if(k&1) a=(a%p*b%p)%p;
b=(b%p*b%p)%p;
k>>=1;
}
return a%p;
}
main(){
#ifndef ONLINE_JUDGE
//freopen("A.in","r",stdin);//Ans=1
//freopen("B.in","r",stdin);//Ans=1154
freopen("C.in","r",stdin);//Ans=322173207
#endif
string s;cin>>s;n=s.size();
for(int i=0;i<n;i++) Wht+=s[i]-'0';Blk=n-Wht;
for(int i=0;i<=n;i++) C[i][0]=1;
for(int i=1;i<=n;i++)
for(int j=1;j<=n;j++)
C[i][j]=(C[i-1][j]%p+C[i-1][j-1]%p+p)%p;
g[1][0]=1;for(int i=2;i<=n;i++) g[i][0]=(g[i-1][0]%p*(i-1)%p+p)%p;
for(int i=1;i<=n;i++)
for(int j=1;j<=n-i;j++)
g[i][j]=(g[i][0]%p*i%p*ksm(i+j,j-1)%p+p)%p;
for(int i=0;i<=Wht;i++)
for(int j=0;j<=Blk;j++)
for(int a=0;a<=i;a++)
for(int b=0;b<=j;b+=2)
if((a|b)&&a+b<=i+j)
f[i][j]=(f[i][j]%p+(C[i][a]%p*C[j][b]%p*g[a+b][i+j-a-b]%p)%p+p)%p;
ans[0][0]=1;
for(int i=0;i<=Wht;i++)
for(int j=0;j<=Blk;j++)
if(i|j)
{
if(i) for(int x=1;x<=i;x++)
for(int y=0;y<=j;y++)
ans[i][j]=(ans[i][j]%p+ans[i-x][j-y]%p*C[i-1][x-1]%p*C[j][y]%p*f[x][y]%p+p)%p;
else for(int y=1;y<=j;y++)
ans[i][j]=(ans[i][j]%p+ans[i][j-y]%p*C[j-1][y-1]%p*f[i][y]%p+p)%p;
}
printf("%lld",ans[Wht][Blk]);
}