結果
| 問題 |
No.2786 RMQ on Grid Path
|
| コンテスト | |
| ユーザー |
|
| 提出日時 | 2024-06-17 18:45:28 |
| 言語 | C++23 (gcc 13.3.0 + boost 1.87.0) |
| 結果 |
AC
|
| 実行時間 | 1,018 ms / 6,000 ms |
| コード長 | 3,713 bytes |
| コンパイル時間 | 4,332 ms |
| コンパイル使用メモリ | 270,216 KB |
| 実行使用メモリ | 41,988 KB |
| 最終ジャッジ日時 | 2024-06-17 18:46:01 |
| 合計ジャッジ時間 | 31,615 ms |
|
ジャッジサーバーID (参考情報) |
judge4 / judge2 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 2 |
| other | AC * 35 |
ソースコード
#include <bits/stdc++.h>
using namespace std;
#include <atcoder/modint>
using namespace atcoder;
#ifdef DEFINED_ONLY_IN_LOCAL
#include <dump.hpp>
#define dump(...) cpp_dump(__VA_ARGS__)
#else
#undef dump
#define dump(...)
#endif
#define rep1(i, a) for (int i = 0; i < (int)(a); i++)
#define rep2(i, a, b) for (int i = (int)(a); i < (int)(b); i++)
#define rep3(i, a, b, c) for (int i = (int)(a); i < (int)(b); i += (int)(c))
#define overloadRep(a, b, c, d, e, ...) e
#define rep(...) overloadRep(__VA_ARGS__, rep3, rep2, rep1)(__VA_ARGS__)
#define rrep(i, a, b) for (int i = (int)(a); i <= (int)(b); i++)
#define drep(i, a, b) for (int i = (int)(a); i >= (int)(b); i--)
#define all(a) a.begin(), a.end()
#define rall(a) a.rbegin(), a.rend()
using ll = long long;
using ull = unsigned long long;
using P = pair<ll, ll>;
using T = tuple<ll, ll, ll>;
const int inf = 1e9;
const ll INF = 1e18;
const int dx[4] = {0, 1, 0, -1};
const int dy[4] = {1, 0, -1, 0};
struct cincout {
cincout() {
ios_base::sync_with_stdio(false);
cin.tie(nullptr);
cout << fixed << setprecision(15);
}
} init;
template <class T>
inline bool chmax(T &a, T b) {
if (a < b) {
a = b;
return true;
}
return false;
}
template <class T>
inline bool chmin(T &a, T b) {
if (a > b) {
a = b;
return true;
}
return false;
}
template <class T>
istream &operator>>(istream &is, vector<T> &v) {
for (T &in : v) {
is >> in;
}
return is;
}
template <class T>
ostream &operator<<(ostream &os, const vector<T> &v) {
rep(i, (int)v.size()) {
os << v[i] << " \n"[i + 1 == (int)v.size()];
}
return os;
}
template <class T>
istream &operator>>(istream &is, vector<vector<T>> &vv) {
for (vector<T> &v : vv) {
is >> v;
}
return is;
}
template <class T>
ostream &operator<<(ostream &os, vector<vector<T>> &vv) {
for (vector<T> &v : vv) {
os << v;
}
return os;
}
template <class T1, class T2>
istream &operator>>(istream &is, pair<T1, T2> &p) {
is >> p.first >> p.second;
return is;
}
template <class T1, class T2>
ostream &operator<<(ostream &os, const pair<T1, T2> &p) {
os << p.first << " " << p.second << '\n';
return os;
}
bool bit(ll x, int p) {
return (x >> p) & 1;
}
bool out(int ni, int nj, int h, int w) {
return (ni < 0 or ni >= h or nj < 0 or nj >= w);
}
int pc(ll x) {
return __builtin_popcountll(x);
}
#include <atcoder/dsu>
int main(){
int h,w;cin>>h>>w;
vector a(h,vector<int>(w));cin>>a;
int q;cin>>q;
vector<P>qs(q);
rep(i,q){
int a,b,c,d;cin>>a>>b>>c>>d;
a--;b--;c--;d--;
qs[i]=P(a*w+b,c*w+d);
}
vector<T>es;
rep(i,h)rep(j,w)rep(d,4){
int ni=i+dx[d],nj=j+dy[d];
if(out(ni,nj,h,w))continue;
int x=i*w+j,y=ni*w+nj;
int mx=max(a[i][j],a[ni][nj]);
es.emplace_back(mx,x,y);
}
sort(all(es));
vector<int>l(q,0),r(q,h*w+1);
while(true){
bool exit=true;
vector<vector<int>>m(h*w+1);
rep(i,q){
if(r[i]-l[i]>1){
exit=false;
int mid=(l[i]+r[i])/2;
m[mid].push_back(i);
}
}
if(exit)break;
dsu uf(h*w);
int cur=0;
rep(i,1,h*w+1){
while(cur<es.size() and get<0>(es[cur])<=i){
auto[w,u,v]=es[cur];
uf.merge(u,v);
cur++;
}
for(auto j:m[i]){
auto[u,v]=qs[j];
if(uf.same(u,v))r[j]=i;
else l[j]=i;
}
}
}
rep(i,q)cout<<r[i]<<endl;
return 0;
}