这是道双倍经验题，做完可以做Luogu P4597
01
( exttt{离散化+滚动数组})
Idea
(dp[i][j])表示前i个数以(a[j])为结尾的满足要求的最少的操作，可是题目给的最大数是(10^9)，二维数组的(j)元素不可能开这么大，所以需要离散化一下，改成前i个数以第(j)个数为结尾的满足要求的最少的操作。
(dp[i][j]=min(dp[i][j], ext{第i-1个位置前j个数的最小操作}+abs(b[j]-a[i]))) (b)数组是原来输入的(a)数组排好序且离散化之后的数组
因为(b)数组是从小到大排好序的，所以在第(i)个位置时第(j)个数肯定比(i-1)个位置的前j个数要大。
最后(dp[5000][5000])也不能存下，那么这个二维数组就要滚一下了。
Code
#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<cstdlib>
#include<algorithm>
#include<string>
#include<queue>
#define ll long long
#define maxn 100001
#define mod 998244353
#define eps 1e-6
#define pi acos(-1.0)
#define de(x) ((x)*(x))
using namespace std;
inline int read(){
int x=0,f=1; char ch=getchar();
while(!isdigit(ch)) {if(ch=='-')f=-1;ch=getchar();}
while(isdigit(ch)) {x=x*10+ch-48;ch=getchar();}
return x*f;
}
const ll inf=1ll<<62;
ll dp[2][5007],a[5007],b[5007];
int main(){
ll n=read();
for(ll i=1;i<=n;i++){
a[i]=read();
b[i]=a[i];
}
sort(b+1,b+n+1);
ll tot=unique(b+1,b+n+1)-b-1;//离散化
dp[1][1]=inf;
for(ll i=1;i<=tot;i++){
dp[1&1][i]=fabs(b[i]-a[1]);
dp[2&1][i]=inf;
}
for(ll i=2;i<=n;i++){
ll minn=inf;
for(ll j=1;j<=tot;j++){
minn=min(minn,dp[(i-1)&1][j]);//取上一个位置中前j颗不同高度的树中花费最小
dp[i&1][j]=min(dp[i&1][j],minn+abs(b[j]-a[i]));//当前位置等于上一个位置中前j颗树花费最小+这个位置是第j颗树的花费
}
for(ll j=1;j<=tot;j++) dp[(i+1)&1][j]=inf;
}
ll ans=inf;
for(ll i=1;i<=tot;i++) ans=min(ans,dp[n&1][i]);//取最后一个位置是第j颗树是花费最小
printf("%lld",ans);
return 0;
}
02
考虑(DP)
Idea
设(dp[i][j])表示前(i)个已经符合要求，而且最大数不大于原序列第(j)个元素最小需要的代价。
那么我们可以得出转移方程(dp[i][j]=min(dp[i][j−1],dp[i−1][j]+abs(a[i]−c[j]))
#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<cstdlib>
#include<algorithm>
#include<string>
#include<queue>
#define ll long long
#define maxn 100001
#define mod 998244353
#define eps 1e-6
#define pi acos(-1.0)
#define de(x) ((x)*(x))
using namespace std;
inline ll read(){
ll x=0,f=1; char ch=getchar();
while(!isdigit(ch)) {if(ch=='-')f=-1;ch=getchar();}
while(isdigit(ch)) {x=x*10+ch-48;ch=getchar();}
return x*f;
}
const ll inf=1ll<<62;
ll ans=inf;
ll a[maxn],c[maxn],dp[2][maxn];
int main(){
int n; scanf("%d",&n);
for(int i=1;i<=n;i++) a[i]=read(),c[i]=a[i];
sort(&c[1],&c[n+1]);
memset(dp,127,sizeof dp);
for(int i=1;i<=n;i++) dp[0][i]=0;
for(int i=1;i<=n;i++){
for(int j=1;j<=n;j++)
dp[1][j]=min(dp[1][j-1],dp[0][j]+abs(a[i]-c[j]));
swap(dp[1],dp[0]);
}
for(int i=1;i<=n;i++) ans=min(ans,dp[0][i]);
printf("%I64d",ans);
return 0;
}
(P.S: ext{这个方法过不了CF原网站的#8})
DP优化
排完序后发现，每个数的更换值都是最小的.
那么有(dp)方程如下
[dp[j]=abs(a[i]-a[j]);
]
[dp[j]=min(dp[j],dp[j-1]);
]
#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<cstdlib>
#include<algorithm>
#include<string>
#include<queue>
#define ll long long
#define maxn 5005
#define inf 2147483647
#define mod 998244353
#define eps 1e-6
#define pi acos(-1.0)
#define de(x) ((x)*(x))
using namespace std;
inline ll read(){
ll x=0,f=1; char ch=getchar();
while(!isdigit(ch)) {if(ch=='-')f=-1;ch=getchar();}
while(isdigit(ch)) {x=x*10+ch-48;ch=getchar();}
return x*f;
}
ll a[maxn],c[maxn],dp[maxn];
int main(){
int n; scanf("%d",&n);
for(int i=1;i<=n;i++) a[i]=read(),c[i]=a[i];
sort(c+1,c+n+1);
for(int i=1;i<=n;i++)
for(int j=1;j<=n;j++){
dp[j]+=abs(a[i]-c[j]);
if(j>1) dp[j]=min(dp[j],dp[j-1]);
}
printf("%I64d",dp[n]);
return 0;
}//洛谷上用时4+ s,CF上用时342 ms
03
用堆写
证明请看(我不会严格证明）
Code
#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<cstdlib>
#include<algorithm>
#include<string>
#include<queue>
#define ll long long
#define maxn 500005
#define inf 2147483647
#define mod 998244353
#define eps 1e-6
#define pi acos(-1.0)
#define de(x) ((x)*(x))
using namespace std;
inline int read(){
int x=0,f=1; char ch=getchar();
while(!isdigit(ch)) {if(ch=='-')f=-1;ch=getchar();}
while(isdigit(ch)) {x=x*10+ch-48;ch=getchar();}
return x*f;
}
priority_queue<int>q;
ll ans;
int main(){
int n=read();
for(int i=1;i<=n;i++){
int x=read();
q.push(x);
if(q.top()>x){
ans+=q.top()-x;
q.pop();
q.push(x);
}
}
printf("%lld",ans);
return 0;
}