BZOJ3615 : MSS
将所有点按横坐标排序,那么最长上升子序列和最长下降子序列里必有一个长度不小于$\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的更多相关文章
- 关于TCP中的MSS
MSS 是TCP选项中最经常出现,也是最早出现的选项.MSS选项占4byte.MSS是每一个TCP报文段中数据字段的最大长度,注意:只是数据部分的字段,不包括TCP的头部.TCP在三次握手中,每一方都 ...
- TCP 协议中MSS的理解
在介绍MSS之前我们必须要理解下面的几个重要的概念.MTU: Maxitum Transmission Unit 最大传输单元MSS: Maxitum Segment Size 最大分段大小PPPoE ...
- tcp协议中mss的理解
在介绍MSS之前我们必须要理解下面的几个重要的概念.<blockquote>MTU: Maxitum Transmission Unit 最大传输单元MSS: Maxitum Segmen ...
- IP分片和TCP分片 MTU和MSS(转)
IP分片和TCP分片 MTU和MSS(转) 访问原文:http://blog.csdn.net/keyouan2008/article/details/5843388 1,MTU(Maximum Tr ...
- MTU & MSS 详解记录(转)
先学习理解一下帧的封装格式: 需要注意的是,区别两种帧封装格式:802标准帧和以太网帧 1,在802标准定义的帧格式中,长度字段是指它后续数据的字节长度,但不包括C R C检验 ...
- tcp netstat用法 TIME_WAIT状态解析 MTU以及MSS
带着问题写博客 问题1:使用netstat查看有源TCP连接的状态时,经常会看到established状态,那么还有哪些状态,这些状态是如何变化的呢? 问题2:TIME_WAIT状态存在的必要? 问题 ...
- TCP层的分段和IP层的分片之间的关系 & MTU和MSS之间的关系 (转载)
首先说明:数据报的分段和分片确实发生,分段发生在传输层,分片发生在网络层.但是对于分段来说,这是经常发生在UDP传输层协议上的情况,对于传输层使用TCP协议的通道来说,这种事情很少发生. 1,MTU( ...
- MTU,MSS基本概念
传输层: 对于UDP协议来说,整个包的最大长度为65535,其中包头长度是65535-20=65515: 对于TCP协议来说,整个包的最大长度是由最大传输大小(MSS,Maxitum Segment ...
- 在linux内核中修改TCP MSS值
MTU: Maxitum Transmission Unit 最大传输单元 MSS: Maxitum Segment Size 最大分段大小 MSS最大传输大小的缩写,是TCP协议里面的一个概念.MS ...
随机推荐
- 开启后台 Service 闪退
04-29 15:36:23.395: E/ActivityThread(15275): Performing stop of activity that is not resumed: {com.e ...
- mybatis 学习!
参考链接 http://www.mybatis.org/spring/zh/mappers.html http://www.cnblogs.com/fangjian0423/p/spring-myba ...
- Eclipse的各种快捷键
参考他人总结好的 Eclipse快捷键大全(转载) Ctrl+1 快速修复(最经典的快捷键,就不用多说了) ---->例如:add unimplement methodCtrl+D: ...
- PAL/NTSC 制电视广播技术有关知识--FPGA
1.PAL和NTSC的区别 常见的电视信号制式是PAL和NTSC,另外还有SECAM等. NTSC即正交平衡调幅制,PAL为逐行倒像正交平衡调幅制. (1)PAL电视标准 PAL电视标准,每秒25帧 ...
- localStorage的使用
HTML5中提供了localStorage对象可以将数据长期保存在客户端,直到人为清除. localStorage提供了几个方法: 1.存储:localStorage.setItem(key,valu ...
- php PDO:数据访问抽象层
<!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN" "http://www.w3.org/ ...
- Jquery.Datatables 结合时间段查询,daterangepicker实现Datatables表格带参数查询
参考:http://datatables.club/example/user_share/send_extra_param.html 下载地址:http://pan.baidu.com/s/1 ...
- 17.观察者模式(Observer Pattern)
using System; using System.Collections.Generic; namespace ConsoleApplication10 { /// <summary> ...
- Delphi中的变量作用域简介
1.在函数或过程中声明的变量是局部变量,只能在函数或过程中使用 2.在单元的implementation部分声明的变量是单元级的(又称模块级.窗体级)变量,该单元的所有函数与过程都能使用,其他引用该单 ...
- tomcat安装服务和内存参数设置
第一:安装服务 在dos窗口进入到tomcat的bin目录下,通过如下命令即可将tomcat安装成服务 service.bat install Tomcat2 其中Tomcat2是服务的名称 如果启动 ...