Removed Interval HDU - 5489(LIS+线段树加速)

思路

一般的思路：先预处理得到f[i]以i结尾的最长上升子序列，g[i]以i开始的最长上升子序列。然后枚举去掉的L，从L+1项开始，那么答案应该是max(f[1]....f[i-L])+g[i]。 所以问题的关键转化为求前i-L项中，小于a[i]，并且f[]最大的值。 对于小于a[i[这个条件，我们可以考虑使用离散化来实现（离散化后，只找它前面的项就可以保证），对于找最大值，那直接线段树或树状数组就OK了。 这样做的正确性在于，我们是从i=L,开始一边查找一边更新线段树的，所以必然不会有错（当前查询，只和之前的操作有关）。 我觉得难点在于理解将值映射到线段树的过程。

代码

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int maxn=1e6+7;
const int inf=0x3f3f3f3f;
int n,L;
int a[maxn],b[maxn];
int f[maxn],g[maxn];
struct LIS //模板，可以取等时用upper_bound,不取等时用lower_bound
{
    int n, s[maxn];
    void init () { n = 0; }
    int find (int x) { return lower_bound(s, s + n, x) - s; }
    void modify(int pos, int val) { s[pos] = val; n = max(n, pos+1); }
}in, de;
void init()
{
    in.init();
    for(int i=1;i<=n;i++)
    {
        f[i]=in.find(a[i]);
        in.modify(f[i],a[i]);
        f[i]+=1;
    }
    de.init();
    for(int i=n;i>=1;i--)
    {
        g[i]=de.find(-a[i]);
        de.modify(g[i],-a[i]);
        g[i]+=1;
    }
}
#define Lson l,m,rt<<1
#define Rson m+1,r,rt<<1|1
struct Tree
{
    int l,r,Max;
}tree[maxn*3];
void push_up(int rt)
{
    tree[rt].Max=max(tree[rt<<1].Max,tree[rt<<1|1].Max);
}
void Build(int l,int r,int rt)
{
    tree[rt].l=l,tree[rt].r=r;
    tree[rt].Max=0;
    if(l==r) return;
    int m=(l+r)>>1;
    Build(Lson);
    Build(Rson);
}
void updata(int k,int v,int l,int r,int rt)
{
    if(l==k&&k==r)
    {
        tree[rt].Max=max(tree[rt].Max,v);
        return;
    }
    int m=(l+r)>>1;
    if(k<=m) updata(k,v,Lson);
    else updata(k,v,Rson);
    push_up(rt);
}
int query(int L,int R,int l,int r,int rt)
{
    if(L>R) return 0;
    if(L<=l&&R>=r)
    {
        return tree[rt].Max;
    }
    int m=(l+r)>>1;
    int ans=-1e9+3;
    if(L<=m) ans=max(ans,query(L,R,Lson));
    if(R>m) ans=max(ans,query(L,R,Rson));
    return ans;
}
int main()
{
    #ifndef ONLINE_JUDGE
        freopen("in.txt","r",stdin);
        freopen("out.txt","w",stdout);
    #endif
    int T;
    cin>>T;
    int Case=0;
    while(T--)
    {
        scanf("%d%d",&n,&L);
        for(int i=1;i<=n;i++)
        {
            scanf("%d",&a[i]);
            b[i]=a[i];
        }
        init();
        sort(b+1,b+1+n);
        int m=unique(b+1,b+1+n)-&b[1];
        int ans=0;
        Build(1,m,1);
        a[n+1]=inf,g[n+1]=0;
        for(int i=L+1;i<=n+1;i++)
        {
            int k=lower_bound(b+1,b+1+n,a[i])-b;
            ans=max(ans,query(1,k-1,1,m,1)+g[i]);
            k=lower_bound(b+1,b+1+m,a[i-L])-b;
            if(i<=n) updata(k,f[i-L],1,m,1);
        }
        printf("Case #%d: ",++Case);
        printf("%d\n",ans);
    }
    return 0;
}