BZOJ 2090 [Poi2010]Monotonicity 2 DP+线段树

题意:
给出N个正整数a[1..N]，再给出K个关系符号（>、<或=）s[1..k]。
选出一个长度为L的子序列（不要求连续），要求这个子序列的第i项和第i+1项的的大小关系为s[(i-1)mod K+1]。
求出L的最大值。
解析:
一眼考虑DP，不过这个DP我的确刚开始认为这是没有什么正确性的，然后我就不会辣，但是大家都说这题就是DP..
并且好像自己一顿想要构造也没有构造出来不合法的数据?
那我就只能回到自己刚开始的DP思路了。
设F[i]表示到第i位拿出来的最长序列长度。
所以下一个符号显然是确定的。
然后等于号的话，我们直接开一个数组记录哪一个值上一次出现的可行位置即可。
小于号，丢到线段树里面，线段树对权值开，记录一段权值区间最大的已出现的f值，查询直接拿最大就好了。
大于号同理，再开一个线段树丢进去即可。
$F[i]=max(F[j]+1$并且$i,j$满足$F[j]$对应的大小关系$)$
复杂度，查询$O(1)$,插入$O(log10^6)$
最坏$O(nlogn)$吧
代码:

#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
#define N 500100
#define M 1000100
#define MM 1000000
#define lson l,mid,rt<<1
#define rson mid+1,r,rt<<1|1 
using namespace std;
int n,k;
int a[N];
int next[M];
int app[M];
int f[N],can[M];
char s[N];
char tmp[2];
int ma[M<<2],pos[M<<2],ma2[M<<2],pos2[M<<2];
void pushup(int rt)
{
    ma[rt]=ma[rt<<1],pos[rt]=pos[rt<<1];
    if(ma[rt<<1|1]>ma[rt])
        ma[rt]=ma[rt<<1|1],pos[rt]=pos[rt<<1|1];
}
void update(int p,int v,int l,int r,int rt)
{
    if(l==r)
    {
        ma[rt]=v;
        pos[rt]=l;
        return;
    }
    int mid=(l+r)>>1;
    if(p<=mid)update(p,v,lson);
    else update(p,v,rson);
    pushup(rt);
}
int query(int L,int R,int l,int r,int rt)
{
    if(L<=l&&r<=R)
    {
        return rt;
    }
    int ret=0;
    int mid=(l+r)>>1;
    if(L<=mid)
    {
        int tmp=query(L,R,lson);
        if(ma[tmp]>=ma[ret])ret=tmp; 
    }
    if(R>mid)
    {
        int tmp=query(L,R,rson);
        if(ma[tmp]>ma[ret])ret=tmp;
    }
    pushup(rt);
    return ret;
}
void pushup2(int rt)
{
    ma2[rt]=ma2[rt<<1],pos2[rt]=pos2[rt<<1];
    if(ma2[rt<<1|1]>ma2[rt])
        ma2[rt]=ma2[rt<<1|1],pos2[rt]=pos2[rt<<1|1];
}
void update2(int p,int v,int l,int r,int rt)
{
    if(l==r)
    {
        ma2[rt]=v;
        pos2[rt]=l;
        return;
    }
    int mid=(l+r)>>1;
    if(p<=mid)update2(p,v,lson);
    else update2(p,v,rson);
    pushup2(rt);
}
int query2(int L,int R,int l,int r,int rt)
{
    if(L<=l&&r<=R)
    {
        return rt;
    }
    int ret=0;
    int mid=(l+r)>>1;
    if(L<=mid)
    {
        int tmp=query2(L,R,lson);
        if(ma2[tmp]>=ma2[ret])ret=tmp; 
    }
    if(R>mid)
    {
        int tmp=query2(L,R,rson);
        if(ma2[tmp]>ma2[ret])ret=tmp;
    }
    pushup2(rt);
    return ret;
}
int main()
{
    scanf("%d%d",&n,&k);
    for(int i=1;i<=n;i++)scanf("%d",&a[i]),app[a[i]]=n+1;
    for(int i=n;i>=1;i--)next[a[i]]=app[a[i]],app[a[i]]=i;
    for(int i=1;i<=k;i++)
    {
        scanf("%s",tmp);
        s[i]=tmp[0];
    }
    f[1]=1;
    if(s[1]=='<')update(a[1],1,1,MM,1);
    else if(s[1]=='>')update2(a[1],1,1,MM,1);
    else can[a[1]]=1;
    for(int i=2;i<=n;i++)
    {
        if(can[a[i]])
            f[i]=f[can[a[i]]]+1;
        int tmp=0;
        if(a[i]>1)
            tmp=query(1,a[i]-1,1,MM,1);
        if(ma[tmp]+1>f[i])
            f[i]=ma[tmp]+1;
        tmp=0;
        if(a[i]<1000000)
            tmp=query2(a[i]+1,MM,1,MM,1);
        if(ma2[tmp]+1>f[i])
            f[i]=ma2[tmp]+1;
        if(s[(f[i]-1)%k+1]=='<')update(a[i],f[i],1,MM,1);
        else if(s[(f[i]-1)%k+1]=='>')update2(a[i],f[i],1,MM,1);
        else can[a[i]]=i;
    }
    int ans=0;
    for(int i=1;i<=n;i++)ans=max(ans,f[i]);
    printf("%d\n",ans);
}