ssoj2455有趣的有趣的家庭菜园（线段树）

【题意】有一个n块的线性菜园，每块菜园只有照到阳光（左右两边没有遮挡）才能收获果实卖出去价值为p，也可以除去费用为c，问最大利益是多少

【思路】枚举n块田地i为最高处，ans即为其左边最大利润加右边最大利润。dp o（n^2）会超时。用线段树维护左边（右边）最大值，点更新，与区间值的更改。

【代码】

#include <iostream>
#include <cstdlib>
#include <cstdio>
#include <cstring>
#include <cmath>
#include <algorithm>
#define ll long long
using namespace std;
const int maxn=100005;
const ll inf=10000000000009;
struct data{
    ll mx,val;
}tree[maxn*4];
ll n,h[maxn],H[maxn];
ll f[maxn],g[maxn],ans=0,p[maxn],c[maxn];
inline ll get(){
    char c;while(!isdigit(c=getchar()));
    int v=c-48;while(isdigit(c=getchar()))v=v*10+c-48;
    return v;
}
inline void down(int node){
    tree[node<<1].val+=tree[node].val;
    tree[(node<<1)+1].val+=tree[node].val;
    tree[node<<1].mx+=tree[node].val;
    tree[(node<<1)+1].mx+=tree[node].val;
    tree[node].val=0;
}
inline void mt(int node){
    tree[node].mx=max(tree[node<<1].mx,tree[(node<<1)+1].mx);
}
inline void build(int node,int l,int r){
    if(l==r){
		if(l==0)tree[node].mx=0;
		else tree[node].mx=-inf;
		tree[node].val=0;
	}
    else{
		int mid=(l+r)>>1;
	    build(node<<1,l,mid);
	    build((node<<1)+1,mid+1,r);
	    mt(node);
	}
}
inline ll query(int node,int l,int r,int x,int y){
    if(l>=x && r<=y)return tree[node].mx;
    else{
		down(node);
	    int mid=(l+r)>>1;
	    if(y<=mid)return query(node<<1,l,mid,x,y);
	    if(x>mid)return query((node<<1)+1,mid+1,r,x,y);
	    return max(query(node<<1,l,mid,x,y),query((node<<1)+1,mid+1,r,x,y));
	}
}
inline void add(int node,int l,int r,int p,ll x){
    if(l==r)tree[node].mx=x;
    else{
		down(node);
	    int mid=(l+r)>>1;
	    if(p<=mid)add(node<<1,l,mid,p,x);
	    else add((node<<1)+1,mid+1,r,p,x);
	    mt(node);
	}
}
inline void update(int node,int l,int r,int x,int y,ll v){
    if(l>=x && r<=y)tree[node].val+=v,tree[node].mx+=v;
    else{
	    down(node);
	    int mid=(l+r)>>1;
	    if(x<=mid)update(node<<1,l,mid,x,y,v);
	    if(y>mid)update((node<<1)+1,mid+1,r,x,y,v);
	    mt(node);
	}
}
int main(){
    n=get();
    for(int i=1;i<=n;++i)H[i]=h[i]=get(),p[i]=get(),c[i]=get();
    sort(H+1,H+1+n);
    int t=unique(H+1,H+1+n)-H-1;
    for(int i=1;i<=n;++i)h[i]=lower_bound(H+1,H+1+t,h[i])-H;
    build(1,0,n);
    for(int i=1;i<=n;++i){
	    f[i]=query(1,0,n,0,h[i])+p[i];
	    add(1,0,n,h[i],f[i]);
	    update(1,0,n,0,h[i]-1,-c[i]);
	}
	memset(tree,0,sizeof(tree));build(1,0,n);//如果不memset貌似会出事
	for(int i=n;i>=1;--i){
	    g[i]=query(1,0,n,0,h[i])+p[i];
	    add(1,0,n,h[i],g[i]);
	    update(1,0,n,0,h[i]-1,-c[i]);
	}
	for(int i=1;i<=n;++i)ans=max(ans,f[i]+g[i]-p[i]);
	printf("%lld\n",ans);
	return 0;
}