Problem: Prime Caves

An international expedition discovered abandoned Buddhist cave temples in a giant cliff standing on the middle of a desert. There were many small caves dug into halfway down the vertical cliff, which were aligned on square grids. The archaeologists in the expedition were excited by Buddha's statues in those caves. More excitingly, there were scrolls of Buddhism sutras (holy books) hidden in some of the caves. Those scrolls, which estimated to be written more than thousand years ago, were of immeasurable value.

The leader of the expedition wanted to collect as many scrolls as possible. However, it was not easy to get into the caves as they were located in the middle of the cliff. The only way to enter a cave was to hang you from a helicopter. Once entered and explored a cave, you can climb down to one of the three caves under your cave; i.e., either the cave directly below your cave, the caves one left or right to the cave directly below your cave. This can be repeated for as many times as you want, and then you can descent down to the ground by using a long rope.

So you can explore several caves in one attempt. But which caves you should visit? By examining the results of the preliminary attempts, a mathematician member of the expedition discovered that (1) the caves can be numbered from the central one, spiraling out as shown in Figure D-1; and (2) only those caves with their cave numbers being prime (let's call such caves prime caves), which are circled in the figure, store scrolls. From the next attempt, you would be able to maximize the number of prime caves to explore.

<image>

Figure D-1: Numbering of the caves and the prime caves


Write a program, given the total number of the caves and the cave visited first, that finds the descending route containing the maximum number of prime caves.

Input

The input consists of multiple datasets. Each dataset has an integer m (1 ≤ m ≤ 106) and an integer n (1 ≤ n ≤ m) in one line, separated by a space. m represents the total number of caves. n represents the cave number of the cave from which you will start your exploration. The last dataset is followed by a line with two zeros.

Output

For each dataset, find the path that starts from the cave n and contains the largest number of prime caves, and output the number of the prime caves on the path and the last prime cave number on the path in one line, separated by a space. The cave n, explored first, is included in the caves on the path. If there is more than one such path, output for the path in which the last of the prime caves explored has the largest number among such paths. When there is no such path that explores prime caves, output "`0 0`" (without quotation marks).

Sample Input


49 22
46 37
42 23
945 561
1081 681
1056 452
1042 862
973 677
1000000 1000000
0 0


Output for the Sample Input


0 0
6 43
1 23
20 829
18 947
10 947
13 947
23 947
534 993541






Example

Input

49 22
46 37
42 23
945 561
1081 681
1056 452
1042 862
973 677
1000000 1000000
0 0


Output

0 0
6 43
1 23
20 829
18 947
10 947
13 947
23 947
534 993541

Solution (in CPP):

Solution:
#include <bits/stdc++.h>

using namespace std;

#ifdef DEBUG_MODE
	#define DBG(n) n;
#else
	#define DBG(n) ;
#endif
#define REP(i,n) for(ll (i) = (0);(i) < (n);++i)
#define rep(i,s,g) for(ll (i) = (s);(i) < (g);++i)
#define rrep(i,s,g) for(ll (i) = (s);i >= (g);--(i))
#define PB push_back
#define MP make_pair
#define FI first
#define SE second
#define SHOW1d(v,n) {for(int W = 0;W < (n);W++)cerr << v[W] << ' ';cerr << endl << endl;}
#define SHOW2d(v,i,j) {for(int aaa = 0;aaa < i;aaa++){for(int bbb = 0;bbb < j;bbb++)cerr << v[aaa][bbb] << ' ';cerr << endl;}cerr << endl;}
#define ALL(v) v.begin(),v.end()
#define Decimal fixed<<setprecision(20)
#define INF 1000000000
#define MOD 1000000007
#define M_MAX 1111

typedef long long ll;
typedef pair<ll,ll> P;

int xx[M_MAX*M_MAX];
int yy[M_MAX*M_MAX];
pair<int,int> dp[M_MAX][M_MAX];
bool isp[M_MAX*M_MAX];
int mp[M_MAX][M_MAX];
int dx[] = {-1,0,1,0};
int dy[] = {0,-1,0,1};

void f(){
	int x = M_MAX-1;
	int y = M_MAX-1;
	int now = M_MAX*M_MAX;
	int cou = 0;
	
	while(now > 0){
		mp[y][x] = now;
		xx[now] = x;
		yy[now] = y;
		now--;
		int tmpx = x + dx[cou];
		int tmpy = y + dy[cou];
		
		if(tmpx >= 0 && tmpx < M_MAX && tmpy >= 0 && tmpy < M_MAX && 
		mp[tmpy][tmpx] == 0){
			x = tmpx;
			y = tmpy;
		}
		else{
			cou = (cou + 1) % 4;
			x += dx[cou];
			y += dy[cou];
		}
	}
	
	REP(i,M_MAX*M_MAX)isp[i] = true;
	isp[0] = isp[1] = false;
	for(int i = 2;i < M_MAX*M_MAX;i++){
		if(isp[i]){
			for(int j = 2*i;j < M_MAX*M_MAX;j += i)isp[j] = false;
		}
	}
}

int main(){
	
	f();
	
	int n,m;
	
	while(cin >> m >> n,n|m){
		
		REP(i,M_MAX)REP(j,M_MAX)dp[i][j] = MP(-1,-1);
		
		dp[yy[n]][xx[n]] = MP(0,0);
		
		REP(i,M_MAX){
			REP(j,M_MAX){
				if(dp[i][j] != MP(-1,-1) && mp[i][j] <= m){
					if(isp[mp[i][j]]){
						dp[i][j] = MP(dp[i][j].FI+1,mp[i][j]);
					}
					REP(k,3){
						int x = j + dx[k];
						int y = i + 1;
						DBG(cout << "xy : " << x << ' ' << y << endl;)
						if(x >= 0 && x < M_MAX && y >= 0 && y < M_MAX){
							dp[y][x] = max(dp[i][j],dp[y][x]);
						}
					}
				}
			}
		}
		
		pair<int,int> ans = MP(0,0);
		
		REP(i,M_MAX){
			REP(j,M_MAX){
				if(mp[i][j] <= m)ans = max(ans,dp[i][j]);
			}
		}
		
		cout << ans.FI << ' ' << ans.SE << endl;
	}

	return 0;
}