例题8-8 防线 UVa1471

1.题目描述：点击打开链接

2.解题思路：本题要求在一个长度为n的序列中，删除一段连续子序列，使得剩下的序列有一个长度最大的连续递增子序列。按照以往的经验，可以事先用f[i]计算出从i出发的最长连续递增子序列的长度，用g[i]表示以i结尾的最长连续递增子序列的长度，然后枚举j,i(j<=i)即可，然而这么做的时间复杂度O(N^2)的，无法通过本题，还需要进一步优化。

本题的优化环节就在于已知一个i，如何快速找到一个j。由于我们希望得到的连续递增子序列尽可能的长，那么我们希望在A[j]<A[i]时，还可以让g[j]尽可能的大。考虑这样一种情形，A[j]<A[k]且g[j]>g[k]，这意思就是说：在序列...A[j]...A[k]...A[i]....中，以A[j]结尾的连续递增子序列长度比以A[k]结尾的连续递增子序列长度更长一点，那么此时我们只需要保留A[j]即可。或者第二种情况，g[i]==g[k]但A[j]<A[k]，那么我们保留A[j]一定不会丢失最优解。这就非常像经典的LIS问题的O(NlogN)算法的做法。仿照LIS问题的解法，我们不妨令s[i]表示g==i时的最小的A[j]，这样不难得到s[1]<=s[2]<=...<=s[n]，可以用lower_bound找到s[p]>=A[i]的第一个下标p，那么p-1就是s[p‘]<A[i]的最后一个下标p’，此时的p其实就代表了g值（因为每个g值都对应一个最小的A[j]），所以此时的答案就是p-1+f[i]，如果此时的A[i]<s[g[i]]，那么更新s[g[i]]=A[i]。

注意：千万不要写成s[p]=A[i]，因为我们更新的是g[i]的最小值，p只是为了得到p-1。

3.代码：

//#pragma comment(linker, "/STACK:1024000000,1024000000")
#include<iostream>
#include<algorithm>
#include<cassert>
#include<string>
#include<sstream>
#include<set>
#include<bitset>
#include<vector>
#include<stack>
#include<map>
#include<queue>
#include<deque>
#include<cstdlib>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<ctime>
#include<cctype>
#include<functional>
using namespace std;

#define me(s)  memset(s,0,sizeof(s))
#define rep(i,n) for(int i=0;i<(n);i++)
typedef long long ll;
typedef unsigned int uint;
typedef unsigned long long ull;
typedef pair <int, int> P;


const int N=200110;
int n,a[N],f[N],g[N],s[N];

int main()
{
    int T;
    for(scanf("%d",&T);T--;)
    {
        scanf("%d",&n);
        for(int i=1;i<=n;i++)
        {
            scanf("%d",&a[i]);
            if(i==1)g[1]=1;
            else if(a[i]>a[i-1])g[i]=g[i-1]+1;//计算以i结尾的最长子序列的长度
            else g[i]=1;
        }
        for(int i=n;i>=1;i--)
        {
            if(i==n)f[i]=1;
            else if(a[i]<a[i+1])f[i]=f[i+1]+1;//计算以i开始的最长子序列的长度
            else f[i]=1;
        }
        for(int i=1;i<=n;i++)
            s[i]=(1e9)+10;
        s[1]=a[1];
        int res=1;
        for(int i=2;i<=n;i++)
        {
            int p=lower_bound(s+1,s+n+1,a[i])-s;
            res=max(res,p-1+f[i]);  //p-1就是A[j]<A[i]且g[j]最大的那个最大值
            if(a[i]<s[g[i]])s[g[i]]=a[i]; //更新g[i]对应的最小的A[i]
        }
        printf("%d\n",res);
    }
}