BZOJ系列1049《[HAOI2006]数字序列》题解

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。

保证所有数列是随机的。

分析：

第一问是个LIS，不过需要注意，一般的LIS是i>j，则a[i]>a[j]。

但这道题因为需要改变，所以要给中间的数留出空间，所以LIS应该是：

a[i]-a[j]>=i-j  ————>>a[i]-i>=a[j]=j；

所以让b[i]=a[i]-i处理，求b[i]的最长不上升子序列即可（注意时间问题，LIS要优化）。

第二问其实用第一问的思路来看是很好写的，因为改变需要基于最少改变次数，

也就是说，改变的数已经确定。

比如一些需要改变的数在i与j之间：

可以把他们的b都改成b[i]或改为b[j]，取小即可。（有点说不清，自己想想）。

改变的数为第k个数时，改变为abs(b[k]-b[i])；

其实实质就是中间的这些数改为上升序列要么跟着i要么跟着j。

额，个人能力有限，写了写发现代码能力疲软，所以借助了他人帮助。

代码如下：

#include<iostream>
#include<cmath>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
typedef long long ll;
const ll MAXN=(int)(35010);
const ll inf=(int)(999999999);
int N,a[MAXN],b[MAXN],aa[MAXN],f[MAXN],linkk[MAXN],len,l=0;
ll g[MAXN],s1[MAXN],s2[MAXN];
struct list
{
    int to,next;
}e[35005];
void insert(int u,int v)
{
    e[++l].to=v;
    e[l].next=linkk[u];
    linkk[u]=l;
}
void init()
{
    scanf("%d",&N);
    for(int i=1;i<=N;i++)
    {
        scanf("%d",&a[i]);
        b[i]=a[i]-i;
    }
}
int find(int x)//二分优化。
{
    int l=1,r=len,t=0;
    while(l<=r)
    {
        int mid=(l+r)/2;
        if(aa[mid]<=x)
        {
            t=mid;
            l=mid+1;
        }
        else r=mid-1;
    }
    return t;
}
void DP()//求b[i]的LIS。
{
    b[++N]=1<<30;
    memset(aa,127,sizeof(aa));
    aa[0]=-1<<30;
    for(int i=1;i<=N;i++)
    {
        int t=find(b[i]);
        f[i]=t+1;
        len=max(len,t+1);
        aa[t+1]=min(aa[t+1],b[i]);
    }
    printf("%d\n",N-len);
}
void work()
{
    for(int i=N;i>=0;i--)
    {
        insert(f[i],i);
        g[i]=1LL<<60;
    }
    g[0]=0;b[0]=-1<<30;
    for(int x=1;x<=N;x++)
        for(int i=linkk[f[x]-1];i;i=e[i].next)
        {
            int p=e[i].to;
            if(p>x) break;
            if(b[p]>b[x]) continue;
            for(int j=p;j<=x;j++)
                s1[j]=abs(b[p]-b[j]),s2[j]=abs(b[x]-b[j]);
            for(int j=p+1;j<=x;j++)
                s1[j]+=s1[j-1],s2[j]+=s2[j-1];
            for(int j=p;j<x;j++)
                g[x]=min(g[x],g[p]+s1[j]-s1[p]+s2[x]-s2[j]);
        }
    printf("%d\n",g[N]);
}
int main()
{
    init();
    DP();
    work();
    return 0;
}