結果

問題 No.1141 田グリッド
ユーザー shimarutshimarut
提出日時 2021-03-28 18:54:52
言語 C++14
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 567 ms / 2,000 ms
コード長 4,240 bytes
コンパイル時間 3,576 ms
コンパイル使用メモリ 198,928 KB
実行使用メモリ 6,272 KB
最終ジャッジ日時 2024-11-29 09:20:09
合計ジャッジ時間 8,881 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge3
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
5,248 KB
testcase_01 AC 1 ms
5,248 KB
testcase_02 AC 1 ms
5,248 KB
testcase_03 AC 2 ms
5,248 KB
testcase_04 AC 1 ms
5,248 KB
testcase_05 AC 2 ms
5,248 KB
testcase_06 AC 2 ms
5,248 KB
testcase_07 AC 2 ms
5,248 KB
testcase_08 AC 3 ms
5,248 KB
testcase_09 AC 3 ms
5,248 KB
testcase_10 AC 4 ms
5,248 KB
testcase_11 AC 5 ms
5,248 KB
testcase_12 AC 4 ms
5,248 KB
testcase_13 AC 129 ms
5,248 KB
testcase_14 AC 138 ms
6,272 KB
testcase_15 AC 135 ms
5,248 KB
testcase_16 AC 129 ms
5,248 KB
testcase_17 AC 131 ms
5,248 KB
testcase_18 AC 118 ms
5,248 KB
testcase_19 AC 116 ms
5,248 KB
testcase_20 AC 116 ms
5,248 KB
testcase_21 AC 135 ms
5,248 KB
testcase_22 AC 131 ms
5,248 KB
testcase_23 AC 128 ms
5,248 KB
testcase_24 AC 135 ms
5,248 KB
testcase_25 AC 554 ms
5,248 KB
testcase_26 AC 168 ms
5,248 KB
testcase_27 AC 174 ms
5,248 KB
testcase_28 AC 567 ms
5,248 KB
testcase_29 AC 120 ms
5,248 KB
testcase_30 AC 123 ms
5,248 KB
testcase_31 AC 116 ms
5,248 KB
testcase_32 AC 116 ms
5,248 KB
testcase_33 AC 114 ms
5,248 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#pragma region
#define _USE_MATH_DEFINES
#include <iostream>
#include <string>
#include <algorithm>
#include <cmath>
#include <cstdlib>
#include <vector>
#include <map>
#include <queue>
#include <stack>
#include <set>
#include <list>
#include <iomanip>
#include <cstdint>
#include <bitset>
#include <sstream>
#include <regex>
#include <fstream>
#include <array>
#include <atcoder/all>

using namespace atcoder;
using mint = modint1000000007;
//using mint = modint998244353;
//using mint = modint;

using namespace std;
typedef long long ll;
using vi = vector<int>;
using vvi = vector<vi>;
using vl = vector<ll>;
using vvl = vector<vl>;
using vmint = vector<mint>;
using vvmint = vector<vmint>;
using pint = pair<int, int>;
using pll = pair<ll, ll>;
//#define rep(i, s, e) for (int(i) = (s); (i) < (e); ++(i))
#define rep(i, e) for (int(i) = 0; (i) < (e); ++(i))
#define rrep(i, s) for (int(i) = (s) - 1; (i) >= 0; --(i))
#define all(x) x.begin(),x.end()
#define rall(x) x.rbegin(),x.rend()

#pragma region UnionFind
struct UnionFind
{
	vector<int> par;

	UnionFind(int n) : par(n, -1) {}
	void init(int n) { par.assign(n, -1); }

	int root(int x)
	{
		if (par[x] < 0) return x;
		else return par[x] = root(par[x]);
	}

	bool issame(int x, int y)
	{
		return root(x) == root(y);
	}

	bool merge(int x, int y)
	{
		x = root(x); y = root(y);
		if (x == y) return false;
		if (par[x] > par[y]) swap(x, y);
		par[x] += par[y];
		par[y] = x;
		return true;
	}

	int size(int x)
	{
		return -par[root(x)];
	}
};
#pragma endregion
#pragma region GCD
ll gcd(ll a, ll b)
{
	if (b == 0)return a;
	return gcd(b, a%b);
}
#pragma endregion
#pragma region LCM
ll lcm(ll a, ll b)
{
	return a / gcd(a, b) * b;
}
#pragma endregion
#pragma region chmin
template<class T> inline bool chmin(T& a, T b)
{
	if (a > b)
	{
		a = b;
		return true;
	}
	return false;
}
#pragma endregion
#pragma region chmax
template<class T> inline bool chmax(T& a, T b)
{
	if (a < b)
	{
		a = b;
		return true;
	}
	return false;
}
#pragma endregion
#pragma region グリッド内チェック
bool out(int x, int y, int h, int w)
{
	if (x < 0 || h <= x || y < 0 || w <= y)return true;
	else return false;
}
#pragma endregion
#pragma region Dijkstra
vl dijkstra(vector<vector<pair<int, ll>>> v, int s)
{
	ll INF = 1e18;
	int MAX = 1e6;
	vl res(MAX, INF);
	priority_queue<pair<ll, int>, vector<pair<ll, int>>, greater<pair<ll, int>>> q;
	q.push({ 0,s });
	while (!q.empty())
	{
		int now;
		ll d;
		tie(d, now) = q.top();
		q.pop();
		if (!chmin(res[now], d))continue;
		for (auto p : v[now])
		{
			int next;
			ll c;
			tie(next, c) = p;
			if (res[next] <= res[now] + c)continue;
			q.push({ res[now] + c,next });
		}
	}
	return res;
}
#pragma endregion
#pragma endregion


int main()
{
	int h, w; cin >> h >> w;
	bool includes_zero = false;
	vvl a(h, vl(w));
	rep(i, h)rep(j, w)
	{
		cin >> a[i][j];
		if (a[i][j] == 0)includes_zero = true;
	}	
	if (!includes_zero)
	{
		vmint rs(h), cs(w);
		mint s = 1;
		rep(i, h)rep(j, w)s *= a[i][j];
		rep(i, h)
		{
			mint t = 1;
			rep(j, w)t *= a[i][j];
			rs[i] = t;
		}
		rep(i, w)
		{
			mint t = 1;
			rep(j, h)t *= a[j][i];
			cs[i] = t;
		}
		int q; cin >> q;
		while (q--)
		{
			int r, c; cin >> r >> c;
			--r, --c;
			mint res = s;
			res *= a[r][c];
			res /= rs[r];
			res /= cs[c];
			cout << res.val() << endl;
		}
	}
	else
	{
		vvi zero(h, vi(w));
		int zero_cnt = 0;
		rep(i, h)rep(j, w)if (a[i][j] == 0)
		{
			zero[i][j] = 1;
			++zero_cnt;
		}
		vi r_zero(h), c_zero(w);
		rep(i, h)
		{
			int cnt = 0;
			rep(j, w)if (a[i][j] == 0)++cnt;
			r_zero[i] = cnt;
		}
		rep(i, w)
		{
			int cnt = 0;
			rep(j, h)if (a[j][i] == 0)++cnt;
			c_zero[i] = cnt;
		}
		mint s = 1;
		rep(i, h)rep(j, w)if (a[i][j] != 0)s *= a[i][j];
		map<pint, mint> mp;
		int q; cin >> q;
		while (q--)
		{
			int r, c; cin >> r >> c;
			--r, --c;
			int z = r_zero[r] + c_zero[c] - zero[r][c];
			if (z == zero_cnt)
			{
				if (mp[{r, c}].val())cout << mp[{r, c}].val() << endl;
				else
				{
					mint res = s;
					rep(i, h)if (a[i][c] != 0)res /= a[i][c];
					rep(i, w)if (a[r][i] != 0)res /= a[r][i];
					if (a[r][c] != 0)res *= a[r][c];
					mp[{r, c}] = res;
					cout << res.val() << endl;
				}
			}
			else cout << 0 << endl;
		}
	}
}
0