最优贸易 trade Dijkstra解法

#include <stdio.h>
#define MAXINT 2139062143
struct edge{
int dest;//下一结点
int value;//权值
struct edge *link;//指针
};
struct edge *G;
int Initial()//初始化图
{
int n,m;//n个结点,m条边
FILE *fp;
int *node;
int i,u,v,flag;
int value;
struct edge *l;
fp=fopen("trade.in","r");
fscanf(fp,"%d %d",&n,&m);
G=(struct edge *)malloc(sizeof(struct edge)*(n+1));//分配n个结点存储空间 node=(int *)malloc(sizeof(int)*(n+1));
memset(G,0,sizeof(struct edge)*(n+1));
if(G==NULL)
return 0;
for(i=1;i<=n;i++)
fscanf(fp,"%d",&node[i]);
for(i=1;i<=m;i++)
{
l=(struct edge *)malloc(sizeof(struct edge));
fscanf(fp,"%d %d %d",&u,&v,&flag);
if(flag==1)
value=node[v]-node[u]>0?(node[v]-node[u]):0;
if(flag==2)
value=(int)abs(node[u]-node[v]);
l->dest=v;
l->value=value;
l->link=G[u].link;
G[u].link=l;
}
fclose(fp);
return n;
}
void Dijkstra(int n,int s)
{
int *S,*D,*P;//D记录了从s到其他顶点的最短长度,P记录了所求到的最短路 struct edge *l;
int i,j;
int min,w;
FILE *fp;
fp=fopen("trade.out","w");
S=(int *)malloc(sizeof(int)*(n+1));
D=(int *)malloc(sizeof(int)*(n+1));
P=(int *)malloc(sizeof(int)*(n+1));
memset(D,-1,sizeof(int)*(n+1));
memset(S,0,sizeof(int)*(n+1));//解集顶点
memset(P,-1,sizeof(int)*(n+1));
S[s]=1;//将顶点s加入到路径中
D[s]=0;//顶点到达自身为0
l=G[s].link;//以点s为始点
//初始化D和P
while(l)
{
D[l->dest]=l->value;//到相邻结点的路径长度为权值
P[l->dest]=s;//s为l->dset的前驱
l=l->link;
}
for(i=1;i<n;i++)//每次找出一个不在S中的顶点最短路
{
min=-1;
for(j=1;j<=n;j++)//找出不在S中的D值可达且最小的顶点,记为w
{
if(S[j]==0&&D[j]>=0)//结点j不在S中且可达
{
if(min<0||D[j]>min)
{
min=D[j];
w=j;
}
}
}
if(min<0) break;
S[w]=1;//将w加入S中,并更新所有由它可直接到达的顶点的D值
l=G[w].link;
while(l)
{
if(D[l->dest]<0||D[w]+l->value>D[l->dest])
{
if(D[w]>l->value)
D[l->dest]=D[w];
else
D[l->dest]=l->value;
P[l->dest]=w;
}
l=l->link;
}
}
for(i=1;i<=n;i++)
{
if(D[i]<0)
{
printf("Can't reach %d from %d\n",i,s);
continue;
}
printf("The shortest path from %d to %i is %d\n",s,i,D[i]); j=i;
do{
printf("%d ",j);
j=P[j];
}while(j>=0);
printf("\n");
}
fprintf(fp,"%d",D[n]);
fclose(fp);
}
int main()
{
int n;
n=Initial();
Dijkstra(n,1);
return 0;
}。