题目链接在这里
问题描述
Given a n*n matrix Cij (1<=i,j<=n),We want to find a n*n matrix Xij (1<=i,j<=n),which is 0 or 1.
Besides,Xij meets the following conditions:
1.X12+X13+...X1n=1
2.X1n+X2n+...Xn-1n=1
3.for each i (1<i<n), satisfies ∑Xki (1<=k<=n)=∑Xij (1<=j<=n).
For example, if n=4,we can get the following equality:
X12+X13+X14=1
X14+X24+X34=1
X12+X22+X32+X42=X21+X22+X23+X24
X13+X23+X33+X43=X31+X32+X33+X34
Now ,we want to know the minimum of ∑Cij*Xij(1<=i,j<=n) you can get.
Hint:
For sample, X12=X24=1,all other Xij is 0.
Input
The input consists of multiple test cases (less than 35 case).
For each test case ,the first line contains one integer n (1<n<=300).
The next n lines, for each lines, each of which contains n integers, illustrating the matrix C, The j-th integer on i-th line is Cij(0<=Cij<=100000).
Output
For each case, output the minimum of ∑Cij*Xij you can get.
Sample Input
4
1 2 4 10
2 0 1 1
2 2 0 5
6 3 1 2
Sample Output
3
题目大意是给出一个n x n的矩阵C,让构建一个n x n的矩阵X,如果Xij为1的话,那么ans += Cij。说白了就是在矩阵C中选出几个点,求它们的和。但是这几个点不是随便找的,需要满足上面的三个条件。
这道题怎么也想不到是用最短路的方法解决的,开始看了邝斌巨犇的题解才懂。题解链接在这里。
我们把C[i][j]当做从i点到j点的边的权值,即有一条边,它的起点是i,终点是j,权值是C[i][j],那么这个题目就可以转换成最短路问题了。
我们看题目给出的第一个条件
X12+X13+X14=1
这个的意思就是,从第1个点出发,只能到2、3、4……这几个点中的其中一个,换句话说第一个点的出度为1.
第二个条件也很明了,第n个点的入度为1
第三个条件的意思就是每个点的入度等于出度
这就很好做了,找一条从点1出发到达点n的最短路不就行了呗。
当然最短路只是情况1,还有另外一种情况:
点1出发走了一圈(未经过点n)又回来了,并且,点n出发走了一圈(未经过点1)又回来了,即有两个闭环。
两种情况所花费的权值取一个最小值输出就行了。
对于闭环的情况,我们可以在初始化时设dis[start] = INF,而start点所连接的点的权值dis[v] = C[start][v],这样就可以找闭环了。
代码如下:
#include <stdio.h>
#include <string.h>
#include <algorithm>
#include <iostream>
#include <queue>
using namespace std;
const int INF = 0x3f3f3f3f;
const int MaxN = 310;
int dis[MaxN];
int cost[MaxN][MaxN];
int n;
void dijkstra(int st){
priority_queue<pair<int, int>, vector<pair<int, int> >, greater<pair<int, int> > > que;
for(int v = 1; v <= n; ++v){
if(v == st){ //这个是找闭环用的
dis[v] = INF;
}
else{
dis[v] = cost[st][v];
que.push(make_pair(dis[v], v));
}
}
while(!que.empty()){
pair<int, int> p = que.top();
que.pop();
int u = p.second;
if(dis[u] < p.first) continue;
for(int v = 1; v <= n; ++v){
if(dis[v] > dis[u] + cost[u][v]){
dis[v] = dis[u] + cost[u][v];
que.push(make_pair(dis[v], v));
}
}
}
}
int main(){
while(~scanf("%d", &n)){
for(int i = 1; i <= n; ++i)
for(int j = 1; j <= n; ++j)
scanf("%d", &cost[i][j]);
dijkstra(1);
int ans = dis[n];
int loop1 = dis[1];
dijkstra(n);
int loopn = dis[n];
printf("%d\n", min(ans, loop1 + loopn));
}
return 0;
}