【bzoj1853】 SCOI2010幸运数字 容斥原理

跪烂晨神，一上来先想数位dp，结果根本没法做，然后看了看晨神的题解，发现幸运数字一共就2000来个，然后加一些剪枝，用容斥原理筛一下就好了，其实容斥就是选了奇数个就加，选了偶数个就减，要从大到小搜，这样容易爆，然后一定要判断如果爆了longlong，或者大于r的时候就退出，如果一个数是另一个数的倍数，那么就不用考虑这个数了。

#include<cstdio>#include<cstring>#include<cstdlib>#include<cmath>#include<iostream>#include<algorithm>using namespace std;long long q[5010];int tot;long long l,r;bool vis[5010];long long ans;void dfs(long long x){if (x>r) return;if (x) q[++tot]=x;dfs(x*10+6);dfs(x*10+8);}bool cmp(long long a,long long b){return a>b;}long long gcd(long long a,long long b){if (b==0) return a;else return gcd(b,a%b);}void dfs1(int x,int d,long long lcm){if (x==tot) return;x++;dfs1(x,d,lcm);d++;if (vis[x]) return;lcm=lcm*q[x]/gcd(lcm,q[x]);if (lcm<=0 ||lcm>r) return;if (d&1) ans+=r/lcm-(l-1)/lcm;else ans-=r/lcm-(l-1)/lcm;dfs1(x,d,lcm);}int main(){scanf("%lld%lld",&l,&r);dfs(0);sort(q+1,q+tot+1,cmp);for (int i=1;i<=tot;i++)   for (int j=i+1;j<=tot;j++)    if (q[i]%q[j]==0) vis[i]=1;dfs1(0,0,1);printf("%lld\n",ans);return 0;}