bzoj2393

女友问我，就去补习了一下容斥原理
要求1到n有多少个数能被5或8整除，5和8的最小公倍数为40,，那么答案就是能被5整除的个数加能被8整除的个数减去能被40整除的个数
同理，在tol个数中，求出任意m（1<=m<=tol）个数的最小公倍数，如果m是奇数,如{5}、{8}，那把所有能被这个公倍数的数的个数加到答案里，否则（如{5,8}）减去
显然可以用二进制枚举枚举所有组合，但是这样会超时或爆空间，考虑到求n里有多少个数能被c整除，是直接算n/c，当c>n时加减n/c都无意义，故可以写成dfs再进行剪枝

#include <iostream>
#include <algorithm>
#include <cstdio>
#include <cstring>
#include <math.h>
#include <map>
#include <set>
#include <queue>
using namespace std;
typedef long long LL;
LL L,R;
LL a[2000],ans,b[2000];
int tol;
LL gcd(LL a,LL b){return b?gcd(b,a%b):a;}
LL lcm(LL a,LL b){return a/gcd(a,b)*b;}
void get29()
{
    int l,r;
    l=1,r=2;
    a[1]=2,a[2]=9;
    tol=2;
    //找到这些数
    while(a[tol]<R)//通过层次遍历二叉树来生成每一个节点的左右节点，这样写有点丑，但是算出来的数组是有序的
    {
        for(int i=l;i<=r;i++)
        {
            a[++tol]=a[i]*10+2;
            a[++tol]=a[i]*10+9;
        }
        l=r+1;
        r=tol;
    }
    while(a[tol]>R)tol--;

    //删掉倍数
    int dtol=tol;
    tol=0;
    for(int i=1;i<=dtol;i++)
    {
        b[++tol]=a[i];
        for(int j=i-1;j>=1;j--)
        {
            if(!(a[i]%a[j]))
            {
                tol--;
                break;
            }
        }
    }
}
void dfs(int idx,int sign,LL z)//idx表示数组下标，sign表示+或-，z是最小公倍数
{
    if(idx>tol)
    {
        if(z!=1)ans=ans+sign*(R/z-(L-1)/z);
        return ;
    }
    dfs(idx+1,sign,z);
    LL zz=lcm(z,b[idx]);
    if(zz>R)return ;             //此处为剪枝
    dfs(idx+1,-sign,zz);
}
int main()
{
    cin>>L>>R;
    get29();
    dfs(1,-1,1);
    cout<<ans;
    return 0;
}