老早以前的坑了
题意大概是,给定平面内 n(≤2000) 个点,每个点有个可正可负的点权,求两条平行的直线,使得所在在这两条直线之间的点权和最大
首先,如果我们知道这两条直线的斜率,假设这条线没有把这些点分成两个点集,那么可以把点按照距离直线的距离排序,记排序后的数列为A。
转动直线,那么会有一些点对交换在A中的位置,总共有 O(n2) 次交换,把这些交换按极角排序,依次处理就好。
然后点权的最大值就是经典的最大字段和,因为在线,所以可以用线段树维护
#include <cstdio>
#include <iostream>
#include <algorithm>
#include <cstring>
#include <cmath>
#include <assert.h>
using namespace std;
typedef long long ll;
const int N=5010;
int n,cnt,b[N],c[N];
ll ans;
struct node{
ll x,y,g,w;
node(ll _x=0,ll _y=0){x=_x,y=_y;}
friend bool operator <(node a,node b){
return a.x<b.x || (a.x==b.x && a.y<b.y);
}
friend ll operator *(node a,node b){
return (ll)a.x*b.y-(ll)a.y*b.x;
}
friend node operator -(node a,node b){
return node(b.x-a.x,b.y-a.y);
}
}p[N],a[N];
struct SWAP{
int x,y,p;
SWAP(int _x=0,int _y=0,int _p=0):x(_x),y(_y),p(_p){}
friend bool operator <(SWAP A,SWAP B){
if(A==B) return A.p>B.p;
return (a[A.y]-a[A.x])*(a[B.y]-a[B.x])>0;
}
friend bool operator ==(SWAP A,SWAP B){
return (a[A.y]-a[A.x])*(a[B.y]-a[B.x])==0;
}
}s[N*N];
ll mx[N<<2],mn[N<<2],tag[N<<2],tmx[N<<2];
inline char nc(){
static char buf[100000],*p1=buf,*p2=buf;
return p1==p2&&(p2=(p1=buf)+fread(buf,1,100000,stdin),p1==p2)?EOF:*p1++;
}
inline void rea(ll &x){
char c=nc(); x=0; int f=1;
for(;c>'9'||c<'0';f=c=='-'?-f:f,c=nc());for(;c>='0'&&c<='9';x=x*10+f*(c-'0'),c=nc());
}
inline void rea(int &x){
char c=nc(); x=0; int f=1;
for(;c>'9'||c<'0';f=c=='-'?-f:f,c=nc());for(;c>='0'&&c<='9';x=x*10+f*(c-'0'),c=nc());
}
inline void add(int g,ll x){
mx[g]+=x; mn[g]+=x; tag[g]+=x;
}
inline void Push(int g){
if(!tag[g]) return ;
add(g<<1,tag[g]); add(g<<1|1,tag[g]);
tag[g]=0;
}
inline void Up(ll g){
tmx[g]=max(max(tmx[g<<1],tmx[g<<1|1]),mx[g<<1|1]-mn[g<<1]);
mx[g]=max(mx[g<<1],mx[g<<1|1]);
mn[g]=min(mn[g<<1],mn[g<<1|1]);
}
void Add(int g,int L,int R,int l,int r,ll x){
if(L==l && r==R) return add(g,x);
int mid=L+R>>1; Push(g);
if(r<=mid) Add(g<<1,L,mid,l,r,x);
else if(l>mid) Add(g<<1|1,mid+1,R,l,r,x);
else Add(g<<1,L,mid,l,mid,x),Add(g<<1|1,mid+1,R,mid+1,r,x);
Up(g);
}
int main(){
rea(n);
for(int i=1;i<=n;i++)
rea(p[i].x),rea(p[i].y),rea(p[i].w),p[i].g=i,a[i]=p[i];
sort(p+1,p+1+n);
memset(tmx,-0x7f,sizeof(tmx));
for(int i=1;i<=n;i++)
b[p[i].g]=i,c[i]=p[i].g,Add(1,0,n,i,n,p[i].w);
for(int i=1;i<=n;i++)
for(int j=i+1;j<=n;j++){
s[++cnt]=SWAP(p[i].g,p[j].g,max(i,j));
}
sort(s+1,s+1+cnt);
ll mn=0,pre=0;
for(int i=1;i<=n;i++){
pre+=p[i].w;
if(p[i].x==p[i+1].x && i!=n) continue;
mn=min(mn,pre); ans=max(pre-mn,ans);
}
for(int i=1;i<=cnt;i++){
if(b[s[i].x]>b[s[i].y]) continue;
swap(c[b[s[i].x]],c[b[s[i].y]]);
Add(1,0,n,b[s[i].x],b[s[i].y]-1,a[s[i].y].w-a[s[i].x].w);
swap(b[s[i].x],b[s[i].y]);
if(s[i]==s[i+1] && i!=cnt) continue;
ans=max(ans,tmx[1]);
}
printf("%lld\n",max(ans,0LL));
return 0;
}