BZOJ 2933([Poi1999]地图-区间Dp)
区间Dp
w[i][j]表示i到j划为一段的代价。。。
f[i][j]表示前i个用j个颜色的mincost
#include<cstdio>#include<cstdlib>#include<cstring>#include<iostream>#include<algorithm>#include<functional>#include<cmath>#include<cctype>#include<cassert>#include<climits>using namespace std;#define For(i,n) for(int i=1;i<=n;i++)#define Rep(i,n) for(int i=0;i<n;i++)#define Fork(i,k,n) for(int i=k;i<=n;i++)#define ForD(i,n) for(int i=n;i;i--)#define Forp(x) for(int p=pre[x];p;p=next[p])#define RepD(i,n) for(int i=n;i>=0;i--)#define MEM(a) memset(a,0,sizeof(a))#define MEMI(a) memset(a,127,sizeof(a))#define MEMi(a) memset(a,128,sizeof(a))#define INF (2139062143)#define F (1000000009)#define MAXN (3000+10)#define MAXM (10+10)typedef long long ll;int n,m;ll a[MAXN],f[MAXN][MAXM]={0},w[MAXN][MAXN]={0};int main(){// freopen("bzoj2933.in","r",stdin); scanf("%d%d",&n,&m); For(i,n) scanf("%lld",&a[i]); sort(a+1,a+1+n); For(j,n) { ForD(i,j-1) { w[i][j]=w[i+1][j]+a[((i+1)+j)/2]-a[i]; } } Rep(i,n+1) Rep(j,m+1) f[i][j]=INF; f[0][0]=0; For(i,n) For(j,m) { f[i][j]=INF; Rep(k,i) { f[i][j]=min(f[i][j],f[k][j-1]+w[k+1][i]); } } cout<<f[n][m]<<endl;// while(1); return 0;}