Hdu4370-0or1(最短路!)

题意:

在这里插入图片描述

思路:

这个题是最短路!!看到题是一脸迷茫啊。。
请看大神的题解:https://www.cnblogs.com/quasar/p/5140558.html
显然,题目给的是一个0/1规划模型。

解题的关键在于如何看出这个模型的本质。

3个条件明显在刻画未知数之间的关系,从图论的角度思考问题,容易得到下面3个结论:

1.X12+X13+…X1n=1 于是1号节点的出度为1

2…X1n+X2n+…Xn-1n=1 于是n号节点的入度为1

3.∑Xki =∑Xij 于是2~n-1号节点的入度必须等于出度

于是3个条件等价于一条从1号节点到n号节点的路径,故Xij=1表示需要经过边(i,j),代价为Cij。Xij=0表示不经过边(i,j)。注意到Cij非负且题目要求总代价最小,因此最优答案的路径一定可以对应一条简单路径。

最终,我们直接读入边权的邻接矩阵,跑一次1到n的最短路即可,记最短路为path。

以上情况设为A

非常非常非常非常非常非常非常非常抱歉,简单路径只是充分条件,但不必要。(对造成困扰的队伍深表歉意)

漏了如下的情况B:

从1出发,走一个环(至少经过1个点,即不能是自环),回到1;从n出发,走一个环(同理),回到n。

容易验证,这是符合题目条件的。且A || B为该题要求的充要条件。

由于边权非负,于是两个环对应着两个简单环。

因此我们可以从1出发,找一个最小花费环,记代价为c1,再从n出发,找一个最小花费环,记代价为c2。(只需在最短路算法更新权值时多加一条记录即可:if(i==S) cir=min(cir,dis[u]+g[u][i]))

故最终答案为min(path,c1+c2)

#include <cstdio>
#include <cstring>
#include <queue>
#include <algorithm>
#pragma comment(linker, "/STACK:102400000,102400000")
typedef long long LL;
const int INF = 0x3f3f3f3f;
const int maxn = 300+5;
using namespace std;
int n,m;
int G[maxn][maxn];

int head[maxn], cur;
struct Edge{
	int v,w,next;
};
Edge edges[100000];
void init(int a){
	for(int i = 0; i <= a; ++i) head[i] = -1;
	cur = 0;
}
void addEdge(int u, int v, int w){
	edges[cur].v = v; edges[cur].w = w;
	edges[cur].next = head[u]; head[u] = cur++;
}

int d[maxn];
bool vis[maxn];
void BellmanFord(int s){
	queue<int> Q;
	for(int i = 1; i <= n; ++i) { 
		if(i == s){
			d[i] = INF; vis[i] = 0;
		}
		else if(G[s][i] != INF){
			d[i] = G[s][i]; vis[i] = 1; Q.push(i);
		}
		else{
			d[i] = INF; vis[i] = 0;
		}
	}
	
	while(!Q.empty()){
		int u = Q.front(); Q.pop();
		vis[u] = 0;
		for(int v = 1; v <= n; ++v){
			if(d[v] > d[u] + G[u][v]){
				d[v] = d[u] + G[u][v];
				if(!vis[v]){ Q.push(v); vis[v] = 1; }
			}
		}
		
	}
}

int main()
{
    //freopen("in.txt","r",stdin);
    
	while(scanf("%d",&n) == 1){
		init(n);
		for(int i = 1; i <= n; ++i)
		for(int j = 1; j <= n; ++j) scanf("%d", &G[i][j]);
		
		BellmanFord(1);
		int ans = d[n];
		int loop1 = d[1];
		
		BellmanFord(n);
		int loopn = d[n];
		
		ans = min(ans, loop1+loopn);
		printf("%d\n", ans);
	}
    fclose(stdin);
	return 0;
}


  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值