bzoj1049: [HAOI2006]数字序列 最长不降子序列的特殊操作

bzoj1049: [HAOI2006]数字序列

Description

现在我们有一个长度为n的整数序列A。但是它太不好看了，于是我们希望把它变成一个单调严格上升的序列。
但是不希望改变过多的数，也不希望改变的幅度太大。

Input

第一行包含一个数n，接下来n个整数按顺序描述每一项的键值。n<=35000,保证所有数列是随机的

Output

第一行一个整数表示最少需要改变多少个数。 第二行一个整数，表示在改变的数最少的情况下，每个数改变
的绝对值之和的最小值。

Sample Input

4
5 2 3 5

Sample Output

1
4

分析

第一问补集转化一下，变成求$a_j<a_i,a_i-a_j\ge i-j$的一个序列。
令$b_i=a_i-i$变成求$b$的最长不降子序列。
对于第二问有一个神仙操作。
首先求出合法的转移，一种办法是把$i$挂到$f_i$上面，那么对于一个$i$，所有合法的最优转移一定是$f_i-1$中比$i$小的且$b$的值小于$b_i$的。
这样的做法是$O(n^2)$的，但是随机情况下实际上答案很少。
考虑一组合法的转移。如果是从$i->j$的，那么这之间的数应该怎么改变？
有一个简单的想法就是肯定是把在$i$到$j$之间的数都挪到一条水平线上。

由于是一组合法的转移，所以一定不存在在$b_i$和$b_j$之间的数。
所以对于比$i$小的数，肯定最少要把它们挪到$i$上，对于比$j$大的数，至少要把他们挪到$j$上。
可是有可能出现下面这种情况。

这个时候发现，对于不合法的那两个节点，只能挪到一边去。
至于挪到哪边，或者是挪到某个中间点，代价是不变的。
因为一个少往上挪多少，另一个就要多往下挪多少，最后挪到一个相等的水平线上。
所以可以枚举一个分界点$k$，$k$之前的放在$b_i$上，$k$之后的放在$b_j$上的即可。
这个过程实际上是$O(n^3)$的，但是在数据随机意义下常数很小，可以通过。

代码

#include<cstdio>
#include<cmath>
#include<algorithm>
const int N = 35500, tf = 0x3f3f3f3f; const long long lf = 0x3f3f3f3f3f3f3f3f;
int ri() {
    char c = getchar(); int x = 0, f = 1; for(;c < '0' || c > '9'; c = getchar()) if(c == '-') f  = -1;
    for(;c >= '0' && c <= '9'; c = getchar()) x = (x << 1) + (x << 3) - '0' + c; return x * f;
}
int a[N], b[N], f[N], pr[N], nx[N], to[N], l, C; long long g[N], s1[N], s2[N];
void add(int u, int v) {to[++C] = v; nx[C] = pr[u]; pr[u] = C;}
int F(int x) {
    int L = 1, R = l, m; 
    for(;L != R; a[b[m = L + R >> 1]] > x ? R = m : L = m + 1) ;
    return L;
}
int Cal(int L, int R, int x) {
    int r = 0;
    for(int i = L;i <= R; ++i) r += abs(a[i] - x);
    return r;
}
int main() {
    int n = ri();
    for(int i = 1;i <= n; ++i) a[i] = ri() - i;
    b[l = 1] = 1; f[1] = 1; a[++n] = tf;
    for(int i = 2;i <= n; ++i) {
        if(a[i] >= a[b[l]]) b[f[i] = ++l] = i;
        else b[f[i] = F(a[i])] = i;
    }
    printf("%d\n", n - l); 
    for(int i = n; ~i; --i) add(f[i], i), g[i] = lf;
    g[0] = 0; a[0] = -tf;
    for(int x = 1;x <= n; ++x) 
        for(int i = pr[f[x] - 1], v; i; i = nx[i]) {
            if((v = to[i]) > x) break;
            if(a[v] > a[x]) continue;
            for(int j = v; j <= x; ++j) s1[j] = abs(a[v] - a[j]), s2[j] = abs(a[x] - a[j]);
            for(int j = v + 1; j <= x; ++j) s1[j] += s1[j - 1], s2[j] += s2[j - 1];
            for(int j = v; j < x; ++j) g[x] = std::min(g[x], g[v] + s1[j] - s1[v] + s2[x] - s2[j]);
        }
    printf("%lld\n", g[n]);
    return 0;
}