容斥原理及其应用.docVIP

容斥原理及其应用 福建省上杭职业中专学校 赖锦德 先来看一道很简单的数学题：求不大于20能被2或3整除的正整数个数？ 先算出不大于20能被2整除的正整数有10个：2、4、6、8、10、12、14、16、18、20； 再算出不大于20能被3整除的正整数有6个：3、6、9、12、15、18。 但是不大于20能被2或3整除的正整数却并不是16个，而是13个。原因很简单：因为6、12、18这三个数都计算了两次，须扣除重复计算的次数。 上面这道看似简单的数学题其实就蕴涵了“容斥原理”的精髓。为了描述容斥原理，通常都要借助集合思想： 假定A和B是集合U的子集，当A∩B=时|A∪B|=|A| + |B|；当A∩B≠时|A| + |B|将|A∩B|计算了两次，所以|A∪B|=|A| + |B| - |A∩B|。这一公式便是容斥原理。 打开百度，用“什么是容斥原理”搜索一下，会得到如下通俗易懂的解答：在计数时，必须注意无一重复，无一遗漏。为了使重叠部分不被重复计算，人们研究出一种新的计数方法，这种方法的基本思想是：先不考虑重叠的情况，把包含于某内容中的所有对象的数目先计算出来，然后再把计数时重复计算的数目排斥出去，使得计算的结果既无遗漏又无重复，这种计数的方法称为容斥原理。 关于容斥原理，不得不提起大多数教材都涉及到的两个公式： [公式1] 设A=A1∪A2∪…Am，则 |A|=+…++…+ [公式2] 设AI，A=A1∪2A∪…Am，则 ||=| =|I|+…++…+ 公式1是容斥原理的一般性描述，公式2是一个推论称为淘汰原则。 容斥原理的两个公式看似繁杂，但只要对公式理解透彻，在应用中突破问题的关键往往还是很轻松的。下面通过“雕塑安置”例题讲述容斥原理的应用，题意描述如下： 为了迎接百年校庆，美化校园，请设计师设计了N座精美的雕塑，准备安置在校园里。整个校园可以抽象为一个N*N 的大网格，并且为了平均分布这些雕塑，学校决定网格的同一行、同一列必须有且只有一座雕塑，还规定不能出现1 座以上的雕塑出现在同一个1*1 网格里的情况。然而某些1*1 的网格恰巧是一片湖或者是食堂，这些网格就不能安置雕塑了。每个雕塑的造型是相同的，这样同一种安置方案中交换排列都算一种。 学校想知道有多少种安置方案，你能解决这个问题吗？ 输入说明： 第一行，两个整数 n,m 用空格分开 n表示n*n的大网格，m 表示不能安置雕塑的位置个数。 第二行至第m+1行，每行两个整数x,y，用空格分开，表示坐标（x,y）的1*1 网格上不能安置雕塑。 输出说明： 仅一行，方案的个数。 样例输入： 6 7 1 1 2 1 2 2 3 3 3 4 4 3 4 4 样例输出： 184 数据范围： 对于50%的数据， 保证n≤10，m≤10 对于100%的数据，保证n≤20，m≤10，方案总数≤263-1。 解题分析： 学生看到这个题，马上联想到经典的“n皇后”。仔细看看题目，与n皇后问题还是有一些区别的： 第1个区别是皇后连同一对角线都不能放，而雕塑没有这个规定； 第2个区别是安置雕塑时还有一些特殊的“禁区”。 所以，如果把题目中“某些1*1 的网格恰巧是一片湖或者是食堂，这些网格就不能安置雕塑了。”这些文字去掉，那么本题就完完全全是n皇后问题了。可是现在题目却加上了这些条件，除了雕塑之间不能“互相攻击”之外，有些网格是不能放置的。怎么办？ 很朴素也很简单的解决方法是：用数组记下这些不能放置的位置，搜索时避开它们， 参考程序如下： var a:array[1..20,1..20]of boolean; han:array[1..20]of boolean; ans:int64; n,i,j,m,aa,bb:integer; procedure search(dep:integer); var i:integer; begin if depn then begin inc(ans); exit; end; for i:=1 to n do if (not a[dep,i]) and (not han[i]) then begin han[i]:=true; search (dep+1); han[i]:=false; end; end; begin assign(input,arrange.in);assign(output,arrange.out); reset(input);rewrite(output); readln(n,m); fillchar(a,siz