UPC-排课表+玉米田(容斥原理+组合数学公式)

简介: UPC-排课表+玉米田(容斥原理+组合数学公式)

排课表

时间限制: 1 Sec 内存限制: 128 MB

[提交] [状态]

题目描述

新学期伊始,作为玉米高中的教务主任W某,又要安排学生们的课程表了。


W某想要知道所有可能的排课表方案,于是他开始在纸上列举所有方案,然而在写满了一摞A4纸后,他发现可能的方案太多了——用尽玉米高中所有的A4纸都写不完。


W某最终放弃了列举所有方案的想法,但他对排课表的方案数产生了兴趣。他的组合数学不太好,所以他找到了正在玉米高中就读的你,请你帮帮TA。


简单地说,玉米高中共有T个班级。


对于其中一个班级i,这个班级每天要上mi节互不相同的课,一共有ni节课可供选择,但这ni节课不能随便安排,其中也有一些限制:

·有ai节课不能安排在第一节上

·有bi节课不能安排在最后一节上

·没有任何一节课既不能在第一节上又不能在最后一节上

你需要求出每个班级排课表的方案数除以998244353的余数。

输入

第1行包含一个正整数T,表示玉米高中的班级数。


第2行到第T+1行,每行包含四个整数,第i+1行的四个整数ni,mi,ai,bi,分别表示班级i可选的课程数,一天的课程数,不能在第一节上的课程数,不能在最后一节上的课程数。

输出

输出T行,第i行表示班级i的排课表方案数除以998244353的余数。

样例输入 Copy

【样例1】

1

3 2 0 1

【样例2】

1

5 3 1 1

样例输出 Copy

【样例1】

4

【样例2】

39

提示

样例1解释

设3节可选的课为a,b,c,其中c不能排在最后一节

4种排课表的方案分别为:ab,ba,ca,cb


所有测试数据满足

·1≤T≤104

·2≤mi≤ni≤105

·ai+bi≤ni

思路:

高中排列组合? 划掉,忘得一干二净

今天刚学了容斥原理,但是依稀记得高中好像也是这么做的

直接计算方案数并不好算,因为要分类讨论很多种情况,我们不妨逆向思维。

用总的方案数 - a放在第一位的方案数 - b放在第一位的方案数 就可以

真的可以吗?

我们在计算a放在第一位的方案数时,其中有不合法的情况是b在最后一位,同理,我们在计算b放在第一位的方案数时,有不合法的情况是a在第一位,好像减多了呢。

只需要把这一部分加上就可以了~

To be short : 总的方案数-第一节是a的方案数-最后一节是b的方案数+第一节是a而且最后一节是b的方案数

其他的就是高中的知识了

化简过程就不放了~

image.jpeg

解释一下第二个式子,其他的都类似:

因为a的位置是确定的,我们首先从a里选出1个放在首位,此时需要选的还有m-1个,能够选的还有n-1个,所以就是从n-1里选m-1个,全排列。

另外,数据较大,需要预处理阶乘和逆元

模板:传送门

代码:

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N = 100010;
const int mod = 998244353;
ll fact[N];//阶乘
ll infact[N];//逆元
ll ksm(ll a,ll b,ll p){
    ll res=1;
    while(b){
//&运算当相应位上的数都是1时,该位取1,否则该为0。
        if(b&1)
            res=1ll*res*a%p;//转换为ll型
        a=1ll*a*a%p;
        b>>=1;//十进制下每除10整数位就退一位
    }
    return res;
}
void init(){
    fact[0]=1;
    infact[0]=1;
    for(int i=1;i<N;i++){
        fact[i]=fact[i-1]*i%mod;
        infact[i]=infact[i-1]*ksm(i,mod-2,mod)%mod;
    }
}
int main(){
    init();
    int t;
    cin>>t;
    int n,m,a,b;
    while(t--){
        cin>>n>>m>>a>>b;
        ll res=(fact[n]%mod*infact[n-m]%mod-fact[n-1]%mod*infact[n-m]%mod*(a+b)%mod+fact[n-2]%mod*infact[n-m]%mod*a*b%mod)%mod;
        if(res<0) res+=mod;
        printf("%lld\n",res);
    }
    return 0;
}

玉米田

时间限制: 1 Sec 内存限制: 128 MB

[提交] [状态]

题目描述

玉米中学的学生社会实践的内容是去玉米田中种玉米。


玉米中学有n块不同的玉米田,这些玉米田编号从1到n,且第i号玉米田与第i+1号玉米田相邻,特殊地,第n号玉米田与第1号玉米田相邻。


现在玉米中学购置了k种不同的玉米,为了美观,学校要求相邻的玉米田中不能种植同一种玉米,现在W某想要知道种植玉米的方案总数。


由于W某耐心有限,因此只需要你求出对20011021取模后的结果即可。

输入

一行两个整数n,k,表示玉米田的数量和玉米的种类数。

输出

加粗样式一行一个整数,表示种植玉米的方案数对20011021取模后的结果。

样例输入 Copy

【样例1】

4 2

【样例2】

4 3

样例输出 Copy

【样例1】

2

【样例2】

18

提示

样例1解释

设2种玉米为a,b

2种种植玉米的方案为:abab,baba


所有数据满足:n,k≤109


直接一波公式:

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N = 100010;
const ll mod = 20011021;
ll fact[N];//阶乘
ll infact[N];//逆元
ll ksm(ll a,ll b,ll p){
    ll res=1;
    while(b){
//&运算当相应位上的数都是1时,该位取1,否则该为0。
        if(b&1)
            res=1ll*res*a%p;//转换为ll型
        a=1ll*a*a%p;
        b>>=1;//十进制下每除10整数位就退一位
    }
    return res;
}
void init(){
    fact[0]=1;
    infact[0]=1;
    for(int i=1;i<N;i++){
        fact[i]=fact[i-1]*i%mod;
        infact[i]=infact[i-1]*ksm(i,mod-2,mod)%mod;
    }
}
int main(){
    ll n,k;
    cin>>n>>k;
    ll res=((k-1)*ksm(-1,n,mod)%mod+ksm(k-1,n,mod))%mod;
    cout<<res;
    return 0;
}

证明见大佬博客

类似的题:题目链接 题解链接

好像并不类似(无奈)

目录
相关文章
|
6月前
|
算法 安全 测试技术
[组合数学]LeetCode:2954:统计感冒序列的数目
[组合数学]LeetCode:2954:统计感冒序列的数目
|
6月前
|
算法 vr&ar 图形学
☆打卡算法☆LeetCode 166. 分数到小数 算法解析
☆打卡算法☆LeetCode 166. 分数到小数 算法解析
|
算法
【算法挨揍日记】day12——153. 寻找旋转排序数组中的最小值、LCR 173. 点名
已知一个长度为 n 的数组,预先按照升序排列,经由 1 到 n 次 旋转 后,得到输入数组。例如,原数组 nums = [0,1,2,4,5,6,7] 在变化后可能得到:
112 1
Visio如何插入高数公式
Visio如何插入高数公式
80 0
|
人工智能
P1208 [USACO1.3]混合牛奶 Mixing Milk(贪心思想加一点特别的循环)
P1208 [USACO1.3]混合牛奶 Mixing Milk(贪心思想加一点特别的循环)
69 0
786. 第 K 个最小的素数分数【我亦无他唯手熟尔】
786. 第 K 个最小的素数分数【我亦无他唯手熟尔】
48 0
Leecode 914.卡牌分组
Leecode 914.卡牌分组
80 0
UPC-趾压板矩阵(强行找规律)
UPC-趾压板矩阵(强行找规律)
102 0
UPC-趾压板矩阵(强行找规律)
UPC-喜爱(打表+二分)
UPC-喜爱(打表+二分)
90 0
UPC-喜爱(打表+二分)
|
机器学习/深度学习 算法 安全
【Python 百练成钢】时间调整、二进制数、回文素数、字母距离、CTF、Huffuman树、抽奖、前后缀最值、纯质数求解、花园灌溉
【Python 百练成钢】时间调整、二进制数、回文素数、字母距离、CTF、Huffuman树、抽奖、前后缀最值、纯质数求解、花园灌溉
【Python 百练成钢】时间调整、二进制数、回文素数、字母距离、CTF、Huffuman树、抽奖、前后缀最值、纯质数求解、花园灌溉