题意:中文题。。。
分析:在10^10以内baka数大概只有几百个设为k个,预处理出来,然后将能互相被整除的两个中去掉大的那个,然后只要容斥就好了,奇数个减偶数个加,因为数字增长得比较快所以其实不会达到2^k,因为当当前这个要处理的lcm>n时我们直接return掉即可,其实在最大的范围内最多就只能同时选7个左右的样子而且还是最小的那几个。
代码:
#include<map>
#include<set>
#include<cmath>
#include<queue>
#include<math.h>
#include<cstdio>
#include<vector>
#include<string>
#include<cstring>
#include<iostream>
#include<algorithm>
#pragma comment(linker, "/STACK:102400000,102400000")
using namespace std;
const int N=1020;
const int MAX=151;
const int MOD1=100000007;
const int MOD2=100000009;
const double EPS=0.00000001;
typedef long long ll;
const ll MOD=9999991;
const ll INF=10000000010;
typedef unsigned long long ull;
int k,q[1050],d[1050];
void deal() {
int i,j,l=1,r=1,L,R;
memset(q,0,sizeof(q));
d[1]=0;q[1]=1;k=0;
for (i=1;i<10;i++) {
L=l;R=r;
for (j=L;j<=R;j++) {
d[++r]=d[j]*10+2;d[++r]=d[j]*10+9;
}
l=R+1;
}
for (i=2;i<=r;i++)
for (j=i+1;j<=r;j++)
if (d[j]%d[i]==0) q[j]=1;
for (i=2;i<=r;i++)
if (!q[i]) d[++k]=d[i];
}
ll ret;
ll gcd(ll a,ll b) {
return b ? a:gcd(b,a%b);
}
void dfs(int a,int b,ll x,ll n) {
if (x>n) return ;
if (a==k+1) {
if (b&1) ret-=n/x;
else ret+=n/x;
return ;
}
dfs(a+1,b,x,n);
dfs(a+1,b+1,x/gcd(x,d[a])*d[a],n);
}
ll get(ll n) {
ret=0;
dfs(1,0,1,n);
return ret;
}
int main()
{
ll l,r;
deal();
scanf("%lld%lld", &l, &r);
printf("%lld\n", r-get(r)-l+1+get(l-1));
return 0;
}