#include using namespace std; using ll = long long; constexpr char newl = '\n'; // https://noshi91.hatenablog.com/entry/2019/03/31/174006 template struct ModInt { using u64 = std::uint_fast64_t; static constexpr u64 MOD = Modulus; u64 val; constexpr ModInt(const u64 x = 0) noexcept : val(x % MOD) {} constexpr ModInt operator+() const noexcept { return ModInt(*this); } constexpr ModInt operator-() const noexcept { ModInt res(*this); if (res.val != 0) res.val = MOD - res.val; return res; } constexpr bool operator==(const ModInt& rhs) const noexcept { return val == rhs.val; } constexpr bool operator!=(const ModInt& rhs) const noexcept { return val != rhs.val; } // prefix increment/decrement constexpr ModInt& operator++() noexcept { return *this += ModInt(1); } constexpr ModInt& operator--() noexcept { return *this -= ModInt(1); } // postfix increment/decrement constexpr ModInt operator++(int) noexcept { ModInt tmp(*this); ++*this; return tmp; } constexpr ModInt operator--(int) noexcept { ModInt tmp(*this); --*this; return tmp; } constexpr ModInt operator+(const ModInt& rhs) const noexcept { return ModInt(*this) += rhs; } constexpr ModInt operator-(const ModInt& rhs) const noexcept { return ModInt(*this) -= rhs; } constexpr ModInt operator*(const ModInt& rhs) const noexcept { return ModInt(*this) *= rhs; } constexpr ModInt operator/(const ModInt& rhs) const noexcept { return ModInt(*this) /= rhs; } constexpr ModInt& operator+=(const ModInt& rhs) noexcept { val += rhs.val; if (val >= MOD) val -= MOD; return *this; } constexpr ModInt& operator-=(const ModInt& rhs) noexcept { if (val < rhs.val) val += MOD; val -= rhs.val; return *this; } constexpr ModInt& operator*=(const ModInt& rhs) noexcept { val = val * rhs.val % MOD; return *this; } // prime Modulus only constexpr ModInt& operator/=(const ModInt& rhs) noexcept { return *this *= rhs.inv(); } // prime Modulus only constexpr ModInt inv() const noexcept { return pow(*this, MOD - 2); } }; template constexpr ModInt pow(ModInt x, std::uint_fast64_t n) { ModInt res(1); while (n) { if (n & 1) res *= x; x *= x; n >>= 1; } return res; } template istream& operator>>(istream& is, ModInt& x) { std::uint_fast64_t val; is >> val; x = ModInt(val); return is; } template ostream& operator<<(ostream& os, const ModInt& x) { return os << x.val; } using mint = ModInt<1000000007>; int main() { cin.tie(nullptr); ios::sync_with_stdio(false); int H, W; cin >> H >> W; vector< vector > a(H, vector(W)); vector rs(H, 1), cs(W, 1); vector r0(H, 0), c0(W, 0); mint all_mul = 1; int cnt0 = 0; for (int i = 0; i < H; i++) { for (int j = 0; j < W; j++) { cin >> a[i][j]; if (a[i][j] == 0) { r0[i]++; c0[j]++; cnt0++; } else { rs[i] *= a[i][j]; cs[j] *= a[i][j]; all_mul *= a[i][j]; } } } int q; cin >> q; for (int i = 0; i < q; i++) { int r, c; cin >> r >> c; --r; --c; if (cnt0 - r0[r] - c0[c] + (a[r][c] == 0) > 0) { cout << 0 << newl; } else { mint ans = all_mul / rs[r] / cs[c]; if (a[r][c] != 0) ans *= a[r][c]; cout << ans << newl; } } return 0; }