Points Division（线段树+DP）2019牛客暑期多校训练营（第一场）

题意：https://ac.nowcoder.com/acm/contest/881/I

给你n个平面上的点，每个点有a、b两个权值，现在让你划分成两个区域（要求所有A集合里的点不能在任何B集合里的点的右下方）。

求MAX（Sigma ai+Sigma bi）。

思路：

dp+线段树。

https://blog.csdn.net/qq_41194925/article/details/97079075

就是有一个要注意的地方：对于枚举的点我们算的是B集合所以dp【i】=max（1~i）+bi，这里必须是bi。

#define IOS ios_base::sync_with_stdio(0); cin.tie(0);
#include <cstdio>//sprintf islower isupper
#include <cstdlib>//malloc  exit strcat itoa system("cls")
#include <iostream>//pair
#include <fstream>//freopen("C:\Users\13606\Desktop\草稿.txt","r",stdin);
#include <bitset>
//#include <map>
//#include<unordered_map>
#include <vector>
#include <stack>
#include <set>
#include <string.h>//strstr substr
#include <string>
#include <time.h>//srand(((unsigned)time(NULL))); Seed n=rand()%10 - 0~9;
#include <cmath>
#include <deque>
#include <queue>//priority_queue<int, vector<int>, greater<int> > q;//less
#include <vector>//emplace_back
//#include <math.h>
//#include <windows.h>//reverse(a,a+len);// ~ ! ~ ! floor
#include <algorithm>//sort + unique : sz=unique(b+1,b+n+1)-(b+1);+nth_element(first, nth, last, compare)
using namespace std;//next_permutation(a+1,a+1+n);//prev_permutation
//******************
int abss(int a);
int lowbit(int n);
int Del_bit_1(int n);
int maxx(int a,int b);
int minn(int a,int b);
double fabss(double a);
void swapp(int &a,int &b);
clock_t __STRAT,__END;
double __TOTALTIME;
void _MS(){__STRAT=clock();}
void _ME(){__END=clock();__TOTALTIME=(double)(__END-__STRAT)/CLOCKS_PER_SEC;cout<<"Time: "<<__TOTALTIME<<" s"<<endl;}
//***********************
#define rint register int
#define fo(a,b,c) for(rint a=b;a<=c;++a)
#define fr(a,b,c) for(rint a=b;a>=c;--a)
#define mem(a,b) memset(a,b,sizeof(a))
#define pr printf
#define sc scanf
#define ls rt<<1
#define rs rt<<1|1
typedef long long ll;
//#define int ll
const double E=2.718281828;
const double PI=acos(-1.0);
//const ll INF=(1LL<<60);
const int inf=(1<<30);
const double ESP=1e-9;
const int mod=(int)1e9+7;
const int N=(int)1e6+10;

struct node
{
    int x,y,va,vb;
    friend bool operator<(node a,node b)
    {
        if(a.x==b.x)
            return a.y>b.y;
        return a.x<b.x;
    }
}a[N];
int temp[N];
int LS(int n)
{
    int m=0;
    for(int i=1;i<=n;++i)
        temp[++m]=a[i].y;
    sort(temp+1,temp+1+m);
    m=unique(temp+1,temp+1+m)-temp-1;
    for(int i=1;i<=n;++i)
        a[i].y=lower_bound(temp+1,temp+1+m,a[i].y)-temp;
    return m;
}

ll add[N<<2],max_[N<<2];

void up(int rt,int l,int r)
{
    max_[rt]=max(max_[rt<<1],max_[rt<<1|1]);
}

void dn(int rt)
{
    if(add[rt]!=0)
    {
        add[ls]+=add[rt];
        add[rs]+=add[rt];
        max_[rt<<1]+=add[rt];
        max_[rt<<1|1]+=add[rt];
        add[rt]=0;
    }
}

void Build(int l,int r,int rt)
{
    max_[rt]=0;
    add[rt]=0;
    if(l==r)
    {
        return;
    }
    int mid=(l+r)>>1;

    Build(l,mid,rt<<1);
    Build(mid+1,r,rt<<1|1);
    up(rt,l,r);
}

void update_dot(int pos,ll V,int l,int r,int rt)
{
    if(l==r)
    {
        max_[rt]=V;
        return;
    }

    int mid=(l+r)>>1;
    dn(rt);
    if(pos<=mid)
        update_dot(pos,V,l,mid,rt<<1);
    else
        update_dot(pos,V,mid+1,r,rt<<1|1);
    up(rt,l,r);
}
void update_qu(int L,int R,int V,int l,int r,int rt)
{
    if(L>R)return;
    if(L<=l&&r<=R)
    {
        max_[rt]+=V;
        add[rt]+=V;
        return;
    }

    int mid=(l+r)>>1;
    dn(rt);
    if(L<=mid)
        update_qu(L,R,V,l,mid,rt<<1);
    if(R>mid)
        update_qu(L,R,V,mid+1,r,rt<<1|1);
    up(rt,l,r);
}
ll Query(int L,int R,int l,int r,int rt)
{
    if(L<=l&&r<=R)    return max_[rt];
    int mid=(l+r)>>1;
    ll res=0;
    dn(rt);
    if(L<=mid)
        res=max(res,Query(L,R,l,mid,rt<<1));
    if(R>mid)
        res=max(res,Query(L,R,mid+1,r,rt<<1|1));
    return res;
}
void check(int pos,int l,int r,int rt)
{
    if(l==r)
    {
        cout<<max_[rt]<<' ';
        return ;
    }
    int mid=(l+r)>>1;

    dn(rt);
    if(pos<=mid)
        check(pos,l,mid,rt<<1);
    else
        check(pos,mid+1,r,rt<<1|1);
}

signed main()
{
    int n;
    while(~sc("%d",&n))
    {
        for(int i=1;i<=n;++i)
        {
            int q,w,e,r;
            sc("%d%d%d%d",&q,&w,&e,&r);
            a[i]={q,w,e,r};
        }
        int tot=LS(n);
        sort(a+1,a+1+n);
        Build(0,tot,1);
        for(int i=1;i<=n;++i)
        {
        //    a[i].y++;
            ll max__=Query(0,a[i].y+1,0,tot,1);
            update_dot(a[i].y,max__+a[i].vb,0,tot,1);
            update_qu(0,a[i].y-1,a[i].va,0,tot,1);
            update_qu(a[i].y+1,tot,a[i].vb,0,tot,1);
    //        for(int j=0;j<=tot;++j)
    //            check(j,0,tot,1);
    //        cout<<endl;
        }
        pr("%lld
",max_[1]);
    }
    return 0;
}

/**************************************************************************************/

int maxx(int a,int b)
{
    return a>b?a:b;
}

void swapp(int &a,int &b)
{
    a^=b^=a^=b;
}

int lowbit(int n)
{
    return n&(-n);
}

int Del_bit_1(int n)
{
    return n&(n-1);
}

int abss(int a)
{
    return a>0?a:-a;
}

double fabss(double a)
{
    return a>0?a:-a;
}

int minn(int a,int b)
{
    return a<b?a:b;
}