UVA10534 Wavio Sequence

链接

https://uva.onlinejudge.org/index.php?option=com_onlinejudge&Itemid=8&page=show_problem&problem=1475

题解

正反各做一遍LIS
设$f_i$为以第$i$个元素结尾的最长严格上升子序列
用$list_i$存储$f$值为$i$的最小的数字，那么这个$list$是严格单调递增的
反证法：
假设$list$有两个元素，分别对应原数组中的$a_i,a_j$，$i<j$，
假设第$i$个元素是用第$k$个转移的，且$a_j<a_i$，如果用$k$去转移$j$的话，就能使得$f_j=f_i$，根据$list$的定义，$i$会被替换成$j$
若$a_j=a_i$，那么必定存在$i<k<j$，使得$a_k<a_j$，且$f_i\le f_k$，那么$i$就应该被替换为$k$
综上，这个$list$是严格单调递增的
所以就可以二分一下，在$O(log_2n)$时间内求得$f_i$

代码

//DP、二分
#include <cstdio>
#include <algorithm>
#include <cstring>
#include <cctype>
#define maxn 10010
#define cl(x) memset(x,0,sizeof(x))
#define inf 0x3f3f3f3f
using namespace std;
int n, f1[maxn], f2[maxn], a[maxn], list[maxn];
int read(int x=0)
{
    char c, f=1;
    for(c=getchar();!isdigit(c) and c^-1;c=getchar())if(c=='-')f=-1;
    for(;isdigit(c);c=getchar())x=(x<<1)+(x<<3)+c-48;
    return x;
}
void calc(int f[])
{
    int i, j;
    *list=0;
    for(i=f==f1?1:n;f==f1?i<=n:i;f==f1?i++:i--)
    {
        f[i]=lower_bound(list+1,list+*list+1,a[i])-list;
        if(f[i]>*list)list[++*list]=a[i];
        else list[f[i]]=min(list[f[i]],a[i]);
    }
}
void work()
{
    int i, ans=-inf;
    for(i=1;i<=n;i++)a[i]=read();
    calc(f1);
    calc(f2);
    for(i=1;i<=n;i++)ans=max(ans,min(f1[i],f2[i])*2-1);
    printf("%lld\n",ans);
}
int main()
{
    for(n=read();n;n=read())work();
    return 0;
}