码迷,mamicode.com
首页 > 其他好文 > 详细

CODEVS_2800 送外卖 状态压缩+动态规划

时间:2014-12-26 20:22:01      阅读:260      评论:0      收藏:0      [点我收藏+]

标签:动态规划   dp   状态压缩   

原题链接:http://codevs.cn/problem/2800/

题目描述 Description

有一个送外卖的,他手上有n份订单,他要把n份东西,分别送达n个不同的客户的手上。n个不同的客户分别在1~n个编号的城市中。送外卖的从0号城市出发,然后n个城市都要走一次(一个城市可以走多次),最后还要回到0点(他的单位),请问最短时间是多少。现在已知任意两个城市的直接通路的时间。

输入描述 Input Description

第一行一个正整数n (1<=n<=15)

接下来是一个(n+1)*(n+1)的矩阵,矩阵中的数均为不超过10000的正整数。矩阵的i行j列表示第i-1号城市和j-1号城市之间直接通路的时间。当然城市a到城市b的直接通路时间和城市b到城市a的直接通路时间不一定相同,也就是说道路都是单向的。

输出描述 Output Description

一个正整数表示最少花费的时间

样例输入 Sample Input
3
0 1 10 10
1 0 1 2
10 1 0 10
10 2 10 0
样例输出 Sample Output

8

数据范围及提示 Data Size & Hint

1<=n<=15

这道题很容易想到状压dp,首先跑一遍floyd求出各个点之间的最短路,然后dp。

状态是dp[i][j]表示状态为i(二进制,表示是否访问过每个点),在位置j时的最短路。

转移就是:dp[i][j]=min(dp[i-(1<<v)][u]+grid[u][v],dp[i][j]),其中v是当前位置,u是上一个状态的位置。

需要注意的是dp的顺序应该是由含1的个数少的二进制到1的个数高的二进制;由于最开始就在0位置,所以dp[(1<<n)-1][0]不可能被转移到,所以最后的答案应该是ans=min(ans,dp[(1<<n)-1][i]),其中0<=i<n。

详见代码:

#include<iostream>
#include<cstring>
#include<vector>
#include<queue>
#include<algorithm>
#include<cstdio>
#include<climits>
#define INF INT_MAX
#define MAX_S 1<<16
#define MAX_N 17
using namespace std;

struct node
{
    int val,num;
    node(int v,int n)
        {val=v;num=n;}
    node(){}
};

node one[MAX_S];
int grid[MAX_N][MAX_N];
int n;
int dp[MAX_S][MAX_N];
int numOfOne(int x)
{
    int res=0;
    for(int i=0;i<n;i++)
        if((x>>i)&1)res++;
    return res;
}

void floyd()
{
    for(int i=0;i<n;i++)
        for(int j=0;j<n;j++)
            for(int k=0;k<n;k++)
                grid[i][j]=min(grid[i][j],grid[i][k]+grid[k][j]);
}

bool cmp(node a,node b)
{
    return a.num<b.num;
}

int main()
{
    scanf("%d",&n);
    n++;
    for(int i=0;i<n;i++)
        for(int j=0;j<n;j++)
            scanf("%d",&grid[i][j]);
    int totS=1<<n;
    for(int i=0;i<totS;i++)one[i]=node(i,numOfOne(i));

    for(int i=0;i<totS;i++)
        for(int j=0;j<n;j++)
            dp[i][j]=INF;

    floyd();
    sort(one,one+totS,cmp);
    dp[1][0]=0;
    for(int i=0;i<totS;i++)
    {
        int s=one[i].val;
        for(int v=0;v<n;v++)
        {
            if(!((s>>v)&1))continue;
            int t=s-(1<<v);
            int tmp=INF;
            for(int u=0;u<n;u++)
            {
                if((!((t>>u)&1))||dp[t][u]==INF)continue;
                tmp=min(tmp,dp[t][u]+grid[u][v]);
            }
            if(tmp!=INF)
                dp[s][v]=tmp;
        }
    }
    int ans=INF;
    for(int i=0;i<n;i++)
        ans=dp[totS-1][i]==INF?ans:min(ans,dp[totS-1][i]+grid[i][0]);
    cout<<ans<<endl;
    return 0;
}


CODEVS_2800 送外卖 状态压缩+动态规划

标签:动态规划   dp   状态压缩   

原文地址:http://blog.csdn.net/harryguo2012/article/details/42175559

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!