【POJ】3259 Wormholes bellman-ford | SPFA

诶，真的是越活越回去了，现在打的程序还不如以前了，慢的爆炸……（先测的是bellman-ford，后测的是SPFA）

题目大意：FJ想要进行时光旅行，现在有F个农场，对于每个农场，已知该农场有N片农田，M条小路，W个虫洞。其中小路是双向的，虫洞是单向的。对于每条小路，有S,E,T三个参数：S和E表示小咯两端农田的序号，通过这条小咯需要的T的时间。对于每个虫洞，同样有S,E,T三个参数：S表示虫洞起点所在的农田的序号，E表示虫洞终点所在的农田的序号，通过这个虫洞需要-T的时间。现在FJ想知道他在每个农场是否能回到从前。

我们从题目中可以清晰的看出，题目就是想让我们求每个农场是否存在负权环。一道大水题。

方法一：bellman-ford，用side数组来存边。若松弛次数超过n次，就说明存在负权环。

#include <cstdio>#include <vector>#include <cstring>using namespace std;struct note{	int x,y,w;}side[2710];int ti,n,m,p,x,y,w,dis[510],g;bool b;int main(void){	scanf("%d",&ti);	while (ti--){		scanf("%d%d%d",&n,&m,&p);		for (int i=1; i<=m; ++i)			scanf("%d%d%d",&side[i].x,&side[i].y,&side[i].w);		for (int i=m+1; i<=m+p; ++i)			scanf("%d%d%d",&side[i].x,&side[i].y,&side[i].w),side[i].w*=-1;		memset(dis,0x3f,sizeof(dis));		b=1,g=0;dis[1]=0;		while (b){			b=0;			if (++g>n) break;			for (int i=1; i<=m; ++i){				if (dis[side[i].y]>dis[side[i].x]+side[i].w)					dis[side[i].y]=dis[side[i].x]+side[i].w,b=1;				if (dis[side[i].x]>dis[side[i].y]+side[i].w)					dis[side[i].x]=dis[side[i].y]+side[i].w,b=1;			}			for (int i=m+1; i<=m+p; ++i)			if (dis[side[i].y]>dis[side[i].x]+side[i].w)					dis[side[i].y]=dis[side[i].x]+side[i].w,b=1;		}		if (g>n) PRintf("YES/n");			else printf("NO/n");	}	return 0;}方法二：SPFA——bellman-Ford的队列优化升级版，可惜在这题中慢的起飞……（好像是因为SPFA计算了许多不必要的路线）若一个节点入队超过n次就说明存在负权环。
#include <cstdio>#include <cstring>#include <queue>using namespace std;queue <int> que;int t,n,m1,m2,x,y,w,map[510][510],c[510],dis[510];bool b,bo[510];bool spfa(int s){	while (!que.empty()) que.pop();	memset(c,0,sizeof c);	memset(bo,0,sizeof bo);	memset(dis,0x3f,sizeof dis);	que.push(s);c[s]=1;bo[s]=1;	while (!que.empty()){		int p=que.front();que.pop();		bo[p]=0;		for (int i=1; i<=n; ++i)			if (dis[i]>dis[p]+map[p][i]){				dis[i]=dis[p]+map[p][i];				if (!bo[i]){					if (++c[i]>n) return 1;					bo[i]=1;					que.push(i);				}			}	}	return 0;}int main(void){	scanf("%d",&t);	while (t--){		scanf("%d%d%d",&n,&m1,&m2);		memset(map,0x3f,sizeof map);		for (int i=1; i<=m1; ++i){			scanf("%d%d%d",&x,&y,&w);			if (w<map[x][y]) map[x][y]=map[y][x]=w;		}		for (int i=1; i<=m2; ++i){			scanf("%d%d%d",&x,&y,&w);			map[x][y]=-w;		}		b=0;		for (int i=1; i<=n; ++i)			if (spfa(i)){				b=1;				break;			}		if (b) printf("YES/n");			else printf("NO/n");	}	return 0;}写在最后：虽然SPFA的运行速度过慢可能跟计算了许多不必要的路线有关，但是和我用的是邻接矩阵有不可分割的关系。在以后的做题过程中，尽量用高级一些的数据结构，狮子搏兔亦用全力。