Codeforces Round #371 (Div. 1) C. Sonya and Problem Wihtout a Legend

最新推荐文章于 2020-11-01 21:12:05 发布

zqh_wz

最新推荐文章于 2020-11-01 21:12:05 发布

阅读量690

点赞数

分类专栏： Codeforces dp 文章标签： Codeforces

本文链接：https://blog.csdn.net/zqh_wz/article/details/52648340

版权

dp 同时被 2 个专栏收录

16 篇文章 0 订阅

订阅专栏

Codeforces

2 篇文章 0 订阅

订阅专栏

Codeforces Round #371 (Div. 1) C. Sonya and Problem Wihtout a Legend

Sonya was unable to think of a story for this problem, so here comes the formal description.

You are given the array containing n positive integers. At one turn you can pick any element and increase or decrease it by 1. The goal is the make the array strictly increasing by making the minimum possible number of operations. You are allowed to change elements in any way, they can become negative or equal to 0.

Input

The first line of the input contains a single integer n (1 ≤ n ≤ 3000) — the length of the array.

Next line contains n integer ai $(1 ≤ ai ≤ 10^9)$ .

Output

Print the minimum number of operation required to make the array strictly increasing.

Examples

Input

7

2 1 5 11 5 9 11

Output

9

Input

5

5 4 3 2 1

Output
12

Note

In the first sample, the array is going to look as follows:

2 3 5 6 7 9 11

|2 - 2| + |1 - 3| + |5 - 5| + |11 - 6| + |5 - 7| + |9 - 9| + |11 - 11| = 9

And for the second sample:

1 2 3 4 5

|5 - 1| + |4 - 2| + |3 - 3| + |2 - 4| + |1 - 5| = 12

Task：

给定一个序列的 $n$ 个数，每次可以对其中的一个数加一或减一。求最终使得整个序列单调递增的最小操作数。

Solution：

CF上这题挂了一个flow标签，并不懂是什么意思…~~（其实是自己功力尚浅）~~

直接这么做似乎有点困难，我们可以用一个小技巧将这个要求转化一下：

A [i] < A [i + 1]

$A[i]<A[i+1]$

A [i] - i < = A [i + 1] - (i + 1)

$A[i]-i<=A[i+1]-(i+1)$
既然将小于变成了小于等于，那么就可以转化为dp问题做了，令

dp[i][j] $dp[i][j]$ 为前i个数，最后一个结尾的数大小为

A[j] $A[j]$ （已减

j $j$ ）所需要的代价。复杂度

O(n2) $O(n^2)$ 。

注意：因为最后的所有值如果在 $A$ 数组中一定不会差于不在 $A$ 数组中的情况，所以我们的第二维只需要开到 $n$ 即可，而不是无上限地增加。

#include<stdio.h>
#include<iostream>
#include<algorithm>
#define ll long long
#define INF (1LL<<62)
#define M 3005
using namespace std;
int A[M],B[M];
ll dp[2][M];
int main(){
    int n;
    scanf("%d",&n);
    for(int i=1;i<=n;i++){
        scanf("%d",&A[i]);
        A[i]-=i;
        B[i-1]=A[i];
    }
    sort(B,B+n);
    int m=unique(B,B+n)-B;
    for(int i=0;i<m;i++)
        dp[1][i]=abs(A[1]-B[i]);
    for(int i=2;i<=n;i++){
        ll mi=INF;
        int cur=i&1;
        for(int j=0;j<m;j++){
            if(dp[!cur][j]<mi)mi=dp[!cur][j];
            dp[cur][j]=mi+abs(A[i]-B[j]);
        }
    }
    ll mi=INF;
    for(int i=0;i<m;i++)
        if(dp[n&1][i]<mi)mi=dp[n&1][i];
    cout<<mi<<endl;
    return 0;
}

然而，CF上的各路大神有一种 $O(nlogn)$ 的解法。Orz…

具体就是每次把当前这个值加入到堆中，从堆中取出最大值。如果最大值大于当前的数：将答案加上这个最大值与当前数的差值，将最大值弹出，而再次加入一个当前数。即相当于将最大值替换为当前数。

说起来有点玄学，本蒟蒻目前也无力证明，只能暂且先当是玄学好了。

#include<stdio.h>
#include<iostream>
#include<queue>
#define ll long long
#define M 3005
using namespace std;
priority_queue<int>Q;
int main(){
    int n;ll ans=0;
    scanf("%d",&n);
    for(int i=1;i<=n;i++){
        int x;
        scanf("%d",&x);
        x-=i;
        Q.push(x);
        int mx=Q.top();
        if(mx>x){
            ans+=mx-x;
            Q.pop();
            Q.push(x);
        }
    }
    cout<<ans<<endl;
    return 0;
}