题目:https://www.acwing.com/activity/content/11/可能需要报课才能做
有下面这样的一个网格棋盘,a,b,c,d 表示了对应边长度,也就是对应格子数。
当 a=b=c=d=2 时,对应下面这样一个棋盘:
要在这个棋盘上放 k 个相互不攻击的车,也就是这 k 个车没有两个车在同一行,也没有两个车在同一列,问有多少种方案。
只需要输出答案 mod100003 后的结果。
输入格式:
共一行,五个非负整数 a,b,c,d,k。
输出格式:
包括一个正整数,为答案 mod100003 后的结果。
数据范围:
0≤a,b,c,d,k≤1000,
保证至少有一种可行方案。
题意:给一个这种表格,然后在里面填入k个车,任何俩个车不能在同一行或者同一列,问有多少种填法。
题解:1. 如果这个表格是一个n*m的长方形,那么填法就是现在n行里面任选k行就是C(n,k),然后每一行里面选一列,每选完一列其他这个列就不能再填了,即P(m,k),所以方案数就是C(n,k)*P(m,k)。
2.现在这个表格是这种形状的,那么就将这个表格切分为上下俩部分,分别是俩个矩形,先填上面那部分,然后枚举在上面填0-k个,下面对应就是填k-0个,方案数就是每种情况的和,当上面那个矩形用了i列以后,下面那个矩形的可选列数就要减去i。
#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int mod=100003;
const int N=1e6+7;
int fact[N],infact[N];
int qmi(int a,int b){//求逆元
int res=1;
while(b){
if(b&1)res=(ll)res*a%mod;
a=(ll)a*a%mod;
b>>=1;
}
return res;
}
int P(int a,int b){
if(a<b)return 0; //控制枚举的时候,i比高度要大的情况
return (ll)fact[a]*infact[a-b]%mod;
}
int C(int a,int b){
if(a<b)return 0;
return (ll)fact[a]*infact[a-b]%mod*infact[b]%mod;
}
int main()
{
fact[0]=infact[0]=1; //求出每一个阶乘和其逆元
for(int i=1;i<N;i++){
fact[i]=(ll)fact[i-1]*i%mod;
infact[i]=qmi(fact[i],mod-2)%mod;
}
int a,b,c,d,k;
cin>>a>>b>>c>>d>>k;
int res=0;
for(int i=0;i<=k;i++){ //根据公式求出答案
res=(res+(ll)C(b,i)*P(a,i)%mod*C(d,k-i)%mod*
P(a+c-i,k-i)%mod)%mod;
}
cout<<res<<endl;
return 0;
}