Problem:
Score : 400 points

Problem Statement

We say that a odd number N is similar to 2017 when both N and (N+1)/2 are prime.

You are given Q queries.

In the i-th query, given two odd numbers l_i and r_i, find the number of odd numbers x similar to 2017 such that l_i ≤ x ≤ r_i.

Constraints

1≤Q≤10^5

1≤l_i≤r_i≤10^5

l_i and r_i are odd.

All input values are integers.

Input

Input is given from Standard Input in the following format:

Q
l_1 r_1
:
l_Q r_Q

Output

Print Q lines. The i-th line (1≤i≤Q) should contain the response to the i-th query.

Sample Input 1

1
3 7

Sample Output 1

2

3 is similar to 2017, since both 3 and (3+1)/2=2 are prime.

5 is similar to 2017, since both 5 and (5+1)/2=3 are prime.

7 is not similar to 2017, since (7+1)/2=4 is not prime, although 7 is prime.

Thus, the response to the first query should be 2.

Sample Input 2

4
13 13
7 11
7 11
2017 2017

Sample Output 2

1
0
0
1

Note that 2017 is also similar to 2017.

Sample Input 3

6
1 53
13 91
37 55
19 51
73 91
13 49

Sample Output 3

4
4
1
1
1
2

Solution in C (GCC 5.4.1) (judge verdict: Time Limit Exceeded)

#include <limits.h>
#include <math.h>
#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <time.h>

// 内部定数
#define D_VAL_MAX		100000									// 最大値

// 内部変数
static int si1Prime[D_VAL_MAX];									// 素数

// 内部変数 - テスト用
#ifdef D_TEST
	static int siRes;
	static FILE *szpFpT, *szpFpA;
#endif

// ルート
int
fSqrt(
	int piVal					// <I> 値 1～
)
{
	int liVal;

	// 正方形(辺の２乗)を１ずつ削る
	for (liVal = 1; piVal > 0; liVal++) {
		piVal -= liVal * 2 + 1;
	}

	return liVal - 1;
}

// 素数チェック
// 戻り値：[1]素数 [-1]素数ではない
int
fIsPrime(
	int piVal					// <I> 値
)
{
	int i;

	// 対象外チェック
	if (piVal < 2) {
		return -1;
	}

	// チェック済判定
	if (si1Prime[piVal] != 0) {
		return si1Prime[piVal];
	}

	// チェック上限
	int liLimit = fSqrt(piVal);

	// チェック
	for (i = 2; i <= liLimit; i++) {
		if (piVal % i == 0) {
			si1Prime[piVal] = -1;
			return -1;
		}
	}

	si1Prime[piVal] = 1;
	return 1;
}

// 実行メイン
int
fMain(
	int piTNo					// <I> テスト番号 1～
)
{
	int i, j, liRet;
	char lc1Buf[1024], lc1Out[1024];

	// テストファイルオープン
#ifdef D_TEST
	sprintf(lc1Buf, ".\\Test\\T%d.txt", piTNo);
	szpFpT = fopen(lc1Buf, "r");
	sprintf(lc1Buf, ".\\Test\\A%d.txt", piTNo);
	szpFpA = fopen(lc1Buf, "r");
	siRes = 0;
#endif

	// クエリ数取得
	int liQCnt;
#ifdef D_TEST
	fgets(lc1Buf, sizeof(lc1Buf), szpFpT);
#else
	fgets(lc1Buf, sizeof(lc1Buf), stdin);
#endif
	sscanf(lc1Buf, "%d", &liQCnt);

	// クエリ数でループ
	for (i = 0; i < liQCnt; i++) {

		// 範囲取得
		int liSNo, liENo;
#ifdef D_TEST
		fgets(lc1Buf, sizeof(lc1Buf), szpFpT);
#else
		fgets(lc1Buf, sizeof(lc1Buf), stdin);
#endif
		sscanf(lc1Buf, "%d%d", &liSNo, &liENo);

		// 対象の個数取得
		int liCnt = 0;
		for (j = liSNo; j <= liENo; j += 2) {

			// 素数チェック
			liRet = fIsPrime(j);
			if (liRet != 1) {
				continue;
			}

			// 素数チェック
			liRet = fIsPrime((j + 1) / 2);
			if (liRet != 1) {
				continue;
			}

			// 対象
			liCnt++;
		}

		// 結果セット
		sprintf(lc1Out, "%d\n", liCnt);

		// 結果表示
#ifdef D_TEST
		fgets(lc1Buf, sizeof(lc1Buf), szpFpA);
		if (strcmp(lc1Buf, lc1Out)) {
			siRes = -1;
		}
#else
		printf("%s", lc1Out);
#endif
	}

	// テストファイルクローズ
#ifdef D_TEST
	fclose(szpFpT);
	fclose(szpFpA);
#endif

	// テスト結果
#ifdef D_TEST
	if (siRes == 0) {
		printf("OK %d\n", piTNo);
	}
	else {
		printf("NG %d\n", piTNo);
	}
#endif

	return 0;
}

int
main()
{

#ifdef D_TEST
	int i;
	for (i = D_TEST_SNO; i <= D_TEST_ENO; i++) {
		fMain(i);
	}
#else
	fMain(0);
#endif

	return 0;
}