【数论】 再数17 容斥原理

学校做了一套题，大家虽然看不到数据，但这几道题都很有意思，我还是写在这里。

再数17

问题描述

之前，小cot发邮件问小cos一个区间所有整数里有多少个“17”，小cos在你的帮助下回答了这个奇怪的问题，然而得到的回信只有“你对了:-)”这几个字……小cos耿耿于怀，始终弄不明白小cot的意图。
终于，小cos决定主动问小cot一些问题。他先写了封邮件，问“17”的含义，回信写道：“我喜欢……17这个数字~”。小cos当然不能接受这样的回答，继续追问她为什么要自己数17，可回信又是简单的一句话：“我喜欢你……数17~”。
小cos急了，狠敲键盘，发了一连串问号过去。几分钟后，对方回道：“想要知道为什么，你再回答我一个问题——
“给你一个有n个正整数的数列{an}。一个正整数x若满足在数列{an}中存在一个正整数ai，使x≡17（mod ai），那么x就是一个‘cost数’。请问1到m的正整数中，有多少个‘cost数’？”
小cos看得一头雾水，只好又来请教你……

输入

输入文件名为c17.in，共2行。
第一行两个正整数n和m，意义见问题描述。
第二行n个正整数，分别为数列{an}中的n个数。

输出

输出文件名为c17.out。
输出一个整数，表示1到m中“cost数”的个数。

样例输入

3 100
18 22 23

样例输出

11

数据范围

对于20%的数据，有n≤20，m≤100,000；
另有40%的数据，n≤3，m＜2^31；
对于100%的数据，有n≤30，m＜2^31，17＜ai＜2^31。
（注：“≡”为同余符号，a≡b (mod k) 即a mod k = b mod k）

限制

2s

这道题数论，没什么好说的，容斥原理，没什么好说的，就是再一看，n≤30，再一看样例，ai不互质，一下就懵逼了，怎么做？
还是容斥公式，见下：

仔细看一看，复杂度是O（$2^nlogm$) 肯定T，要优化。
然后就把lcm（最小公倍数）大于m的return，深搜就好了。
这道题其实不难，写起来也很简单，就是怕dfs栈溢出，不敢写。
下面是代码，供参考

#include<iostream>
#include<cstring>
#include<cstdio>
#include<queue>
#include<cmath>
#include<algorithm>
#define M 32
#define ll long long
using namespace std;

ll n,m,ans,a[M];

ll gcd(ll m,ll n)
{
    return (n==0)?m:gcd(n,m%n);
}

ll lcm(ll m,ll n)
{
    return m/gcd(m,n)*n;
}

bool cmp(ll a,ll b)
{
    return a>b;
}

void dfs(ll i,ll s,ll b)
{
    if(s>m)return ;
    if(i==n+1)return ;
    ll l=lcm(s,a[i]);
    if(b&1)ans+=m/l;
    else ans-=m/l;
    dfs(i+1,s,b);
    dfs(i+1,l,b+1);
}

int main()
{
//  freopen("c17.in","r",stdin);
//  freopen("c17.out","w",stdout);
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++)scanf("%d",&a[i]);
    if(m>17)m-=17;
    else
    {
        cout<<0;
        return 0;
    }
    sort(a+1,a+n+1,cmp);
    dfs(1,1,1);
    cout<<ans+1;
}

大概就这样，如果有什么问题，或错误，请在评论区提出，谢谢。
最后吐槽一句，c17这个名字，让我想起了c85，c86，c87。。。