【bzoj 1049】 数字序列 【HAOI2006】

Description

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

Input

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

Output

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

Sample Input

4
5 2 3 5

Sample Output

1
4

这道题第一问显然是裸的LIS,第二问设g[i]表示到第i位为止的最小值,则g[i]=min(g[j]+w[j][i])(f[j]+1=f[i]),其中f[i]为第一问DP数组,w[i][j]ij的最小修改,有一个结论即w[i][j]=min(a[i]*(k-i+1)+a[j]*(j-k+1))(i\leq k\leq j),然后直接DP就可以了(注意long long),下面是程序:

#include<stdio.h>
#include<iostream>
#include<string.h>
#include<algorithm>
#define ll long long
using namespace std;
const int N=35005;
struct node{
	int j,next;
}e[N];
int head[N],f[N],a[N],q[N],n,len,k;
ll g[N];
ll abs(ll n){
	return n>0?n:-n;
}
void add(int u,int v){
	e[++k]=(node){v,head[u]};
	head[u]=k;
}
int find(int x){
	int l=1,r=len,m,s=0;
	while(l<=r){
		m=l+r>>1;
		if(q[m]<=x){
			s=m;
			l=m+1;
		}
		else{
			r=m-1;
		}
	}
	return s;
}
void LIS(){
	int i,k;
	memset(q,127,sizeof(q));
	q[0]=-1<<30;
	for(i=1;i<=n;i++){
		k=find(a[i]);
		f[i]=k+1;
		len=max(len,k+1);
		q[k+1]=min(q[k+1],a[i]);
	}
	printf("%d\n",n-len);
}
int main(){
	int i,j,k,now;
	ll s1,s2;
	scanf("%d",&n);
	for(i=1;i<=n;i++){
		scanf("%d",&a[i]);
		a[i]-=i;
	}
	a[++n]=1<<30;
	LIS();
	for(i=n;i>=0;i--){
		add(f[i],i);
		g[i]=1ll<<60;
	}
	a[0]=-1<<30;
	g[0]=0;
	for(i=1;i<=n;i++){
		for(now=head[f[i]-1];now;now=e[now].next){
			j=e[now].j;
			s1=s2=0;
			if(j>i){
				break;
			}
			if(a[j]>a[i]){
				continue;
			}
			for(k=j;k<=i;k++){
				s2+=abs(a[k]-a[i]);
			}
			for(k=j;k<=i;k++){
				g[i]=min(g[i],s1+s2+g[j]);
				s1+=abs(a[j]-a[k]);
				s2-=abs(a[i]-a[k]);
			}
		}
	}
	printf("%lld\n",g[n]);
	return 0;
}

 

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值