【BZOJ1853/Scoi2010】幸运数字

题目：幸运数字

解析：

  简单容斥原理。
  需要注意中间过程可能爆$longlong$，所以可以先强转为$double$，详见代码。
  还有一个小优化：搜索顺序从大往小搜会快很多。

代码：

#include <bits/stdc++.h>
#define int long long
using namespace std;

const int Max=12025;
int l,r,s1,s2,ans;
int a[Max],b[Max],vis[Max];

inline void search(int now)
{
	if(now>r) return;
	if(now) a[++s1]=now;
	search(now*10+6),search(now*10+8);
}

inline void dfs(int now,int sum,int lcm) //now当前数sum奇偶
{
	if(now>s2)
	{
	  if(sum&1) ans+=r/lcm-(l-1)/lcm;
	  else if(sum) ans-=r/lcm-(l-1)/lcm;
	  return;
	}
	dfs(now+1,sum,lcm);
	int s=lcm/__gcd(a[now],lcm);
	if((double)a[now]*s>r) return;   //题解中的强转
	dfs(now+1,sum+1,s*a[now]);
}

signed main()
{
	cin>>l>>r;
	search(0);
	sort(a+1,a+s1+1);
	for(int i=1;i<=s1;i++)
	{
	  if(vis[i]) continue;
	  b[++s2]=a[i];
	  for(int j=i+1;j<=s1;j++) if(!(a[j]%a[i])) vis[j]=1;
	}
	for(int i=1;i<=s2;i++) a[s2-i+1]=b[i];
	dfs(1,0,1);
	cout<<ans;
	return 0;
}