码迷,mamicode.com
首页 > 其他好文 > 详细

bzoj1853【SCOI2010】幸运数字

时间:2016-05-30 15:41:34      阅读:180      评论:0      收藏:0      [点我收藏+]

标签:

1853: [Scoi2010]幸运数字

Time Limit: 2 Sec  Memory Limit: 64 MB
Submit: 1837  Solved: 673
[Submit][Status][Discuss]

Description

在中国,很多人都把6和8视为是幸运数字!lxhgww也这样认为,于是他定义自己的“幸运号码”是十进制表示中只包含数字6和8的那些号码,比如68,666,888都是“幸运号码”!但是这种“幸运号码”总是太少了,比如在[1,100]的区间内就只有6个(6,8,66,68,86,88),于是他又定义了一种“近似幸运号码”。lxhgww规定,凡是“幸运号码”的倍数都是“近似幸运号码”,当然,任何的“幸运号码”也都是“近似幸运号码”,比如12,16,666都是“近似幸运号码”。 现在lxhgww想知道在一段闭区间[a, b]内,“近似幸运号码”的个数。

Input

输入数据是一行,包括2个数字a和b

Output

输出数据是一行,包括1个数字,表示在闭区间[a, b]内“近似幸运号码”的个数

Sample Input

【样例输入1】
1 10
【样例输入2】
1234 4321

Sample Output

【样例输出1】
2
【样例输出2】
809

HINT

【数据范围】
对于30%的数据,保证1 < =a < =b < =1000000
对于100%的数据,保证1 < =a < =b < =10000000000

Source




容斥原理+搜索剪枝

方法同bzoj2393,但是这道题数据范围略大,所以DFS得时候要从后往前DFS,这样剪枝优化更明显。




#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cmath>
#include<cstring>
#include<algorithm>
#define F(i,j,n) for(int i=j;i<=n;i++)
#define D(i,j,n) for(int i=j;i>=n;i--)
#define ll long long
#define maxn 3005
using namespace std;
int cnt,tot;
ll l,r,ans,a[maxn],b[maxn];
bool tag[maxn];
void find(ll x)
{
	if (x>r) return;
	a[++cnt]=x;
	find(x*10+6);find(x*10+8);
}
ll gcd(ll a,ll b)
{
	return !b?a:gcd(b,a%b);
}
void dfs(int x,int flag,ll t)
{
	if (!x)
	{
		if (t!=1) ans+=(r/t-(l-1)/t)*flag;
		return;
	}
	dfs(x-1,flag,t);
	ll tmp=gcd(t,a[x]);
	if ((double)t*a[x]/tmp>(double)r) return;
	dfs(x-1,-flag,t/tmp*a[x]);
}
int main()
{
	scanf("%lld%lld",&l,&r);
	find(6);find(8);
	sort(a+1,a+cnt+1);
	F(i,1,cnt) F(j,1,i-1) if (a[i]%a[j]==0){tag[i]=1;break;}
	for(int i=1;a[i];i++) if (!tag[i]) b[++tot]=a[i];
	cnt=tot;
	F(i,1,cnt) a[i]=b[i];
	dfs(cnt,-1,1);//要从大往小搜索,不然会TLE 
	printf("%lld\n",ans);
}


bzoj1853【SCOI2010】幸运数字

标签:

原文地址:http://blog.csdn.net/aarongzk/article/details/51520376

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!