bzoj2393: Cirno的完美算数教室

链接：http://www.lydsy.com/JudgeOnline/problem.php?id=2393

题意：中文题。。。

分析：在10^10以内baka数大概只有几百个设为k个，预处理出来，然后将能互相被整除的两个中去掉大的那个，然后只要容斥就好了，奇数个减偶数个加，因为数字增长得比较快所以其实不会达到2^k，因为当当前这个要处理的lcm>n时我们直接return掉即可，其实在最大的范围内最多就只能同时选7个左右的样子而且还是最小的那几个。

代码：

#include<map>
#include<set>
#include<cmath>
#include<queue>
#include<math.h>
#include<cstdio>
#include<vector>
#include<string>
#include<cstring>
#include<iostream>
#include<algorithm>
#pragma comment(linker, "/STACK:102400000,102400000")
using namespace std;
const int N=1020;
const int MAX=151;
const int MOD1=100000007;
const int MOD2=100000009;
const double EPS=0.00000001;
typedef long long ll;
const ll MOD=9999991;
const ll INF=10000000010;
typedef unsigned long long ull;
int k,q[1050],d[1050];
void deal() {
    int i,j,l=1,r=1,L,R;
    memset(q,0,sizeof(q));
    d[1]=0;q[1]=1;k=0;
    for (i=1;i<10;i++) {
        L=l;R=r;
        for (j=L;j<=R;j++) {
            d[++r]=d[j]*10+2;d[++r]=d[j]*10+9;
        }
        l=R+1;
    }
    for (i=2;i<=r;i++)
        for (j=i+1;j<=r;j++)
        if (d[j]%d[i]==0) q[j]=1;
    for (i=2;i<=r;i++)
    if (!q[i]) d[++k]=d[i];
}
ll ret;
ll gcd(ll a,ll b) {
    return b ? a:gcd(b,a%b);
}
void dfs(int a,int b,ll x,ll n) {
    if (x>n) return ;
    if (a==k+1) {
        if (b&1) ret-=n/x;
        else ret+=n/x;
        return ;
    }
    dfs(a+1,b,x,n);
    dfs(a+1,b+1,x/gcd(x,d[a])*d[a],n);
}
ll get(ll n) {
    ret=0;
    dfs(1,0,1,n);
    return ret;
}
int main()
{
    ll l,r;
    deal();
    scanf("%lld%lld", &l, &r);
    printf("%lld\n", r-get(r)-l+1+get(l-1));
    return 0;
}