結果
問題 |
No.3214 small square
|
ユーザー |
|
提出日時 | 2025-07-28 21:58:20 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 953 ms / 3,000 ms |
コード長 | 4,441 bytes |
コンパイル時間 | 3,252 ms |
コンパイル使用メモリ | 224,360 KB |
実行使用メモリ | 53,568 KB |
最終ジャッジ日時 | 2025-09-05 07:32:02 |
合計ジャッジ時間 | 30,356 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 40 |
ソースコード
#include<bits/stdc++.h> using namespace std; //#pragma GCC optimize("Ofast") #define rep(i,n) for(ll i=0;i<n;i++) #define repl(i,l,r) for(ll i=(l);i<(r);i++) #define per(i,n) for(ll i=(n)-1;i>=0;i--) #define perl(i,r,l) for(ll i=r-1;i>=l;i--) #define fi first #define se second #define pb push_back #define ins insert #define pqueue(x) priority_queue<x,vector<x>,greater<x>> #define all(x) (x).begin(),(x).end() #define CST(x) cout<<fixed<<setprecision(x) #define vtpl(x,y,z) vector<tuple<x,y,z>> #define rev(x) reverse(x); using ll=long long; using vl=vector<ll>; using vvl=vector<vector<ll>>; using pl=pair<ll,ll>; using vpl=vector<pl>; using vvpl=vector<vpl>; const ll MOD=1000000007; const ll MOD9=998244353; const int inf=1e9+10; const ll INF=4e18; //const ll dy[9]={1,0,-1,0,1,1,-1,-1,0}; //const ll dx[9]={0,1,0,-1,1,-1,1,-1,0}; template<class T> inline bool chmin(T& a, T b) { if (a > b) { a = b; return true; } return false; } template<class T> inline bool chmax(T& a, T b) { if (a < b) { a = b; return true; } return false; } template <typename X, typename M>//faがint lenになってる struct SegTreeLazy {//遅延セグ木 単位元に注意(updateなら選ばれない数、affineなら(1,0)) using FX = function<X(X, X)>; using FA = function<X(X, M, int)>; using FM = function<M(M, M)>; int n; FX fx; FA fa; FM fm; const X ex; const M em; vector<X> dat; vector<M> lazy; SegTreeLazy(int n_, FX fx_, FA fa_, FM fm_, X ex_, M em_) : n(), fx(fx_), fa(fa_), fm(fm_), ex(ex_), em(em_), dat(n_ * 4, ex), lazy(n_ * 4, em) { int x = 1; while (n_ > x) x *= 2; n = x; } void set(int i, X x) { dat[i + n - 1] = x; } void build() { for (int k = n - 2; k >= 0; k--) dat[k] = fx(dat[2 * k + 1], dat[2 * k + 2]); } /* lazy eval */ void eval(int k, int len) { if (lazy[k] == em) return; // 更新するものが無ければ終了 if (k < n - 1) { // 葉でなければ子に伝搬 lazy[k * 2 + 1] = fm(lazy[k * 2 + 1], lazy[k]); lazy[k * 2 + 2] = fm(lazy[k * 2 + 2], lazy[k]); } // 自身を更新 dat[k] = fa(dat[k],lazy[k],len);//fa(dat[k], fp(lazy[k], len)); lazy[k] = em; } void update(int a, int b, M x, int k, int l, int r) { eval(k, r - l); if (a <= l && r <= b) { // 完全に内側の時 lazy[k] = fm(lazy[k], x); eval(k, r - l); } else if (a < r && l < b) { // 一部区間が被る時 update(a, b, x, k * 2 + 1, l, (l + r) / 2); // 左の子 update(a, b, x, k * 2 + 2, (l + r) / 2, r); // 右の子 dat[k] = fx(dat[k * 2 + 1], dat[k * 2 + 2]); } } void update(int a, int b, M x) { update(a, b, x, 0, 0, n); } X query_sub(int a, int b, int k, int l, int r) { eval(k, r - l); if (r <= a || b <= l) { // 完全に外側の時 return ex; } else if (a <= l && r <= b) { // 完全に内側の時 return dat[k]; } else { // 一部区間が被る時 X vl = query_sub(a, b, k * 2 + 1, l, (l + r) / 2); X vr = query_sub(a, b, k * 2 + 2, (l + r) / 2, r); return fx(vl, vr); } } X query(int a, int b) { return query_sub(a, b, 0, 0, n); } X operator[](int i){ return query(i,i+1); } }; int main(){ ll n,m;cin >> n >> m;m*=2; vector<vl> g; rep(i,n){ ll a,b,c;cin >> a >> b >> c;a*=2;b*=2; g.push_back((vl){a-m,b,c}); g.push_back((vl){a+1,b,-c}); } vl z; for(auto p:g){ z.emplace_back(p[1]); z.emplace_back(p[1]+m); } sort(all(z)); z.erase(unique(all(z)),z.end()); SegTreeLazy<ll,ll> st(n*2,[](ll a,ll b){return max(a,b);},[](ll a,ll b,int len){return a+b;},[](ll a,ll b){return a+b;},-INF,0); rep(i,n*2)st.set(i,0);st.build(); //sort(all(g)); map<ll,vpl> mp; for(auto p:g){ mp[p[0]].push_back({p[1],p[2]}); } ll ans=0; for(auto [p,q]:mp){ for(auto r:q){ ll lx=lower_bound(all(z),r.first)-z.begin(); ll rx=lower_bound(all(z),r.first+m)-z.begin(); st.update(lx,rx+1,r.second); } chmax(ans,st.query(0,n*2)); } cout << ans << endl; }