题意:
思路:
这个题是最短路!!看到题是一脸迷茫啊。。
请看大神的题解: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;
}