Description
问至少有 A A 行列全是 1 1 的的 01 01 矩阵个数
Input
多组用例,每组用例输入四个整数 n,m,A,B(1≤n,m,A,B≤3000) n , m , A , B ( 1 ≤ n , m , A , B ≤ 3000 )
Output
输出方案数,结果模 998244353 998244353
Sample Input
3 4 1 2
Sample Output
169
Solution
先考虑行,选出 i i 行为黑、至少列为黑的方案数为 Cin⋅ C n i ⋅ 所选 i i 行为黑、至少列为黑的方案数,但在选 j<i j < i 行为黑时会重复计算,考虑容斥,假设每部分实际应该计算的次数为 fi f i 次,那么有 fi=1−∑k=ai−1fkCki f i = 1 − ∑ k = a i − 1 f k C i k ,同理可以求出选出 j j 列为黑实际被计算的次数,而所选 i i 行列为黑的方案数即为 2(n−i)(m−j) 2 ( n − i ) ( m − j ) ,进而答案即为 ans=∑i=an∑j=bmfigjCinCjm2(n−i)(m−j) a n s = ∑ i = a n ∑ j = b m f i g j C n i C m j 2 ( n − i ) ( m − j ) ,预处理组合数和 2 2 的幂次即可,方案数
Code
#include<cstdio>
using namespace std;
typedef long long ll;
#define maxn 3005
#define mod 998244353
int mul(int x,int y)
{
ll z=1ll*x*y;
return z-z/mod*mod;
}
int add(int x,int y)
{
x+=y;
if(x>=mod)x-=mod;
return x;
}
int C[maxn][maxn],f2[maxn*maxn],f[maxn],g[maxn];
void init(int n=3000)
{
f2[0]=1;
for(int i=1;i<=n*n;i++)f2[i]=add(f2[i-1],f2[i-1]);
C[0][0]=1;
for(int i=1;i<=n;i++)
{
C[i][0]=C[i][i]=1;
for(int j=1;j<i;j++)C[i][j]=add(C[i-1][j-1],C[i-1][j]);
}
}
int main()
{
init();
int n,m,a,b;
while(~scanf("%d%d%d%d",&n,&m,&a,&b))
{
for(int i=a;i<=n;i++)
{
f[i]=1;
for(int j=a;j<i;j++)f[i]=add(f[i],mod-mul(f[j],C[i][j]));
}
for(int i=b;i<=m;i++)
{
g[i]=1;
for(int j=b;j<i;j++)g[i]=add(g[i],mod-mul(g[j],C[i][j]));
}
int ans=0;
for(int i=a;i<=n;i++)
for(int j=b;j<=m;j++)
{
int temp=mul(mul(f[i],C[n][i]),mul(g[j],C[m][j]));
temp=mul(temp,f2[(n-i)*(m-j)]);
ans=add(ans,temp);
}
printf("%d\n",ans);
}
return 0;
}