将所有点按横坐标排序,那么最长上升子序列和最长下降子序列里必有一个长度不小于$\sqrt{n}$。

因为如果最长上升子序列很短,那么根据Dilworth定理,它的反链长度与它成反比。

因此可以将$n$个点划分为不超过$\sqrt{n}$条链或者反链,每一个部分里都相当于只有一维坐标。

对于每一个部分,建立线段树,支持合并、分裂以及打标记即可。

时间复杂度$O(q\sqrt{n}\log n)$。

#include<cstdio>
#include<algorithm>
using namespace std;
typedef long long ll;
const int N=50010,M=100010,E=5000000;
const ll inf=1LL<<60;
int n,m,cnt,cset,mx,i,j,b[N],d[N],q[M][4],rem,dir,del[N],f[N],bit[N];ll ans[M][3];
struct P{int x,y,w,p;}a[N],c[N];
inline bool cmp(const P&a,const P&b){return a.x<b.x;}
inline int lower(int x){
int l=1,r=n,mid,t;
while(l<=r)if(b[mid=(l+r)>>1]<=x)l=(t=mid)+1;else r=mid-1;
return t;
}
inline void ins(int x,int p){for(;x<=n;x+=x&-x)if(bit[x]<p)bit[x]=p;}
inline int ask(int x){int t=0;for(;x;x-=x&-x)if(t<bit[x])t=bit[x];return t;}
int T[M*2+N],tot,l[E],r[E],v[E];ll mi[E],ma[E],sum[E],tag[E];
inline void tag1(int x,ll p){if(x)mi[x]+=p,ma[x]+=p,sum[x]+=p*v[x],tag[x]+=p;}
inline void pb(int x){if(tag[x])tag1(l[x],tag[x]),tag1(r[x],tag[x]),tag[x]=0;}
inline void up(int x){
v[x]=v[l[x]]+v[r[x]];
mi[x]=min(mi[l[x]],mi[r[x]]);
ma[x]=max(ma[l[x]],ma[r[x]]);
sum[x]=sum[l[x]]+sum[r[x]];
}
int build(int a,int b,int c,int p){
int x=++tot;
v[x]=1,mi[x]=ma[x]=sum[x]=p;
if(a==b)return x;
int mid=(a+b)>>1;
if(c<=mid)l[x]=build(a,mid,c,p);else r[x]=build(mid+1,b,c,p);
return x;
}
int merge(int x,int y,int a,int b){
if(!x||!y)return x+y;
if(a==b){
v[x]+=v[y];
sum[x]+=sum[y];
mi[x]=min(mi[x],mi[y]);
ma[x]=max(ma[x],ma[y]);
return x;
}
int mid=(a+b)>>1;
pb(x),pb(y);
l[x]=merge(l[x],l[y],a,mid);
r[x]=merge(r[x],r[y],mid+1,b);
up(x);
return x;
}
int split(int x,int a,int b,int c,int d){
if(!x)return 0;
if(c<=a&&b<=d)return x;
pb(x);
int y=++tot,mid=(a+b)>>1;
if(c<=mid)l[y]=split(l[x],a,mid,c,d);
if(d>mid)r[y]=split(r[x],mid+1,b,c,d);
up(x),up(y);
return y;
}
inline int findx(int x){
int l=1,r=cnt,mid,t=0;
while(l<=r)if(c[mid=(l+r)>>1].x<=x)l=(t=mid)+1;else r=mid-1;
return t;
}
inline int findy0(int x){
int l=1,r=cnt,mid,t=0;
while(l<=r)if(c[mid=(l+r)>>1].y<=x)l=(t=mid)+1;else r=mid-1;
return t;
}
inline int findy1(int x){
int l=1,r=cnt,mid,t=cnt+1;
while(l<=r)if(c[mid=(l+r)>>1].y<=x)r=(t=mid)-1;else l=mid+1;
return t;
}
inline void solve(){
int i;
for(cset=n,i=1;i<=n;i++)T[i]=0;
for(i=1;i<=cnt;i++)T[c[i].p]=build(1,cnt,i,c[i].w);
for(i=0;i<m;i++){
int x=q[i][1],y=q[i][2],z=q[i][3];
if(!q[i][0]){
T[x]=merge(T[x],T[y],1,cnt);
T[y]=0;
}
if(q[i][0]==1){
cset+=2;
T[cset-1]=T[cset]=0;
if(!y){
int o=findx(z);
if(o)T[cset-1]=split(T[x],1,cnt,1,o);
if(o<cnt)T[cset]=split(T[x],1,cnt,o+1,cnt);
}else if(!dir){
int o=findy0(z);
if(o)T[cset-1]=split(T[x],1,cnt,1,o);
if(o<cnt)T[cset]=split(T[x],1,cnt,o+1,cnt);
}else{
int o=findy1(z);
if(o<=cnt)T[cset-1]=split(T[x],1,cnt,o,cnt);
if(o>1)T[cset]=split(T[x],1,cnt,1,o-1);
}
T[x]=0;
}
if(q[i][0]==2){
ans[i][0]=max(ans[i][0],ma[T[x]]);
ans[i][1]=min(ans[i][1],mi[T[x]]);
ans[i][2]+=sum[T[x]];
}
if(q[i][0]==3)tag1(T[x],y);
}
for(i=1;i<=tot;i++)l[i]=r[i]=v[i]=mi[i]=ma[i]=sum[i]=tag[i]=0;
tot=0;
}
int main(){
mi[0]=inf,ma[0]=-inf;
scanf("%d",&n);
for(i=1;i<=n;i++)scanf("%d%d%d",&a[i].x,&a[i].y,&a[i].w),a[i].p=i;
scanf("%d",&m);
for(i=0;i<m;i++){
char op[9];
scanf("%s",op);
if(op[0]=='M')scanf("%d%d",&q[i][1],&q[i][2]);
if(op[0]=='S')q[i][0]=1,scanf("%d%d%d",&q[i][1],&q[i][2],&q[i][3]);
if(op[0]=='Q'){
q[i][0]=2,scanf("%d",&q[i][1]);
ans[i][0]=-inf,ans[i][1]=inf;
}
if(op[0]=='A')q[i][0]=3,scanf("%d%d",&q[i][1],&q[i][2]);
}
rem=n;
for(i=1;i<=n;i++)b[i]=a[i].y;
sort(b+1,b+n+1);
sort(a+1,a+n+1,cmp);
for(i=1;i<=n;i++)d[i]=lower(a[i].y);
while(rem){
for(i=1;i<=n;i++)bit[i]=0;
for(i=1;i<=n;i++)if(!del[i])ins(d[i],f[i]=ask(d[i])+1);
mx=ask(n);
if(mx>=rem/mx){
rem-=mx;
dir=0;
for(i=n;i;i--)if(!del[i]&&f[i]==mx)break;
j=i;
c[cnt=mx]=a[j],del[j]=1;
for(i=j-1;i;i--)if(!del[i]&&f[i]+1==f[j]&&d[i]<d[j])c[--mx]=a[i],del[j=i]=1;
solve();
continue;
}
for(i=1;i<=n;i++)bit[i]=0;
for(i=1;i<=n;i++)if(!del[i])ins(n-d[i]+1,f[i]=ask(n-d[i]+1)+1);
mx=ask(n);
rem-=mx;
dir=1;
for(i=n;i;i--)if(!del[i]&&f[i]==mx)break;
j=i;
c[cnt=mx]=a[j],del[j]=1;
for(i=j-1;i;i--)if(!del[i]&&f[i]+1==f[j]&&d[i]>d[j])c[--mx]=a[i],del[j=i]=1;
solve();
}
for(i=0;i<m;i++)if(q[i][0]==2)printf("%lld %lld %lld\n",ans[i][0],ans[i][1],ans[i][2]);
return 0;
}

  

BZOJ3615 : MSS的更多相关文章

  1. 关于TCP中的MSS

    MSS 是TCP选项中最经常出现,也是最早出现的选项.MSS选项占4byte.MSS是每一个TCP报文段中数据字段的最大长度,注意:只是数据部分的字段,不包括TCP的头部.TCP在三次握手中,每一方都 ...

  2. TCP 协议中MSS的理解

    在介绍MSS之前我们必须要理解下面的几个重要的概念.MTU: Maxitum Transmission Unit 最大传输单元MSS: Maxitum Segment Size 最大分段大小PPPoE ...

  3. tcp协议中mss的理解

    在介绍MSS之前我们必须要理解下面的几个重要的概念.<blockquote>MTU: Maxitum Transmission Unit 最大传输单元MSS: Maxitum Segmen ...

  4. IP分片和TCP分片 MTU和MSS(转)

    IP分片和TCP分片 MTU和MSS(转) 访问原文:http://blog.csdn.net/keyouan2008/article/details/5843388 1,MTU(Maximum Tr ...

  5. MTU & MSS 详解记录(转)

              先学习理解一下帧的封装格式:   需要注意的是,区别两种帧封装格式:802标准帧和以太网帧 1,在802标准定义的帧格式中,长度字段是指它后续数据的字节长度,但不包括C R C检验 ...

  6. tcp netstat用法 TIME_WAIT状态解析 MTU以及MSS

    带着问题写博客 问题1:使用netstat查看有源TCP连接的状态时,经常会看到established状态,那么还有哪些状态,这些状态是如何变化的呢? 问题2:TIME_WAIT状态存在的必要? 问题 ...

  7. TCP层的分段和IP层的分片之间的关系 & MTU和MSS之间的关系 (转载)

    首先说明:数据报的分段和分片确实发生,分段发生在传输层,分片发生在网络层.但是对于分段来说,这是经常发生在UDP传输层协议上的情况,对于传输层使用TCP协议的通道来说,这种事情很少发生. 1,MTU( ...

  8. MTU,MSS基本概念

    传输层: 对于UDP协议来说,整个包的最大长度为65535,其中包头长度是65535-20=65515: 对于TCP协议来说,整个包的最大长度是由最大传输大小(MSS,Maxitum Segment ...

  9. 在linux内核中修改TCP MSS值

    MTU: Maxitum Transmission Unit 最大传输单元 MSS: Maxitum Segment Size 最大分段大小 MSS最大传输大小的缩写,是TCP协议里面的一个概念.MS ...

随机推荐

  1. 001课-java_web开发入门

    一.Tomcat服务器常见启动问题:(1).Java_home环境变量,由于tomcat服务器的bin目录中的一些jar文件必须使用到java类库,所以必须先配置Java_home环境变量.(2).端 ...

  2. Eclipse的各种快捷键

     参考他人总结好的 Eclipse快捷键大全(转载) Ctrl+1 快速修复(最经典的快捷键,就不用多说了)     ---->例如:add unimplement methodCtrl+D:  ...

  3. 三、jQuery--jQuery基础--jQuery基础课程--第1章 初识jQuery

    环境搭建 搭建一个jQuery的开发环境非常方便,可以通过下列几个步骤进行. 下载jQuery文件库 在jQuery的官方网站(http://jquery.com)中,下载最新版本的jQuery文件库 ...

  4. python以指定宽度及格式输出字符串

  5. java中常用的工具类(一)

    我们java程序员在开发项目的是常常会用到一些工具类.今天我汇总了一下java中常用的工具方法.大家可以在项目中使用.可以收藏!加入IT江湖官方群:383126909 我们一起成长 一.String工 ...

  6. 在MAVEN仓库中添加ORACLE JDBC驱动

    本文转载自 http://www.cnblogs.com/leiOOlei/archive/2013/10/21/3380568.html 因为已经是第二次遇到,所以COPY过来,怕以后别人的BLOG ...

  7. LSM-Tree (BigTable 的理论模型)(转)

    Google的BigTable架构在分布式结构化存储方面大名鼎鼎,其中的MergeDump模型在读写之间找到了一个较好的平衡点,很好的解决了web scale数据的读写问题. MergeDump的理论 ...

  8. PHP不同域名cookie共享(单点登录实现原理)

    PHP使用P3P完成COOKIE跨域操作实际实用中,类似的需求有,比如说我们有两个域名,我们想实现在一个域名登录后,能自动完成另一个域名的登录,也就是单点登录(SSO)功能.为了测试的方便,先编辑ho ...

  9. web.xml中同一servlet/filter配置多个url-pattern

    转自:http://blog.sina.com.cn/s/blog_4c2c2a0c0100dh67.html 若你的servlet要多个地址,或你的filter需要过滤不同的url如有*.jsp,* ...

  10. APP设计尺寸规范大全,APP界面设计新手教程【官方版】(转)

    正值25学堂一周年之际,同时站长和APP设计同仁们在群里(APP界面设计 UI设计交流群,APP界面设计⑥群 APPUI设计③群58946771 APP设计资源⑤群 386032923欢迎大家加入交流 ...