結果
問題 | 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; }