【bzoj1853】 SCOI2010幸运数字容斥原理

最新推荐文章于 2019-02-25 19:05:00 发布

qingdaobaibai

最新推荐文章于 2019-02-25 19:05:00 发布

阅读量906

点赞数

分类专栏：组合数学

本文链接：https://blog.csdn.net/u012288458/article/details/47274555

版权

组合数学专栏收录该内容

9 篇文章 0 订阅

订阅专栏

跪烂晨神，一上来先想数位dp，结果根本没法做，然后看了看晨神的题解，发现幸运数字一共就2000来个，然后加一些剪枝，用容斥原理筛一下就好了，其实容斥就是选了奇数个就加，选了偶数个就减，要从大到小搜，这样容易爆，然后一定要判断如果爆了longlong，或者大于r的时候就退出，如果一个数是另一个数的倍数，那么就不用考虑这个数了。

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<cmath>
#include<iostream>
#include<algorithm>

using namespace std;

long long q[5010];
int tot;
long long l,r;
bool vis[5010];
long long ans;

void dfs(long long x)
{
	if (x>r) return;
	if (x) q[++tot]=x;
	dfs(x*10+6);dfs(x*10+8);
}

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

long long gcd(long long a,long long b)
{
	if (b==0) return a;
	else return gcd(b,a%b);
}

void dfs1(int x,int d,long long lcm)
{
	if (x==tot) return;
	x++;
	dfs1(x,d,lcm);
	d++;
	if (vis[x]) return;
	lcm=lcm*q[x]/gcd(lcm,q[x]);
	if (lcm<=0 ||lcm>r) return;
	if (d&1) ans+=r/lcm-(l-1)/lcm;
	else ans-=r/lcm-(l-1)/lcm;
	dfs1(x,d,lcm);
}

int main()
{
	scanf("%lld%lld",&l,&r);
	dfs(0);
	sort(q+1,q+tot+1,cmp);
	for (int i=1;i<=tot;i++) 
	  for (int j=i+1;j<=tot;j++)
	    if (q[i]%q[j]==0) vis[i]=1;
	dfs1(0,0,1);
	printf("%lld\n",ans);
	return 0;
}