3042

/*
DP

这是个很有趣的DP问题

状态的设计是这样子的 f[i][j][2];
f[i][j][0] 表示吃完i到j之间所有的草,并且最后吃i,整体腐烂度的最小值
f[i][j][1] 表示吃完i到j之间所有的草,并且最后吃j, 整体腐烂度的最小值

我们要求 f[1][n][0] 和 f[1][n][1],那么如何得到呢
f[i][j][0] = min( f[i+1][j][0]+dst[i][i+1]*(i+n-j),f[i+1][j][1]+dst[i][j]*(i+n-j) );
f[i][j][1] = min( f[i][j-1][0]+dst[i][j]*(i+n-j),f[i][j-1][1]+dst[j-1][j]*(i+n-j) );
这样的状态转移方程是好理解的,因为可以左右走
那么f[i][j][0] 可以是(i+1,j)两个端点过来,所以就有了这个方程。方程的构造还是不大容易想到的


*/

// include file
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <cmath>
#include <cctype>
#include <ctime>

#include <iostream>
#include <sstream>
#include <fstream>
#include <iomanip>
#include <bitset>

#include <algorithm>
#include <string>
#include <vector>
#include <queue>
#include <set>
#include <list>
#include <functional>

using namespace std;

// typedef
typedef long long LL;
typedef unsigned long long ULL;

// 
#define read freopen("in.txt","r",stdin)
#define write freopen("out.txt","w",stdout)
#define FORi(a,b,c) for(int i=(a);i<(b);i+=c)
#define FORj(a,b,c) for(int j=(a);j<(b);j+=c)
#define FORk(a,b,c) for(int k=(a);k<(b);k+=c)
#define FORp(a,b,c) for(int p=(a);p<(b);p+=c)
#define FORii(a,b,c) for(int ii=(a);ii<(b);ii+=c)
#define FORjj(a,b,c) for(int jj=(a);jj<(b);jj+=c)
#define FORkk(a,b,c) for(int kk=(a);kk<(b);kk+=c)

#define FF(i,a)    for(int i=0;i<(a);i++)
#define FFD(i,a)   for(int i=(a)-1;i>=0;i--)

#define Z(a) (a<<1)
#define Y(a) (a>>1)

const double eps = 1e-6;
const double INFf = 1e100;
const int INFi = 1000000000;
const LL INFll = (LL)1<<62;
const double Pi = acos(-1.0);

template<class T> inline T sqr(T a){return a*a;}
template<class T> inline T TMAX(T x,T y)
{
	if(x>y) return x;
	return y;
}
template<class T> inline T TMIN(T x,T y)
{
	if(x<y) return x;
	return y;
}
template<class T> inline void SWAP(T &x,T &y)
{
	T t = x;
	x = y;
	y = t;
}
template<class T> inline T MMAX(T x,T y,T z)
{
	return TMAX(TMAX(x,y),z);
}
template<class T> inline T MMIN(T x,T y,T z)
{
	return TMIN(TMIN(x,y),z);
}


// code begin
//
//
int N,L;
int data[1010];
int dst[1010][1010];
int DP[1010][1010][2];
bool cmp(int a,int b)
{
	return a<b;
}
int main()
{
	read;
	write;
	int l,r;
	while(scanf("%d %d",&N,&L)!=-1)
	{
		FORi(1,N+1,1)
		{
			scanf("%d",data+i);
		}
		data[++N] = L;

		sort(data+1,data+N+1);

		FORi(1,N+1,1)
		{
			FORj(i+1,N+1,1)
			{
				dst[j][i] = dst[i][j] = abs(data[i]-data[j]);
			}
		}

		//
		FORi(1,N+1,1)
		{
			DP[i][i][0] = DP[i][i][1] = INFi;
			if(data[i]==L)
			{
				DP[i][i][0] = DP[i][i][1] = 0;
			}
			//printf("%d ",data[i]);
		}
		//printf("\n");
		

		FORi(2,N+1,1)
		{
			FORj(1,N+1,1)
			{
				// j j+i-1;
				if(j+i-1>N) break;
				l = j;
				r = j+i-1;
				DP[l][r][0] = TMIN( DP[l+1][r][0]+dst[l][l+1]*(l+N-r), DP[l+1][r][1]+dst[l][r]*(l+N-r));
				DP[l][r][1] = TMIN( DP[l][r-1][0]+dst[l][r]*(l+N-r), DP[l][r-1][1]+dst[r-1][r]*(l+N-r));
			}
		}

		int ans = TMIN( DP[1][N][0],DP[1][N][1]);
		printf("%d\n",ans);
	}
	return 0;
}

posted @ 2011-04-11 23:44  AC2012  阅读(233)  评论(0)    收藏  举报