【2018/11/01搜索大赛T1】cost 数

【题目】

传送门

【分析】

20 pts：

直接枚举。时间复杂度：O$(nm)$。

另外 40 pts：

过这 $40pts$，我们需要知道一点小知识：

• $1$ 到 $n$ 的正整数中，能被 $x$ 整除的数的个数 $nu{m}_x=\frac{n}{x}$（这里是整除，下同）；
• $1$ 到 $n$ 的正整数中，同时被 $x$、$y$ 整除的数个数 $nu{m}_{xandy}=\frac{n}{lcm(x,y)}$（$lcm$ 为最小公倍数）。
• 依此类推，同时被 $x_1,x_2,...,x_i$ 整除的数的个数为 $\frac{n}{lcm(x_1,x_2,...,x_i)}$

现在回到这道题来

将题目中的式子转化一下：$x\equiv 17(mod{a}_i)\to xmod{a}_i=17mod{a}_i\to (x-17)mod{a}_i=0$（注意已保证 $a_i>17$，所以 $17mod{a}_i=17$），于是问题变为求 $0$ 到 $m-17$ 整数中至少能被 $a$ 中一个数整除的有多少个。

根据容斥原理，$|A\cup B|=|A|+|B|-|A\cap B|$（$|A|$ 表示集合 $A$ 元素个数），可得 $1$ 到 $n$ 的正整数中，能被 $x$ 或 $y$ 整除的数个数 $nu{m}_{xory}=nu{m}_x+nu{m}_y-nu{m}_{xandy}=\frac{n}{x}+\frac{n}{y}-\frac{n}{lcm(x,y)}$。对于三个数 $x$、$y$、$z$ 的情况，只要类比一下，画画图就可以知道：

$nu{m}_{xoryorz}=nu{m}_x+nu{m}_y+nu{m}_z-nu{m}_{xandy}-nu{m}_{xandz}-nu{m}_{yandz}+nu{m}_{xandyandz}=\frac{n}{x}+\frac{n}{y}+\frac{n}{z}-\frac{n}{lcm(x,y)}-\frac{n}{lcm(x,z)}-\frac{n}{lcm(y,z)}+\frac{n}{lcm(x,y,z)}$

时间复杂度：O$(logm)$

100 pts：

容斥原理是可以扩展到 $n$ 个集合的，具体就是口诀：奇加偶减

这时我们需要用个 $dfs$ 枚举每个 $a_i$ 选与不选，若选出 $k$ 个数，且 $k$ 为奇数，则 $ans$ 加上 $\frac{m}{lcm}$（这里的 $lcm$ 就是选出的数的最小公倍数，下同），若为偶数 $ans$ 减去 $\frac{m}{lcm}$。

当然，这个算法复杂度还比较大，还要加上优化，包括可行性剪枝（选出的数的 $lcm$ 要小于等于 $m$。由于 $a_i>17$，该剪枝十分有效）、优化搜索顺序（先搜大数）。

时间复杂度：O$(2^n\ast logm)$（实际小得多）

【代码】

#include<cstdio>
#include<cstring>
#include<algorithm>
#define N 105
#define ll long long
using namespace std;
int n,m,ans,a[N];
bool comp(int a,int b)  {return a>b;}
ll Lcm(ll a,ll b)  {return a*b/__gcd(a,b);}
void dfs(int x,int num,ll lcm)
{
	if(lcm>m)  return;
	if(x>n&&lcm==1)  {ans++;return;}
	if(x>n)  {num&1?ans+=m/lcm:ans-=m/lcm;return;}
	dfs(x+1,num,lcm),dfs(x+1,num+1,Lcm(lcm,a[x]));
}
int main()
{
//	freopen("count.in","r",stdin);
//	freopen("count.out","w",stdout);
	int i;
	scanf("%d%d",&n,&m);
	for(i=1;i<=n;++i)
	  scanf("%d",&a[i]);
	sort(a+1,a+n+1,comp);
	m-=17,dfs(1,0,1);
	printf("%d",ans);
//	fclose(stdin);
//	fclose(stdout);
	return 0;
}