能被整除的数

能被整除的数

题目描述

---

样例解释

1到10的元素有{1，2，3，4，5，6，7，8，9，10}，能被2整除的元素集合为{2，4，6，8，10}为5个，能被3整除的元素集合为{3、6、9}为3个，由于6重复被多算了一次，因此总的次数就是5+3-1=7.

---

前置知识

容斥原理

容斥原理它总共有$2^n-1$项。即从n个集合中选出若干个集合，对于每个集合，有两种状态，选它或者不选它，因此有n个集合，于是就有$2^n$种状态，排除一个集合动不选的这个状态，还剩下$2^n-1$个状态。 $C_k^1-C_k^2+C_k^3-C_k^4\cdots +(-1{)}^{k+1}{C}_k^k$

由二项式定理：$C_n^0+C_n^1+C_n^2+\cdots +C_n^n=2^n$，于是有$C_n^1+C_n^2+\cdots +C_n^n=2^n-1$

如何理解$|\bigcup _{i=1}^n{S}_i|=\sum _{i=1}^n|S_i|-\sum _{1\le i<j\le n}|S_i\bigcap S_j|+\sum _{1\le i<j<k\le n}|S_i\bigcap S_j\bigcap S_k|$$\cdots +(-1{)}^{n+1}|S_1\bigcap S_2\cdots \bigcap S_n|$呢？

---

核心思路

设$S_i$表示1到n中能整除$p_i$的集合，那么根据容斥原理, 所有数的个数为各个集合的并集。我们设$|S|$表示集合$S$所含有的元素个数

以题目样例解释：

设集合$S_1$表示能被质数2整除的集合，则$S_1$={2，4，6，8，10}

设集合$S_2$表示能被质数3整除的集合，则$S_2$={3，6，9}

得到的结果就是$S_1$和$S_2$的并集，即$S=S_1\bigcup S_2$={2，3，4，6，8，9，10}，那么答案就是集合$S$所含有的元素个数，即为7.

每个集合实际上并不需要知道具体元素是什么，只要知道这个集合的大小即可。比如集合$S_1$表示能被质数2整除的集合，那么我们怎么知道1到n中能被$p_i$=2整除的元素的个数呢？其实有个公式，就是$n/p_i$，只要算出了$n/p_i$也就是算出了1到n中能被$p_i$整除的元素的个数了，也就是知道了这个集合$S_i$的大小了。即$|S_i|=n/p_i$，向下取整。比如题目中$|S_1|=10/2=5$，$|S_2|=10/3=3$

交集的大小如何确定呢？

比如说$p_i=2$，$p_i=3$，那么1到n中可能有同时能被2和3整除的，比如元素6。也就是如何确定$S_1\bigcap S_2$这个集合的大小呢？这里因为题目已经说了$p_i$是质数，那么这些质数的乘积就是他们的最小公倍数，n除这个最小公倍数就是交集的大小。举个例子吧：比如n=24，有质数$p_i$是2,3。那么1到24中能同时被2和3整除的元素有6，12，18，24共有4个，用公式计算就是：先算出2和3的最小公倍数，是6。那么交集的大小为24/6=4。在这个题目中就是$|S_1\bigcap S_2|=\frac{n}{p_1\ast p_2}=\frac{10}{2\ast 3}=1$

依次类推，$p_1,p_2,\cdots ,p_x$的交集大小就是$\frac{n}{p_1\ast p_2\ast \cdots \ast p_x}$

如何用代码表示每个集合的状态？

有m个质数，每个质数都确定了一种集合。这里可以考虑用二进制。以m = 4为例，每个二进制位可以表示选和不选该集合的状态，用1表示选了该集合，用0表示不选该集合，m=4因此需要4个二进制位。比如1011则表示选择了$S_1,S_2,S_4$这三个集合。其中$S_1$集合表示1到n中能整除$p_1$的所有元素的集合，$S_2$集合表示1到n中能整除$p_2$的所有元素的集合，$S_4$集合表示1到n中能整除$p_4$的所有元素的集合，故这个集合中元素的个数为$\frac{n}{p_1\ast p_2\ast p_4}$。 因为集合个数是3个，根据公式，前面的系数为$(-1{)}^{3-1}=1$。所以到当前这个状态时，应该是$res+=\frac{n}{p_1\ast p_2\ast p_4}$。这样用二进制0000到1111就可以表示m=4时的集合状态了。

由于至少要包含一个集合，因此枚举时是从0001到$2^m-1$，即区间是$[1,2^m-1]$。

如何理解$t\ast =p[j]$呢？

比如$p_1=2,p_2=3,p_3=5$，那么$t=p_1\ast p_2\ast p_3=2\ast 3\ast 5=30$。

如何理解$t\ast p[j]>n$呢？

比如$p_1=2,p_2=3,p_3=5$，n=24，那么$t\ast p[j]=2\ast 3\ast 5=30$，但是24小于30，1到24中并不存在能够被30整除的元素，因为30已经超过n这个数本身了，因此需要提前break。

由容斥原理可知，总共有$2^m-1$项，比如m=2，那么就有3项，则$|S_1\bigcup S_2|=|S_1|-|S_1\bigcap S_2|+|S_2|$，当i=1时，说明在枚举第一项，即$|S_1|$，当i=2时，说明在枚举第二项，即$|S_1\bigcap S_2|$，当i=3时，说明在枚举第三项，即$|S_3|$。

对于$|S_1|+|S_2|+\cdots +|S_n|$拆开来看的话，就是$|S_1|$,只选1个，奇数，所以前面符号是正+；对于$|S_2|$，只选1个，奇数，所以前面符号是正+；对于$|S_n|$，只选1个，奇数，所以前面符号是正+

也就是说我们也统计了$|S_1|、|S_2|、\cdots 、|S_n|$了。这正对应了容斥原理中等式右边的第一项。

---

代码

#include<iostream>
using namespace std;
typedef long long LL;
const int N=20;
int n,m;
int p[N];
int main()
{
    scanf("%d%d",&n,&m);
    for(int i=0;i<m;i++)
        scanf("%d",&p[i]);
    int res=0;  //表示满足条件的整数的个数。
    //由于至少要包含一个集合,因此是从00...001枚举到2^m-1
    //即枚举区间是[1,2^m-1]
    for(int i=1;i<(1<<m);i++)
    {
        int t=1;     //选中集合对应质数的乘积
        int cnt=0;  //选中的集合数量
        //枚举当前状态的每一位
        for(int j=0;j<m;j++)    //判断这m个二进制位
        {
            if((i>>j&1))//如果该二进制位是1
            {
                //乘积大于n, 则n/t = 0, 跳出这轮循环
                if((LL)t*p[j]>n)
                {
                    t=-1;
                    break;
                }
                cnt++;       //有一个1，集合数量+1
                t*=p[j];
            }
        }
        if(t==-1)
            continue;
        //选中奇数个集合, 由容斥原理公式可知,前面符号是正号+, n/t为当前这种状态的集合数量
        if(cnt&1)
            res+=n/t;
        //选中偶数个集合, 由容斥原理公式可知,前面符号是负号-, n/t为当前这种状态的集合数量
        else
            res-=n/t;
    }
    printf("%d\n",res);
    return 0;
}

---