#include using namespace std; #define DEBUG(x) cerr<<#x<<": "< #define vl vector #define vii vector< vector > #define vll vector< vector > #define vs vector #define pii pair #define pis pair #define psi pair #define pll pair template pair operator+(const pair &s, const pair &t) { return pair(s.first + t.first, s.second + t.second); } template pair operator-(const pair &s, const pair &t) { return pair(s.first - t.first, s.second - t.second); } template ostream& operator<<(ostream& os, pair p) { os << "(" << p.first << ", " << p.second << ")"; return os; } #define X first #define Y second #define rep(i,n) for(int i=0;i<(n);i++) #define rep1(i,n) for(int i=1;i<=(n);i++) #define rrep(i,n) for(int i=(n)-1;i>=0;i--) #define rrep1(i,n) for(int i=(n);i>0;i--) #define REP(i,a,b) for(int i=a;i bool chmax(T &a, const T &b) { if (a bool chmin(T &a, const T &b) { if (a>b) { a = b; return 1; } return 0; } #define UNIQUE(v) v.erase(std::unique(v.begin(), v.end()), v.end()); const ll inf = 1000000001; const ll INF = 1e18; const ll MOD = 1000000007; //const ll MOD = 998244353; const double pi = 3.14159265358979323846; #define Sp(p) cout<> n >> m; const ll ADD = 100000000000000000; vll a(m, vl(m)); rep (i, m) { rep (j, m) { cin >> a[i][j]; } } vl x(n), y(n); rep (i, n) { cin >> y[i] >> x[i]; y[i]--; x[i]--; } rep (q, n) { a[y[q]][x[q]] += ADD; vll sum(m + 1, vl (m + 1)); rep (i, m) { rep (j, m) { sum[i + 1][j + 1] = sum[i + 1][j] + a[i][j]; } } rep (i, m) { rep (j, m + 1) { sum[i + 1][j] += sum[i][j]; } } vector diff(m, vll(m + 1)); rep (i0, m) { for (int i1 = i0 + 1; i1 <= m; i1++) { rep (j, m + 1) { diff[i0][i1].push_back(sum[i1][j] - sum[i0][j]); } sort(all(diff[i0][i1])); } } int ans = 0; rep (i0, m) { for (int i1 = i0 + 1; i1 <= m; i1++) { rep (j, m + 1) { ll now = -sum[i1][j] + sum[i0][j]; if (now >= ADD / 2 || now <= -ADD / 2) { break; } ll nokori = ADD - now; ans += upper_bound(all(diff[i0][i1]), nokori) - lower_bound(all(diff[i0][i1]), nokori); } } } cout << ans << endl; a[y[q]][x[q]] -= ADD; } }