2393: Cirno的完美算数教室 容斥原理

这题好像做过一次。。想复习一下容斥原理。。
qwq其实自己随便yy一下就好了。。没什么好复习的。。

#include<iostream>#include<cstdio>#include<algorithm>#define ll long long using namespace std;int l,r,cnt,tot;ll ans;ll a[2005],b[2005];bool vis[2005];inline void pre(ll x){    if (x>r) return;    a[++cnt]=x;    pre(x*10+2); pre(x*10+9);}ll gcd(ll a,ll b){    return b==0?a:gcd(b,a%b);}void dfs(int x,int y,ll z){    if (x>tot)    {        if (y&1) ans+=r/z-(l-1)/z;        else if (y) ans-=r/z-(l-1)/z;        return;    }    dfs(x+1,y,z);    ll tmp=b[x]*z/gcd(b[x],z);    if (tmp<=r) dfs(x+1,y+1,tmp);}int main(){    scanf("%d%d",&l,&r);    pre(2); pre(9);     sort(a+1,a+cnt+1);    for (int i=1;i<=cnt;i++)        if (!vis[i])        {            b[++tot]=a[i];            for (int j=i+1;j<=cnt;j++)                if (!(a[j]%a[i])) vis[j]=1;        }    dfs(1,0,1);    cout << ans << endl;    return 0;}