POJ3666Making the Grade[DP 离散化 LIS相关]

Making the Grade
Time Limit: 1000MS Memory Limit: 65536K
Total Submissions: 6445 Accepted: 2994

Description

A straight dirt road connects two fields on FJ's farm, but it changes elevation more than FJ would like. His cows do not mind climbing up or down a single slope, but they are not fond of an alternating succession of hills and valleys. FJ would like to add and remove dirt from the road so that it becomes one monotonic slope (either sloping up or down).

You are given N integers A1, ... , AN (1 ≤ N ≤ 2,000) describing the elevation (0 ≤ Ai ≤ 1,000,000,000) at each of N equally-spaced positions along the road, starting at the first field and ending at the other. FJ would like to adjust these elevations to a new sequence B1, . ... , BN that is either nonincreasing or nondecreasing. Since it costs the same amount of money to add or remove dirt at any position along the road, the total cost of modifying the road is

| A B 1| + | A B 2| + ... + | AN -  BN |

Please compute the minimum cost of grading his road so it becomes a continuous slope. FJ happily informs you that signed 32-bit integers can certainly be used to compute the answer.

Input

* Line 1: A single integer: N
* Lines 2..N+1: Line i+1 contains a single integer elevation: Ai

Output

* Line 1: A single integer that is the minimum cost for FJ to grade his dirt road so it becomes nonincreasing or nondecreasing in elevation.

Sample Input

7
1
3
2
4
5
3
9

Sample Output

3

Source


和那道 CF#371(div.2)E一样,只是不严格单增单减各一遍
//
//  main.cpp
//  poj3666
//
//  Created by Candy on 9/22/16.
//  Copyright © 2016 Candy. All rights reserved.
//

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
using namespace std;
const int N=2005,INF=1e9+5;
int read(){
    char c=getchar();int x=0,f=1;
    while(c<'0'||c>'9'){if(c=='-')f=-1; c=getchar();}
    while(c>='0'&&c<='9'){x=x*10+c-'0'; c=getchar();}
    return x*f;
}
int n,a[N],mp[N],k,ans=INF;
int f[N][N];
void dp(){
    sort(mp+1,mp+1+n);
    for(int i=1;i<=n;i++){
        int mn=INF;
        for(int j=1;j<=k;j++){
            mn=min(mn,f[i-1][j]);
            f[i][j]=mn+abs(a[i]-mp[j]);
        }
    }
    for(int j=1;j<=k;j++) ans=min(ans,f[n][j]);
    
    for(int i=1;i<=n;i++){
        int mn=INF;
        for(int j=k;j>=1;j--){
            mn=min(mn,f[i-1][j]);
            f[i][j]=mn+abs(a[i]-mp[j]);
        }
    }
    for(int j=1;j<=k;j++) ans=min(ans,f[n][j]);
}
int main(int argc, const char * argv[]) {
    n=read();
    for(int i=1;i<=n;i++){
        a[i]=mp[i]=read();
    }
    k=unique(mp+1,mp+1+n)-mp-1;
    dp();
    printf("%d",ans);
    return 0;
}

 

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值