容斥原理简析

一、初步了解

韦恩图

简单的容斥原理的例子， 求出那图形的面积

n个圆相交的情况，面积：
单个圆的面积 - 两个圆（相交的面积） + 三个（相交的面积） - 四个（圆相交的面积）+五个（圆相交的面积）…+(-1) ^(n-1) * n（n个圆相交的面积）

容斥原理

1. 也可以考虑某些集合相交后的元素个数。
2. 时间复杂度：2 **n - 1
   在计数时，必须注意没有重复，没有遗漏 。
   为了使重叠部分不被重复计算，人们研究出一种新的计数方法，这种方法的基本思想是：先不考虑重叠的情况，把包含于某内容中的所有对象的数目先计算出来，然后再把计数时重复计算的数目排斥出去，使得计算的结果既无遗漏又无重复，这种计数的方法称为容斥原理。
   

二、能被整除的数

题目链接
题目描述：
给定一个整数n和m个不同的质数p1,p2,…,pm。
请你求出1~n中能被p1,p2,…,pm中的至少一个数整除的整数有多少个。
输入格式
第一行包含整数n和m。
第二行包含m个质数。
输出格式
输出一个整数，表示满足条件的整数的个数。
数据范围
1≤m≤16,
1≤n,pi≤109

分析

1. 1 - n 中 p 的 个数？ = n / p 下取整
   能同时被 k 个质数整除的数 = n / (p1 * p2 * .....pk)
2. 一共 2 **n - 1种
3. 选一个集合出来，符号为正， 选两个出来，符号为负的， 选奇数个集合，符号为正， 偶数个集合符号为负
4. 枚举 每种情况， int i 看成 n 位二进制数， 该位是1看成这集合被选，0 表示这集合不被选。
   如何确定每一位是不是1
   比如看 i 的第 k 位是不是1：
   i >> k & 1

代码

#include <iostream>
#include <algorithm>
using namespace std;
typedef long long LL;
const int N = 20;
int n, m;
int p[20];
int main(void){
    scanf("%d%d", &n, &m);
    for(int i = 0; i < m; i ++) scanf("%d", &p[i]);
    int res = 0;
    for(int i = 1; i < 1 << m; i++){//i左移m位，2 ** m次方,枚举每种方案
        int t = 1;//当前所有质数的乘积
        int cnt = 0; //表示当前i有几个1，即当前有几个集合
        for(int j = 0; j < m; j ++){
            if(i >> j & 1){ //如果这一位是1
                if((LL) t * p[j] > n){
                    t = -1;
                    break; //几个质数的乘积超过了n的范围，不符合题意，自然要舍去
                }
                t *= p[j];
                cnt ++;
            }
        }
        if(t != -1){
            if(cnt % 2){
                res += n/t;
            }
            else res -= n/t;
        }
    }
    cout << res << endl;
    return 0;
}