快来踩爆这个蒟蒻吧|

Little_corn

园龄:1年1个月粉丝:10关注:17

2024-04-25 12:57阅读: 41评论: 0推荐: 0

容斥原理

容斥原理:

容斥原理是一种在知道所有集合之间的交,求集合之间的并的数学方法。(注:交即为两个集合之间相同的部分,记作 |A||B|

problem:

U 中元素有 n 种不同的属性,而第 i 种属性称为 Pi,拥有属性 Pi 的元素构成集合 Si,现在请求出 U 中有哪些元素。

易知:

也就是:

(摘自 OI-Wiki)

不定方程非负整数解计数:

problem:

给定 n 个变量,第 i 个记作 xi。对于 xi 有限制 xiai 。且对于所有 x ,有 x1+x2+...+xn=m ,问有多少组解。

solution:(by OI-Wiki)

对于这个问题,我们可以抽象出一个容斥模型:

1.全集:i=1nxi=m 解的个数。

2.集合:Si=xi

3.属性:即为限制条件。

解即为 |i=1nSi|

我们可以通过求这个解集的补集,并用全集减去补集即可。

易知解集的补集为:|i=1n!Si|(不会补集怎么打) 使用插板法+容斥原理即可。

例题

贴个代码:

#include<bits/stdc++.h>
#define int long long
using namespace std;
const int mod=1e9+7;
int n,m,a[50],ans;
int inv(int x);
int c_zh(int x,int y);
signed main(){
cin>>n>>m;
for(int i=1;i<=n;i++)cin>>a[i];
for(int s=1;s<=(1<<n)-1;s++){
int cnt=0,sum=0;
for(int i=0;i<n;i++){
if(s&(1<<i)){
sum+=a[n-i]+1;
cnt++;
}
}
int f=-1;
if(cnt&1)f=1;
ans=(ans+f*c_zh(n-1,m-sum-1+n)%mod+mod)%mod;//这里是因为将限制减掉后变成了 0 ,所以还要加上n。
}
cout<<(c_zh(n-1,m-1+n)%mod-ans%mod+mod)%mod<<endl;
return 0;
}
int fast_pow(int x,int y){
int ret=x%mod,res=1;
while(y){
if(y&1)res*=ret;
ret*=ret;
res%=mod;ret%=mod;
y>>=1;
}
return res;
}
int inv(int x){
return fast_pow(x,mod-2);
}
int c_zh(int x,int y){
if(y<0)return 0;
int res=1;
for(int i=1;i<=x;i++){
res=(res*((y-i+1)%mod))%mod;// 注意这个可能溢出
}
for(int i=1;i<=x;i++){
res=(res*inv(i))%mod;
}
//cout<<"x: "<<x<<" y: "<<y<<" ans: "<<res<<endl;
return res;
}

错位排列:

problem:对于一个有 n 个元素的排列 P ,如果对于所有 in , 都有 pii,则说排列 P 是一个错位排列。问对于有 n 个元素的序列有多少个错排列。

solution:

还是考虑使用全集 - 补集的方法去求解。

模型:

1.全集:n 个元素的排列:n!

2.集合:Pi

3.属性:Pii

那么很容易发现 Pi 的补集 Ui 即为满足 Pi=i 的一种解集。

所以答案即为:

n!|i=1nUi|

通过容斥原理以及组合数的定义,我们可以知道答案为(这里不再推式子):

n!i=0n(1)ii!

本文作者:Little_corn

本文链接:https://www.cnblogs.com/little-corn/p/18157425

版权声明:本作品采用知识共享署名-非商业性使用-禁止演绎 2.5 中国大陆许可协议进行许可。

posted @   Little_corn  阅读(41)  评论(0编辑  收藏  举报
点击右上角即可分享
微信分享提示
评论
收藏
关注
推荐
深色
回顶
收起