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=(<<);
 const double ESP=1e-;
 const int mod=(int)1e9+;
 const int N=(int)1e6+;

 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=;
     for(int i=;i<=n;++i)
         temp[++m]=a[i].y;
     sort(temp+,temp++m);
     m=unique(temp+,temp++m)-temp-;
     for(int i=;i<=n;++i)
         a[i].y=lower_bound(temp+,temp++m,a[i].y)-temp;
     return m;
 }

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

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

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

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

     Build(l,mid,rt<<);
     Build(mid+,r,rt<<|);
     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)>>;
     dn(rt);
     if(pos<=mid)
         update_dot(pos,V,l,mid,rt<<);
     else
         update_dot(pos,V,mid+,r,rt<<|);
     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)>>;
     dn(rt);
     if(L<=mid)
         update_qu(L,R,V,l,mid,rt<<);
     if(R>mid)
         update_qu(L,R,V,mid+,r,rt<<|);
     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)>>;
     ll res=;
     dn(rt);
     if(L<=mid)
         res=max(res,Query(L,R,l,mid,rt<<));
     if(R>mid)
         res=max(res,Query(L,R,mid+,r,rt<<|));
     return res;
 }
 void check(int pos,int l,int r,int rt)
 {
     if(l==r)
     {
         cout<<max_[rt]<<' ';
         return ;
     }
     int mid=(l+r)>>;

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

 signed main()
 {
     int n;
     while(~sc("%d",&n))
     {
         for(int i=;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+,a++n);
         Build(,tot,);
         for(int i=;i<=n;++i)
         {
         //    a[i].y++;
             ll max__=Query(,a[i].y+,,tot,);
             update_dot(a[i].y,max__+a[i].vb,,tot,);
             update_qu(,a[i].y-,a[i].va,,tot,);
             update_qu(a[i].y+,tot,a[i].vb,,tot,);
     //        for(int j=0;j<=tot;++j)
     //            check(j,0,tot,1);
     //        cout<<endl;
         }
         pr("%lld\n",max_[]);
     }
     return ;
 }

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

 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-);
 }

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

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

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