【最短路】hdu2544 SPFA以及SPFA的优化

hdu2544 题意：给出N个点和M条有向边，求1到N的最短距离。

其实在此题中这个优化有可能更慢。

首先是SPFA的模版 hdu2544

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<queue>
#include<iostream>
using namespace std;

const int N=110,M=11000,INF=(int)1e6;
int n,m,len;
int first[N],d[N];
bool vis[N];

struct node{
	int x,y,d,next;
}a[M];

void ins(int x,int y,int d)
{
	len++;
	a[len].x=x;a[len].y=y;a[len].d=d;
	a[len].next=first[x];first[x]=len;
}

int main()
{
	int i,j;
	while(1)
	{    
	    scanf("%d%d",&n,&m);
	    if(!n && !m) return 0;
		len=0;
		memset(first,0,sizeof(first));
		for(i=1;i<=m;i++)
		{
		   int x,y,d;
		   scanf("%d%d%d",&x,&y,&d);
		   ins(x,y,d);
		   ins(y,x,d);			
		}
		//for(i=1;i<=len;i++)
		// printf("%d %d %d\n",a[i].x,a[i].y,a[i].d);
		d[1]=0;
		for(i=2;i<=n;i++) d[i]=INF;
        memset(vis,0,sizeof(vis));
        queue<int> q;
		q.push(1);
        while(!q.empty())
        {
        	int x=q.front();q.pop();
        	vis[x]=1;
        	for(i=first[x];i;i=a[i].next)
        	{
        		int y=a[i].y;
	        	if(d[x]+a[i].d<d[y])
	        	{
	        		d[y]=d[x]+a[i].d;
	        		if(!vis[y])
	        		{
		        		vis[y]=1;
		        		q.push(y);
		        	}
	        	}
	        }
	        vis[x]=0;
        }
        printf("%d\n",d[n]);
	}
}

SPFA有个经典的优化SLF，对于专门卡SPFA的题目可以用此优化：

用双向队列，对于每一个要加入队列的y，先判断是否比队首元素z的d值小，若d[y]<d[z]，则将y插入队首，否则插入队尾。注意要判断队列是否为空，否则没有队首元素。

//SPFA+SLF优化  hdu2544 
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<deque>
#include<iostream>
using namespace std;

const int N=110,M=11000,INF=(int)1e6;
int n,m,len;
int first[N],d[N];
bool vis[N];

struct node{
	int x,y,d,next;
}a[M];

void ins(int x,int y,int d)
{
	len++;
	a[len].x=x;a[len].y=y;a[len].d=d;
	a[len].next=first[x];first[x]=len;
}

int main()
{
	int i,j;
	while(1)
	{    
	    scanf("%d%d",&n,&m);
	    if(!n && !m) return 0;
		len=0;
		memset(first,0,sizeof(first));
		for(i=1;i<=m;i++)
		{
		   int x,y,d;
		   scanf("%d%d%d",&x,&y,&d);
		   ins(x,y,d);
		   ins(y,x,d);			
		}
		//for(i=1;i<=len;i++)
		// printf("%d %d %d\n",a[i].x,a[i].y,a[i].d);
		d[1]=0;
		for(i=2;i<=n;i++) d[i]=INF;
        memset(vis,0,sizeof(vis));
        deque<int> q;
		q.push_back(1);
        while(!q.empty())
        {
        	int x=q.front();q.pop_front();
        	vis[x]=1;
        	for(i=first[x];i;i=a[i].next)
        	{
        		int y=a[i].y;
	        	if(d[x]+a[i].d<d[y])
	        	{
	        		d[y]=d[x]+a[i].d;
	        		if(!vis[y])
	        		{
		        		vis[y]=1;
		        		if(!q.empty())
		        		{
		        			if(d[y]<d[q.front()]) q.push_front(y);
						    else q.push_back(y);
		        		}
						else q.push_back(y);
		        	}
	        	}
	        }
	        vis[x]=0;
        }
        printf("%d\n",d[n]);
	}
}