HDU 1227

Fast Food

动态规划

sum【i】【j】表示懂第i个超市到第j个超市建立一个仓库的最短路程，显然

sum【i】【j】 = sum【i】【j】+=abs(a【k】-a【mid】);(i<=k<=j)

dp【i】【j】表示前i个超市简历j个仓库的最短路程

状态转移方程dp【i】【j】 = min{ dp【k】【j-1】+sum【k+1】【i】}(1<=k<=i-1)(其实k从j-1开始就ok，因为j-1个仓库至少要给j-1个超市供货，才能保证路程最短除了仓库很多)

#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#include <iostream>
#include<string>
#include <queue>
#include <stack>
#include <vector>
#define MAXN 100005
#define max(a,b) a>b?a:b
#define min(a,b) a<b?a:b
#define MAXX 99999
using namespace std;
int a[201];
int dp[201][31];
int sum[201][201];
void init(int n,int m)
{
    for(int i=1;i<=n;i++)
        for(int j=i;j<=n;j++)
        {
            if(i==j) sum[i][j] = 0;
            else
            {
                int mid = (i+j)/2;
                sum[i][j] = 0;
                for(int k=i;k<=j;k++)
                {
                    sum[i][j]+=abs(a[k]-a[mid]);
                }
            }
        }

    //memset(dp,99,sizeof(dp));
    //cout<<dp[3][1]<<endl;

    for(int i=1;i<=n;i++)
        for(int j=1;j<=m;j++)
        {
            if(j>=i) dp[i][j] = 0;
            else if(j==1) dp[i][j] = sum[1][i];
            else
            {


            int minn=dp[i-1][j-1];
            for(int k=1;k<=i-1;k++)
            {
                if(minn>dp[k][j-1]+sum[k+1][i])
                    minn = dp[k][j-1]+sum[k+1][i];
            }
            dp[i][j] = minn;
            }
        }
}
int main()
{
    int n,m;
    int casee = 1;
    while(cin>>n>>m)
    {
        if(!(n||m)) break;
        for(int i=1;i<=n;i++)
            cin>>a[i];
        sort(a+1,a+n+1);///为了保险起见最好先排序
        init(n,m);
        printf("Chain %d\n",casee++);
        printf("Total distance sum = %d\n\n",dp[n][m]);
    }
    return 0;
}