#include using namespace std; #define ALL(x) (x).begin(),(x).end() #define REP(i, n) for(ll i=0; i<(ll)(n); i++) #define PER(i, n) for(ll i=(ll)(n)-1; i>=0; i--) template int LB(const vector& v, T x) { return lower_bound(ALL(v),x)-v.begin(); } template int UQ(T& v) { sort(ALL(v)); v.erase(unique(ALL(v)),v.end()); return v.size(); } template bool chmax(T &a, T b) { return a bool chmin(T &a, T b) { return a>b ? a=b, true : false; } template using rpriority_queue = priority_queue,greater>; using ll=long long; const int INF=1e9+10; const ll INFL=4e18; using ld=long double; using lll=__int128_t; using ull=unsigned long long; using VI=vector; using VVI=vector; using VL=vector; using VVL=vector; using PL=pair; using VP=vector; using WG=vector>>; #ifdef LOCAL #include "./debug.hpp" #else #define debug(...) #define print_line #endif /// @brief セグメント木 template struct SegTree { using Type=typename Monoid::Type; SegTree()=default; /// @brief 要素数 n のセグ木を構築する SegTree(int n) { this->n=n; dat.assign(n<<1,Monoid::id()); cand.reserve(100),cand_l.reserve(100),cand_r.reserve(100); } /// @brief 配列 v からセグ木を構築する /// @note O(N) SegTree(const vector& v) { this->n=v.size(); dat.assign(n<<1,Monoid::id()); for(int i=0; i0; i--) dat[i]=Monoid::op(dat[i<<1],dat[i<<1|1]); cand.reserve(100),cand_l.reserve(100),cand_r.reserve(100); } /// @brief i 番目の要素を x に変更する /// @note O(log(N)) void set(int i, Type x) { i+=n; dat[i]=x; while(i>>=1) dat[i]=Monoid::op(dat[i<<1],dat[i<<1|1]); } /// @brief 区間 [l, r) のモノイド積を返す /// @note O(log(N)) Type fold(int l, int r) { Type retl=Monoid::id(),retr=Monoid::id(); l+=n,r+=n; while(l>=1,r>>=1; } return Monoid::op(retl,retr); } /// @brief 区間 [l, x) のモノイド積が f を満たすような最大の x >= l を返す /// @attention `f(Monoid::id())=true` が成り立つ必要がある /// @note O(log(N)) template int find_right(int l, F f) { assert(f(Monoid::id())); if(l==n) return n; l+=n; int r=n+n; cand_l.clear(),cand_r.clear(); while(l>=1,r>>=1; } cand=cand_l; reverse(cand_r.begin(),cand_r.end()); cand.insert(cand.end(),cand_r.begin(),cand_r.end()); Type val=Monoid::id(); for(int i:cand) { if(f(Monoid::op(val,dat[i]))) { val=Monoid::op(val,dat[i]); } else { while(i int find_left(int r,F f) { assert(f(Monoid::id())); if(r==0) return 0; r+=n; int l=n; cand_l.clear(),cand_r.clear(); while(l>=1,r>>=1; } cand=cand_r; reverse(cand_l.begin(),cand_l.end()); cand.insert(cand.end(),cand_l.begin(),cand_l.end()); Type val=Monoid::id(); for(int i:cand) { if(f(Monoid::op(dat[i],val))) { val=Monoid::op(dat[i],val); } else { while(i dat; vector cand,cand_l,cand_r; }; /// @brief モノイド namespace Monoid { /// @brief Minモノイド /// @tparam max_value 単位元 template struct Min { using Type=T; static Type id() { return max_value; } static Type op(const Type& a, const Type& b) { return min(a,b); } }; /// @brief Maxモノイド /// @tparam min_value 単位元 template struct Max { using Type=T; static Type id() { return min_value; } static Type op(const Type& a, const Type& b) { return max(a,b); } }; /// @brief 和 template struct Sum { using Type=T; static Type id() { return 0; } static Type op(const Type& a, const Type& b) { return a+b; } }; /// @brief (和,区間の長さ) template struct SumPair { using Type=pair; static Type id() { return make_pair(T(0),0); } static Type op(const Type& a, const Type& b) { return {a.first+b.first,a.second+b.second}; } }; } /// @brief 区間クエリ namespace RangeQuery { /// @brief 1点変更 / 区間 min template struct Min { using Type=struct SegTree>; }; /// @brief 1点変更 / 区間 max template struct Max { using Type=struct SegTree>; }; /// @brief 1点変更 / 区間和 template struct Sum { using Type=struct SegTree>; }; } //---------------------------------------------------------- void solve() { ll N; cin>>N; VL A(N); REP(i,N) cin>>A[i]; VL P(N+1); REP(i,N) P[i+1]=P[i]+A[i]; auto Q=P; UQ(Q); debug(P,Q); using segt=RangeQuery::Max::Type; vector small(3,segt(N+1)); vector large(3,segt(N+1)); //small[0].set(LB(Q,0ll),0); large[0].set(LB(Q,0ll),0); for(ll i=1; i<=N; i++) { ll idx=LB(Q,P[i]); for(int t=1; t<=2; t++) { ll res=-INFL; chmax(res,small[t-1].fold(0,idx+1)+P[i]); chmax(res,large[t-1].fold(idx,N+1)-P[i]); debug(i,t,res); if(i==N && t==2) { cout<>T; while(T--) solve(); }