您的位置:首页 > 其它

poj3259( spfa 和 bellman_ford)

2014-03-26 20:57 316 查看
一拿到这一题什么复杂度都没有分析,直接用floyd做,之后悲催的没有过才想起来有复杂度这回事。之后就考虑用BF来做,wa了几次,主要是因为存边的时候忘写另一半了,也忘更新了,改了之后就过了。时间还可以。之后用spfa来做,是看学姐博客上的,再加一个num[ maxn ]来判断,是否加入队列超过n次,若超过则有负环,但这种用我的代码来实现的话,比较慢。下面给出代码:

bellman_ford:

#include<cstring>
#include<iostream>
#include<cstdio>
#include<algorithm>
using namespace std;
const int INF =0x3f3f3f3f;
const int maxl=5200;
const int maxn=505;

int u[maxl],v[maxl],w[maxl];
int dis[maxn];

void init()
{
memset(u,0,sizeof(u));
memset(v,0,sizeof(v));
memset(w,0,sizeof(w));
}

bool bellman_ford(int n,int m)
{
memset(dis,INF,sizeof(dis));
dis[1]=0;
for(int k=0;k<n;k++){
bool finish=true;
for(int i=0;i<m;i++){
int x=u[i],y=v[i],len=w[i];
if(w[i]<INF&&dis[x]<INF&&dis[y]>dis[x]+w[i]){
dis[y]=dis[x]+w[i];
finish=false;
}
}
if(finish)break;
}
for(int i=0;i<m;i++){
int x=u[i],y=v[i];
if(dis[y]>dis[x]+w[i])
return true;
}
return false;
}

int main()
{
//freopen("int.txt","r",stdin);
int F;
scanf("%d",&F);
while(F--){
init();
int N,M,W;
scanf("%d%d%d",&N,&M,&W);
int from ,to, len,cnt=0;
while(M--){
scanf("%d%d%d",&from ,&to ,&len);
u[cnt]=from ,v[cnt]=to, w[cnt]=len;
cnt++;
u[cnt]=to ,v[cnt]=from, w[cnt]=len;
cnt++;
}
while(W--){
scanf("%d%d%d",&from,&to ,&len);
len=-len;
u[cnt]=from ,v[cnt] =to,w[cnt] =len;
cnt++;
}
if(bellman_ford(N,cnt))printf("YES\n");
else printf("NO\n");
}
return 0;
}


spfa  :

#include<cstring>
#include<cstdio>
#include<iostream>
#include<queue>
#include<algorithm>
using namespace std;
const int maxn=505;
const int INF=0x3f3f3f3f;

int dis[maxn],vis[maxn],num[maxn];
int g[maxn][maxn];

queue<int > que;

void init( )
{
while( !que.empty() ) que.pop();
memset(g,INF,sizeof(g));
for(int i=0;i<maxn;i++)g[i][i]=0;
}

bool spfa(int n)
{
memset(dis,INF,sizeof(dis));
memset(vis,0,sizeof(dis));
memset(num,0,sizeof(num));
dis[1]=0;
vis[1]=1;
num[1]++;
que.push(1);
while(!que.empty()){
int x=que.front();
vis[x]=0; //注意反向标记
que.pop();
for(int y=1;y<=n;y++){
if(dis[x]<INF&&g[x][y]<INF&&dis[y]>dis[x]+g[x][y]){
dis[y]=dis[x]+g[x][y];
if(!vis[y]){
que.push(y);
vis[y]=1;
num[y]++;
if(num[y]>n)return true;
}
}
}
}
return false;
}

int main()
{
//freopen("int.txt","r",stdin);
int F;
scanf("%d",&F);
while(F--){
init();
int N,M,W;
scanf("%d%d%d",&N,&M,&W);
int from ,to, len;
while(M--){
scanf("%d%d%d",&from ,&to ,&len);
if(g[from][to]>len)
g[from][to]=g[to][from]=len;
}
while(W--){
scanf("%d%d%d",&from,&to ,&len);
len=-len;
if(g[from][to]>len)
g[from][to]=len;
}
if(spfa(N))printf("YES\n");
else printf("NO\n");
}
return 0;
}
内容来自用户分享和网络整理,不保证内容的准确性,如有侵权内容,可联系管理员处理 点击这里给我发消息