2019杭电多校第六场hdu6638 Snowy Smile （最大权值和矩形）

传送门

题意：
$$给出二维坐标轴上的n个点，每个点有整数权值。用一个矩形包含若干个点使得总权值和最大。$$

思路
（扯淡
拖到现在才补这题。。。
赛中尝试使用二维离散化前缀和（$wa$了并且找到了错误样例想法有误），又尝试了扫描线但是$T$了(姿势不对？)

甚至还尝试了线段树但是就是没往最大子段和上面想）

正解：先对坐标离散化。 按照$x$从小到大排序保证矩形每次都是从左往右去圈点，我们枚举每一个$x$, 线段树维护对应的$y$
坐标的权值。

这样子我们可以$O(n^2)$遍历矩形的长（其实就是$x$坐标），$O(logn)$更新$y$坐标的最大子段和,最后更新答案。

#include<iostream>
#include<algorithm>
#include<cstdio>
#include<stdio.h>
#include<string.h>
#include<queue>
#include<cmath>
#include<map>
#include<set>
#include<vector> 
using namespace std;
#define inf 0x3f3f3f3f
#define lson l,mid,rt<<1
#define rson mid+1,r,rt<<1|1
#define mem(a,b) memset(a,b,sizeof(a));
#define lowbit(x)  x&-x;  
#define debugint(name,x) printf("%s: %d\n",name,x);
#define debugstring(name,x) printf("%s: %s\n",name,x);
typedef long long ll;
typedef unsigned long long ull;
const double eps = 1e-6;
const int maxn = 1e5+5;
const int mod = 1e9+7;
inline int read()
{
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
    return x*f;
}
struct node{
    int x,y;
    ll w;
    bool operator<(const node C)const{
        return x < C.x;
    }
}p[maxn];
vector<ll>xx,yy;

ll Maxsum[maxn<<2], lsum[maxn<<2], rsum[maxn<<2], sum[maxn<<2];
void pushup(int rt){
    sum[rt] = sum[rt<<1] + sum[rt<<1|1];
    Maxsum[rt] = max(max(Maxsum[rt<<1],Maxsum[rt<<1|1]), rsum[rt<<1]+lsum[rt<<1|1]);
    lsum[rt] = max(lsum[rt<<1], sum[rt<<1] + lsum[rt<<1|1]);
    rsum[rt] = max(rsum[rt<<1|1], sum[rt<<1|1] + rsum[rt<<1]); 
}
void build(int l,int r,int rt){
    Maxsum[rt] = lsum[rt] = rsum[rt] = sum[rt] = 0;
    if(l == r) return;

    int mid = (l+r)>>1;
    build(lson);
    build(rson);
    pushup(rt);
}

void update(int l,int r,int rt,int pos, ll val){
    if(l == r){
        sum[rt] = lsum[rt] = rsum[rt] = Maxsum[rt] = Maxsum[rt] + val;
        return ;
    }
    int mid = (l+r)>>1;
    if(pos <= mid) update(lson,pos,val);
    else update(rson,pos,val);
    pushup(rt); 
}


int main(){
    int T;
    scanf("%d", &T);
    while(T--){
        xx.clear(), yy.clear();
        int n;
        scanf("%d", &n);
        for(int i = 1; i <= n; i++){
            scanf("%d%d%lld", &p[i].x, &p[i].y, &p[i].w);
            xx.push_back(p[i].x), yy.push_back(p[i].y);
        }
        sort(xx.begin(), xx.end());
        xx.erase(unique(xx.begin(), xx.end()), xx.end());
        sort(yy.begin(), yy.end());
        yy.erase(unique(yy.begin(), yy.end()), yy.end());

        for(int i = 1; i <= n; i++){
            p[i].x = lower_bound(xx.begin(), xx.end(), p[i].x) - xx.begin() + 1;
            p[i].y = lower_bound(yy.begin(), yy.end(), p[i].y) - yy.begin() + 1;
            //printf("%d %d\n",p[i].x,p[i].y);
        }

        sort(p+1, p+n+1);
        ll ans = 0;
        int m = yy.size();
        for(int i = 1; i <= n; i++){  // x
            if(i > 1 && (p[i].x == p[i-1].x)) continue;
            build(1,m,1);
            update(1,m,1,p[i].y,p[i].w);
            for(int j = i+1; j <= n; j++){ // x
                if(p[j].x != p[j-1].x) ans = max(ans, Maxsum[1]);
                update(1,m,1,p[j].y,p[j].w);
            }
            ans = max(ans, Maxsum[1]);
        }

        printf("%lld\n",ans);
    }
    return 0;
}