/*
_/_/_/_/    _/_/_/_/_/  _/_/_/
_/      _/      _/    _/      _/
_/      _/      _/    _/      _/
_/      _/      _/    _/      _/
_/      _/      _/    _/  _/  _/
_/      _/  _/  _/    _/    _/_/
_/_/_/_/      _/_/     _/_/_/_/_/

_/_/_/_/    _/    _/  _/      _/
_/      _/   _/  _/   _/_/  _/_/
_/      _/    _/_/    _/ _/_/ _/
_/      _/     _/     _/  _/  _/
_/      _/    _/_/    _/      _/
_/      _/   _/  _/   _/      _/
_/_/_/_/    _/    _/  _/      _/

_/_/_/_/_/ _/_/_/_/_/ _/_/_/_/_/
    _/         _/     _/
    _/         _/     _/
    _/         _/     _/_/_/_/
    _/         _/     _/
    _/         _/     _/
    _/     _/_/_/_/_/ _/_/_/_/_/

_/_/_/_/_/ _/_/_/_/_/ _/_/_/_/_/
    _/         _/     _/
    _/         _/     _/
    _/         _/     _/_/_/_/
    _/         _/     _/
    _/         _/     _/
    _/     _/_/_/_/_/ _/_/_/_/_/
*/
#include <cassert>
#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <cmath>
#include <ctime>
#include <cctype>

#include <algorithm>
#include <random>
#include <bitset>
#include <queue>
#include <functional>
#include <set>
#include <map>
#include <vector>
#include <chrono>
#include <iostream>
#include <iomanip>
#include <limits>
#include <numeric>

#define LOG(FMT...) fprintf(stderr, FMT)

using namespace std;

typedef long long ll;
typedef unsigned long long ull;

// mt19937 rng(chrono::steady_clock::now().time_since_epoch().count());

template<class T>
istream &operator>>(istream &is, vector<T> &v) {
  for (T &x : v)
    is >> x;
  return is;
}

template<class T>
ostream &operator<<(ostream &os, const vector<T> &v) {
  if (!v.empty()) {
    os << v.front();
    for (int i = 1; i < v.size(); ++i)
      os << ' ' << v[i];
  }
  return os;
}

const int P = 1000000007;

int norm(int x) { return x >= P ? (x - P) : x; }

void add(int &x, int y) { if ((x += y) >= P) x -= P; }

void sub(int &x, int y) { if ((x -= y) < 0) x += P; }

void exGcd(int a, int b, int &x, int &y) {
  if (!b) {
    x = 1;
    y = 0;
    return;
  }
  exGcd(b, a % b, y, x);
  y -= a / b * x;
}

int inv(int a) {
  int x, y;
  exGcd(a, P, x, y);
  return norm(x + P);
}

const int N = 65;

int n, k;
char g[N][N];
int f[N][N][N], sum[N][N][N];
int l[N][N][N], r[N][N][N];

int main() {
#ifdef ELEGIA
  freopen("test.in", "r", stdin);
  int nol_cl = clock();
#endif
  ios::sync_with_stdio(false);
  cin.tie(nullptr);

  int q;
  cin >> n >> k >> q;
  for (int i = 1; i <= n; ++i) {
    cin >> (g[i] + 1);
    for (int j = 1; j <= n; ++j) g[i][j] -= '0';
  }
  for (int i = 1; i <= k; ++i) {
    for (int u = 1; u <= n; ++u) {
      l[i][u][u] = r[i][u][u] = 1;
      for (int v = 1; v <= n; ++v)
        for (int w = 1; w <= n; ++w)
          if (g[w][u]) // v - w - [u]
            add(l[i][u][v], sum[i - 1][v][w]);
      for (int v = 1; v <= n; ++v)
        for (int w = 1; w <= n; ++w)
          if (g[u][v]) // [u] - v - w
            add(r[i][u][w], sum[i - 1][v][w]);
      for (int v = 1; v <= n; ++v)
        for (int w = 1; w <= n; ++w)
          f[i][v][w] = (f[i][v][w] + l[i][u][v] * (ull) r[i][u][w]) % P;
    }
    for (int j = 1; j <= n; ++j)
      for (int t = 1; t <= n; ++t)
        sum[i][j][t] = norm(sum[i - 1][j][t] + f[i][j][t]);
  }
  while (q--) {
    int a, s, b, t;
    cin >> a >> s >> b >> t;
    static int x[N][N], y[N][N], xs[N][N], ys[N][N];
    memset(x, 0, sizeof(x));
    memset(y, 0, sizeof(y));
    memset(xs, 0, sizeof(xs));
    memset(ys, 0, sizeof(ys));
    x[a][s] = 1;
    for (int i = a; i <= k; ++i) {
      for (int u = 1; u <= n; ++u)
        for (int v = 1; v <= n; ++v)
          if (g[u][v])
            add(x[i][v], xs[i - 1][u]);
      for (int u = 1; u <= n; ++u)
        for (int v = 1; v <= n; ++v)
          xs[i][v] = (xs[i][v] + x[i][u] * (ull) r[i][u][v]) % P;
      for (int u = 1; u <= n; ++u) add(xs[i][u], xs[i - 1][u]);
    }
    y[b][t] = 1;
    for (int i = b; i <= k; ++i) {
      for (int u = 1; u <= n; ++u)
        for (int v = 1; v <= n; ++v)
          if (g[u][v])
            add(y[i][u], ys[i - 1][v]);
      for (int u = 1; u <= n; ++u)
        for (int v = 1; v <= n; ++v)
          ys[i][u] = (ys[i][u] + y[i][v] * (ull) l[i][v][u]) % P;
      for (int u = 1; u <= n; ++u) add(ys[i][u], ys[i - 1][u]);
    }

    int ans = 0;
    for (int i = 1; i <= k; ++i)
      for (int j = 1; j <= n; ++j) {
        ans = (ans + x[i][j] * (ull) y[i][j]) % P;
      }
    cout << ans << '\n';
  }

#ifdef ELEGIA
  LOG("Time: %dms\n", int ((clock()
          -nol_cl) / (double)CLOCKS_PER_SEC * 1000));
#endif
  return 0;
}