19 hdu多校 Snowy Smile //线段树维护最大子矩阵和

题意：n(<2000) 个点，xi(1e9),yi(1e9) ，wi有正有负，用矩形框出一个最大的sigma（wi）。

思路：枚举长的左端点，不停向线段树加入点，离散化后用线段树维护最大字段和。

单点动态最大字段和(线段树维护）

维护：

struct Tno{
    LL sum,pa,suf,mx; //和,最大前缀,最大后缀,最大子段
}T[MX];
void pushUp(int rt){
    T[rt].sum=T[rt<<1].sum+T[rt<<1|1].sum;
    T[rt].pa=max(T[rt<<1].pa,T[rt<<1].sum+T[rt<<1|1].pa);
    T[rt].suf=max(T[rt<<1|1].suf,T[rt<<1|1].sum+T[rt<<1].suf);
    T[rt].mx=max(max(T[rt<<1].mx,T[rt<<1|1].mx),T[rt<<1].suf+T[rt<<1|1].pa);
}

#include<bits/stdc++.h>
using namespace std;
#define LL long long
#define FI first
#define SE second
#define MP make_pair
#define PII pair<int,int>
const LL mod = 1e9+7;
const int LG = 50;
const int MX = 1e6+10;
struct no{int x,y;LL w;}a[2005];
bool cmpX(no x,no y){return x.x<y.x;}
vector<int>v;
struct Tno{
    LL sum,pa,suf,mx; //和,最大前缀,最大后缀,最大子段
}T[MX];
void pushUp(int rt){
    T[rt].sum=T[rt<<1].sum+T[rt<<1|1].sum;
    T[rt].pa=max(T[rt<<1].pa,T[rt<<1].sum+T[rt<<1|1].pa);
    T[rt].suf=max(T[rt<<1|1].suf,T[rt<<1|1].sum+T[rt<<1].suf);
    T[rt].mx=max(max(T[rt<<1].mx,T[rt<<1|1].mx),T[rt<<1].suf+T[rt<<1|1].pa);
}
void build(int l,int r,int rt){
    if(l==r){
        T[rt].sum=T[rt].pa=T[rt].suf=T[rt].mx=0;
        return;
    }
    int m=(l+r)>>1;
    build(l,m,rt<<1);
    build(m+1,r,rt<<1|1);
    pushUp(rt);
}
void update(int l,int r,int rt,int p,LL c){
    if(l==r){
        T[rt].mx+=c;T[rt].pa+=c;
        T[rt].sum+=c;T[rt].suf+=c;
        return;
    }
    int m=(l+r)>>1;
    if(m>=p) update(l,m,rt<<1,p,c);
    else update(m+1,r,rt<<1|1,p,c);
    pushUp(rt);
}
int main(){
    int TT;cin>>TT;while(TT--){
        int n;scanf("%d",&n);v.clear();
        for(int i=1;i<=n;i++)scanf("%d%d%lld",&a[i].x,&a[i].y,&a[i].w);
        sort(a+1,a+1+n,cmpX);
        for(int i=1;i<=n;i++){
            v.push_back(a[i].y);
        }
        sort(v.begin(),v.end());v.erase(unique(v.begin(),v.end()),v.end());
        LL ans=0;
        for(int i=1;i<=n;i++) {
            if(a[i].x==a[i-1].x&&i!=1)continue;
            build(1,n,1);
            for(int j=i;j<=n;j++){
                int id=lower_bound(v.begin(),v.end(),a[j].y)-v.begin()+1;
                update(1,n,1,id,a[j].w);
                if(a[j].x==a[j+1].x&&j!=n) continue;
                ans=max(ans,T[1].mx);
            }
        }
        cout<<ans<<'\n';
    }
    return 0;
}
/*
10
4
0 1 10
1 0 10
-1 0 10
0 -1 10
6
0 1 123456789
1 0 123456789
-1 0 123456789
0 -1 123456789
0 0 -1000000000000
1 1 1
*/