紫薯例题+1。

题意:给你一个长度为n(n<=200000)的序列a[n],求删除一个连续子序列后的可能的最长连续上升子序列的长度。

首先对序列进行分段,每一段连续的子序列的元素递增,设L[i]为下标i对应的元素向左能延伸到的最大长度(在同一段内),R[i]为向右能延伸到的最大长度,则问题转化成了对于每个下标i,找到在它前面的下标j中a[j]<a[i]且L[j]最大的j,然后用R[i]+L[j]去更新ans。

第一种方法是用一个二元组(x,y)表示大小为x的元素所对应的L的值,用一个set保存所有这样的二元组,则任意两个二元组不存在包含关系,即满足x递增且y也递增。于是对原序列从左往右扫一遍,每扫到一个元素,就可以从set中二分找到最大的x使得a[i]>x,用R[i]+y去更新ans即可。然后向set中插入(a[i],L[i]),并根据set中左右两边的二元组的值来确定是否需要保留它,并清掉一些无用的二元组。

这里用到了一点小技巧,就是向set中插入两个虚的元素作为左右边界,这样就不用特判了,方便了很多。

 #include<bits/stdc++.h>

 using namespace std;
typedef long long ll;
const int N=2e5+,inf=0x3f3f3f3f;
int a[N],L[N],R[N],n,ans;
struct nd {
int x,y;
bool operator<(const nd& b)const {return x!=b.x?x<b.x:y<b.y;}
};
set<nd> st; int main() {
int T;
for(scanf("%d",&T); T--;) {
st.clear();
st.insert({,});
st.insert({inf,inf});
ans=;
scanf("%d",&n);
for(int i=; i<n; ++i)scanf("%d",&a[i]);
for(int l=,r; l<n; l=r) {
for(r=l+; r<n&&a[r]>a[r-]; ++r);
for(int j=l; j<r; ++j)L[j]=j-l+,R[j]=r-j;
}
for(int i=; i<n; ++i) {
set<nd>::iterator it,it2;
it=st.lower_bound({a[i],-}),it--;
ans=max(ans,it->y+R[i]);
st.insert({a[i],L[i]});
it=st.find({a[i],L[i]});
while() {
it2=it,it2++;
if(it2->y<=it->y)st.erase(it2);
else break;
}
it2=it,it2--;
if(it2->y>=it->y)st.erase(it);
}
printf("%d\n",ans);
}
return ;
}

第二种方法是每扫到一个i,直接查询a[j]<a[i]的元素中最大的L,怎么查呢?用树状数组(BIT)呗!树状数组不仅能动态维护前缀和,还能动态维护所有小于k的元素中的最大值哦!(当然,前提是没有将元素减小的操作)将原序列离散化后的值作为下标插入BIT中就行了,注意下标要从1开始。

 #include<bits/stdc++.h>

 using namespace std;
typedef long long ll;
const int N=2e5+,inf=0x3f3f3f3f;
int a[N],b[N],c[N],L[N],R[N],n,nn,ans;
int lowbit(int x) {return x&-x;}
void add(int u,int x) {for(; u<=nn; u+=lowbit(u))c[u]=max(c[u],x);}
int get(int u) {int ret=; for(; u; u-=lowbit(u))ret=max(ret,c[u]); return ret;} int main() {
int T;
for(scanf("%d",&T); T--;) {
ans=;
scanf("%d",&n);
for(int i=; i<n; ++i)scanf("%d",&a[i]),b[i]=a[i];
sort(b,b+n);
nn=unique(b,b+n)-b;
for(int i=; i<n; ++i)a[i]=lower_bound(b,b+nn,a[i])-b+;
for(int l=,r; l<n; l=r) {
for(r=l+; r<n&&a[r]>a[r-]; ++r);
for(int j=l; j<r; ++j)L[j]=j-l+,R[j]=r-j;
}
memset(c,,sizeof c);
for(int i=; i<n; ++i) {
ans=max(ans,R[i]+get(a[i]-));
add(a[i],L[i]);
}
printf("%d\n",ans);
}
return ;
}

第三种方法是最简洁最高效但也是最抽象的一种方法,就是仿照求LIS的算法,维护一个单调的数组c,记录所有L值为i的元素中最小的a[i](和BIT正好相反),每扫到一个元素a[i],在c中通过二分查找到比a[i]小的最大的元素,它的下标值就是最大的L,然后令c[L[i]]=min(c[L[i]],a[i])即可。其实这种方法和方法一原理是相同的,由于这道题比较特殊——L值是连续的!因为假设出现了某个L值,则L-1一定要出现在它前面,以此类推,因此不会出现下标空缺的现象,完全可以用数组来代替set。

 #include<bits/stdc++.h>

 using namespace std;
typedef long long ll;
const int N=2e5+,inf=0x3f3f3f3f;
int a[N],c[N],L[N],R[N],n,ans; int main() {
int T;
for(scanf("%d",&T); T--;) {
ans=;
scanf("%d",&n);
for(int i=; i<n; ++i)scanf("%d",&a[i]);
for(int l=,r; l<n; l=r) {
for(r=l+; r<n&&a[r]>a[r-]; ++r);
for(int j=l; j<r; ++j)L[j]=j-l+,R[j]=r-j;
}
memset(c,inf,sizeof c);
c[]=;
for(int i=; i<n; ++i) {
ans=max(ans,int(R[i]+lower_bound(c,c+n,a[i])-c-));
c[L[i]]=min(c[L[i]],a[i]);
}
printf("%d\n",ans);
}
return ;
}

虽然这种题已经做过n遍了,但做完之后还是有很大的感悟,里面蕴含的思想真的不容小觑啊。

UVA - 1471 Defense Lines (set/bit/lis)的更多相关文章

  1. UVA 1471 Defense Lines 防线 (LIS变形)

    给一个长度为n的序列,要求删除一个连续子序列,使剩下的序列有一个长度最大的连续递增子序列. 最简单的想法是枚举起点j和终点i,然后数一数,分别向前或向后能延伸的最长长度,记为g(i)和f(i).可以先 ...

  2. UVA - 1471 Defense Lines 树状数组/二分

                                  Defense Lines After the last war devastated your country, you - as the ...

  3. Uva 1471 Defense Lines(LIS变形)

    题意: 给你一个数组,让你删除一个连续的子序列,使得剩下的序列中有最长上升子序列, 求出这个长度. 题解: 预处理:先求一个last[i],以a[i]为开始的合法最长上升子序列的长度.再求一个pre[ ...

  4. UVa 1471 Defense Lines - 线段树 - 离散化

    题意是说给一个序列,删掉其中一段连续的子序列(貌似可以为空),使得新的序列中最长的连续递增子序列最长. 网上似乎最多的做法是二分查找优化,然而不会,只会值域线段树和离散化... 先预处理出所有的点所能 ...

  5. uva 1471 Defense Lines

    题意: 给一个长度为n(n <= 200000) 的序列,你删除一段连续的子序列,使得剩下的序列拼接起来,有一个最长的连续递增子序列 分析: 就是最长上升子序列的变形.需要加一个类似二分搜索就好 ...

  6. UVa 1471 Defense Lines (二分+set优化)

    题意:给定一个序列,然后让你删除一段连续的序列,使得剩下的序列中连续递增子序列最长. 析:如果暴力枚举那么时间复杂度肯定受不了,我们可以先进行预处理,f[i] 表示以 i 结尾的连续最长序列,g[i] ...

  7. uva 1471 defence lines——yhx

    After the last war devastated your country, you - as the king of the land of Ardenia - decided it wa ...

  8. 1471 - Defense Lines

    After the last war devastated your country, you - as the king of the land of Ardenia - decided it wa ...

  9. 【二分】Defense Lines

    [UVa1471] Defense Lines 算法入门经典第8章8-8 (P242) 题目大意:将一个序列删去一个连续子序列,问最长的严格上升子序列 (N<=200000) 试题分析:算法1: ...

随机推荐

  1. C++中引用编译过的C代码为什么要用“extern c”

    函数经过编译系统的翻译成汇编,函数名对应着汇编标号.  因为C编译函数名与得到的汇编代号基本一样,如:fun()=>_fun, main=>_main  但是C++中函数名与得到的汇编代号 ...

  2. first application

    <!DOCTYPE html> <html> <head> <title>Create a Map</title> <meta htt ...

  3. [原创]spring及springmvc精简版--AOP

    接上一篇:[原创]spring及springmvc精简版--IOC 理解AOP.java是一种面向对象的语言.而AOP是面向切面,在我看来是面向逻辑或者业务编程,它是对一组逻辑的抽象和分配. 经典例子 ...

  4. Python自然语言处理-系列一

    一:python基础,自然语言概念 from nltk.book import * 1,text1.concordance("monstrous")      用语索引 2,tex ...

  5. Adding Flexcan driver support on Kernel

    Adding Flexcan driver support on Kernel On kernel menuconfig, add the following items: [*] Networkin ...

  6. require.js和sea.js的区别

    下面为大家讲解一下require.js和sea.js的区别.纯属个人意见,不喜勿喷. 首先原理上的区别 sea.js遵循CMD规范.书写方式类似node.js的书写模板代码.依赖的自动加载,配置的简洁 ...

  7. Hadoop相关知识整理系列之一:HBase基本架构及原理

    1. HBase框架简单介绍 HBase是一个分布式的.面向列的开源数据库,它不同于一般的关系数据库,是一个适合于非结构化数据存储的数据库.另一个不同的是HBase基于列的而不是基于行的模式.HBas ...

  8. Linux自定义别名alias重启失效问题

    Linux上的别名功能非常方便,例如ll可以显示文件列表的长信息,但是却不是以human能读懂的方式显示,所以我尝试直接在命令行中自定义一个别名: alisa lk='ls -lh' 然后lk就能正常 ...

  9. IDEA配置文件的保存目录

    IntelliJ IDEA 的实时代码模板保存在 /templates 目录下,其他系统目录位置如下:(因为目录名在各个系统上是一致的,建议用硬盘搜索工具搜索即可) Windows: . Linux: ...

  10. 剑指Offer——圆圈中最后剩下的数(约瑟夫环)

    Question 每年六一儿童节,牛客都会准备一些小礼物去看望孤儿院的小朋友,今年亦是如此.HF作为牛客的资深元老,自然也准备了一些小游戏.其中,有个游戏是这样的:首先,让小朋友们围成一个大圈.然后, ...