组合数学 - 车的放置 - AcWing 1309

组合数学 - 车的放置 - AcWing 1309

有下面这样的一个网格棋盘，a,b,c,d 表示了对应边长度，也就是对应格子数。

当 a=b=c=d=2 时，对应下面这样一个棋盘：

要在这个棋盘上放 k 个相互不攻击的车，也就是这 k 个车没有两个车在同一行，也没有两个车在同一列，问有多少种方案。

只需要输出答案 mod 100003 后的结果。

输入格式

共一行，五个非负整数 a,b,c,d,k。

输出格式

包括一个正整数，为答案 mod 100003 后的结果。

数据范围

$1\le a,b,c,d,k\le 1000，保证至少有一种可行方案。$

输入样例：

2 2 2 2 2

输出样例：

38

---

分析：

先考虑规则图形:

$在下图中放入i个车。$

$先考虑行:同一行不能放，我们需要在b行中选择i行放，方案数C_b^i$

$再考虑列：同一列不能放，我们需要在a列中选择i列放，需要考虑不同的顺序，方案总数A_a^i。$

$故上图的方案总数为C_b^i\times A_a^i。$

注： $行与列有一个考虑顺序即可，即C_b^i\times A_a^i=C_a^i\times A_b^i。$

接着考虑本题：

$我们将这{个}^{\prime \prime }{L}^{\prime \prime }型以红线一份为二，分别计算上下两部分的方案总数，再利用乘法原理求整体的方案总数。$

$我们枚举上半部分放i个车，$

$因为上半部分对下半部分的影响是确定的，而放在下半部分，可能会影响上半部分，也可能不影响。$

$①、上半部分方案总数：C_b^i\times A_a^i$

$②、下半部分方案总数：C_d^{k-i}\times A_{a+c-i}^{k-i}$

$整体的方案总数：C_b^i\times A_a^i\times C_d^{k-i}\times A_{a+c-i}^{k-i}$

本题直接定义法求组合数。

代码：

#include<iostream>

#define ll long long

using namespace std;

const int N=2010, mod=100003;

int fact[N],infact[N];
int a,b,c,d,k;

int quick_pow(int a,int b,int mod)
{
    int res=1;
    while(b)
    {
        if(b&1) res=(ll)res*a%mod;
        a=(ll)a*a%mod;
        b>>=1;
    }
    return res;
}

void Init()
{
    fact[0]=infact[0]=1;
    for(int i=1;i<N;i++)
    {
        fact[i]=(ll)fact[i-1]*i%mod;
        infact[i]=quick_pow(fact[i],mod-2,mod);  //等价于infact[i-1]*i^(mod-2)
    }
}

int C(int a,int b)
{
    if(a<b) return 0;
    return (ll)fact[a]*infact[a-b]%mod*infact[b]%mod;
}

int A(int a,int b)
{
    if(a<b) return 0;
    return (ll)fact[a]*infact[a-b]%mod;
}

int main()
{
    Init();
    
    cin>>a>>b>>c>>d>>k;
    
    int res=0;
    for(int i=0;i<=k;i++)
        res=(res+(ll)C(b,i)*A(a,i)%mod*C(d,k-i)*A(a+c-i,k-i)%mod)%mod;
        
    cout<<res<<endl;
    
    return 0;
}