Problem: F - Point Distance

Problem Statement

You work for invention center as a part time programmer. This center researches movement of protein molecules. It needs how molecules make clusters, so it will calculate distance of all pair molecules and will make a histogram. Molecules’ positions are given by a N x N grid map. Value C_{xy} of cell (x, y) in a grid map means that C_{xy} molecules are in the position (x, y).

You are given a grid map, please calculate histogram of all pair molecules.

Input

Input is formatted as follows.


N
C_{11} C_{12} ... C_{1N}
C_{21} C_{22} ... C_{2N}
...
C_{N1} C_{N2} ... C_{NN}


First line contains the grid map size N (1 \leq N \leq 1024). Each next N line contains N numbers. Each numbers C_{xy} (0 \leq C_{xy} \leq 9) means the number of molecule in position (x,y). There are at least 2 molecules.

Output

Output is formatted as follows.


D_{ave}
d_{1} c_{1}
d_{2} c_{2}
...
d_{m} c_{m}


Print D_{ave} which an average distance of all pair molecules on first line. Next, print histogram of all pair distance. Each line contains a distance d_i and the number of pair molecules c_i(0 \lt c_i). The distance d_i should be calculated by squared Euclidean distance and sorted as increasing order. If the number of different distance is more than 10,000, please show the first 10,000 lines. The answer may be printed with an arbitrary number of decimal digits, but may not contain an absolute or relative error greater than or equal to 10^{-8}.

Sample Input 1


2
1 0
0 1


Output for the Sample Input 1


1.4142135624
2 1


Sample Input 2


3
1 1 1
1 1 1
1 1 1


Output for the Sample Input 2


1.6349751825
1 12
2 8
4 6
5 8
8 2


Sample Input 3


5
0 1 2 3 4
5 6 7 8 9
1 2 3 2 1
0 0 0 0 0
0 0 0 0 1


Output for the Sample Input 3


1.8589994382
0 125
1 379
2 232
4 186
5 200
8 27
9 111
10 98
13 21
16 50
17 37
18 6
20 7
25 6






Example

Input

2
1 0
0 1


Output

1.4142135624
2 1

Solution (in JAVA):

Solution:
import java.util.Arrays;
import java.util.Comparator;
import java.util.HashMap;
import java.util.Map.Entry;
import java.util.Scanner;

public class Main {
	public static void main(String[] args) {
		new Main().run();
	}

	void run() {
		Scanner sc = new Scanner(System.in);
		int n = sc.nextInt();
		long[][] a = new long[n][n];
		long[][] b = new long[n][n];
		int[][] C = new int[n][n];
		for (int i = 0; i < n; ++i) {
			for (int j = 0; j < n; ++j) {
				int v = sc.nextInt();
				C[i][j] = v;
				a[i][j] = v;
				b[n - 1 - i][n - 1 - j] = v;
			}
		}
		a = mul2D(a, b);
		long[][] cnt = new long[n][n];
		for (int i = 0; i < 2 * n; ++i) {
			for (int j = 0; j < 2 * n; ++j) {
				if (a[i][j] == 0)
					continue;
				int di = Math.abs(n - 1 - i);
				int dj = Math.abs(n - 1 - j);
				cnt[di][dj] += (int) a[i][j];
			}
		}
		double ave = 0;
		HashMap<Integer, Long> map = new HashMap<>();
		for (int i = 0; i < n; ++i) {
			for (int j = 0; j < n; ++j) {
				if (cnt[i][j] == 0 || (i | j) == 0)
					continue;
				ave += Math.sqrt(i * i + j * j) * cnt[i][j];
				if (map.containsKey(i * i + j * j)) {
					map.put(i * i + j * j, map.get(i * i + j * j) + cnt[i][j]);
				} else {
					map.put(i * i + j * j, cnt[i][j]);
				}
			}
		}
		{
			long s = 0;
			for (int i = 0; i < n; ++i) {
				for (int j = 0; j < n; ++j) {
					if (C[i][j] > 1)
						s += C[i][j] * (C[i][j] - 1);
				}
			}
			if (s > 0)
				map.put(0, s);
		}

		long sz = 0;
		for (int i = 0; i < n; ++i) {
			for (int j = 0; j < n; ++j) {
				sz += C[i][j];
			}
		}
		System.out.println(ave / (sz * (sz - 1)));
		long[][] ans = new long[map.size()][2];
		{
			int p = 0;
			for (Entry<Integer, Long> es : map.entrySet()) {
				ans[p][0] = es.getKey();
				ans[p][1] = es.getValue() / 2;
				++p;
			}
		}
		Arrays.sort(ans, new Comparator<long[]>() {
			@Override
			public int compare(long[] o1, long[] o2) {
				return Long.compare(o1[0], o2[0]);
			}
		});
		for (int i = 0; i < Math.min(10000, ans.length); ++i) {
			System.out.println(ans[i][0] + " " + ans[i][1]);
		}
	}

	final long MODULO = 1012924417L;
	final long root = 5;

	long[][] mul2D(long[][] a, long[][] b) {
		int n1 = Integer.highestOneBit(a.length + b.length) << 1;
		int n2 = Integer.highestOneBit(a[0].length + b[0].length) << 1;
		long[][] a2 = new long[n1][n2];
		long[][] b2 = new long[n1][n2];
		for (int i = 0; i < a.length; ++i)
			for (int j = 0; j < a[i].length; ++j)
				a2[i][j] = a[i][j];
		for (int i = 0; i < b.length; ++i)
			for (int j = 0; j < b[i].length; ++j)
				b2[i][j] = b[i][j];
		a2 = fft2D(a2, false);
		b2 = fft2D(b2, false);
		for (int i = 0; i < a2.length; ++i)
			for (int j = 0; j < a2[i].length; ++j)
				a2[i][j] = a2[i][j] * b2[i][j] % MODULO;
		a2 = fft2D(a2, true);
		long invn1 = pow(n1, MODULO - 2);
		long invn2 = pow(n2, MODULO - 2);
		for (int i = 0; i < a2.length; ++i)
			for (int j = 0; j < a2[i].length; ++j)
				a2[i][j] = a2[i][j] * invn1 % MODULO * invn2 % MODULO;
		return a2;
	}

	long[][] fft2D(long[][] a, boolean inv) {
		int n = a.length;
		for (int i = 0; i < n; ++i) {
			a[i] = fft(a[i], inv);
		}
		for (int i = 0; i < n; ++i) {
			for (int j = i; j < n; ++j) {
				long d = a[i][j];
				a[i][j] = a[j][i];
				a[j][i] = d;
			}
		}
		for (int i = 0; i < n; ++i) {
			a[i] = fft(a[i], inv);
		}
		for (int i = 0; i < n; ++i) {
			for (int j = i; j < n; ++j) {
				long d = a[i][j];
				a[i][j] = a[j][i];
				a[j][i] = d;
			}
		}
		return a;
	}

	long[] mul(long[] a, long[] b) {
		int n = Integer.highestOneBit(a.length + b.length) << 1;
		a = Arrays.copyOf(a, n);
		b = Arrays.copyOf(b, n);
		a = fft(a, false);
		b = fft(b, false);
		long ninv = pow(n, MODULO - 2);
		for (int i = 0; i < n; ++i) {
			a[i] = a[i] * b[i] % MODULO;
		}
		a = fft(a, true);
		for (int i = 0; i < n; ++i)
			a[i] = a[i] * ninv % MODULO;
		return a;
	}

	long[] fft(long[] a, boolean inv) {
		int n = a.length;
		int c = 0;
		for (int i = 1; i < n; ++i) {
			for (int j = n >> 1; j > (c ^= j); j >>= 1)
				;
			if (c > i) {
				long d = a[c];
				a[c] = a[i];
				a[i] = d;
			}
		}

		for (int i = 1; i < n; i <<= 1) {
			long w = pow(root, (MODULO - 1) / (2 * i));
			if (inv)
				w = pow(w, MODULO - 2);
			for (int j = 0; j < n; j += 2 * i) {
				long wn = 1;
				for (int k = 0; k < i; ++k) {
					long u = a[k + j];
					long v = a[k + j + i] * wn % MODULO;
					a[k + j] = (u + v) % MODULO;
					a[k + j + i] = (u - v + MODULO) % MODULO;
					wn = wn * w % MODULO;
				}
			}
		}
		return a;
	}

	long pow(long a, long n) {
		long ret = 1;
		for (; n > 0; n >>= 1, a = a * a % MODULO) {
			if (n % 2 == 1)
				ret = ret * a % MODULO;
		}
		return ret;
	}

	static void tr(Object... objects) {
		System.out.println(Arrays.deepToString(objects));
	}

}