1853: [Scoi2010]幸运数字

Time Limit: 2 Sec Memory Limit: 64 MB
Submit: 2117 Solved: 779
[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

Day1

[Submit][Status][Discuss]

分析

不难想到用容斥原理来统计。

先预处理出所有小于1e10的幸运数字，并不是很多。但是发现枚举所有的组合还是会爆炸的，需要一些剪枝。

1. 对于两个幸运数字，x<y，如果有y为x的倍数，则y可以忽略，因为x可以完全覆盖y的倍数。

2. 对于一种组合，如果目前的积已经大于N，即再进行下去得到的都是加减0的无意义操作，可以直接跳出。

3. 可以把GCD函数写成非递归的形式，但貌似没多大用，跑出来的结果差距不是很大，也许是我写得不好。

4. 枚举的时候从大往小枚举，据说有奇效，因为懒癌晚期，我并没有对比验证。

代码

 #include <bits/stdc++.h>

 using namespace std;

 typedef long long LL;

 const int N = ;

 const LL lim = 10000000000LL;

 template <class T>

 __inline void read(T &x)

 {

     x = ; char c = getchar();

     while (c < '')

         c = getchar();

     while (c >= '')

     {

         x = x* + c - '';

         c = getchar();

     }

 }

 LL gcd(LL a, LL b)

 {

     if (a < b)

     {

         a ^= b;

         b ^= a;

         a ^= b;

     }

     while (b)

     {

         a %= b;

         a ^= b;

         b ^= a;

         a ^= b;

     }

     return a;

 }

 LL num[N]; int tot = ;

 __inline void prework(void)

 {

     int t, tail = ;

     for (t = ; num[t] <= lim; ++t)

     {

         num[++tail] = num[t] *  + ;

         num[++tail] = num[t] *  + ;

     }

     for (int i = ; i <= t; ++i)

     {

         bool flag = true;

         for (int j = ; j <= tot; ++j)

             if (num[i] % num[j] == )

                 { flag = false; break; }

         if (flag)num[++tot] = num[i];

     }

 }

 LL answer, limit;

 void search(int t, bool f, LL sum)

 {

     if (t)

     {

         search(t - , f, sum);

         LL GCD = gcd(num[t], sum);

         if (sum / GCD <= limit / num[t])

         {

             LL LCM = sum / GCD * num[t];

             if (f)

                 answer -= limit / LCM;

             else

                 answer += limit / LCM;

             search(t - , !f, LCM);

         }

     }

 }

 __inline LL count(LL n)

 {

     limit = n;

     answer = ;

     int pos = ; 

     while (pos <= tot

     && num[pos] <= n)++pos;

     search(pos - , , );

     return answer;

 }

 signed main(void)

 {

     prework();

     LL a; read(a);

     LL b; read(b);

     printf("%lld\n",

         count(b)

     -    count(a - )

     );

 }

BZOJ_1853.cpp

@Author: YouSiki

BZOJ 1853: [Scoi2010]幸运数字的更多相关文章

Bzoj 1853: [Scoi2010]幸运数字容斥原理,深搜
1853: [Scoi2010]幸运数字 Time Limit: 2 Sec Memory Limit: 64 MBSubmit: 1774 Solved: 644[Submit][Status] ...
bzoj 1853: [Scoi2010]幸运数字容斥
1853: [Scoi2010]幸运数字 Time Limit: 2 Sec Memory Limit: 64 MBSubmit: 1170 Solved: 406[Submit][Status] ...
BZOJ 1853: [Scoi2010]幸运数字（容斥原理）
http://www.lydsy.com/JudgeOnline/problem.php?id=1853 题意: 在中国,很多人都把6和8视为是幸运数字!lxhgww也这样认为,于是他定义自己的“幸运 ...
bzoj 1853: [Scoi2010]幸运数字&&2393: Cirno的完美算数教室【容斥原理】
翻了一些blog,只有我用状压预处理嘛2333,.把二进制位的0当成6,1当成8就行啦.(2393是2和9 然后\( dfs \)容斥,加上一个数的\( lcm \),减去两个数的\( lcm \), ...
●BZOJ 1853 [Scoi2010]幸运数字
题链: http://www.lydsy.com/JudgeOnline/problem.php?id=1853 题解: 容斥原理,暴力搜索,剪枝(这剪枝剪得真玄学) 首先容易发现,幸运号码不超过 2 ...
【BZOJ 1853】 1853: [Scoi2010]幸运数字（容斥原理）
1853: [Scoi2010]幸运数字 Time Limit: 2 Sec Memory Limit: 64 MBSubmit: 2472 Solved: 911 Description 在中国 ...
1853: [Scoi2010]幸运数字[容斥原理]
1853: [Scoi2010]幸运数字 Time Limit: 2 Sec Memory Limit: 64 MBSubmit: 2405 Solved: 887[Submit][Status] ...
BZOJ2393 & 1853 [Scoi2010]幸运数字【搜索 + 容斥】
题目在中国,很多人都把6和8视为是幸运数字!lxhgww也这样认为,于是他定义自己的"幸运号码"是十进制表示中只包含数字6和8的那些号码,比如68,666,888都是" ...
AC日记——[SCOI2010]幸运数字 bzoj 1853
1853: [Scoi2010]幸运数字 Time Limit: 2 Sec Memory Limit: 64 MBSubmit: 2405 Solved: 887[Submit][Status] ...

随机推荐

js 增删改查方法
push() 向数组的末尾添加一个或多个元素 pop() 删除数组内部并返回数组的最后一个元素 shift() 把数组内部的第一个元素从其中删除,并返回第一个元素的值 unshift() 向数组外部的 ...
codevs 1033 蚯蚓的游戏问题
Description 在一块梯形田地上,一群蚯蚓在做收集食物游戏.蚯蚓们把梯形田地上的食物堆积整理如下: a(1,1) a(1,2)…a(1,m) a(2,1) a(2,2) a(2,3)…a ...
java反射详解（转）
本篇文章依旧采用小例子来说明,因为我始终觉的,案例驱动是最好的,要不然只看理论的话,看了也不懂,不过建议大家在看完文章之后,在回过头去看看理论,会有更好的理解. 下面开始正文. [案例1]通过一个对象 ...
redis 学习笔记(2)-client端示例代码
redis提供了几乎所有主流语言的client,java中主要使用二种:Jedis与Redisson 一.Jedis的使用 <dependency> <groupId>redi ...
VMware Fusion DHCP方式下如何指定虚拟机IP地址
默认情况下,vmware fusion中的虚拟机,网卡设置成dhcp(动态分配 )时,会分配一个IP地址,但这个IP通常很难记,如果我们想为某台虚拟机挑一个好记的IP地址,可以按如下步骤操作: 命令行 ...
无法连接windows虚拟机oracle的解决办法
在mac机上玩基于oracle db的开发真心不容易,oracle公司死活不出oracle express edition for mac OS,曾经发布过的oracle 10 for mac下载地址 ...
React入门 (1)—使用指南（包括ES5和ES6对比）
前言本篇会简明扼要的介绍一下React的使用方法.代码会用JSX+ES5和JSX+ES6两种方式实现. React简介 React来自Facebook,于2013年开源.至今不断修改完善,现在已经到 ...
大学回顾和C与PHP之路
我去年毕业,从事PHP学习和开发一年多. background:medical muti-media electric web; 先讲一下我的背景吧,我大学的学校是一个医科学校,然而专业是计算机动漫设 ...
shell note
1 输出重定向:ll > aaa 将输出内容添加到aaa文件中 ll >> aaa将输出内容追加到aaa中 ll &>> abc 将输出内容不论正确或错误都保存 ...
把php上传sae问题要使用IO
应用移植指南一,为什么要移植应用 SAE禁止IO写操作,代码目录不能写入.这意味着普通程序的上传图片.生成缓存等操作都不能在SAE上正常运行,这时候你需要对这些代码进行修改后才能让你的程序运行在SA ...

BZOJ 1853: [Scoi2010]幸运数字