题目链接
题目大意:
给一个N*N的矩阵,让你自己找一个矩阵(矩阵的数非1即0),并求出∑C ij*X ij(1<=i,j<=n)的最小值
刚开始在最短路的知识部分接触到这个题根本没有头绪,而且也没注意到所求矩阵的特点(非0即1,很重要)
下面分析一下所要求矩阵的三个特点:
我们可以把行理解为图中出度,列理解为入度
1.X 12+X 13+…X 1n=1 第一行挑一个为1,转化为在图中点1的出度为1
2.X 1n+X 2n+…X n-1n=1 第N列挑一个为1,转化为在图中点N的入度为1
3.∑Xki =∑Xij,转化为编号在2-N-1的点的入度==出度
便于理解,题目给了个示例:
For example, if n=4,we can get the following equality:
X 12+X 13+X 14=1
X 14+X 24+X 34=1
X 12+X 22+X 32+X 42=X 21+X 22+X 23+X 24
X 13+X 23+X 33+X 43=X 31+X 32+X 33+X 34
即为中间编号在2-n-1之间的点的入度和出度相等
题目给的矩阵我们可以理解为每个点都和所有的点(包括本身)有且只有一条权值为C(i,j)的单向边,
所以为了满足以上三个条件,这个图可以有以下两种情况:
A. 从1到N有一条最短边
B.从1出去回到1有最短边,从N出去回到N有最短边
比较一下A B两种情况哪个短选哪个
区别于普通的最短路,因为要求从1到1至少经过其他一个点的最短路,所以在SPFA算法初始化的时候,不能将点1入队列,而是将除了1之外的所有其他点入队,这样点1再跑最短路的时候必须经过其他点
# include <iostream>
# include <stdio.h>
# include <algorithm>
# include <queue>
using namespace std;
# define LL long long
# define MAXN 310///标记点的最大数目
# define MAX 0x3f3f3f3f
LL data[MAXN][MAXN];///记录所给矩阵的值,即各点之间边的权值
LL dis[MAXN];///记录各点到始点的最短路
int N;
LL SPFA(int start, int end)
{
queue<int>Q;
bool flag[MAXN];///标记点是否已在队列中
for(int i=1;i<=N;i++)
{
if(i==start)
{
dis[i]=MAX;
flag[i]=false;
}
else
{
Q.push(i);
dis[i]=data[start][i];
flag[i]=true;
}
}
while(!Q.empty())
{
int t=Q.front();
Q.pop();
flag[t]=false;
for(int i=1;i<=N;i++)
{
if(dis[i]>dis[t]+data[t][i])
{
dis[i]=dis[t]+data[t][i];
if(!flag[i])
{
Q.push(i);
flag[i]=true;
}
}
}
}
return dis[N];
}
int main()
{
while(cin>>N)
{
for(int i=1;i<=N;i++)
{
for(int j=1;j<=N;j++)
scanf("%I64d",&data[i][j]);
}
SPFA(1,N);
LL a=dis[N];
LL b=dis[1];
SPFA(N,N);
LL c=dis[N];
cout<<min(a,b+c)<<endl;
}
return 0;
}