车的放置acwing(放置东西,不能在同一行或同一列&&逆元预处理求组合计数)

题目:https://www.acwing.com/activity/content/11/可能需要报课才能做
有下面这样的一个网格棋盘,a,b,c,d 表示了对应边长度,也就是对应格子数。
1.png
当 a=b=c=d=2 时,对应下面这样一个棋盘:
2.png
要在这个棋盘上放 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;
}

  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值