HDU 5489 Removed Interval (LIS变形)

题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=5489

给你n个数，要删去其中连续的L个，问你删去之后的LIS最大是多少？

我们先预处理出以i下标为开头的LIS，存到数组中。

然后可以枚举长为L的区间，每次移动，左边增加一个，右边删除一个。

最长上升子序列长度 = 窗口右边以右边第一个元素开头的最长上升子序列 + 窗口左边最大元素小于窗口右边第一个元素的最长上升子序列。

比如 1 2 [4 3] 2 3 5  ,  LIS = 3 + 1 = 4

求以i开头的LIS 只要倒着求最长下降子序列 or 倒着a[i]变负求最长上升子序列即可。

代码写的有点乱。

 //#pragma comment(linker, "/STACK:102400000, 102400000")
 #include <algorithm>
 #include <iostream>
 #include <cstdlib>
 #include <cstring>
 #include <cstdio>
 #include <vector>
 #include <cmath>
 #include <ctime>
 #include <list>
 #include <set>
 #include <map>
 using namespace std;
 typedef long long LL;
 typedef pair <int, int> P;
 const int N = 1e5 + ;
 int dp1[N], dp2[N], a[N], inf = 1e9 + ;
 int y[N], x[N]; // y[i]表示以i下标为开头的LIS，x[i]表示以i为结尾的LIS

 int main()
 {
     int t, n, m;
     scanf("%d", &t);
     for(int ca = ; ca <= t; ++ca) {
         scanf("%d %d", &n, &m);
         for(int i = ; i <= n; ++i) {
             scanf("%d", a + i);
             dp1[i] = dp2[i] = inf;
         }
         dp2[] = dp1[] = inf;
         int ans = ;
         for(int i = n; i >= ; --i) { //倒着求变负，求以i开头的LIS
             int pos = lower_bound(dp2, dp2 + n, -a[i]) - dp2;
             y[i] = pos + ;
             //cout << y[i] << endl;
             dp2[pos] = -a[i];
             if(i == m + ) {
                 ans = lower_bound(dp2, dp2 + n, inf) - dp2;
             }
         }
         //cout << ans << endl;
         for(int i = ; i <= n; ++i) {
             int pos = lower_bound(dp1, dp1 + n, a[i]) - dp1;
             x[i] = pos + ;
             dp1[pos] = a[i];
             if(i == n - m) {
                 ans = max(int(lower_bound(dp1, dp1 + n, inf) - dp1), ans);
                 break;
             }
         }
         //cout << ans << endl;
         printf("Case #%d: ", ca);
         dp1[] = inf;
         for(int i = ; i <= n; ++i) {
             dp1[i] = inf;
         }
         for(int i = m + ; i <= n; ++i) {
             int pos = lower_bound(dp1, dp1 + n, a[i] - ) - dp1; //求滑窗左边刚好小于滑窗右边第一个数的LIS
             ans = max(ans, y[i] + (a[i] -  == a[pos] ? pos +  : pos));
             *lower_bound(dp1, dp1 + n, a[i - m]) = a[i - m];
         }
         printf("%d\n", ans);
     }
     return ;
 }