#include<bits/stdc++.h>
using namespace std;
int n,a[1001][1001];
namespace square_sum {
	int b[1001][1001];                     // Prefix sum 2D
	// b_{x,y}=\sum\limits_{i=1}^x\sum\limits_{j=1}^ya_{i,j}
	void Set() {                           // initialization
		for(int i=1;i<=n;i++)
			for(int j=1;j<=n;j++)
				b[i][j]=a[i][j]+b[i-1][j]+b[i][j-1]-b[i-1][j-1];
	}
	int Get(int x1,int y1,int x2,int y2) { // Get \sum\limits_{i=x_1}^{x_2}\sum\limits_{j=y_1}^{y_2}a_{i,j}
		x1=max(0,x1);x2=min(n,x2);
		y1=max(0,y1);y2=min(n,y2);
		if(y1>y2||x1>x2) return 0;
		return b[x2][y2]-b[x2][y1-1]-b[x1-1][y2]+b[x1-1][y1-1];
	}
}
namespace triangle_sum {
	int b[2001][1001];
	// b_{x,y}=\sum\limits_{i=1}^y\sum\limits_{j=\max(1,x-i+1)}^na_{i,j}
	void Set() {
		for(int i=1;i<=n*2;i++)
			for(int j=1;j<=n;j++)
				b[i][j]=b[i][j-1]+square_sum::Get(j,max(1,i-j+1),j,n);
	}
	int Get(int x,int y,int k) {           // Get \sum\limits_{i=1}^y\sum\limits_{j=\max(1,x-i+1)}^na_{i,j}
		return b[max(1,x+y-k-1)][y]-(y<=k+1?0:b[x+y-k-1][y-k-1])-square_sum::Get(x+1,y-k,n,y);
	}
}
int main() {
	cin >> n;
	for(int i=1;i<=n;i++)
		for(int j=1;j<=n;j++)
			cin >> a[i][j];
	square_sum::Set();
	triangle_sum::Set();
}