图论

#图论

图的储存

1.邻接矩阵

2.邻接表

struct tyc{int x,t,l,ne;} e[M];
int v[N];
void put (int x,int y,int l){
    num++;e[num]={x,y,l,v[x]};
    v[x]=num;
}

//枚举点x的所有出边
for(int i=v[x];i;i=e[i].ne) i是出边

---

最短路径

floyd

//复杂度O(n^3)
for(int k=1;k<=n;k++)
    for(int i=1;i<=n;i++)
        for(int j=1;j<=n;j++)
            dis[i][j]=min(dis[i][k]+dis[k][j],dis[i][j]);

dijkstra

需要注意的是：dijkstra无法用于处理带负权值的图

dijkstra的主要思路是分蓝白点，确定最短路长度的是白点，蓝点反之。
复杂度是相对比较优秀的O(n^2)；

	int S=1;
	memset(vis,0,sizeof(vis));
	for(int i=1;i<=n;i++){
	    if(dis[S][i]!=-1) d[i]=dis[S][i];
		else d[i]=inf;
	}
	int k=0;
	vis[S]=true;
	d[0]=inf;
	for(int i=1;i<n;i++){
	    int mn=inf,k=0;
	    for(int j=1;j<=n;j++)
	        if(!vis[j] && d[j]<=mn)  mn=d[j],k=j;//用变量mn筛出此时距离S最近的点
	    vis[k]=1;
	    for(int j=1;j<=n;j++){
	        if(!vis[j] && dis[k][j]!=-1)
	            d[j]=min(d[j],dis[k][j]);
	    }
	}

上述代码，如果加上常数，其复杂度应是O(2n^2)。
值的注意的是，通过堆进行优化的dijkstra复杂度更加优秀，可以达到O(n(logn+n))，当数据很大的时候，会比上述复杂度优秀不少。

//堆优化后的dij
struct tyc {
	int id,v;//v 表示上文中的d[x]
	friend bool operator < (tyc a,tyc b) {
		return a.v>b.v;
	} 
};

priority_queue<tyc> que;

void Dijkstra(){
	int S=1;
	memset(vis,0,sizeof(vis));
	for(int i=1;i<=n;i++){
	    if(dis[S][i]!=-1) d[i]=dis[S][i];
		else d[i]=inf;
	}
	que.push((tyc){1,0});
    rep(i,1,n){
        que.push((tyc){i,d[i]});
    }
	int x,y;
    d[0]=inf;
	while(!que.empty()){
		x = que.top().id, y = que.top().v;
		que.pop();
		if(vis[x]) continue;
		vis[x]=1;
		for(int j=1;j<=n;j++) {
			if(vis[j]) continue;
			if(dis[x][j]!=-1 && dis[x][j] < d[j]) {
				d[j] = dis[x][j];
				que.push((tyc){j,d[j]});
			}
		}
	}
}

---

最小生成树

Prim

基本思路与Dijkstra相同

int Prim(){	
	int ans=0;
	int S=1;
	memset(vis,0,sizeof(vis));
	for(int i=1;i<=n;i++){
	    if(dis[S][i]!=-1) d[i]=dis[S][i];
		else d[i]=inf;
	}
	int k=0;
	vis[S]=true;
	d[0]=inf;
	for(int i=1;i<n;i++){
	    int mn=inf,k=0;
	    for(int j=1;j<=n;j++)
	        if(!vis[j] && d[j]<=mn)  mn=d[j],k=j;//用变量mn筛出此时距离S最近的点
	    vis[k]=1;
	    ans+=d[k];//!!!!!!!!!!!
	    for(int j=1;j<=n;j++){
	        if(!vis[j] && dis[k][j]!=-1)
	            d[j]=min(d[j],dis[k][j]);
	    }
	}
	return ans;
}

堆优化版：

struct tyc {
	int id,v;//v 表示上文中的d[x]
	friend bool operator < (tyc a,tyc b) {
		return a.v>b.v;
	} 
};

priority_queue<tyc> que;



int Prim(){	
	int ans=0;
	int S=1;
	memset(vis,0,sizeof(vis));
	for(int i=1;i<=n;i++){
	    if(dis[S][i]!=-1) d[i]=dis[S][i];
		else d[i]=inf;
	}
	que.push((tyc){1,0});
    rep(i,1,n){
        if(d[i]==inf) continue;
        que.push((tyc){i,d[i]});
    }
	int x,y;
    d[0]=inf;
	while(!que.empty()){
		x = que.top().id, y = que.top().v;
		que.pop();
		if(vis[x]) continue;
		vis[x]=1;
        ans+=y;
		for(int j=1;j<=n;j++) {
			if(vis[j]) continue;
			if(dis[x][j]!=-1 && dis[x][j] < d[j]) {
				d[j] = dis[x][j];
				que.push((tyc){j,d[j]});
			}
		}
	}