10.06 容斥练习T3 简单容斥+去重搜索

Description

~Cirno发现了一种baka数，这种数呢~只含有2和⑨两种数字~~
现在Cirno想知道~一个区间中~~有多少个数能被baka数整除~
但是Cirno这么天才的妖精才不屑去数啦
只能依靠聪明的你咯。

Input

一行正整数L R( 1 < L < R < 10^10)

Output

一个正整数，代表所求的答案

Sample Input

1 100

Sample Output

58

 

 

 

与之前的SCOI的幸运数字其实是一个道理，就是预处理出所有2,9相关的数字，然后吧这些数字去重之后放进去跑（如果不去重会相当的恐怖，18000+个去重后446个）搜索

搜索非常普通，注意往后转移的时候是lcm，其他的参见幸运数字（一次就过了好少见）

code：

#include<iostream>
#include<cstdio>
#include<algorithm>
#define N 100005
using namespace std;
long long ans,a[N],l,r,tot,check[N];
void pre(long long dep,long long val){
//    cout<<val<<endl;
//    system("pause");
    if(dep==9){
        if(!val)return;
        a[++tot]=val;
        return;
    }
    pre(dep+1,val*10+2);
    pre(dep+1,val*10+9);
    pre(dep+1,val);
}
void uni(){
    sort(a+1,a+tot+1);
    tot=unique(a+1,a+tot+1)-a-1;
    long long temp=0;
    long long b[20000];
    for(long long i=1;i<=tot;i++)b[i]=a[i];
    for(long long i=1;i<=tot;i++){
        if(check[i])continue;
        a[++temp]=b[i];
        for(long long j=i+1;j<=tot;j++){
            if(b[j]%b[i])continue;
            check[j]=1;
        }    
    }
    tot=temp;
    return;
}
long long gcd(long long a,long long b){
    if(a<b)swap(a,b);
    while(a=a%b)swap(a,b);
    return b;
}
long long lcm(long long a,long long b){
    return a/gcd(a,b)*b;
}
void dfs(long long dep,long long val,int now){
    if(val>r)return;
    if(dep==tot){
        if(val==1)return;
    //    cout<<val<<endl;
        if(now%2==1)
            ans+=r/val-(l-1)/val;    
        else
            ans-=r/val-(l-1)/val;
        return ;
    }
    dfs(dep+1,lcm(val,a[dep+1]),now+1);
    dfs(dep+1,val,now);
}
int main(){
    pre(0,0);
    uni();
//    for(int i=1;i<=tot;i++){
//        cout<<a[i]<<endl;
//        if(i%100==0)
//        system("pause");
//    }
    cin>>l>>r;
    dfs(0,1,0);
    cout<<ans;
    return 0;
}

over