BZOJ2393: Cirno的完美算数教室

2393: Cirno的完美算数教室

Time Limit: 10 Sec  Memory Limit: 128 MB
Submit: 138  Solved: 83
[Submit][Status]

Description

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

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

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

只能依靠聪明的你咯。

 

 

Input

一行正整数L R

( 1 < L < R < 10^10)

 

Output

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

 

Sample Input

1 100

Sample Output

58

HINT

题解：

同幸运数字，第一次排rank1，哈哈！

代码：

#include<cstdio>

#include<cstdlib>

#include<cmath>

#include<cstring>

#include<algorithm>

#include<iostream>

#include<vector>

#include<map>

#include<set>

#include<queue>

#include<string>

#define inf 1000000000

#define maxn 1200

#define maxm 500+100

#define eps 1e-10

#define ll unsigned long long

#define pa pair<int,int>

#define for0(i,n) for(int i=0;i<=(n);i++)

#define for1(i,n) for(int i=1;i<=(n);i++)

#define for2(i,x,y) for(int i=(x);i<=(y);i++)

#define for3(i,x,y) for(int i=(x);i>=(y);i--)

#define mod 1000000007

using namespace std;

inline ll read()

{

    ll x=0,f=1;char ch=getchar();

    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}

    while(ch>='0'&&ch<='9'){x=10*x+ch-'0';ch=getchar();}

    return x*f;

}
ll l,r,ans,tot,n,a[maxn],b[maxn];
ll v[maxn];
inline void dfs(ll x)
{
    if(x>r)return;
    if(x)a[++tot]=x;
    dfs(10*x+2);
    dfs(10*x+9);
}    
inline ll gcd(ll x,ll y)
{
    return y?gcd(y,x%y):x;
}
inline void calc(ll x,int y,int z)
{
    if(y>n)
    {
            if(z&1)ans+=r/x-(l-1)/x;
            else if(z)ans-=r/x-(l-1)/x;
            return;    
    }
    calc(x,y+1,z);
    ll t=(x/gcd(x,a[y]))*a[y];
    if(t<=r)calc(t,y+1,z+1);
}    

int main()

{

    freopen("input.txt","r",stdin);

    freopen("output.txt","w",stdout);

    l=read();r=read();
    dfs(0);
    sort(a+1,a+tot+1);
    for1(i,tot)
    if(!v[i])
    {
        b[++n]=a[i];
        for2(j,i+1,tot)if(a[j]%a[i]==0)v[j]=1;
    }
    for1(i,n)a[n+1-i]=b[i];
    calc(1,1,0);
    printf("%lld\n",ans);

    return 0;

}