UESTC 1544 当咸鱼也要按照基本法 组合数学

题目链接
参考博客
此题还出自

玲珑学院 1138 震惊，99%+的中国人都会算错的问题 【容斥】【技巧】

DESCRIPTION
众所周知zhu是一个大厨，zhu一直有自己独特的咸鱼制作技巧.
tang是一个咸鱼供应商，他告诉zhu在他那里面有N条咸鱼（标号从1到N）可以被用来制作.每条咸鱼都有一个咸鱼值Ki，初始时所有Ki都是0
zhu是一个特别的人，他有M个咸数（咸鱼数字）, 对于每个咸数x，他都会让所有满足标号是x倍数的咸鱼的咸鱼值异或上1.zhu现在想知道经过了这M个咸数的筛选之后，最终有多少条的咸鱼的咸鱼值是1?
INPUT
输入的第一行包含一个整数T(1≤T≤1000)，表示有T组数据.
对于每组数据：
输入第一行只有两个整数N(1≤N≤10^9),M(1≤M≤15).
接下来一行有
M
M个整数，依次对应zhu的每个咸数(1≤咸数≤2∗105).
OUTPUT
对于每组数据，输出答案.
SAMPLE INPUT
2
10 1
3
10 1
1
SAMPLE OUTPUT
3
10
题解：考虑容斥，答案为一个数的倍数-两个数lcm的倍数+三个数lcm的倍数
具体计算每一类个数的时候，有如下计算公式

系数(该类的个数)(该类向下递推时每一项包含下一类的个数)
以2 3 5 6 为例
第一类 一个数的{2,3,5,6}
第二类 两个数的{(2,3),(2,5),(2,6)(3,5)(3,6)(5,6)}
第三类 三个数的{(2,3,5)(2,3,6)(2,5,6)(3,5,6)}
则由第一类和第二类向第三类递推时的计算式子为
1*4*3-2*6*2+(4*4*1) = 4
那么第三类的系数4即括号内的第一个4是如何算出来的呢?
很明显最后第三类有4个，第一类每一个（例如2包含了 2 3 5,2 3 6,2 5 6）都包含第三类的三个，乘上第一类个数4*第一类系数1，第二类同理，因为这样只有一个未知数所以就可以算出第三类的系数是4，这样找规律可以看出第k类的系数正好是(1<<(k-1))，
那么为什么是1<<(k-1)呢？

因为从当n一个一个增加时，其内部的计算是2倍2倍的增长的，具体看上面的第一个从一推到二的时候因为跟平时的容斥不一样此题的规则是偶消奇不消,所以求的是阴影部分的面积即需要减去2AB (记为原理1),以后维度增加时原理也是一样的,比如从二推到三,此时需要先看图三,然后从图三推到图二,这样当n增加时,其实是原理1的叠加,因此第k类的系数正好是(1<<(k-1))
下面代码中的n/cnt 计算的则是该类向下递推时每一项包含下一类的个数 ,cnt为当前 前面所有数的最小公倍数

#include<cstdio>
#include<algorithm>
#include<iostream>
#include<cstdlib>
#include<string.h>
#include<vector>
using namespace std;

int m;
long long n;
const int maxn  = 200006;
int vis[maxn];
long long  a[299990];
long long ans = 0;
long long  lcm(long long a,long long b)
{
    return a*b/(__gcd(a,b));
}
void dfs(int o,int x,long long cnt)
{
    cnt = lcm(cnt,a[o]);
    if(x&1)ans+=n/cnt*(1LL<<(x-1));
    else ans-=n/cnt*(1LL<<(x-1));
    for(int i=o+1;i<m;i++)
    {
        dfs(i,x+1,cnt);
    }
}
int main()
{
    int t;
    cin>>t;

    while(t--)
    {
        cin>>n>>m;
        for(int i=0;i<m;i++)
            scanf("%lld",&a[i]);
        ans =0 ;
        for(int i=0;i<m;i++)
        dfs(i,1LL,1LL);
        cout<<ans<<endl;
    }
    return 0;
}