#pragma GCC optimize("-O2")
#pragma GCC optimize("-Ofast")

#undef _GLIBCXX_DEBUG
 
#include <bits/stdc++.h>
 
using namespace std;
 
const int N = 1505;
 
short L[2][N][N];
short R[2][N][N];
int a[N][N];
int row[N][N];
int up[N][N];
string board[N];
 
int main() {
  ios::sync_with_stdio(false);
  cin.tie(0);
	cout.tie(0);
  int n;
  cin >> n;
  vector<string> board(n);
  for (int i = 0; i < n; i++) {
    cin >> board[i];
  }
  for (int i = 0; i < n; i++) {
    for (int j = 0; j < n; j++) {
      a[i][j] = (board[i][j] == '#' ? 0 : (int) (board[i][j] - '0'));
    }
  }
  for (int i = 0; i < n; i++) {
    row[i][0] = 0;
    for (int j = 0; j < n; j++) {
      row[i][j + 1] = row[i][j] + a[i][j];
    }
  }
  long long ans = 0;
  function<void(int,int)> put = [&](int i, int j) {
    if (0 <= i && i < n && 0 <= j && j < n && a[i][j] > 0) {
      if (i == 0 || a[i - 1][j] == 0) {
        if (j == 0 || a[i][j - 1] == 0) {
          for (int k = j; k < n; k++) {
            row[i][k + 1] -= a[i][j];
          }
          a[i][j] = 0;
          put(i + 1, j);
          put(i, j + 1);
        }
      }
    }
  };
  for (int i = n - 1; i >= 0; i--) {
    for (int j = n - 1; j >= 0; j--) {
      if (a[i][j] == 0) {
        continue;
      }
      int sum = 0;
      if (i < n - 1 && a[i + 1][j] != 0 && j < n - 1 && a[i][j + 1] != 0) {
        L[i & 1][j][i] = j;
        R[i & 1][j][i] = R[i & 1][j + 1][i];
        sum += row[i][R[i & 1][j][i] + 1] - row[i][L[i & 1][j][i]];
        up[i][j] = max(up[i + 1][j], up[i][j + 1]);
        int to = min(up[i + 1][j], up[i][j + 1]);
        memcpy(L[i & 1][j] + (i + 1), L[(i + 1) & 1][j] + (i + 1), sizeof(short) * (to - (i + 1) + 1));
        memcpy(R[i & 1][j] + (i + 1), R[i & 1][j + 1] + (i + 1), sizeof(short) * (to - (i + 1) + 1));
        for (int k = i + 1; k <= to; k++) {
          sum += row[k][R[i & 1][j][k] + 1] - row[k][L[i & 1][j][k]];
        }
        if (up[i + 1][j] < up[i][j + 1]) {
          memcpy(L[i & 1][j] + (to + 1), L[i & 1][j + 1] + (to + 1), sizeof(short) * (up[i][j] - (to + 1) + 1));
          memcpy(R[i & 1][j] + (to + 1), R[i & 1][j + 1] + (to + 1), sizeof(short) * (up[i][j] - (to + 1) + 1));
          for (int k = to + 1; k <= up[i][j]; k++) {
            sum += row[k][R[i & 1][j][k] + 1] - row[k][L[i & 1][j][k]];
          }
        }
        if (up[i + 1][j] > up[i][j + 1]) {
          memcpy(L[i & 1][j] + (to + 1), L[(i + 1) & 1][j] + (to + 1), sizeof(short) * (up[i][j] - (to + 1) + 1));
          memcpy(R[i & 1][j] + (to + 1), R[(i + 1) & 1][j] + (to + 1), sizeof(short) * (up[i][j] - (to + 1) + 1));
          for (int k = to + 1; k <= up[i][j]; k++) {
            sum += row[k][R[i & 1][j][k] + 1] - row[k][L[i & 1][j][k]];
          }
        }
      } else {
        if (i < n - 1 && a[i + 1][j] != 0) {
          L[i & 1][j][i] = j;
          R[i & 1][j][i] = j;
          up[i][j] = up[i + 1][j];
          memcpy(L[i & 1][j] + (i + 1), L[(i + 1) & 1][j] + (i + 1), sizeof(short) * (up[i][j] - (i + 1) + 1));
          memcpy(R[i & 1][j] + (i + 1), R[(i + 1) & 1][j] + (i + 1), sizeof(short) * (up[i][j] - (i + 1) + 1));
        } else {
          if (j < n - 1 && a[i][j + 1] != 0) {
            L[i & 1][j][i] = j;
            R[i & 1][j][i] = R[i & 1][j + 1][i];
            up[i][j] = up[i][j + 1];
            memcpy(L[i & 1][j] + (i + 1), L[i & 1][j + 1] + (i + 1), sizeof(short) * (up[i][j] - (i + 1) + 1));
            memcpy(R[i & 1][j] + (i + 1), R[i & 1][j + 1] + (i + 1), sizeof(short) * (up[i][j] - (i + 1) + 1));
          } else {
            L[i & 1][j][i] = j;
            R[i & 1][j][i] = j;
            up[i][j] = i;
          }
        }
        for (int k = i; k <= up[i][j]; k++) {
          sum += row[k][R[i & 1][j][k] + 1] - row[k][L[i & 1][j][k]];
        }
      }
      ans += a[i][j] * (sum - a[i][j]);
      put(i, j);
    }
  }
  cout << ans << '\n';
  cerr << "time = " << clock() << " ms" << "\n";
  return 0;
}

Measured memory_bytes: 46976