容斥原理(道理很简单,重在实践)

2024-02-17 22:38

本文主要是介绍容斥原理(道理很简单,重在实践),希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

前言

在组合数学中如何求2或2个以上集合的并集中元素的个数,若这几个集合两两无关,也就是没有交集,可以用加法原理,这个比较好理解,把所有集合的大小加起来就行了。若出现交集,就需要使用容斥原理。

理解

我认为容斥原理就是一个把多算的那块去掉,多去掉的那块补上的过程,可能这么说比较抽象,但确实如此。

这里举个例子:

计数1到600之间不能被6整除的整数个数。

这里600个是包含能被6整除的个数,所以要把多算的那块去掉,从1开始每6个数都有一个可以被6整除,所以600/6 = 100答案就是600-100=500

这个例子是减法原理,但也是容斥原理最简单的例子。

这个例子抽象一下是对于有限集合S,有一个性质p不满足某个性质,那么求满足性质的集合大小。

将这个问题拓展一下,p1,p2,p3…pm不满足某个性质,这些不满足性质的集合对应的是 A i A_i Ai,我们要求满足的集合大小。所求的对象个数可以由下面的交错表达式给出:
∣ A 1 ‾ ⋂ A 2 ‾ ⋂ ⋯ ⋂ A m ‾ ∣ = ∣ S ∣ − ∑ ∣ A i ∣ + ∑ ∣ A i ⋂ A j ∣ − ∑ ∣ A i ⋂ A j ⋂ A k ∣ + ⋯ + ( − 1 ) m ∣ A 1 ⋂ A 2 ⋂ ⋯ ⋂ A m ∣ \mid \overline{A_1}\bigcap\overline{A_2}\bigcap\cdots\bigcap\overline{A_m} \mid =\mid S \mid - \sum{\mid A_i \mid} + \sum{\mid A_i\bigcap A_j \mid} - \sum{\mid A_i\bigcap A_j \bigcap A_k\mid} + \cdots+(-1)^{m}\mid {A_1}\bigcap{A_2}\bigcap\cdots\bigcap{A_m}\mid A1A2Am=SAi+AiAjAiAjAk++(1)mA1A2Am
根据上述的表达式我们可以解决

计数1到1000之间不能被5,6,8整除的整数个数

很明显|S| = 1000

|A1| = 1000/5 = 200

|A2| = 1000/6 = 166

|A3| = 1000/8=125

然后我们要求的交集部分

(A1,A2):1000/lcm(5, 6)=33

(A1,A3):1000/lcm(5, 8)=25

(A1,A2):1000/lcm(6, 8)=33

(A1,A2,A3):1000/lcm(5, 6, 8) = 8

答案:1000-(200+166+125)+(33+25+41)-8=600

题目

一班容斥原理的题目不能直接套那个公式,需要我们真正明白把多算的那块去掉,多去掉的那块补上的过程的含义

The Lottery

UVA - 10325

这道题和上面第二例子一样,只不过不能被N个数整除,只要枚举被1~N个数整除,把所有情况都找出来,套公式即可,这里值得注意的是lcm和gcd都具有结合律,也就是说lcm(a,b,c) = lcm(lcm(a,b),c)

代码

#include <bits/stdc++.h>
using namespace std;typedef long long ll;
const int maxm = 20;
ll num[maxm], n;
ll gcd(ll a, ll b) {return b == 0 ? a : gcd(b, a % b);
}ll lcm(ll a, ll b) {return a / gcd(a, b) * b;
}
int m;
ll ans;
void dfs(int c, int cur, int i, ll now) {if(cur == c) {ans += ((c & 1) ? -1 : 1) * (n/now);return;}while(i < m) {dfs(c, cur+1, i+1, lcm(now, num[i]));i++;}
}
int main()
{// freopen("/Users/maoxiangsun/MyRepertory/acm/i.txt", "r", stdin);while(~scanf("%lld%d", &n, &m)) {for(int i = 0; i < m; i++) {scanf("%lld", num+i);}ans  = n;for(int i = 1; i <= m; i++) {dfs(i, 0, 0, 1);}printf("%lld\n", ans);}return 0;
}

Cheerleaders

UVA - 11806

这道题的意思是一个n*m的方格表,每个方格最多放一个东西,我们要放k个东西,问在上下左右那几条边都有东西的前提下有多少种方式。

也是套那个公式,最后可以推出这样的结果:
a n s = C ( n m , k ) − 2 ( C ( ( n − 1 ) m , k ) + C ( ( m − 1 ) n , k ) ) + ( C ( ( n − 2 ) m , k ) + C ( ( m − 2 ) n , k ) + 4 C ( ( m − 1 ) ( n − 1 ) , k ) − 2 C ( m − 2 ) ( n − 1 ) , k ) − 2 C ( C ( m − 1 ) ( n − 2 ) , k ) + C ( ( m − 2 ) ( n − 2 ) , k ) ; ans = C(nm, k) -2(C((n-1)m, k)+C((m-1)n,k))\\ +(C((n-2)m,k)+C((m-2)n,k)+4C((m-1)(n-1), k)\\ -2C(m-2)(n-1),k)-2C(C(m-1)(n-2),k) +C((m-2)(n-2), k); ans=C(nm,k)2(C((n1)m,k)+C((m1)n,k))+(C((n2)m,k)+C((m2)n,k)+4C((m1)(n1),k)2C(m2)(n1),k)2C(C(m1)(n2),k)+C((m2)(n2),k);
还有就是组合数取模问题,利用杨辉三角那个样子打表处理这种n*m不是很大,比较合适,再不行就用Lucas定理……

#include <bits/stdc++.h>
using namespace std;typedef long long ll;
const ll maxn = 420;
const ll mod = 1000007;
ll CC[maxn][maxn] = {0};
void init() {CC[0][0] = 1;for(ll i = 1; i < maxn; i++) {CC[i][0] = 1;for(ll j = 1; j <= i; j++) {CC[i][j] = (CC[i-1][j-1] + CC[i-1][j])%mod;}}
}
ll C(ll n, ll m) {return CC[n][m];
}
int main()
{ios::sync_with_stdio(false); cin.tie(0);ll n,m,k;ll T; init();cin >> T;int cas=0;while(T--) {cin >> n >> m >> k;ll ans =  (C(n*m,k)-(2*C((n-1)*m, k) + 2*(C(n*(m-1), k)) -C((m-2)*n, k) -C((n-2)*m, k) - 4*C((m-1)*(n-1) , k) + 2*C((m-2)*(n-1), k) + 2*C((m-1)*(n-2),k)-C((m-2)*(n-2), k) ) );while(ans < 0) {ans += mod;}cout << "Case "<<++cas<<": ";cout << ans % mod << endl; }return 0;
}

2018山东省赛F题

https://ac.nowcoder.com/acm/contest/123/F

这道题的意思是给出四个区间左右端点li,ri,从中取出x1,x2,x3,x4,问这四个值相邻元素不等的取法有多少种,可以用容斥想,S-不符合条件的总数。

情况1,只考虑两两相邻相等,可以分出四组来

情况2,在三个相邻相等和x1=x2&&x3=x4与x2=x3&&x1=x4公6组,不难发现在情况1中包含了情况2的可能所以要把多去的情况1加上去,加多少呢?加一倍,为啥?比如情况1在计算x1=x2的数量时已经计算了x3=x4的数量,而x3=x4的数量时又计算了x1=x2所以计算了两遍x1=x2&&x3=x4

情况3,x1=x2=x3=x4

这是最容易错的地方,我已开始做也是因为这wa了,为啥,情况1多算4次,情况2把6次减去了,情况3需要1次所以需要加上3倍的情况3才行!

代码:

#include <bits/stdc++.h>
using namespace std;typedef long long ll;
ll ans;
ll l[5],r[5];
const ll mod = 1e9 + 7;int main()
{int T; scanf("%d", &T);while(T--) {ans = 1;for(int i = 0; i < 4; i++) {// cin >> l[i] >> r[i];scanf("%lld%lld", &l[i], &r[i]);ans = ans * (r[i] - l[i] + 1) % mod;}ll left, right;for(int i = 0; i < 4; i++) {left = max(l[i], l[(i+1)%4]);right = min(r[i], r[(i+1)%4]);if(right >= left) ans=((ans-(right-left+1)*(r[(i+2)%4]-l[(i+2)%4]+1)%mod*(r[(i+3)%4]-l[(i+3)%4]+1)%mod)%mod+mod)%mod;} for(int i = 0; i < 4; i++) {left = max(l[i], max(l[(i+1)%4], l[(i+2)%4]));right = min(r[i], min(r[(i+1)%4], r[(i+2)%4]));if(right >= left) ans = (ans + (right - left + 1) * (r[(i+3)%4] - l[(i+3)%4] + 1)%mod)%mod;}ll l1,l2,r1,r2;for(int i = 0; i < 2; i++) {l1 = max(l[i], l[i+1]);r1 = min(r[i], r[i+1]);l2 = max(l[(i+2)%4], l[(i+3)%4]);r2 = min(r[(i+2)%4], r[(i+3)%4]);if(l1 <= r1 && l2 <= r2) {ans = (ans + (r1-l1+1)*(r2-l2+1)%mod)%mod;}}left = l[0];right = r[0];for(int i = 1; i < 4; i++) {left = max(l[i], left);right = min(r[i], right);}if(right >= left) {ans = (ans - (3*(right - left + 1))%mod + mod ) % mod;}printf("%lld\n", ans);}   return 0;
}

总结

容斥原理这类题难得地方在于理解模型中集合与集合的公共部分,不要盲目的套公式!

这篇关于容斥原理(道理很简单,重在实践)的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/719196

相关文章

Java Spring 中 @PostConstruct 注解使用原理及常见场景

《JavaSpring中@PostConstruct注解使用原理及常见场景》在JavaSpring中,@PostConstruct注解是一个非常实用的功能,它允许开发者在Spring容器完全初... 目录一、@PostConstruct 注解概述二、@PostConstruct 注解的基本使用2.1 基本代

Golang HashMap实现原理解析

《GolangHashMap实现原理解析》HashMap是一种基于哈希表实现的键值对存储结构,它通过哈希函数将键映射到数组的索引位置,支持高效的插入、查找和删除操作,:本文主要介绍GolangH... 目录HashMap是一种基于哈希表实现的键值对存储结构,它通过哈希函数将键映射到数组的索引位置,支持

Spring Boot 整合 SSE的高级实践(Server-Sent Events)

《SpringBoot整合SSE的高级实践(Server-SentEvents)》SSE(Server-SentEvents)是一种基于HTTP协议的单向通信机制,允许服务器向浏览器持续发送实... 目录1、简述2、Spring Boot 中的SSE实现2.1 添加依赖2.2 实现后端接口2.3 配置超时时

Python使用getopt处理命令行参数示例解析(最佳实践)

《Python使用getopt处理命令行参数示例解析(最佳实践)》getopt模块是Python标准库中一个简单但强大的命令行参数处理工具,它特别适合那些需要快速实现基本命令行参数解析的场景,或者需要... 目录为什么需要处理命令行参数?getopt模块基础实际应用示例与其他参数处理方式的比较常见问http

Java Optional的使用技巧与最佳实践

《JavaOptional的使用技巧与最佳实践》在Java中,Optional是用于优雅处理null的容器类,其核心目标是显式提醒开发者处理空值场景,避免NullPointerExce... 目录一、Optional 的核心用途二、使用技巧与最佳实践三、常见误区与反模式四、替代方案与扩展五、总结在 Java

Spring Boot循环依赖原理、解决方案与最佳实践(全解析)

《SpringBoot循环依赖原理、解决方案与最佳实践(全解析)》循环依赖指两个或多个Bean相互直接或间接引用,形成闭环依赖关系,:本文主要介绍SpringBoot循环依赖原理、解决方案与最... 目录一、循环依赖的本质与危害1.1 什么是循环依赖?1.2 核心危害二、Spring的三级缓存机制2.1 三

C#中async await异步关键字用法和异步的底层原理全解析

《C#中asyncawait异步关键字用法和异步的底层原理全解析》:本文主要介绍C#中asyncawait异步关键字用法和异步的底层原理全解析,本文给大家介绍的非常详细,对大家的学习或工作具有一... 目录C#异步编程一、异步编程基础二、异步方法的工作原理三、代码示例四、编译后的底层实现五、总结C#异步编程

Python 中的 with open文件操作的最佳实践

《Python中的withopen文件操作的最佳实践》在Python中,withopen()提供了一个简洁而安全的方式来处理文件操作,它不仅能确保文件在操作完成后自动关闭,还能处理文件操作中的异... 目录什么是 with open()?为什么使用 with open()?使用 with open() 进行

Go 语言中的select语句详解及工作原理

《Go语言中的select语句详解及工作原理》在Go语言中,select语句是用于处理多个通道(channel)操作的一种控制结构,它类似于switch语句,本文给大家介绍Go语言中的select语... 目录Go 语言中的 select 是做什么的基本功能语法工作原理示例示例 1:监听多个通道示例 2:带

鸿蒙中@State的原理使用详解(HarmonyOS 5)

《鸿蒙中@State的原理使用详解(HarmonyOS5)》@State是HarmonyOSArkTS框架中用于管理组件状态的核心装饰器,其核心作用是实现数据驱动UI的响应式编程模式,本文给大家介绍... 目录一、@State在鸿蒙中是做什么的?二、@Spythontate的基本原理1. 依赖关系的收集2.