1.图的表示
对于一个图(graph)G=(V,E)由顶点集V(vertex)和边集E(edges)组成。每一条边就是一个点对(u,w),其中u、w属于V。
边:单向边,双向边
权值:有点权,边权
2.一般做法
建图+算法
3.建边
建边有三种方法:
1.邻接矩阵 a[maxn][maxn]
int main()
{
memset(a,0,sizeof(a));
int u,v,val;
//若u,v之间有权值为val的双向边
a[u][v]=a[v][u]=val;
//u->v单向边,无权值
a[u][v]=1;
}
2.stl:vector存边
没有权值vector存int ,有全权值的话要用结构体
struct node
{
int to,val;
};
vector<node>v1[maxn];//储存双向边
vector<int>v[maxn];//储存单向边
int main()
{
int n,x,y,val;
for(int i=1;i<=n;i++) v[i].clear();
//x->y 单向边没有权值
v[x].push_back(y);
//x<->y 双向边有权值val,要加两次
v1[x].push_back(node{y,val});
v1[y].push_back(node{x,val});
}
3.邻接表(链式向前星)
较为抽象,但是也是用的最多的。
struct node
{
int to,next,val;
}edge[maxn*2];
void init()
{
sign=0;
for(int i=0;i<=n;i++)
head[i]=-1;
}
void add(int u,int v,int val)
{
edge[sign].to=v;
edge[sign].next=head[u];
edge[sign].val=val;
head[u]=sign++;
}
最短路算法:
1.dij 复杂度O(n^2)
#include<bits/stdc++.h>
#define mem(a,b) memset(a,b,sizeof(a))
#define inf 0x3f3f3f3f
using namespace std;
int mp[105][105],vis[105],d[105],n,m;
void init()
{
mem(vis,0);
for(int i=2; i<=n; i++)
d[i]=inf;
d[1]=0;
}
void dij(int start)
{
init();
int i,j,pos,x=0;
for(i=1; i<=n; i++)
{
pos=-1;
for(j=1; j<=n; j++)
{
if(!vis[j])
{
if(pos==-1||d[j]<d[pos])
pos=j;
}
}
if(d[pos]==inf)
break;
vis[pos]=1;
for(j=1; j<=n; j++)
{
if(!vis[j])
{
x=d[pos]+mp[pos][j];
d[j]=min(x,d[j]);
}
}
}
}
int main()
{
while(scanf("%d %d",&n,&m)&&(n+m))
{
int a,b,c,i,j;
for (i=1; i<=n; i++)
{
for (j=1; j<=n; j++)
mp[i][j]=inf;
mp[i][i]=0;
}
for(i=1; i<=m; i++)
{
scanf("%d%d%d",&a,&b,&c);
mp[a][b]=c;
mp[b][a]=c;
}
dij(1);
printf("%d\n",d[n]);
}
return 0;
}
vector版本的
#include<bits/stdc++.h>
#define mem(a,b) memset(a,b,sizeof(a))
#define inf 0x3f3f3f3f
using namespace std;
int vis[105],d[105],n,m;
struct node
{
int pre,val;
};
vector<node>v[105];
void init()
{
mem(vis,0);
for(int i=2; i<=n; i++)
d[i]=inf;
d[1]=0;
}
void dij(int start)
{
init();
int i,j,pos;
for(i=1; i<=n; i++)
{
pos=-1;
for(j=1; j<=n; j++)
{
if(!vis[j])
{
if(pos==-1||d[j]<d[pos])
pos=j;
}
}
if(d[pos]==inf)
break;
vis[pos]=1;
for(j=0; j<v[pos].size(); j++)
{
int pre=v[pos][j].pre;
int val=v[pos][j].val;
if(d[pos]+val<d[pre])
d[pre]=d[pos]+val;
}
}
}
int main()
{
while(scanf("%d %d",&n,&m)&&(n+m))
{
int a,b,c,i;
for (i=1; i<=n; i++)
v[i].clear();
for(i=1; i<=m; i++)
{
scanf("%d%d%d",&a,&b,&c);
v[a].push_back(node{b,c});
v[b].push_back(node{a,c});
}
dij(1);
printf("%d\n",d[n]);
}
return 0;
}
复杂度O(mlogn)优先队列+dij
#include<bits/stdc++.h>
#define mem(a,b) memset(a,b,sizeof(a))
#define inf 0x3f3f3f3f
using namespace std;
int D[105][105],d[105],n,m;//D用来储存任意两点之间的权值 d储存的是i点到起点的最短距离
struct node
{
int num,val;
} nod;
priority_queue<node>qq;
bool operator<(node a,node b)//使该优先队列从小到大进行输出
{
if(a.val==b.val)
return a.num>b.num;
return a.val>b.val;//a的优先级更高 先输出
}
int main()
{
while(scanf("%d %d",&n,&m)&&(n+m))
{
int a,b,c,i;
mem(D,-1);
for(i=2; i<=n; i++)
d[i]=inf;
d[1]=0;
while(!qq.empty())
qq.pop();//将队列清空进行初始化
nod.num=1;
nod.val=0;
qq.push(nod);//先将起点放入队列之中
while(m--)
{
scanf("%d%d%d",&a,&b,&c);
D[a][b]=D[b][a]=c;
}
while(!qq.empty())
{
for(i=2; i<=n; i++) //遍历所有节点,找和起点距离最短的节点并将其加入队列
{
if(D[qq.top().num][i]!=-1&&d[i]>D[qq.top().num][i]+d[qq.top().num])//将栈顶的点弹出来寻找与它距离最短的节点
{
d[i]=D[qq.top().num][i]+d[qq.top().num];
nod.num=i;
nod.val=d[i];
qq.push(nod);//找到这个节点将该点加入栈顶
}
}
qq.pop();
}
printf("%d\n",d[n]);
}
return 0;
}
2 .SPFA 期望复杂度O(E)
可以解决边权有负值的,最长路等等dij无法处理的。
#include<bits/stdc++.h>
using namespace std;
const int inff=0x3f3f3f3f;
const int maxn=105;
int d[maxn];//d[i]表示i点到起点的最短距离
int vis[maxn];//标记是否走过
int n,m;
struct node
{
int to,val;
};
vector<node>v[maxn];
void spfa(int start)
{
for(int i=1;i<=n;i++) d[i]=inff;
memset(vis,0,sizeof(vis));
d[start]=0;
queue<int> q;
q.push(start);
vis[start]=1;
while(!q.empty())
{
int u=q.front();
q.pop();
vis[u]=0;
for(int i=0;i<v[u].size();i++)
{
int e=v[u][i].to;
int val=v[u][i].val;
if(d[e]>d[u]+val)
{
d[e]=d[u]+val;
if(!vis[e])
{
q.push(e);
vis[e]=1;
}
}
}
}
}
int main()
{
int x,y,val;
while(cin>>n>>m,n+m)
{
for(int i=1;i<=n;i++) v[i].clear();
while(m--)
{
scanf("%d %d %d",&x,&y,&val);
v[x].push_back(node{y,val});
v[y].push_back(node{x,val});
}
spfa(1);
printf("%d\n",d[n]);
}
return 0;
}
邻接表版:
const int maxn=105;
const int maxm=2e5+4;
struct node
{
int to,p,val;
}edge[maxm];
int d[maxn],vis[maxn],head[maxn];
int n,m,sign;
void add(int u,int v,int val)
{
edge[sign]=node{v,head[u],val};
head[u]=sign++;
}
void init()
{
memset(head,-1,sizeof(head));
memset(d,inff,sizeof(d));
}
void spfa()
{
d[1]=0;
queue<int>q;
q.push(1);
while(!q.empty())
{
int u=q.front();
q.pop();
vis[u]=0;
for(int i=head[u];~i;i=edge[i].p)
{
int v=edge[i].to;
int val=edge[i].val;
if(d[v]>d[u]+val)
{
d[v]=d[u]+val;
if(!vis[v])
{
q.push(v);
vis[v]=1;
}
}
}
}
}
int main()
{
int x,y,z;
while(cin>>n>>m,n+m)
{
init();
while(m--)
{
scanf("%d%d%d",&x,&y,&z);
add(x,y,z);
add(y,x,z);
}
spfa();
printf("%d\n",d[n]);
}
return 0;
}
3.flyod 复杂度O(n^3)
#include<bits/stdc++.h>
using namespace std;
const int inff=0x3f3f3f3f;
const int maxn=105;
int a[maxn][maxn];
int n,m;
int main()
{
int x,y,z;
while(cin>>n>>m,n+m)
{
for(int i=1;i<=n;i++)
{
a[i][i]=0;
for(int j=i+1;j<=n;j++)
a[i][j]=a[j][i]=inff;
}
for(int i=1;i<=m;i++)
scanf("%d%d%d",&x,&y,&z),a[x][y]=a[y][x]=z;
for(int k=1;k<=n;k++)
for(int i=1;i<=n;i++)
for(int j=1;j<=n;j++)
a[i][j]=min(a[i][j],a[i][k]+a[k][j]);
cout<<a[1][n]<<endl;
}
return 0;
}