#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);
	}
}

bool usedc[2000][2000] = {0};

void dfs_c(point & p, int &count, point &c)
{
	usedc[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] || usedc[t.x][t.y])
			continue;
		dfs_c(t, 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)
{
	usedc[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] || usedc[t.x][t.y])
			continue;
		dfs_s(t, c, yes, no, rad);
		break;
	}
}

void f(const point &p)
{
	memset(usedc, 0, sizeof usedc);
	int count = 0;
	point c(0, 0);
	dfs_c(point(p), count, c);
	c.x /= count;
	c.y /= count;
	memset(usedc, 0, sizeof usedc);
	int yes = 0,
		no = 0;
	int rad = abs(p.x - c.x);
	dfs_s(point(p), c, yes, no, rad);
	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};
bool useder[2000][2000] = {0};
bool usedch[2000][2000] = {0};
void dfs_er(point & p, int &count)
{
	useder[p.x][p.y] = 1;
	++count;
	if(count > 100) return;
	for(int i = 0; i < 4; ++i)
	{
		int tx = p.x + dx[i],
			ty = p.y + dy[i];
		if(tx < 0 || ty < 0 || tx >= n || ty >= n ||
			!a[tx][ty] || useder[tx][ty]) continue;
		dfs_er(point(tx, ty), count);
	}
}

void dfs_change(const point & p, int c)
{
	usedch[p.x][p.y] = 1;
	a[p.x][p.y] = c;
	for(int i = 0; i < 4; ++i)
	{
		int tx = p.x + dx[i],
			ty = p.y + dy[i];
		if(tx < 0 || ty < 0 || tx >= n || ty >= n ||
			a[tx][ty] != c || usedch[tx][ty]) continue;
		dfs_change(point(tx, ty), 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 i = 0; i < n; ++i)
	{
		for(int j = 0; j < n; ++j)
		{
			if(useder[i][j] || !a[i][j]) continue;
			int count = 0;
			point cur(i, j);
			dfs_er(cur, count);
			if(count > 30)
				continue;
			memset(usedch, 0, sizeof usedch);
			dfs_change(cur, 1);
			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;
}
### Prompt

Please provide a Cpp coded solution to the problem described below:

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 .

### Response

```cpp
#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);
	}
}

bool usedc[2000][2000] = {0};

void dfs_c(point & p, int &count, point &c)
{
	usedc[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] || usedc[t.x][t.y])
			continue;
		dfs_c(t, 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)
{
	usedc[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] || usedc[t.x][t.y])
			continue;
		dfs_s(t, c, yes, no, rad);
		break;
	}
}

void f(const point &p)
{
	memset(usedc, 0, sizeof usedc);
	int count = 0;
	point c(0, 0);
	dfs_c(point(p), count, c);
	c.x /= count;
	c.y /= count;
	memset(usedc, 0, sizeof usedc);
	int yes = 0,
		no = 0;
	int rad = abs(p.x - c.x);
	dfs_s(point(p), c, yes, no, rad);
	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};
bool useder[2000][2000] = {0};
bool usedch[2000][2000] = {0};
void dfs_er(point & p, int &count)
{
	useder[p.x][p.y] = 1;
	++count;
	if(count > 100) return;
	for(int i = 0; i < 4; ++i)
	{
		int tx = p.x + dx[i],
			ty = p.y + dy[i];
		if(tx < 0 || ty < 0 || tx >= n || ty >= n ||
			!a[tx][ty] || useder[tx][ty]) continue;
		dfs_er(point(tx, ty), count);
	}
}

void dfs_change(const point & p, int c)
{
	usedch[p.x][p.y] = 1;
	a[p.x][p.y] = c;
	for(int i = 0; i < 4; ++i)
	{
		int tx = p.x + dx[i],
			ty = p.y + dy[i];
		if(tx < 0 || ty < 0 || tx >= n || ty >= n ||
			a[tx][ty] != c || usedch[tx][ty]) continue;
		dfs_change(point(tx, ty), 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 i = 0; i < n; ++i)
	{
		for(int j = 0; j < n; ++j)
		{
			if(useder[i][j] || !a[i][j]) continue;
			int count = 0;
			point cur(i, j);
			dfs_er(cur, count);
			if(count > 30)
				continue;
			memset(usedch, 0, sizeof usedch);
			dfs_change(cur, 1);
			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;
}
```