bzoj1853 bzoj2393

两题是类似的，这里说一下bzoj1853

首先我们求出所有的幸运号码，注意如果存在x是y的倍数则x不算在内，避免之后重复计算

下面我们就要统计幸运号码的倍数了，这显然是要用到容斥原理的

但是幸运号码很多，如果直接暴力找几个幸运号码的公倍数做容斥原理弄会TLE的；

因此我们想到在搜索中剪枝，如果几个幸运号码的公倍数已经大于r，

那么我们一定不会再用这几个幸运号码和别的幸运号码求公倍数了

为了体现这个剪枝的威力，我们穷举幸运号码的时候应该从大往小搜索；

var b:array[0..10010] of int64;
    a:array[0..20] of longint;
    m,t,k,i:longint;
    s,x,l,r,ans:int64;
    f:boolean;

function gcd(a,b:int64):int64;
  begin
    if b=0 then exit(a)
    else exit(gcd(b,a mod b));
  end;

procedure dfs(j,t:longint;x:int64);
  var y:int64;
  begin
    if j=0 then
    begin
      if t mod 2=1 then ans:=ans+r div x-(l-1) div x  //容斥原理
      else if t>0 then ans:=ans-r div x+(l-1) div x
    end
    else begin
      dfs(j-1,t,x);
      y:=x div gcd(b[j],x);
      if double(b[j])*double(y)<=r then  //double是防止爆int64
        dfs(j-1,t+1,b[j]*y);
    end;
  end;

begin
  readln(l,r);
  t:=0;
  x:=r;
  while x<>0 do
  begin
    inc(t);
    x:=x div 10;
  end;
  fillchar(a,sizeof(a),255);
  m:=0;
  while a[0]=-1 do
  begin
    s:=0;
    for i:=1 to t do
      if a[i]=0 then s:=s*10+6
      else if a[i]=1 then s:=s*10+8;
    if s>r then break
    else if s<>0 then
    begin
      f:=true;
      for i:=1 to m do
        if s mod b[i]=0 then
        begin
          f:=false;
          break;
        end;
      if f then
      begin
        inc(m);
        b[m]:=s;
      end;
    end;
    k:=t;
    while a[k]=1 do
    begin
      a[k]:=0;
      dec(k);
    end;
    inc(a[k]);
  end;
  dfs(m,0,1);
  writeln(ans);
end.