dfs版容斥原理+剪枝——bzoj1853

学了一种爆搜版+剪枝的容斥方法，即类似数位dp时按位进行容斥，同时需要在搜索过程中进行剪枝

/*
容斥原理，先在打出的表里筛掉所有倍数，然后用容斥原理+1个的倍数-2个lcm的倍数+3个lcm的倍数...
注意剪枝，判断防止爆long long
*/
#include<bits/stdc++.h>
using namespace std;
#define ll long long
#define maxn 200005
ll l,r,ans;
ll a[maxn],m,b[maxn],n;
void init(ll x){
    if(x>r)return;
    a[++m]=x;
    init(x*+);
    init(x*+);
}

void dfs(int pos,int num,ll now){
    if(pos>n){
        if(num%)
            ans+=r/now-(l-)/now;
        else if(num)ans-=r/now-(l-)/now;
        return;
    }
    dfs(pos+,num,now);//不取第pos位
    ll tmp=now/__gcd(b[pos],now);//防止爆ll
    if((double)b[pos]*tmp<=(double)r)
        dfs(pos+,num+,b[pos]*tmp);
}

int vis[maxn];
int main(){
    cin>>l>>r;
    init();init();
    sort(a+,a++m);//排成有序数列 

    for(int i=;i<=m;i++)
        if(vis[i]==){
            for(int j=i+;j<=m;j++)
                if(a[j]%a[i]==)
                    vis[j]=;
        }
    for(int i=m;i>=;i--)
        if(!vis[i])b[++n]=a[i];
    dfs(,,);
    cout<<ans<<endl;
}