[bzoj1856][Scoi2010]幸运数字

1853: [Scoi2010]幸运数字

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

sol:
我们发现幸运号码显然很少，最多有2^10*2个。
那么对于这些幸运号码，某些是另一些的倍数，那这些倍数可以先去掉。
那么对于剩下的幸运号码，他们互相的lcm可能会重复，那我们容斥原理。+至少奇数个幸运号码的lcm-至少偶数个幸运号码的lcm。

这题因为要卡常所以我写了一些糟糕的东西。看起来很gay的东西就是来卡常的。

#include<iostream>
#include<algorithm>
#include<string>
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<cmath>
#define ll long long
using namespace std;
const int N=300000;
ll n,l,r;

inline int read()
{
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
    return x*f;
}
ll number[N],cnt=-1;
void get_num(ll num)
{
    if(num>r) return;
    number[++cnt]=num;
    get_num(num*10+6);
    get_num(num*10+8);
}
ll gcd(ll m,ll n)
{
    return m%n?gcd(n,m%n):n;
}
ll ans;
void dfs(int now,int num,ll lcm)
{
    if(now>n)
    {
        if(!num) return;
        if(num&1) ans+=r/lcm-(l-1)/lcm;
        else ans-=r/lcm-(l-1)/lcm;
        return;
    }
    dfs(now+1,num,lcm);
    ll tmp=lcm/gcd(number[now],lcm);
    if((double)tmp*number[now]<=r)
    dfs(now+1,num+1,tmp*number[now]);
}
inline bool cmp(const ll a,const ll b)
{
    return a>b;
}
bool vis[N];
int main()
{
//  freopen("1853.txt","r",stdin);
//  freopen(".out","w",stdout);
    int i,j;
    scanf("%lld%lld",&l,&r);
    get_num(0);
    sort(number+1,number+1+cnt);
    for(int i=1;i<=cnt;++i)
    if(!vis[i])
    for(int j=i+1;j<=cnt;++j)
    if(number[j]%number[i]==0) vis[j]=1;
    n=0;
    for(int i=1;i<=cnt;++i)
    if(!vis[i]) number[++n]=number[i];
    sort(number+1,number+1+n,cmp);
    dfs(1,0,1);
    printf("%lld",ans);
}