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 = ; 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 = ;i <= n;++i) scanf("%d",a+i),a[i] -= i,bac[++tot] = a[i];
     a[++n] = inf; bac[++tot] = a[n];
     sort(bac+,bac+tot+); tot = unique(bac+,bac+tot+) - bac - ;
     for (int i = ;i <= n;++i)
     {
         int pos = lower_bound(bac+,bac+tot+,a[i]) - bac;
         f[i] = calc(pos) + ;
         change(pos,f[i]);
     }
     printf("%d\n",n-f[n]);
     for (int i = ;i <= n;++i) vec[f[i]].push_back(i);
     a[] = -inf;
     memset(g,0x7,sizeof(g)); g[] = ;
     for (int i = ;i <= n;++i)
     {
         int nn = vec[f[i] - ].size();
         for (int p = ;p < nn;++p)
         {
             int j = vec[f[i]-][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+;k <= i;++k) s1[k] += s1[k-],s2[k] += s2[k-];
             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 ;
 }