传送门
文章目录
- 题意:
- 思路:
题意:
给你nnn个点(xi,yi)(x_i,y_i)(xi,yi),每个点有个价值cic_ici,现在你可以框一个正方形,要求左下角和右上角的坐标(x,y)(x,y)(x,y)必须x=yx=yx=y,也就是说必须在x=yx=yx=y这条直线上。现在你可以框一个正方形,被正方形框到的点算入总价值,求能得到的最大的总价值。总价值=被框的点的价值−-−正方形的长度。
思路:
在二维平面画了很多都没找到规律,看了题解才知道,这个题有一个很巧妙的转换,设正方形左下角坐标(x,x)(x,x)(x,x),右上角坐标(y,y)(y,y)(y,y),对于(xi,yi)(x_i,y_i)(xi,yi)如果他能被这个正方形包含在内部,那么必须满足x<=min(xi,yi)<=max(xi,yi)<=yx<=min(x_i,y_i)<=max(x_i,y_i)<=yx<=min(xi,yi)<=max(xi,yi)<=y。看到这个式子,很明显它可以将每个点转换成一个区间的形式,即[min(xi,yi),max(xi,yi)][min(x_i,y_i),max(x_i,y_i)][min(xi,yi),max(xi,yi)]。
现在我们重新描述以下这个问题,给你若干个区间,每个区间有一个价值,你需要覆盖一段区间,使得覆盖到的区间的价值−-−区间长度最大。当然这里覆盖是必须完全覆盖。
转换成这个问题,就比较容易做了,这里介绍一种线段树的做法。
首先需要发现一个性质,那就是我们选的区间起点一定是某个区间的起点,终点一定是某个区间的终点。这个比较显然。
那么我们将左端点按从小到大排序,让后倒着扫,每次都将当前左端点相等的区间都加入,即将[yi,1e9][y_i,1e9][yi,1e9]都加上cic_ici。这样可保证以左端点为起点能求解出最佳答案。
以上操作显然可仍线段树里面,现在我们问题就转换成了如何求出来tree(xi,1e9)max−(r−xi)tree(x_i,1e9)_{max}-(r-x_i)tree(xi,1e9)max−(r−xi),前面一部分就是线段树区间最大值的查询,对于后面,由于我们枚举的xix_ixi,所以xix_ixi已经知道了,我们可将其移动一下变成xi+tree(xi,1e9)max−rx_i+tree(x_i,1e9)_{max}-rxi+tree(xi,1e9)max−r,那么rrr怎么办呢?我们可以在建线段树的时候将初始值置为其rrr即可,查询的时候需要加上左端点。
让后可以离散化一下比较好些,当然也可以动态开点?
复杂度O(nlogn)O(nlogn)O(nlogn)
// Problem: F. Choose a Square
// Contest: Codeforces - Educational Codeforces Round 73 (Rated for Div. 2)
// URL: https://codeforces.com/contest/1221/problem/F
// Memory Limit: 256 MB
// Time Limit: 6000 ms
//
// Powered by CP Editor (https://cpeditor.org)//#pragma GCC optimize("Ofast,no-stack-protector,unroll-loops,fast-math")
//#pragma GCC target("sse,sse2,sse3,ssse3,sse4.1,sse4.2,avx,avx2,popcnt,tune=native")
//#pragma GCC optimize(2)
#include<cstdio>
#include<iostream>
#include<string>
#include<cstring>
#include<map>
#include<cmath>
#include<cctype>
#include<vector>
#include<set>
#include<queue>
#include<algorithm>
#include<sstream>
#include<ctime>
#include<cstdlib>
#include<random>
#include<cassert>
#define X first
#define Y second
#define L (u<<1)
#define R (u<<1|1)
#define pb push_back
#define mk make_pair
#define Mid ((tr[u].l+tr[u].r)>>1)
#define Len(u) (tr[u].r-tr[u].l+1)
#define random(a,b) ((a)+rand()%((b)-(a)+1))
#define db puts("---")
using namespace std;//void rd_cre() { freopen("d://dp//data.txt","w",stdout); srand(time(NULL)); }
//void rd_ac() { freopen("d://dp//data.txt","r",stdin); freopen("d://dp//AC.txt","w",stdout); }
//void rd_wa() { freopen("d://dp//data.txt","r",stdin); freopen("d://dp//WA.txt","w",stdout); }typedef long long LL;
typedef unsigned long long ULL;
typedef pair<int,int> PII;const int N=1000010,mod=1e9+7,INF=0x3f3f3f3f;
const LL inf=0x3f3f3f3f3f3f3f3f;
const double eps=1e-6;int n,se;
struct node {int X,Y,c;bool operator < (const node &x) const {return x.X>X;}
}p[N];
vector<int>v;
struct Node {int l,r;LL now,mx,lazy,id;
}tr[N<<2];int find(int x) {return lower_bound(v.begin(),v.end(),x)-v.begin()+1;
}void pushup(int u) {if(tr[L].mx>tr[R].mx) {tr[u].mx=tr[L].mx;tr[u].id=tr[L].id;} else {tr[u].mx=tr[R].mx;tr[u].id=tr[R].id;}
}void pushdown(int u) {LL lazy=tr[u].lazy; tr[u].lazy=0;tr[L].lazy+=lazy; tr[L].mx+=lazy;tr[R].lazy+=lazy; tr[R].mx+=lazy;
}void build(int u,int l,int r) {tr[u]={l,r,0,0,0,0};if(l==r) {tr[u].mx=-v[l-1];tr[u].id=v[l-1];return ;}build(L,l,Mid); build(R,Mid+1,r);pushup(u);
}void change(int u,int l,int r,int c) {if(tr[u].l>=l&&tr[u].r<=r) {tr[u].mx+=c; tr[u].lazy+=c;return;}pushdown(u);if(l<=Mid) change(L,l,r,c);if(r>Mid) change(R,l,r,c);pushup(u);
}pair<LL,int> query(int u,int l,int r) {if(tr[u].l>=l&&tr[u].r<=r) return {tr[u].mx,tr[u].id};pushdown(u);pair<LL,int> ans={-inf,0};if(l<=Mid) ans=max(ans,query(L,l,r));if(r>Mid) ans=max(ans,query(R,l,r));return ans;
}int main()
{
// ios::sync_with_stdio(false);
// cin.tie(0);LL ans=0;int l=mod,r=mod;scanf("%d",&n);for(int i=1;i<=n;i++) {scanf("%d%d%d",&p[i].X,&p[i].Y,&p[i].c);if(p[i].X>p[i].Y) swap(p[i].X,p[i].Y);v.pb(p[i].X); v.pb(p[i].Y);}sort(v.begin(),v.end()); v.erase(unique(v.begin(),v.end()),v.end());sort(p+1,p+1+n); se=v.size();for(int i=1;i<=n;i++) p[i].X=find(p[i].X),p[i].Y=find(p[i].Y);build(1,1,se);for(int i=n;i>=1;i--) {int xx=p[i].X;while(i>=1&&p[i].X==xx) change(1,p[i].Y,se,p[i].c),i--; i++;pair<LL,int> now=query(1,p[i].X,se);now.X+=v[p[i].X-1];if(now.X>ans) {ans=now.X;l=v[p[i].X-1],r=now.Y;}}printf("%lld\n%d %d %d %d\n",ans,l,l,r,r);return 0;
}
/**/