Problem: Offering the ABBYY Cup participants a problem written by the Smart Beaver is becoming a tradition. He proposed the following problem.

You are given a monochrome image, that is, an image that is composed of two colors (black and white). The image is given in raster form, that is, as a matrix of pixels' colors, and the matrix's size coincides with the size of the image.

The white color on the given image corresponds to the background. Also, the image contains several black geometric shapes. It is known that the image can contain only two types of shapes: squares and circles. Your task is to count the number of circles and the number of squares which the given image contains.

The squares on the image can be rotated arbitrarily. In addition, the image can possibly contain some noise arranged as follows: each pixel of the original image can change its color to the opposite with the probability of 20%.

<image> An example of an image that has no noise and the sides of the squares are parallel to the coordinate axes (two circles and three squares).  <image> An example of an image that has no noise and the squares are rotated arbitrarily (two circles and three squares).  <image> An example of an image that has noise and the squares are rotated arbitrarily (one circle and three squares). 

Input

The first input line contains a single integer n (1000 ≤ n ≤ 2000), which is the length and the width of the original image. 

Next n lines describe the matrix of colors of the image pixels. The i-th line contains exactly n integers aij (0 ≤ aij ≤ 1), separated by spaces. Value of aij = 0 corresponds to a white pixel and aij = 1 corresponds to a black one. 

It is guaranteed that the lengths of the sides of the squares and the diameters of the circles in the image are at least 15 pixels, and the distance between any two figures is at least 10 pixels. It is also guaranteed that a human can easily calculate the number of circles and squares in the original image. The total number of figures in the image doesn't exceed 50.

The input limitations for getting 20 points are: 

  * These test cases have no noise and the sides of the squares are parallel to the coordinate axes. 



The input limitations for getting 50 points are: 

  * These test cases have no noise, but the squares are rotated arbitrarily. 



The input limitations for getting 100 points are: 

  * These test cases have noise and the squares are rotated arbitrarily. 

Output

Print exactly two integers, separated by a single space — the number of circles and the number of squares in the given image, correspondingly.

Examples

Note

You are given a sample of original data for each difficulty level. The samples are available at http://codeforces.ru/static/materials/contests/178/e-samples.zip .

Solution (in CPP):

Solution:
#define _USE_MATH_DEFINES
#define _CRT_SECURE_NO_DEPRECATE
#include <cmath>
#include <cstdio>
#include <string>
#include <algorithm>
#include <map>
#include <queue>
#include <vector>
#include <set>
#include <iomanip>
#include <ctime>
#include <iostream>
#include <sstream>
#pragma comment(linker, "/STACK:256000000")

using namespace std;

typedef long long int64;
typedef unsigned long long uint64;
template<class T> inline T sqr(T a) {return a * a;}
#define INF 123456789
#define MOD 1000000007
#define PRIME 1103
#define TEST "test"
#define eps 10e-1

int a[2000][2000];
bool was[2000][2000];
struct point 
{
	int x, y;
	point(){}
	point(int _x, int _y)
	{
		x = _x;
		y = _y;
	}
};
int r = 0, q = 0;
int n;
int dx[4] = {0, 0, 1, -1},
	dy[4] = {-1, 1, 0, 0};


void dfs(const point &p)
{
	was[p.x][p.y] = 1;
	for(int i = 0; i < 4; ++i)
	{
		point t(p.x + dx[i], p.y + dy[i]);
		if(t.x < 0 || t.x >= n || 
			t.y < 0 || t.y >= n || was[t.x][t.y] ||
			!a[t.x][t.y]) continue;
		dfs(t);
	}
}


void dfs_c(point & p, bool used[2000][2000], int &count, point &c)
{
	was[p.x][p.y] = 1;
	++count;
	c.x += p.x;
	c.y += p.y;
	for(int i = 0; i < 4; ++i)
	{
		point t (p.x + dx[i], p.y + dy[i]);
		if(t.x < 0 || t.y < 0 ||
			t.x >= n || t.y >= n || 
			!a[t.x][t.y] || was[t.x][t.y])
			continue;
		dfs_c(t, used, count, c);
	}
}

int dx1[] = {0, 1, 1, 1, 0},
	dy1[] = {1, 1, 0, -1, -1};
void dfs_s(point & p, point & c, int &yes, int &no,
	double rad, bool used[2000][2000])
{
	used[p.x][p.y] = 1;
	double rad1 = sqrt(sqr(p.x - c.x) + sqr(p.y - c.y) + .0);
	if(abs(rad - rad1) < 1.)
		++yes;
	else ++no;
	for(int i = 0; i < 5; ++i)
	{
		point t (p.x + dx1[i], p.y + dy1[i]);
		if(t.x < 0 || t.y < 0 ||
			t.x >= n || t.y >= n ||
			!a[t.x][t.y] || used[t.x][t.y])
			continue;
		dfs_s(t, c, yes, no, rad, used);
		break;
	}
}

void f(const point &p)
{
	bool used[2000][2000] = {0};
	int count = 0;
	point c(0, 0);
	dfs_c(point(p), used, count, c);
	c.x /= count;
	c.y /= count;
	memset(used, 0, sizeof used);
	int yes = 0,
		no = 0;
	int rad = abs(p.x - c.x);
	dfs_s(point(p), c, yes, no, rad, used);
	if(!no)
		++r;
	else ++q;
	dfs(p);
}

int dx2[] = {-1, 0, 1, 1, 1, 0, -1, -1},
	dy2[] = {0, 1, 1, 1, 0, -1, -1, -1};
/*
void change(int x, int y, int c)
{
	int oth = 0;
	for(int i = 0; i < 8; ++i)
	{
		int tx = x + dx2[i],
			ty = y + dy2[i];
		if(tx < 0 || ty < 0 || tx >= n || ty >= n ||
			a[tx][ty] != c) continue;
		++oth;
	}
	if(oth > 4) a[x][y] = c;
}
void calc()
{
	for(int i = 0; i < n; ++i)
	{
		a[0][i] = 0;
		a[i][0] = 0;
		a[n - 1][i] = 0;
		a[i][n - 1] = 0;
	}
	for(int it = 0; it < 5; ++it)
	{
		for(int i = 0; i < n; ++i)
		{
			for(int j = 0; j < n; ++j)
			{
				if(!a[i][j]) continue;
				change(i, j, 0);
			}
		}
	}
	for(int it = 0; it < 10; ++it)
	{
		for(int i = 0; i < n; ++i)
		{
			for(int j = 0; j < n; ++j)
			{
				if(a[i][j]) continue;
				change(i, j, 1);
			}
		}
	}
}*/

void dfs_er(point & p, bool used[2000][2000], int count)
{
	used[p.x][p.y] = 1;
	++count;
	for(int i = 0; i < n; ++i)
	{
		int tx = p.x + dx2[i],
			ty = p.y+ dy2[i];
		if(tx < 0 || ty < 0 || tx >= n || ty >= n ||
			!a[tx][ty] || used[tx][ty]) continue;
		dfs_er(point(tx, ty), used, count);
	}
}

void dfs_change(point & p, bool used[2000][2000], int c)
{
	used[p.x][p.y] = 1;
	a[p.x][p.y] = c;
	for(int i = 0; i < n; ++i)
	{
		int tx = p.x + dx2[i],
			ty = p.y + dy2[i];
		if(tx < 0 || ty < 0 || tx >= n || ty >= n ||
			a[tx][ty] != c || used[tx][ty]) continue;
		dfs_change(point(tx, ty), used, c);
	}
}

void calc()
{
	for(int i = 0; i < n; ++i)
	{
		a[0][i] = 0;
		a[i][0] = 0;
		a[n - 1][i] = 0;
		a[i][n - 1] = 0;
	}
	bool used[2000][2000] = {0};
	for(int i = 0; i < n; ++i)
	{
		for(int j = 0; j < n; ++j)
		{
			if(used[i][j] || !a[i][j]) continue;
			int count = 0;
			dfs_er(point(i, j), used, count);
			if(count < 30)
			{
				bool used1[2000][2000] = {0};
				dfs_change(point(i, j), used1, 0);
				continue;
			}
		}
	}


}
int main()
{
#ifndef ONLINE_JUDGE
	freopen(TEST ".in", "r", stdin);
	freopen(TEST ".out", "w", stdout);
#endif
	scanf("%d", &n);
	for(int i = 0; i < n; ++i)
		for(int j = 0; j < n; ++j)
			scanf("%d", &a[i][j]);
	//calc();
	/*for(int i = 0; i < n; ++i)
	{
		for(int j = 0; j < n; ++j)
			cout << a[i][j];
		cout << endl;
	}*/
	for(int i = 0; i < n; ++i)
	{
		for(int j = 0; j < n; ++j)
		{
			if(!a[i][j] || was[i][j]) continue;
			f(point(i, j));
		}
	}
	cout << r << " "<< q;
	return 0;
}