景观美化

景观美化

题目

农夫约翰最近决定来美化他的花园，他需要运输很多的泥土。花园是由N块花圃组成的。第i块花圃初始的时候有Ai数量的泥土。为了达到美化的目的，必须使得第i块花圃的泥土数量Ai变成Bi。

约翰有三个选择：第一，他可以买一个单位的泥土放进任意花圃中，代价是X；第二，他可以将一个单位的泥土从某一个花圃中除去，代价是Y；第三，他可以将第i块花圃中的一个单位的泥土搬运到第j块花圃中，大家是Z*|i-j|。

请帮助约翰计算为了达到目的最小需要花费的代价。

输入

第一行四个整数，分别是N，X，Y，Z。
接下来N行，每行两个整数，分别表示Ai和Bi。

输出

只有一行一个整数，表示最小的代价。

输入样例

4 100 200 1 
1 4 
2 3 
3 2 
4 0 

输出样例

210

样例说明

从第4个花圃中所有的土必须被除去，其中1个单位的土被直接除去，代价是200，剩下3个单位的土从第4个花圃到第1个花圃

数据范围

1<=N<=100
0<=Ai,Bi<=10
0<=X,Y,Z<=1000

思路

我们将每个花圃的泥土都分开来，拿样例来说，就变成了：
1223334444
111122233
然后就变成经典DP“编辑字符距离”了
设 f[i][j] 为原前i堆土变成目标前j堆土所需的最小代价
有三种情况：

1. 填土，则动态转移方程为：$f[i][j-1]+x$
2. 删土，则动态转移方程为：$f[i-1][j]+y$
3. 转移土，则动态转移方程为：$f[i-1][j-1]+z\ast abs(a[i]-b[j])$

每一次求三种情况哪一种更小，最后输出 f[l][ll] 就可以了。
（l为第一个字符串的长度，ll为第二个字符串的长度）

代码

#include<iostream>
#include<cstdio>
#include<cstring>
#define abs(aa) ((aa)<0?-(aa):(aa))//求绝对值
using namespace std;
int x,y,z,n,xx,yy,a[1001],b[1001],f[1001][1001],l,ll;
int main()
{
	scanf("%d%d%d%d",&n,&x,&y,&z);//读入
	for (int i=1;i<=n;i++)
	{
		scanf("%d%d",&xx,&yy);//读入
		for (int j=1;j<=xx;j++)
		a[++l]=i;//把花园的泥土分开
		for (int j=1;j<=yy;j++)
		b[++ll]=i;//把花园的泥土分开（美化后）
	}
	memset(f,0x7f,sizeof(f));//初始化
	for (int i=0;i<=l;i++)
	f[i][0]=y*i;//初始化
	for (int i=0;i<=ll;i++)
	f[0][ll]=x*i;//初始化
	for (int i=1;i<=l;i++)
	for (int j=1;j<=ll;j++)
	f[i][j]=min(f[i-1][j-1]+z*abs(a[i]-b[j]),min(f[i-1][j]+y,f[i][j-1]+x));//动态转移方程
	printf("%d",f[l][ll]);//输出
}