Problem:


宝くじ

ある国の王様は素数とギャンブルが大好きです。この国の通貨の単位はプライムといいます。2007 年 11 月 1 日現在のプライムのクロス円レートは 9973 とちょうど素数になっていたので、王様は大喜びです。この国では 1/101 プライムを1サブプライムとする補助貨幣が使われています。

この国の政府は、国家財政の安定、国民の娯楽、王様の趣味を同時に満足させることを目的に宝くじ振興公社を設立し、宝くじ事業を始めることにしました。素数が大好きな王様は、素数が話題になればそれだけで満足して、賞金をあげたくなります。そこで振興公社は次のような宝くじを考えました。

くじには 0 から MP までの一連の番号がつけられている。MP はこの国で知られている最大の素数で、毎月一日に官報で告示される。同時に、MP 以下のすべての素数も発表される。それ以上大きな素数が発見されても、その月の間は、宝くじ振興公社を含む全ての公的な機関では、一日に発表されたMP を最大の素数として扱う。2007 年 11 月 1 日 にはMP = 999983 (1000000 以下の最大の素数) が発表された。

宝くじの販売価格は 1 サブプライム。

宝くじの抽選では、当たりくじの番号 p と賞金算出のための数 m の組(p, m) が何組か選ばれる。p, m はそれぞれ0 以上 MP 以下の整数。

この国で知られている素数の中で p - m 以上 p + m 以下のものが X 個あるとすると、抽選結果(p, m) の賞金は、X プライムとなる。

抽選結果 (p, m) の賞金 X プライムは番号 p の宝くじを持っている当選者に支払われるが、X = 0 の場合もあり、この場合には当選者には賞金は支払われない。

賞金のうち 1 プライムは宝くじの売り上げから支出され、X - 1 プライムは王様の内廷費から支出される（王様が払ってくれる）。X = 0 ならば宝くじの売り上げから 1 プライムが内廷費に繰り入れられる。（王様に対し支払われる）

ひとつの番号 p が複数の当たりくじになることもある。この場合にはそれぞれの抽選結果(p, m) から算出される賞金の合計が当選者に支払われる。

このくじでは当たりくじを買った人は抽選結果から関係する素数を探し、その個数を数えるので、国民の話題に素数がよく出てきて王様はとてもご満悦です。宝くじ振興公社としては当たりくじ1本当たり公社負担が 1 プライム、販売価格が 1 サブプライムだから、当たりくじの本数 n を、販売した宝くじ 101 本あたり 1 本以下となるようにすれば (すなわち、n ≤ (販売した本数)/101 とすれば) 赤字にはなりません。

問題は内廷費からの支出額です。あなたの仕事は、抽選結果を入力として、2007 年 11 月における宝くじ振興公社が王様に請求する賞金の額を出力するプログラムを作成することです。ただし、請求する賞金の額が負になることはないものとします。

注意

この国における素数の定義も日本の学校教育で学習する内容と同じです。即ち、素数とは 1 と自分自身以外の約数を持たない自然数をいいます（なお、1 は素数ではありません）。

我々は 1000003 が素数であることを知っていますが、この国では 2007 年 11 月段階では知られていません。そのため、999963 以上 1000003 以下の範囲にあるこの国で知られている素数は 999983 と 999979 の 2 個しかありません。

Input

複数のデータセットの並びが入力として与えられます。入力の終わりはゼロひとつの行で示されます。
各データセットは以下の形式で与えられます。

n
p1 m1
p2 m2
:
pn mn

１行目に抽選結果の数 n (1 ≤ n ≤ 1000000)、続く n 行に i 番目の抽選結果の情報 pi,  mi が空白区切りで与えられます。

データセットの数は 20 を超えません。

Output

データセットごとに王様への請求額をプライム単位（整数）で１行に出力します。

Sample Input

4
5 0
9 1
3 10
11 3
1
999983 20
0

Output for the Sample Input

5
1



Sample input:
4
5 0
9 1
3 10
11 3
1
999983 20
0


Sample output:
5
1

Submission status: Time Limit Exceeded

Solution code:
n=1000000
p=[1]*(n+1)
p[0],p[1]=0,0
for i in range(2,int(n**0.5)+1):
    if p[i]:
        for j in range(i*i,n,i):
            p[j]=0
while 1:
    l=int(input())
    if l==0:break
    ans=-l
    while l:
        l-=1
        a,m=map(int,input().split())
        ans+=sum(p[max(1,a-m)-1:min(n,a+m)])
    print(ans)

Submission status: Time Limit Exceeded

Solution code:
n=1000000
p=[1]*(n+1)
pr=[]
p[0],p[1]=0,0
for i in range(2,int(n**0.5)+1):
    if p[i]:        
        for j in range(i*i,n,i):
            p[j]=0
for i in range(2,n):
    if p[i]:
        pr.append(int(i))
while 1:
    l=int(input())
    if l==0:break
    ans=-l
    while l:
        l-=1
        a,m=map(int,input().split())
        j=0
        while j<len(pr) and pr[j]<=a+m:
            if(pr[j]>=a-m):ans+=1
            j+=1
    print(ans)

Submission status: Time Limit Exceeded

Solution code:
import math
N=1000000
prime=[1]*(N+1)
prime[0]=0
prime[1]=0
for i in range(2,int(math.sqrt(N))):
    if prime[i]==1:
        j=i
        while i*j<=N:
            prime[i*j]=0
            j+=1

while True:
    n=input()
    if n==0:break
    cnt=0
    for i in range(n):
        p,m=map(int,raw_input().split())
        for j in range(max(p-m,0),min(p+m+1,N+1)):
            if prime[j]==1:
                cnt+=1
    print cnt-n

Submission status: Time Limit Exceeded

Solution code:
import math
N=1000000
prime=[1]*(N+1)
prime[0]=0
prime[1]=0
for i in range(2,int(math.sqrt(N))):
    if prime[i]==1:
        j=i
        while i*j<=N:
            prime[i*j]=0
            j+=1
prime=[i for i in range(N+1) if prime[i]==1]

while True:
    n=input()
    if n==0:break
    cnt=0
    for i in range(n):
        p,m=map(int,raw_input().split())
        for j in range(len(prime)):
            if p-m<=prime[j]<=p+m:
                cnt+=1
            if p+m<=prime[j]:
                break
    print cnt-n

Submission status: Accepted

Solution code:
import bisect as bs
def prime(m):
  N = range(1,m+1, 2)
  N[0] = 2
  for i in range(1, len(N)):
    x = N[i]
    if x: N[i+x::x] = [0]*len(N[i+x::x])
  return filter(None, N)

P = prime(1000000)
while 1:
  n = input()
  if n==0: break
  s = 0
  for _ in [0]*n:
    p, m = map(int, raw_input().split())
    s += bs.bisect_right(P, p+m) - bs.bisect_left(P, p-m)-1
  print s

Submission status: Accepted

Solution code:
import bisect
n=1000000
p=[1]*(n+1)
p[0],p[1]=0,0
for i in range(2,int(n**0.5)+1):
    if p[i]:        
        for j in range(i*i,n+1,i):
            p[j]=0
p=[i for i in range(n+1) if p[i]==1]
while 1:
    n=int(input())
    ans=-n
    if n==0:break
    while n:
        n-=1
        a,m=map(int,input().split())
        ans+=bisect.bisect_left(p,a+m+1)-bisect.bisect_left(p,a-m)
    print(ans)

Submission status: Accepted

Solution code:
import bisect

LENGTH = 1000000


def eratosthenes(length):
    from math import sqrt

    def is_prime_number(num):
        limit = int(sqrt(num)) + 1

        for lp in range(3, limit, 2):
            if num % lp == 0:
                return False

        return True

    is_prime_number_list = [True] * (LENGTH + 1)
    is_prime_number_list[0] = False
    is_prime_number_list[1] = False
    is_prime_number_list[2] = True

    for index in range(4, length + 1, 2):
        is_prime_number_list[index] = False

    limit = int(sqrt(length)) + 1

    for lp in range(3, limit, 2):
        if is_prime_number(lp):

            for index in range(lp * 2, length + 1, lp):
                is_prime_number_list[index] = False

    return is_prime_number_list


is_prime_number_list = eratosthenes(LENGTH)
prime_number_list = [index for index, item in enumerate(is_prime_number_list) if item]

while True:
    input_count = int(input())

    if input_count == 0:
        break

    pay = 0

    for _ in range(input_count):
        p, m = [int(item) for item in input().split(" ")]

        lower = p - m if 0 < p - m else 0
        upper = p + m if p + m < LENGTH else LENGTH

        lower_index = bisect.bisect_left(prime_number_list, lower)
        upper_index = bisect.bisect_right(prime_number_list, upper)

        prize = upper_index - lower_index
        pay += prize - 1

    print(pay)

Submission status: Accepted

Solution code:
# AOJ 0165 Lottery
# Python3 2018.6.23 bal4u

#define MAX  1000000
#define SQRT 1000     // sqrt(MAX)



MAX = 1000000
SQRT = 1000     # sqrt(MAX)
comp = [0]*(MAX+2)

def sieve():
	for i in range(3, SQRT, 2):
		if comp[i] == 0:
			for j in range(i*i, MAX, i): comp[j] = 1
sieve()
tbl = [0]*(MAX+2)
tbl[2] = k = 1
for i in range(3, MAX, 2):
	if comp[i] == 0: k += 1
	tbl[i+1] = tbl[i] = k

while 1:
	n = int(input())
	if n == 0: break
	ans = 0;
	for i in range(n):
		p, m = map(int, input().split())
		a, b = p-m, p+m
		if a < 2: a = 2
		if b > 1000000: b = 1000000
		ans += tbl[b]-tbl[a-1] - 1;
	if ans < 0: ans = 0
	print(ans)

Submission status: Accepted

Solution code:
import math
import bisect
N=1000000
prime=[1]*(N+1)
prime[0]=0
prime[1]=0
for i in range(2,int(math.sqrt(N))):
    if prime[i]==1:
        j=i
        while i*j<=N:
            prime[i*j]=0
            j+=1
prime=[i for i in range(N+1) if prime[i]==1]

while True:
    n=input()
    if n==0:break
    cnt=0
    for i in range(n):
        p,m=map(int,raw_input().split())
        l=bisect.bisect_left(prime,p-m)
        u=bisect.bisect_left(prime,p+m+1)
        cnt+=u-l
    print cnt-n

Submission status: Accepted

Solution code:
M = 1000001

p = [1]*M
p[0] = p[1] = 0
for i in range(2, int(M**.5)+1):
    if p[i]:
        for j in range(i*i, M, i):
            p[j] = 0
cs = [0]*M
cur = 0
for i in range(M):
    if p[i]:
        cur += 1
    cs[i] = cur

while 1:
    N = int(input())
    if N == 0:
        break
    ans = 0
    for i in range(N):
        p, m = map(int, input().split())
        x = cs[min(p+m, M-1)] - cs[max(p-m-1, 0)]
        if x:
            ans += x-1
        else:
            ans -= 1
    print(ans)

Submission status: Accepted

Solution code:
from itertools import accumulate
MAX = 1000000
is_prime = [1] * MAX
is_prime[0] = is_prime[1] = 0
for i in range(2, int(MAX ** (1 / 2)) + 1):
  if is_prime[i]:
    for j in range(i * i, MAX, i):
      is_prime[j] = 0
acc_prime = list(accumulate(is_prime))

while True:
  n = int(input())
  if n == 0:
    break
  ans = 0
  for _ in range(n):
    p, m = map(int, input().split())
    ans += acc_prime[min(p + m, MAX - 1)] - acc_prime[max(0, p - m - 1)] - 1
  print(ans)