bzoj 1853 容斥 + 搜索

思路：先把所有幸运数字找出来， 把没有用的去掉，然后爆搜容斥，因为最多只会搜十几个就超过限制了，

所以是可行的。

#include<bits/stdc++.h>
#define LL long long
#define fi first
#define se second
#define mk make_pair
#define PII pair<int, int>
#define y1 skldjfskldjg
#define y2 skldfjsklejg
using namespace std;

const int N =  + ;
const int M = 5e5 + ;
const int inf = 0x3f3f3f3f;
const LL INF = 0x3f3f3f3f3f3f3f3f;
const int mod = ;

LL a[N], b[N], L, R, ans;
int tot, top;

void init(LL val) {
    if(val > R) return;
    if(val) a[++tot] = val;
    init(val *  + );
    init(val *  + );
}

void dfs(int p, int cnt, LL lcm) {
    if(p > top) {
        if(cnt & ) ans += R / lcm - (L - ) / lcm;
        else if(cnt) ans -= R / lcm - (L - ) / lcm;
        return;
    }
    dfs(p + , cnt, lcm);
    long long tmp = lcm / __gcd(b[p], lcm);
    if(1.0 * b[p] * tmp <= R) dfs(p + , cnt + , b[p] * tmp);
}

int main() {
    scanf("%lld%lld", &L, &R);
    init();
    sort(a + , a +  + tot);
    for(int i = ; i <= tot; i++) {
        bool flag = true;
        for(int j = ; j <= top; j++) {
            if(a[i] % b[j] == ) {
                flag = false;
                break;
            }
        }
        if(flag) b[++top] = a[i];
    }
    reverse(b + , b +  + top);
    dfs(, , );

    printf("%lld\n", ans);
    return ;
}