【HDU】5489 Removed Interval（2015 ACM/ICPC Asia Regional Hefei Online）

Removed Interval

---

题目链接

• Removed Interval

---

题目大意

    变相LIS，从A中间某一位置删去一部分连续的数，让你求剩下部分的最大LIS。

---

题解

DP

    因为被分开了，所以考虑分开来做，用线段树做$O(nlog_2n)$的LIS，有几个数组需要维护：

• r[i]：从i开始的，最大LIS长度。
• l[i]：以i结尾的，最大LIS长度。
• dp[i]:以i结尾的，删去$[i-L,i-1]$之后的最长LIS长度。

    可以看到最终的结果就是$max(dp[i]+r[i]-1)\ \ \ \ (L<i<=n)$注意要特别判断删除最后L个的情况（把循环写到n+1也可以）。

---

代码

#include <iostream>
#include <cstring>
#include <cstdio>
#include <algorithm>
#define maxn 100005

using namespace std;

struct tree
{
    int l,r,val;
};
tree d[maxn<<2],D[maxn<<2];
struct Mine
{
    int val,p;
};
Mine a[maxn];
int T,l[maxn],r[maxn],n,L,dp[maxn];

bool cmp2(const Mine &a,const Mine &b)
{
    if (a.val!=b.val) return a.val<b.val;
    else return a.p>b.p;
}

void buildtree(int n,int l,int r)
{
    d[n].l=l; d[n].r=r; d[n].val=0;
    if (l==r) return ;
    int mid=(l+r)>>1;
    buildtree(n<<1,l,mid);
    buildtree(n<<1|1,mid+1,r);
}

int Count(int n,int l,int r)
{
    if (l>r) return 0;
    if (l==d[n].l && r==d[n].r) return d[n].val;
    int mid=(d[n].l+d[n].r)>>1;
    if (r<=mid) return Count(n<<1,l,r);
    else if (l>mid) return Count(n<<1|1,l,r);
    else return max(Count(n<<1,l,mid),Count(n<<1|1,mid+1,r));
}

void Insert(int n,int p,int val)
{
    if (d[n].r==d[n].l)
    {
        d[n].val=val;
        return ;
    }
    int mid=(d[n].l+d[n].r)>>1;
    if (p<=mid) Insert(n<<1,p,val);
    else Insert(n<<1|1,p,val);
    d[n].val=max(d[n<<1].val,d[n<<1|1].val);
}

int main()
{
    int Case=1;
    scanf("%d",&T);
    while(T--)
    {
        memset(d,0,sizeof(d));
        memset(r,0,sizeof(r));
        memset(dp,0,sizeof(dp));
        scanf("%d%d",&n,&L);
        for (int i=1;i<=n;i++)
        {
            scanf("%d",&a[i].val);
            a[i].p=i;
        }
        int i=n,j=n,p=n;
        sort(a+1,a+1+n,cmp2);

        //solve r[i]
        buildtree(1,1,n);
        for (int i=n;i>=1;i--)
        {
            int c=Count(1,a[i].p+1,n);
            r[a[i].p]=c+1;
            Insert(1,a[i].p,r[a[i].p]);
        }
        buildtree(1,1,n);
        for (int i=1;i<=n;i++)
        {
            int c=Count(1,1,a[i].p-1);
            l[a[i].p]=c+1;
            Insert(1,a[i].p,l[a[i].p]);
            if (a[i].p-L>=1) dp[a[i].p]=Count(1,1,a[i].p-L-1)+1;
        }
        int ans=0;
        for (int i=L+1;i<=n;i++) ans=max(ans,dp[i]+r[i]-1);
        for (int i=1;i<=n-L;i++) ans=max(ans,l[i]);
        printf("Case #%d: %d\n",Case++,ans);
    }
    return 0;
}