上一次讲了两种关于图的求最短路的方法,今天来讲讲剩下的两种:
三、Bellman-Ford算法(搜边法)
此算法是将所有输入的已知边在已知起点和终点的情况下(int w[ v ][ 3 ] (v表示总边数,w[ i ][ 0 ]表示此边的起点,w[ i ][ 1 ]表示此边的终点,w[ i ][ 2 ]表示权值) ),进行点数次遍历,由此来重置从dis[ i ](起点到第 i 个点的长度),样例代码:
for ( i=1 ; i<=n(点数) ; i++ )
for ( j=1 ; j<=v(边数) ; j++ )
dis[ w[ j ][ 1 ] ] = min( dis[ w[ j ][ 1 ] ] , dis[ w[ j ][ 0 ] ] + w[ j ][ 2 ] ) ;
例题还是待会再说吧……
四、SPFA算法(广搜盗版)(不懂广搜的童鞋来看看<( ̄3 ̄)>——>点击打开链接)
相当于广搜的思想,将被重置过的点入队(如果已经在队中,则需使用bool数组),将此点能延伸到的并能重置的点入队,直到队空为止(注意:此算法无法处理负边权!)
下面来看一道例题:
香甜的黄油(butter.cpp)
题目描述
农夫John发现做出全威斯康辛州最甜的黄油的方法:糖。
把糖放在一片牧场上,他知道N(1<=N<=500)只奶牛会过来舔它,这样就能做出能卖好价钱的超甜黄油。当然,他将付出额外的费用在奶牛上。 农夫John很狡猾。像以前的Pavlov,他知道他可以训练这些奶牛,让它们在听到铃声时去一个特定的牧场。他打算将糖放在那里然后下午发出铃声,以至他可以在晚上挤奶。 农夫John知道每只奶牛都在各自喜欢的牧场(一个牧场不一定只有一头牛)。给出各头牛在的牧场和牧场间的路线,找出使所有牛到达的路程和最短的牧场(他将把糖放在那)
输入
第1行: 3个数:奶牛数N,牧场数(2<=P<=800),牧场间道路数C(1<=C<=1450)
第2行到第N+1行: 1到N头奶牛所在的牧场号
第N+2行到第N+C+1行: 每行有三个数:相连的牧场A、B,两牧场间距离D(1<=D<=255),当然,连接是双向的
输出
一行 输出奶牛必须行走的最小的距离和
样例输入
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
|
|
5 7 3
|
| C3
C2 @--5--@
P3 P4
}
样例输出
8
{说明:
放在4号牧场最优
}
不说了,直接丢代码:
Bellman-Ford算法:
#include<cstdio>
#include<cmath>
#include<cstring>
#include<algorithm>
using namespace std;
int n,cow[501],way[1451][3],m,nn;
void scan()
{
int i;
scanf("%d%d%d",&nn,&n,&m);
for(i=1;i<=nn;i++)
scanf("%d",&cow[i]);
for(i=1;i<=m;i++)
scanf("%d%d%d",&way[i][0],&way[i][1],&way[i][2]);
}
int bellman(int b)
{
int i,j,s=0;
int p[801];
memset(p,127,sizeof(p));
p[b]=0;
for(i=1;i<n;i++)
for(j=1;j<=m;j++)
{
p[way[j][1]]=min(p[way[j][1]],p[way[j][0]]+way[j][2]);
p[way[j][0]]=min(p[way[j][0]],p[way[j][1]]+way[j][2]);
}
for(i=1;i<=nn;i++)
s+=p[cow[i]];
return s;
}
int main()
{
int i;
int k=1<<30;
scan();
for(i=1;i<=n;i++)
k=min(k,bellman(i));
printf("%d",k);
}
SPFA算法:
#include<cstdio>
#include<cmath>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
int cow[501],point[801][801],l[1450][1450],num,n,m;
void scan()
{
int i,x,y;
scanf("%d%d%d",&num,&n,&m);
for(i=1;i<=num;i++)
scanf("%d",&cow[i]);
for(i=1;i<=m;i++)
{
scanf("%d%d",&x,&y);
point[x][++point[x][0]]=y;
point[y][++point[y][0]]=x;
scanf("%d",&l[x][y]);
l[y][x]=l[x][y];
}
}
int spfa(int f)
{
int p[801];
bool v[801]={0};
int h[100001],t=0,w=1,i;
memset(p,127,sizeof(p));
p[f]=0;
h[1]=f;
int s=0;
do
{
v[h[t]]=0;
t++;
for(i=1;i<=point[h[t]][0];i++)
{
if(p[h[t]]+l[h[t]][point[h[t]][i]]<p[point[h[t]][i]])
{
p[point[h[t]][i]]=p[h[t]]+l[h[t]][point[h[t]][i]];
if(!v[point[h[t]][i]])
{
v[point[h[t]][i]]=1;
w++;
h[w]=point[h[t]][i];
}
}
}
}while(t<w);
for(i=1;i<=num;i++)
s+=p[cow[i]];
return s;
}
int main()
{
int i;
int k=1<<30;
memset(l,127,sizeof(l));
scan();
for(i=1;i<=n;i++)
k=min(k,spfa(i));
printf("%d",k);
}
仿照代码,自己去试一下吧~(≧▽≦)/~