Problem:
問題 E Fox Number

問題文

きつねのしえるは数字について考えるのが好きである．ある日，しえるはこの世にねこの数字というものがあることを知り，自分もそのような数字が欲しくなってしまった．そこで，次の性質を満たす数字をとりあえず Fox Number と呼ぶことにした．

[性質]: 整数 N は k (≥ 1) 個の素数 p1, ..., pk と正整数 e1, ..., ek で，p1 < ... < pk，e1 ≥ ... ≥ ek を満たすものによって N = p1e1 × ... × pkek と書けるとき，Fox Number であると呼ぶことにする．

ところでしえるはこの性質をその場の思いつきで決めてしまったのでこれがどういう性質を持っているのか全くわかっていない．そこで，この数字の性質を調べるために区間 [A-B, A+B] に含まれる Fox Number の個数を出力して欲しい．

入力形式

1 行目に以下の 2 つの整数が与えられる．

A B

出力形式

1 行目に，A-B 以上 A+B 以下の Fox Number の個数を出力せよ．

制約

1 ≤ A ≤ 1012

0 ≤ B ≤ 5 × 105

入出力例

入力例 1

18 2

出力例 1

4

16 = 24, 17 = 171, 19 = 191, 20 = 22 × 51 は Fox Number である．一方で 18 = 2 × 32 は Fox Number ではない．合計して 4 つの Fox Number がある．

入力例 2

100 10

出力例 2

18

入力例 3

123456 789

出力例 3

1464

Solution in C++ (judge verdict: Accepted)

#include <iostream>
#include <sstream>
#include <string>
#include <algorithm>
#include <vector>
#include <stack>
#include <queue>
#include <set>
#include <map>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <cmath>
#include <cassert>

using namespace std;

#define FOR(i,k,n) for(int i=(k); i<(int)(n); ++i)
#define REP(i,n) FOR(i,0,n)
#define FORIT(i,c) for(__typeof((c).begin())i=(c).begin();i!=(c).end();++i)

template<class T> void debug(T begin, T end){ for(T i = begin; i != end; ++i) cout<<*i<<" "; cout<<endl; }
inline bool valid(int x, int y, int W, int H){ return (x >= 0 && y >= 0 && x < W && y < H); }

typedef long long ll;
const int INF = 100000000;
const double EPS = 1e-8;
const int MOD = 1000000007;
int dx[8] = {1, 0, -1, 0, 1, -1, -1, 1};
int dy[8] = {0, 1, 0, -1, 1, 1, -1, -1};
int count_fox_number(ll a, ll b){
  if(a < 2) a = 2;
  if(!(a < b)) return 0;
  int ans = b - a;
  bool prime_small[1000100];
  for(int i = 2; (ll) i * i < b; i++) prime_small[i] = true;
  int cnt_pow[1000100];
  REP(i, b - a) cnt_pow[i] = INF;
  for(int i = 2; (ll)i * i < b; i++){
    if(prime_small[i]){
      for(int j = 2 * i; (ll) j * j < b; j += i) prime_small[j] = false;
      for(ll j = max(2LL, (a + i - 1)/i) * i; j < b; j += i){
        if(cnt_pow[j - a] == -1) continue;
        ll t = j;
        int cnt = 0;
        while(t % i == 0){
          t /= i;
          cnt++;
        }
        assert(cnt != 0);
        if(cnt_pow[j - a] < cnt){
          cnt_pow[j - a] = -1;
          ans--;
        }else{
          cnt_pow[j - a] = cnt;
        }
      }
    }
  }
  return ans;
}

int main(){
  ll a, b;
  while(cin>>a>>b){
    cout<<count_fox_number(a - b, a + b + 1)<<endl;
  }
  return 0;
}