Time Limit:10000MS Memory Limit:65536K
Total Submit:126 Accepted:55
Case Time Limit:1000MS
Description
农夫John发现做出全威斯康辛州最甜的黄油的方法:糖。把糖放在一片牧场上,他知道N(1<=N<=500)只奶牛会过来舔它,这样就能做出能卖好价钱的超甜黄油。当然,他将付出额外的费用在奶牛上。
农夫John很狡猾。像以前的Pavlov,他知道他可以训练这些奶牛,让它们在听到铃声时去一个特定的牧场。他打算将糖放在那里然后下午发出铃声,以至他可以在晚上挤奶。
农夫John知道每只奶牛都在各自喜欢的牧场(一个牧场不一定只有一头牛)。给出各头牛在的牧场和牧场间的路线,找出使所有牛到达的路程和最短的牧场(他将把糖放在那)
Input
第一行: 三个数:奶牛数N,牧场数(2<=P<=800),牧场间道路数C(1<=C<=1450)
第二行到第N+1行: 1到N头奶牛所在的牧场号
第N+2行到第N+C+1行: 每行有三个数:相连的牧场A、B,两牧场间距离D(1<=D<=255),当然,连接是双向的
Output
一行 输出奶牛必须行走的最小的距离和
Sample Input
3 4 5
2
3
4
1 2 1
1 3 5
2 3 7
2 4 3
3 4 5
样例图形
P2
P1 @--1--@ C1
`b \ |\
\ | \
5 7 3
\ | \
\| \ C3
C2 @--5--@
P3 P4
Sample Output
8
{说明:放在4号牧场最优 }
解题思路
这题可以采用最棒的SPFA算法
算法解析:SPFA 在形式上和广度优先搜索非常类似,不同的是广度优先搜索中一个点出了队列就不可能重新进入队列,但是SPFA中一个点可能在出队列之后再次被放入队列,也就是说一个点修改过其它的点之后,过了一段时间可能会获得更短的路径,于是再次用来修改其它的点,这样反复进行下去。
算法时间复杂度很低,不超过边长的5倍。
代码
#include<iostream>
#include<cstdio>
#include<cmath>
#include<cstring>
#include<iomanip>
using namespace std;
const int INF=0X7fffffff;
int h=0,t=1,v[100000],q[100000],n,m,c,dis[100000],head[100000],k,b[100000],maxn,ans=INF;
//dis[i]记录从起点s到i的最短路径,v记录一个点是否现在存在在队列中,q是队列
int x,y,z;
struct c{//邻接表记录连接i,j的边,让后面不用枚举所有点
int x,w,next;
}a[100000];
void bfs(int s)
{
int x;
memset(v,0,sizeof(v));
memset(dis,0X7f,sizeof(dis));
memset(q,0,sizeof(q));
v[s]=1;
q[1]=s;
dis[s]=0;
h=0;
t=1;
do{
h++;
x=q[h];//头指针向下移一位,取出指向的点u
for(int i=head[x];i;i=a[i].next)//注意不要去枚举所有点,用数组模拟邻接表存储
{
if(dis[x]+a[i].w<dis[a[i].x])//[松弛操作](https://baike.baidu.com/item/%E6%9D%BE%E5%BC%9B%E6%93%8D%E4%BD%9C/7554089?fr=aladdin).
{
dis[a[i].x]=dis[x]+a[i].w;
if(!v[a[i].x])//队列中不存在v点,v入队
{
v[a[i].x]=1;
t++;
q[t]=a[i].x;//尾指针下移一位,v入队
}
}
}
v[x]=0;
}while(h<t);
}
void add(int x,int y,int z){
++k;
a[k].x=y;
a[k].w=z;
a[k].next=head[x];
head[x]=k;
}
int main(){
scanf("%d%d%d",&n,&m,&c);
for(int i=1;i<=n;i++)
scanf("%d",&b[i]);
for(int i=1;i<=c;i++)
{
scanf("%d%d%d",&x,&y,&z);
add(x,y,z);
add(y,x,z);
}
for(int i=1;i<=m;i++)
{
maxn=0;
bfs(i);
for(int j=1;j<=n;j++)
maxn+=dis[b[j]];
if(maxn<ans)
ans=maxn;
}
cout<<ans;
}
STL简化
#include<iostream>
#include<cstdio>
#include<cmath>
#include<cstring>
#include<iomanip>
#include<queue>
using namespace std;
const int INF=0X7fffffff;
int h=0,t=1,v[100000],q[100000],n,m,c,dis[100000],head[100000],k,b[100000],maxn,ans=INF;
int x,y,z;
struct c{
int x,w,next;
}a[100000];
void bfs(int s)
{
int x;
memset(v,0,sizeof(v));
memset(dis,0X7f,sizeof(dis));
v[s]=1;
queue<int>q;
q.push(s);
dis[s]=0;
while(!q.empty())
{
x=q.front();
q.pop();
for(int i=head[x];i;i=a[i].next)
{
if(dis[x]+a[i].w<dis[a[i].x])
{
dis[a[i].x]=dis[x]+a[i].w;
if(!v[a[i].x])
{
v[a[i].x]=1;
q.push(a[i].x);
}
}
v[x]=0;
}
}
}
void add(int x,int y,int z){
++k;
a[k].x=y;
a[k].w=z;
a[k].next=head[x];
head[x]=k;
}
int main(){
scanf("%d%d%d",&n,&m,&c);
for(int i=1;i<=n;i++)
scanf("%d",&b[i]);
for(int i=1;i<=c;i++)
{
scanf("%d%d%d",&x,&y,&z);
add(x,y,z);
add(y,x,z);
}
for(int i=1;i<=m;i++)
{
maxn=0;
bfs(i);
for(int j=1;j<=n;j++)
maxn+=dis[b[j]];
if(maxn<ans)
ans=maxn;
}
cout<<ans;
}