bzoj1853幸运数字——容斥原理

题目：http://www.lydsy.com/JudgeOnline/problem.php?id=1853

dfs实现容斥原理即可。

注意：若在init中写“cnt++”，则出来后需要先cnt--再继续！！

代码如下：

#include<iostream>
#include<cstdio>
#include<algorithm>
using namespace std;
typedef long long ll;
ll l,r,a[],b[],ans,cnt,ct;
bool del[];
void init(ll x)
{
    if(x>r)return;
    a[cnt++]=x;
    init(*x+);
    init(*x+);
}
ll gcd(ll x,ll y)
{
    return x%y==?y:gcd(y,x%y);
}
void dfs(ll x,ll y,ll z)
{
    if(x>ct)
    {
        if(y&)ans+=r/z-(l-)/z;
        else if(y)ans-=r/z-(l-)/z;
        return;
    }
    dfs(x+,y,z);
    ll tmp=z/gcd(a[x],z);
    if((double)a[x]*tmp<=r)dfs(x+,y+,a[x]*tmp);
}
int main()
{
    scanf("%lld%lld",&l,&r);
    init();
    cnt--;//!!!
    sort(a+,a+cnt+);
    for(ll i=;i<=cnt;i++)
    {
        if(del[i])continue;
        for(ll j=i+;j<=cnt;j++)
            if(a[j]%a[i]==)del[j]=;
        b[++ct]=a[i];
    }
    for(ll i=;i<=ct;i++)//倒序！
        a[i]=b[ct-i+];
    dfs(,,);
    printf("%lld",ans);
    return ;
}