【bzoj2393】Cirno的完美算数教室 数论容斥

Description

~Cirno发现了一种baka数，这种数呢~只含有2和⑨两种数字~~

现在Cirno想知道~一个区间中~~有多少个数能被baka数整除~

但是Cirno这么天才的妖精才不屑去数啦

只能依靠聪明的你咯。

Input

一行正整数L R

( 1 < L < R < 10^10)

Output

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

Sample Input

1 100

Sample Output

58

HINT

此题数据范围应该是10^9

题解：

　　处理处所有数然后容斥。

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

 #define ll long long
 using namespace std;
 inline int read()
 {
     int x=,f=;char ch=getchar();
     while(ch>''||ch<''){if (ch=='-') f=-;ch=getchar();}
     while(ch<=''&&ch>=''){x=(x<<)+(x<<)+ch-'';ch=getchar();}
     return x*f;
 }

 int L,R;
 int n,m,ans;
 ll a[],b[];
 bool vis[];

 void pre(ll x,int R)
 {
     if (x>R) return;
     if (x) a[++m]=x;
     pre(x*+,R);pre(x*+,R);
 }
 ll gcd(ll a,ll b){return b?gcd(b,a%b):a;}
 void dfs(int x,int y,ll z)
 {
     if (x>n)
     {
         if (y&) ans+=R/z-(L-)/z;
         else if (y) ans-=R/z-(L-)/z;
         return;
     }
     if (z>R) return;
     dfs(x+,y,z);
     z=z/gcd(a[x],z)*a[x];
     dfs(x+,y+,z);
 }
 int main()
 {
     L=read();R=read();
     pre(0LL,R);
     sort(a+,a+m+);
     for (int i=;i<=m;i++)
         if (!vis[i])
         {
                b[++n]=a[i];
             for (int j=i+;j<=m;j++)if (a[j]%a[i]==) vis[j]=;
         }
     for (int i=;i<=n;i++)
         a[n-i+]=b[i];
     dfs(,,);
     printf("%d",ans);
 }