题意:
给
出
二
维
坐
标
轴
上
的
n
个
点
,
每
个
点
有
整
数
权
值
。
用
一
个
矩
形
包
含
若
干
个
点
使
得
总
权
值
和
最
大
。
给出 二维坐标轴上的n个点,每个点有整数权值。用一个矩形包含若干个点使得总权值和最大。
给出二维坐标轴上的n个点,每个点有整数权值。用一个矩形包含若干个点使得总权值和最大。
思路
(扯淡
拖到现在才补这题。。。
赛中尝试使用二维离散化前缀和(
w
a
wa
wa了并且找到了错误样例想法有误),又尝试了扫描线但是
T
T
T了(姿势不对?)
甚至还尝试了线段树但是就是没往最大子段和上面想)
正解:先对坐标离散化。 按照
x
x
x从小到大排序保证矩形每次都是从左往右去圈点,我们枚举每一个
x
x
x, 线段树维护对应的
y
y
y
坐标的权值。
这样子我们可以 O ( n 2 ) O(n^2) O(n2)遍历矩形的长(其实就是 x x x坐标), O ( l o g n ) O(logn) O(logn)更新 y y 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;
}