2019杭电多校6 hdu6638 Snowy Smile（二维最大矩阵和 线段树）

http://acm.hdu.edu.cn/showproblem.php?pid=6638

题意：给你一些点的权值，让找一个矩形圈住一部分点，问圈住点的最大权值和

分析：由于是稀疏图，明显要先把x,y坐标离散化，暴力是n^3？(枚举边界n^2，求和是n)显然过不了，那可以枚举y的边界，然后对于x就是最大子段和的问题了，用线段树维护，n^2logn可以过。

 #include<bits/stdc++.h>
 using namespace std;
 const int maxn = 4e3+;
 const int inf = 0x3f3f3f;
 typedef pair<int,int> P;
 typedef long long ll;
 int cas,n;
 struct point{
     int x,y,val;
     bool operator<(const point &a) const{
         return y<a.y;
     }
 }a[maxn];
 vector<int> p[maxn];
 int bx[maxn],by[maxn];

 struct node{
     ll sum,lmax,rmax,lrmax;
 }tree[maxn<<];

 inline void pushup(int rt){
     tree[rt].sum = tree[rt<<].sum+tree[rt<<|].sum;
     tree[rt].lmax = max(tree[rt<<].lmax,tree[rt<<].sum+tree[rt<<|].lmax);
     tree[rt].rmax = max(tree[rt<<|].rmax,tree[rt<<].rmax+tree[rt<<|].sum);
     tree[rt].lrmax = max(max(tree[rt<<].lrmax,tree[rt<<|].lrmax),tree[rt<<].rmax+tree[rt<<|].lmax);
 }

 inline void update(int L,int l,int r,int rt,int c){
     if(l==r){
         tree[rt].sum += 1ll*c;
         tree[rt].lrmax = tree[rt].lmax = tree[rt].rmax = tree[rt].sum;
         return;
     }
     int mid = l+r>>;
     if(L<=mid) update(L,l,mid,rt<<,c);
     else update(L,mid+,r,rt<<|,c);
     pushup(rt);
 }

 int main(){
     cin>>cas;
     while(cas--){
         cin>>n;
         for(int i=;i<=n;i++){
             cin>>a[i].x>>a[i].y>>a[i].val;
             bx[i] = a[i].x,by[i] = a[i].y;
         }
         sort(bx+,bx++n);
         int xlen = unique(bx+,bx++n)-bx-;
         sort(by+,by++n);
         int ylen = unique(by+,by++n)-by-;
         for(int i=;i<=n;i++){
             a[i].x = lower_bound(bx+,bx+xlen+,a[i].x)-bx;
             a[i].y = lower_bound(by+,by+ylen+,a[i].y)-by;
         }
         sort(a+,a++n);
         ll ans = ;
         for(int i=;i<=ylen;i++){
             memset(tree,,(xlen*+)*sizeof(node));
             int pos = ;
             while(a[pos].y<i&&pos<=n) pos++;
             for(int j=i;j<=ylen;j++){
                 while(a[pos].y<=j&&pos<=n){
                     update(a[pos].x,,xlen,,a[pos].val);
                     pos++;
                 }
                 ans = max(ans,tree[].lrmax);
             }

         }
         cout<<ans<<endl;
     }
     return ;
 }