BZOJ 1049 数字序列

Description

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

Input

第一行包含一个数n，接下来n个整数按顺序描述每一项的键值。

Output

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

Sample Input

4
5 2 3 5

Sample Output

1
4

HINT

 

【数据范围】

90%的数据n<=6000。

100%的数据n<=35000。

保证所有数列是随机的。

 

Source

 

这道题感觉太神了。

令原数组为a。

对于第一问我们可以反过来思考——要求最少的改动的，不就是求最多的不改动的吗？

f[i]表示前i位最多的不改动的数字个数，f[i]=max(f[j]+1)，j需要满足a[i]-a[j]≥i-j。对于条件，我们移一下项，化为a[i]-i≥a[j]-j，令b[i]=a[i]-i，不就是b[i]≥b[j]。细心的朋友一定看出来了，这不就是最长不下降子序列吗！！！O(nlogn)的求法上起。

第二问稍微麻烦一点，我们要知道一个结论：另g[i]表示前i个数，在改动数最少的前提下，最少改动的值。

那么对于所有合法的转移i，j(i>j，f[i]=f[j]+1，b[i]≥b[j])，最优解一定是在i与j之间某个k，k到j的值全为b[j]，k+1到i的值全为b[i]。因此就可以dp了。证明自己脑补一下就可以了。。。(自己画画图，根据条件想想应该是可以明白的)

 

#include<cstring>
#include<vector>
#include<algorithm>
#include<iostream>
#include<cstdio>
#include<cstdlib>
using namespace std;

#define inf (1<<30)
#define maxn 35010
int tree[maxn],tot,n,f[maxn],a[maxn],bac[maxn];
long long g[maxn],s1[maxn],s2[maxn];
vector <int> vec[maxn];

inline int lowbit(int x) { return x & -x; }

inline void change(int a,int b) { for (;a <= tot;a += lowbit(a)) tree[a] = max(tree[a],b); }

inline int calc(int a) { int ret = 0; for (;a;a -= lowbit(a)) ret = max(ret,tree[a]); return ret; }

int main()
{
    freopen("1049.in","r",stdin);
    freopen("1049.out","w",stdout);
    scanf("%d",&n);
    for (int i = 1;i <= n;++i) scanf("%d",a+i),a[i] -= i,bac[++tot] = a[i];
    a[++n] = inf; bac[++tot] = a[n];
    sort(bac+1,bac+tot+1); tot = unique(bac+1,bac+tot+1) - bac - 1;
    for (int i = 1;i <= n;++i)
    {
        int pos = lower_bound(bac+1,bac+tot+1,a[i]) - bac;
        f[i] = calc(pos) + 1;
        change(pos,f[i]);
    }
    printf("%d\n",n-f[n]);
    for (int i = 0;i <= n;++i) vec[f[i]].push_back(i);
    a[0] = -inf;
    memset(g,0x7,sizeof(g)); g[0] = 0;
    for (int i = 1;i <= n;++i)
    {
        int nn = vec[f[i] - 1].size();
        for (int p = 0;p < nn;++p)
        {
            int j = vec[f[i]-1][p];
            if (j >= i) break; if (a[i] < a[j]) continue;
            for (int k = j;k <= i;++k) s1[k]=abs(a[k]-a[j]),s2[k]=abs(a[k]-a[i]);
            for (int k = j+1;k <= i;++k) s1[k] += s1[k-1],s2[k] += s2[k-1];
            for (int k = j;k < i;++k) g[i] = min(g[i],g[j]+s1[k]-s1[j]+s2[i]-s2[k]);
        }
    }
    printf("%lld",g[n]);
    fclose(stdin); fclose(stdout);
    return 0;
}