[BZOJ2393]Cirno的完美算数教室（容斥原理+DFS）

=== ===

这里放传送门

=== ===

题解

（这题目让ATP想到了菜鸟ATP在刚开始打非想天则的故事模式的时候被电脑easy模式虐死的惨状。。。。。）

容斥的思路比较显然，就是至少是一个数的倍数-至少是两个数的倍数+至少是三个数的倍数……并且如果预处理所有只含有2和9的数字的话它最多只有2000多个。那么直接用dfs来构造数字就很科学了，因为一旦当前维护的数字x大于给定的R就可以剪枝退出了。

甚至当处理出全部的2000多个数以后，如果存在两个数a和b并且b是a的倍数的话就可以把b扔掉，因为a的倍数肯定会全部包含b的倍数，所以b就不用单独计算了。这样的话剩下的数最多只有500个左右。那么继续爆搜，在搜索过程中枚举选哪几个数字，同时维护最小公倍数。因为最小公倍数的倍数一定同时都是选定的这几个数的倍数。同样，当前数字大于R的时候就剪枝，这样就可以很快出解了！

注意的问题是在计算答案的时候要判断选择的数字是不是为0再减，不然它会莫名其妙在答案里多减一个R。

代码

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
int l,r,num[3010],cnt,n;
long long ans,f[3010];
long long gcd(long long a,long long b){
    long long r=a%b;
    while (r!=0){a=b;b=r;r=a%b;}
    return b;
}
void dfs(long long x){
    if (x>r) return;
    num[++cnt]=x;
    dfs(x*10+2);dfs(x*10+9);
}
void dfs_again(int i,int c,long long lcm){
    if (i==n+1){
        if (c&1) ans=ans+(r/lcm-(l-1)/lcm);
        else if (c!=0) ans=ans-(r/lcm-(l-1)/lcm);
        return;
    }
    dfs_again(i+1,c,lcm);
    long long tmp=lcm*f[i]/gcd(lcm,f[i]);//计算最小公倍数
    if (tmp<=r) dfs_again(i+1,c+1,tmp);
}
int main()
{
    scanf("%d%d",&l,&r);
    cnt=-1;dfs(0);
    for (int i=1;i<=cnt;i++)
      if (num[i]!=0){
        f[++n]=num[i];
        for (int j=i+1;j<=cnt;j++)
          if (num[j]%num[i]==0) num[j]=0;
      }
    sort(f+1,f+n+1);
    dfs_again(1,0,1);
    printf("%I64d\n",ans);
    return 0;
}

偏偏在最后出现的补充说明

很多比较简单的容斥题目都和搜索相关联啦。有时候想到搜索的思路还是要去写一写的，说不定加上剪枝以后复杂度还是比较科学的呢。。并且打表也是很重要的解题过程啦，不然这个题不打表的话谁知道它符合条件的数那么少啊。。。