=== ===
这里放传送门
=== ===
题解
(这题目让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;
}
偏偏在最后出现的补充说明
很多比较简单的容斥题目都和搜索相关联啦。有时候想到搜索的思路还是要去写一写的,说不定加上剪枝以后复杂度还是比较科学的呢。。并且打表也是很重要的解题过程啦,不然这个题不打表的话谁知道它符合条件的数那么少啊。。。