0.PTA得分截图

1.本周学习总结

1.1 总结图内容

1.1.1 图的概念:

图形结构是最普遍的一类数据结构,具有广泛的实际应用。
图中数据元素的关系是多对多的关系,(线性结构一对一,树形结构一对多)

1.1.2 图的定义和基本术语:

图:顶点集 V 和顶点间的关系:边集合E组成的数据结构。

  • 图的逻辑结构描述:
    Graph = (V , E)

G=(V1,E1)
V1={A, B, C, D, E}
E1={<A,B>, <A,E>, <B,C>, <C,D>,<D,B>,<D,A>, <E,C> }

  • 有向图和无向图:
    1.有向图:
    有向图称由顶点集和弧(“弧”是有方向的边,弧可以用尖括号表示。)集构成的图。
    eg:
    有向图G1:
    G1 = (V1, E1)

    V1={A, B, C, E, F}
    E1=

2.无向图:
无向图:没方向边(若<v, w>E 必有<w, v>VR。边可以用圆括号表示。)
eg:
无向图G2:
G2 = (V2, E2)

V2={A, B, C, D, E, F}
E2={(A,B), (A,E), (B,E), (C,D), (D,F),(B,F), (C,F) }

  • 图的基本术语:
  1. 端点和邻接点:
    无向图:若存在一条边(i,j),则称顶点i和顶点j互为邻接点。
    有向图:存在一条边<i,j>,则称此边是顶点i的一条出边,同时也是顶点j的一条入边;称顶点i 和顶点j 互为邻接点。

2、顶点的度、入度和出度:
无向图:以顶点i为端点的边数称为该顶点的度。
有向图:以顶点i为终点的入边的数目,称为该顶点的入度。以顶点i为始点的出边的数目,称为该顶点的出度。一个顶点的入度与出度的和为该顶点的度。
若一个图中有n个顶点和e条边,每个顶点的度为di(0≤i≤n-1),则有:e=(d0+d1+d2+...+dn-1)/2

3、完全图:
无向图:每两个顶点之间都存在着一条边,称为完全无向图, 包含有n(n-1)/2条边。
有向图:每两个顶点之间都存在着方向相反的两条边,称为完全有向图,包含有n(n-1)条边。

4、稠密图、稀疏图:
当一个图接近完全图时,则称为稠密图。
相反,当一个图含有较少的边数(即当e<<n(n-1))时,则称为稀疏图。

5、子图:
设有两个图G=(V,E)和G'=(V',E'),若V'是V的子集,即V'V,且E'是E的子集,即E'E,则称G'是G的子图。

  1. 路径和路径长度:
    从顶点i到顶点j的一条路径是一个顶点序列(i,i1,i2,…,im,j),序列中边(i,i1),(i1,i2),…,(im-1,im),(im,j)属于E(G);
    路径长度是指一条路径上经过的边的数目。
    简单路径:一条路径上除开始点和结束点可以相同外,其余顶点均不相同,

  2. 回路或环:
    回路或环:一条路径上的开始点与结束点为同一个顶点。
    简单回路或简单环:开始点与结束点相同的简单路径。

8、连通、连通图和连通分量:
无向图:若从顶点i到顶点j有路径,则称顶点i和j是连通的。
连通图:若图中任意两个顶点都连通,否则称为非连通图。
连通分量:无向图G中的极大连通子图。
任何连通图的连通分量只有一个,即本身而非连通图有多个连通分量。
有向图:若任意两个顶点之间都存在一条有向路径,则称此有向图为强连通图。
否则,其各个强连通子图称作它的强连通分量。
*在一个非强连通中找强连通分量的方法:
1.在图中找有向环。
2.扩展该有向环:如果某个顶点到该环中任一顶点有路径,并且该环中任一顶点到这个顶点也有路径,则加入这个顶点。

  1. 权和网
    权:图中每一条边都可以附有一个对应的数值,这种与边相关的数值称为权。
    网:边上带有权的图称为带权图,也称作网。

1.1.3 图存储结构

  • 图的两种主要存储结构:

1.邻接矩阵(二维数组);

  • 邻接矩阵的主要特点: 
    一个图的邻接矩阵表示是唯一的。
    特别适合于稠密图的存储。(一个图的邻接矩阵表示是唯一的。特别适合于稠密图的存储。)

  • 邻接矩阵存储表示:
    1.顶点信息:记录各个顶点信息的顶点表。
    2.边或弧信息:各个顶点之间关系的邻接矩阵。
    设图 A = (V, E)是一个有 n 个顶点的图
    A.Vex[n]:表示顶点信息集
    二维数组 A.edge[n][n]表示边的关系
    如果<i,j>属于E,或者(i,j)属于E: A.edge[i][j]=1。否则A.edge[i][j]=0.

  • 邻接矩阵结构体定义:

#define  MAXV  <最大顶点个数>	
typedef struct 
{    int no;			//顶点编号
     InfoType info;		//顶点其他信息
} VertexType;//声明顶点的类型

typedef struct  			//图的定义
{    int edges[MAXV][MAXV]; 	//邻接矩阵
     int n,e;  			//顶点数,边数
     VertexType vexs[MAXV];	//存放顶点信息
}  MatGraph;//声明的邻接矩阵类型

 MatGraph g;//声明邻接矩阵存储的图

  • 两种图的邻接矩阵:
    无向图对称:

    有向图非对称:

  • 网络的邻接矩阵:
    A[i][j]=W(i,j),如果i!=j且<i,j>∈E或(i,j)∈E
    A[i][j]=∞,否则,但是i!=j
    A[i][j]=0,对角线i==j
    eg:

  • 借助于邻接矩阵容易求得顶点的度:
    在无向图中,统计第i行(列)1的个数可得顶点i的度。
    即:顶点i的度=第i行或列总度数
    在有向图中:
    统计第i行1的个数可得顶点i的出度OD;
    统计第j列1的个数可得顶点j的入度ID。

2.邻接表

邻接表是对图中每个顶点i建立一个单链表,将顶点i的所有邻接点链起来。
每个单链表上添加一个表头结点(表示顶点信息)。并将所有表头结点构成一个数组,下标为i的元素表示顶点i的表头结点。
图的邻接表存储方法是一种顺序分配与链式分配相结合的存储方法。 

  • 图的邻接表存储类型定义:
typedef struct Vnode
{    Vertex data;			//顶点信息
     ArcNode *firstarc;		//指向第一条边
}  VNode;

typedef struct ANode
{     int adjvex;			//该边的终点编号
      struct ANode *nextarc;	//指向下一条边的指针
      InfoType info;		//该边的权值等信息
}  ArcNode;

typedef struct 
{     VNode adjlist[MAXV] ;	//邻接表
       int n,e;			//图中顶点数n和边数e
} AdjGraph;
AdjGraph *G;//声明一个邻接表存储的图G

  • 两种图的邻接表:
    无向图:

有向图:

  • 带权值网络的邻接表:

1.1.4 图的基本运算

1.创建图并初始化:

void CreateAdj(AdjGraph*& G, int n, int e) //创建图邻接表
{
	int i, j, a, b;
	ArcNode* p;
	G = new AdjGraph;
	for (i = 0; i < n; i++)   G->adjlist[i].firstarc = NULL;
	//给邻接表中所有头结点的指针域置初值
	for (i = 1; i <= e; i++)		 //根据输入边建图      
	{
		cin >> a >> b;
		p = new ArcNode;	//创建一个结点p
		p->adjvex = b;		 //存放邻接点
		p->nextarc = G->adjlist[a].firstarc;  //采用头插法插入结点p
		G->adjlist[a].firstarc = p;
	}
	G->n = n; G->e = n;
}

2.输出图:

void DispAdj(AdjGraph *G)	//输出邻接表G
{      int i;
       ArcNode *p;
       for (i=0;i<G->n;i++)
       {	p=G->adjlist[i].firstarc;//访问第一个节点
	printf("%3d: ",i);
	while (p!=NULL)
	{       printf("%3d[%d]→",p->adjvex,p->weight);
	         p=p->nextarc;
	}
	printf("∧\n");
       }
}

3.销毁图:

void DestroyAdj(AdjGraph *&G)   //销毁邻接表
{      int i; ArcNode *pre,*p;
       for (i=0;i<G->n;i++)		//扫描所有的单链表
       {	pre=G->adjlist[i].firstarc;//p指向第i个单链表的首结点
	if (pre!=NULL)
	{      p=pre->nextarc;
	        while (p!=NULL)	//释放第i个单链表的所有边结点
	        {	free(pre);
		pre=p; p=p->nextarc;
	         }
	         delete pre;
	}
        }
        delete G;			//释放头结点数组
}

  • 邻接表的特点:
    1.邻接表表示不唯一。
    2.特别适合于稀疏图存储。 (邻接表的存储空间为O(n+e))

  • 邻接表和邻接矩阵的相互转换:
    1.邻接矩阵转邻接表:

void MatToList(MGraph g,ALGraph *&G)
//将邻接矩阵g转换成邻接表G
{  int i,j,n=g.n; ArcNode *p; 	//n为顶点数
   G=new AdjGraph;
   for (i=0;i<n;i++)     //给所有头节点的指针域置初值
      G->adjlist[i].firstarc=NULL;
   for (i=0;i<n;i++)	 //检查邻接矩阵中每个元素
      for (j=n-1;j>=0;j--)
         if (g.edges[i][j]!=0) 			
         {  p=new ArcNode; //创建节点*p
	      p->adjvex=j;
	      p->nextarc=G->adjlist[i].firstarc;
                       //将*p链到链表头
	      G->adjlist[i].firstarc=p;
	   }
   G->n=n;G->e=g.e;
}

2.邻接表转邻接矩阵:

void ListToMat(ALGraph *G,MGraph &g)
  {  int i,j,n=G->n;ArcNode *p;
     for (i=0;i<n;i++) 
     {  p=G->adjlist[i].firstarc;
        while (p!=NULL) 
	  {   g.edges[i][p->adjvex]=1;
	      p=p->nextarc;
	  }
     }
     g.n=n;g.e=G->e;
  } 

1.1.6 图遍历及应用:

  • 图的遍历的概念:
    从给定图中任意指定的顶点(称为初始点)出发,按照某种搜索方法沿着图的边访问图中的所有顶点,使每个顶点仅被访问一次。

  • 图的两种遍历方式:
    1.深度优先遍历(DFS):
    *深度优先搜索遍历的过程是:
      (1)从图中某个初始顶点v出发,首先访问初始顶点v。
      (2)选择一个与顶点v相邻且没被访问过的顶点w为初始顶点,再从w出发进行深度优先搜索,直到图中与当前顶点v邻接的所有顶点都被访问过为止。 
     
    *深度优先遍历图的实质:对每个顶点查找其邻接点的过程。

*DFS遍历(访问v节点,遍历v的邻接点w,若w未被访问,递归访问w节点):

void DFS(ALGraph *G,int v)  
{    ArcNode *p;
      visited[v]=1;                   //置已访问标记
      printf("%d  ",v); 		
      p=G->adjlist[v].firstarc;      	
     while (p!=NULL) 
	{
                  if (visited[p->adjvex]==0)  DFS(G,p->adjvex);    
	     p=p->nextarc;              	
	}
}

邻接表:时间复杂度为O(n+e)
邻接矩阵:O(n2)

  • 非连通图的深度优先搜索遍历:
    首先将图中每个顶点的访问标志设为 FALSE, 之后搜索图中每个顶点,如果未被访问,则以该顶点为起始点,进行深度优先搜索遍历,否则继续检查下一顶点。
void DFSTraverse(Graph G) {
   // 对非连通图 G 作深度优先遍历。
  for (v=0; v<G.vexnum; ++v) 
     visited[v] = FALSE; // 访问标志数组初始化
  for (v=0; v<G.vexnum; ++v) 
     if (!visited[v])  DFS(G,v);  // 对尚未访问的顶点调用DFS
}

2.广度优先遍历(BFS):
*广度优先搜索遍历的过程是:
(1)访问初始点v,接着访问v的所有未被访问过的邻接点。
(2)按照次序访问每一个顶点的所有未被访问过的邻接点。  
(3)依次类推,直到图中所有顶点都被访问过为止。

*eg:

BFS序列:2 1 3 4 0
类似树的层次遍历:队列

*BFS搜索思路:

建一个访问队列q
访问v节点,加入队列q
while(队列不空)
    取队头元素w
    遍历w的邻接表    
         取邻接点j
         若j未被访问,则加入队列q,并访问j。
end while

邻接表:时间复杂度为O(n+e)
邻接矩阵:O(n2)

*用BFS遍历非连通图:

void  BFS1(AdjGraph *G)
{      int i;
        for (i=0;i<G->n;i++)     //遍历所有未访问过的顶点
             if (visited[i]==0) 
                  BFS(G,i);
}

*用DFS判断图是否连通:
采用某种遍历方式来判断无向图G是否连通。这里用深度优先遍历方法,先给visited[]数组(为全局变量)置初值0,然后从0顶点开始遍历该图。
在一次遍历之后,若所有顶点i的visited[i]均为1,则该图是连通的;否则不连通。

int  visited[MAXV];
bool Connect(AdjGraph *G) 	//判断无向图G的连通性
{     int i;
      bool flag=true;
      for (i=0;i<G->n;i++)		 //visited数组置初值
	visited[i]=0;
      DFS(G,0); 	//调用前面的中DSF算法,从顶点0开始深度优先遍历
      for (i=0;i<G->n;i++)
            if (visited[i]==0)
           {     flag=false;
	   break;
           }
      return flag;
}

*如何查找图路径:
利用回溯的深度优先遍历方法。
从顶点u开始进行深度优先遍历。增加path和d记录存走过的路径。
若当前扫描的顶点u = v时,表示找到了一条路径,则输出路径path。
当从顶点u出发的路径找完后,置visited[u]=0,即回溯。

void FindAllPath(AGraph *G,int u,int v,int path[],int d)
{ //d表示path中的路径长度,初始为-1
  int w,i;  ArcNode *p;
  d++; path[d]=u;		//路径长度d增1,顶点u加入到路径中
  visited[u]=1;		//置已访问标记
  if (u==v && d>=1)		//找到一条路径则输出
        {	for (i=0;i<=d;i++)
	    printf("%2d",path[i]);
	printf("\n");
        }
        p=G->adjlist[u].firstarc;	//p指向顶点u的第一个相邻点
        while (p!=NULL)
        {	 w=p->adjvex;		//w为顶点u的相邻顶点
	 if (visited[w]==0)	//若w顶点未访问,递归访问它
	     FindAllPath(G,w,v,path,d);
	 p=p->nextarc;		//p指向顶点u的下一个相邻点
        }
       visited[u]=0;//恢复环境,使该顶点可重新使用
}

*如何找最短路径:
定义非循环队列结构体类型:

typedef struct
{      int data;	//顶点编号
       int parent;	//前一个顶点的位置
} QUERE;		

以路径上经过的边数来衡量路径长度:

void ShortPath(AdjGraph *G,int u,int v)
{   //输出从顶点u到顶点v的最短逆路径
       qu[rear].data=u;//第一个顶点u进队
        while (front!=rear)//队不空循环
        {      front++;		//出队顶点w
               w=qu[front].data;
              if (w==v)   根据parent关系输出路径break; 
              while(遍历邻接表)   
                {         rear++;//将w的未访问过的邻接点进队
		 qu[rear].data=p->adjvex;
		 qu[rear].parent=front;
	  }
         }	      
}

广度优先遍历找到的路径一定是最短路径,而深度优先遍历则不一定。
深度优先遍历能找所有路径,而广度优先遍历难以实现。

1.1.6 最小生成树相关算法及应用

  • 生成树和最小生成树:
    一个连通图的生成树是一个极小连通子图,它含有图中全部n个顶点和构成一棵树的(n-1)条边。不能回路。  
    如果在一棵生成树上添加一条边,必定构成一个环。

    一个连通图的生成树不一定是唯一的。

*最小生成树的概念:
对于带权连通图G ,n个顶点,n-1条边
根据深度遍历或广度遍历生成生成树,树不唯一
其中权值之和最小的生成树称为图的最小生成树。

*非连通图和生成树:
非连通图:需多次调用遍历过程。
每个连通分量中的顶点集和遍历时走过的边一起构成一棵生成树。所有连通分量的生成树组成非连通图的生成森林。

  • 寻找最小生成树:
    一.普里姆算法:

普里姆(Prim)算法是一种构造性算法,用于构造最小生成树。过程如下:
(1)初始化U={v}。v到其他顶点的所有边为候选边;
(2)重复以下步骤n-1次,使得其他n-1个顶点被加入到U中:
从候选边中挑选权值最小的边输出,设该边在V-U中的顶点是k,将k加入U中;
考察当前V-U中的所有顶点j,修改候选边:若(j,k)的权值小于原来和顶点k关联的候选边,则用(k,j)取代后者作为候选边

*构造的最小生成树不一定唯一,但最小生成树的权值之和一定是相同的。

*如何选j到U顶点集的最小边?
设置2个辅助数组。
1.closest[i]:最小生成树的边依附在U中顶点编号。
2.lowcost[i]表示顶点i(i ∈ V-U)到U中顶点的边权重,取最小权重的顶点k加入U。并规定lowcost[k]=0表示这个顶点在U中
3.(closest[k],k)构造最小生成树一条边。
最小生成树求解:
1.找最小lowcost
2.加入新顶点,修正lowcost, closest

  • prim算法代码的实现:
#define INF 32767//INF表示∞
void  Prim(Graph G,int v)
{     
      int lowcost[MAXV];
      int min;
      int closest[MAXV], i, j, k;
      for (i=0;i<G.n;i++)//给lowcost[]和closest[]置初值
      {	
        lowcost[i]=G.edges[v][i];
	closest[i]=v;
      }

       for (i=1;i<G.n;i++)//输出(n-1)条边
       {	
          min=INF;
	  for (j=0;j<G.n;j++) //在(V-U)中找出离U最近的顶点k
	     if (lowcost[j]!=0 && lowcost[j]<min)
	     {	
                min=lowcost[j];
		k=j;//k记录最近顶点编号
	     }
	     lowcost[k]=0;//标记k已经加入U
             for (j=0;j<G.n;j++)//修改数组lowcost和closest
	     if (lowcost[j]!=0 && G.edges[k][j]<lowcost[j])
	     {	
                  lowcost[j]=G.edges[k][j];
		  closest[j]=k;
	     }
        }
}
  • 贪心算法:
    算法原理:以当前情况为基础作最优选择,而不考虑各种可能的整体情况,所以贪心法不要回溯。
    算法优点:因为省去了为寻找解而穷尽所有可能所必须耗费的大量时间,因此算法效率高。
    贪婪算法的精神就是“只顾如何获得眼前最大的利益”,有时不一定是最优解。

prim算法的应用:
公路村村通:
现有村落间道路的统计数据表中,列出了有可能建设成标准公路的若干条道路的成本,求使每个村落都有公路连通所需要的最低成本。
解决代码:

#include<stdio.h>
#include <iostream>
using namespace std;
#define INF 32767
#define MAX 1001
typedef struct
{
	int n, e;//定义图的顶点数和边数
	int edge[MAX][MAX];//邻接矩阵表示顶点之间的关系
}MINGraph;
//创建村庄图
void CreateEdge(MINGraph*& G, int n, int e)
{
	int v1, v2;//同一条边的两端顶点
	int weight;//边的权值
	G = new MINGraph;
	G->n = n;
	G->e = e;
	int i, j;
	//邻接矩阵的初始化
	for (i = 1; i <= n; i++)
	{
		for (j = 1; j <= n; j++)
		{
			if (i == j)
			{
				G->edge[i][j] = 0;
			}
			else
			{
				G->edge[i][j] = INF;
			}
		}
	}
	//为邻接矩阵赋值
	for (i = 1; i <= e; i++)
	{
		cin >> v1 >> v2 >> weight;
		G->edge[v1][v2] = G->edge[v2][v1] = weight;
	}
}
int Prim(MINGraph *G, int v)
{
	//边权重    最小边权重  边依附在U中顶点
	int lowcost[MAX], min, closest[MAX], i, j, k=0;
	int fee=0;//最小费用
	for (i = 1; i <=G->n; i++)	//给lowcost[]和closest[]置初值
	{
		lowcost[i] = G->edge[v][i];
		closest[i] = v;
	}
	for (i = 1; i < G->n; i++)	  //找出(n-1)个顶点
	{
		min = INF;
		for (j = 1; j <= G->n; j++) //     在(V-U)中找出离U最近的顶点k
			if (lowcost[j] != 0 && lowcost[j] < min)
			{
				min = lowcost[j];
				k = j; // k记录最近顶点的编号
			}
		lowcost[k] = 0;		//标记k已经加入U
		fee = fee + min;//最小权值的和
		for (j = 1; j <= G->n; j++)	//修改数组lowcost和closest
			if (lowcost[j] != 0 && G->edge[k][j] < lowcost[j])
			{
				lowcost[j] = G->edge[k][j];
				closest[j] = k;
			}
	}
	if (G->e < n - 1)
		return -1;
	for (i = 1; i <= G->n; i++)
	{
		if (lowcost[i] != 0)
			return -1;
	}
	return fee;
}
int main()
{
	MINGraph* G;
	int n, e;
	cin >> n >> e;
	CreateEdge(G, n, e);
	int fee;
	fee = Prim(G,1);
	cout << fee;
	return 0;

}

二.克鲁斯卡尔算法(Kruskal)
克鲁斯卡尔(Kruskal)算法也是一种求带权无向图的最小生成树的构造性算法。是一种按权值的递增次序选择合适的边来构造最小生成树的方法。

  • 克鲁斯卡尔(Kruskal)算法过程(时间复杂度为O(n2)):
    (1)置U的初值等于V(即包含有G中的全部顶点),TE(表示最小生成树的边集)的初值为空集(即图T中每一个顶点都构成一个连通分量)。
    (2)将图G中的边按权值从小到大的顺序依次选取:
    若选取的边未使生成树T形成回路,则加入TE;
    否则舍弃,直到TE中包含(n-1)条边为止。

  • Kruskal算法实现的具体代码:

typedef struct 
{    int u;     //边的起始顶点
     int v;      //边的终止顶点
     int w;     //边的权值
} Edge; 
Edge E[MAXV];
void Kruskal(Graph G)
{    
      int i,j,u1,v1,sn1,sn2,k;
      int vset[MAXV];
      Edge E[MaxSize];	//存放所有边
      k=0;			//E数组的下标从0开始计
      for (i=0;i<G.n;i++)	//由G产生的边集E
	for (j=0;j<G.n;j++)
	      if (G.edges[i][j]!=0 && G.edges[i][j]!=INF)
	      {     
                     E[k].u=i;  E[k].v=j;  E[k].w=G.edges[i][j];
	             k++;
	      }
        InsertSort(E,G.e);	//用直接插入排序对E数组按权值递增排序
        for (i=0;i<g.n;i++) 	//初始化辅助数组
	vset[i]=i;
	k=1;		//k表示当前构造生成树的第几条边,初值为1
	j=0;		//E中边的下标,初值为0
	while (k<G.n)	//生成的边数小于n时循环
	{ 
              u1=E[j].u;v1=E[j].v;	//取一条边的头尾顶点
	      sn1=vset[u1];
	      sn2=vset[v1];		//分别得到两个顶点所属的集合编号
 	      if (sn1!=sn2)  	//两顶点属于不同的集合
	      {
                k++;		   	//生成边数增1
		for (i=0;i<g.n;i++)  	//两个集合统一编号
		    if (vset[i]==sn2) 	//集合编号为sn2的改为sn1
		       vset[i]=sn1;
	      }
	     j++;			   //扫描下一条边
        }
}

并查集改良Kruskal算法(时间复杂度是O(elog2e)):

typedef struct 
{    int u;     //边的起始顶点
     int v;      //边的终止顶点
     int w;     //边的权值
} Edge; 
Edge E[MAXV];
void Kruskal(AdjGraph *g)
{     int i,j,u1,v1,sn1,sn2,k;
      int vset[MAXV]; //集合辅助数组
      Edge E[MaxSize];	//存放所有边
      k=0;			//E数组的下标从0开始计
    for (i=0;i<g.n;i++)	//由g产生的边集E,邻接表
	{   p=g->adjlist[i].firstarc;
        while(p!=NULL)    
	  {	E[k].u=i;E[k].v=p->adjvex;
            E[k].w=p->weight;
		k++; p=p->nextarc;
	  }
     }
     Sort(E,g.e);	//用快排对E数组按权值递增排序
      for (i=0;i<g.n;i++) 	//初始化集合
	vset[i]=i;
        k=1;		//k表示当前构造生成树的第几条边,初值为1
	j=0;		//E中边的下标,初值为0
	while (k<g.n)	//生成的顶点数小于n时循环
	{ 
                    u1=E[j].u;v1=E[j].v;	//取一条边的头尾顶点
	      sn1=vset[u1];
	      sn2=vset[v1];	//分别得到两个顶点所属的集合编号
 	      if (sn1!=sn2)  	//两顶点属于不同的集合
	      {	printf("  (%d,%d):%d\n",u1,v1,E[j].w);
		k++;		   	//生成边数增1
		for (i=0;i<g.n;i++)  	//两个集合统一编号
		      if (vset[i]==sn2) 	//集合编号为sn2的改为sn1
			vset[i]=sn1;
	     }
	     j++;			   //扫描下一条边
            }
}
  • 2种算法比较;
    普里姆算法:O(n2)、适用于稠密图
    克鲁斯卡尔算法:O(eloge)、适用于稀疏图
    实现Prim算法,选择图的邻接矩阵存储结构
    实现克鲁斯卡尔算法,选择图的邻接表存储结构

1.1.7 最短路径

  • 最短路径问题是指:
    如果从图中某一顶点(源点)到达另一顶点(终点)的路径可能不止一条,如何找到一条路径使得沿此路径上各边的权值总和(称为路径长度)达到最小。

  • 两种常见的最短路径问题:
    1、 单源最短路径—用Dijkstra(迪杰斯特拉)算法(一顶点到其余各顶点)
    2、所有顶点间的最短路径—用Floyd(弗洛伊德)算法(任意两顶点之间)

Dijkstra算法:

  • 过程:
初始化顶点合集S和未选顶点合集T;
1.S={入选顶点集合,初值V0},T={未选顶点集合}。
  若存在<V0,Vi>,距离值为<V0,Vi>弧上的权值
  若不存在<V0,Vi>,距离值为∞
  从T中选取一个其距离值为最小的顶点W, 加入S
2.从T中选取一个其距离值为最小的顶点W, 加入S
3.S中加入顶点w后,对T中顶点的距离值进行修改:
  若加进W作中间顶点,从V0到Vj的距离值比不加W的路径要短,则修改此距离值;
4.重复上述步骤1,直到S中包含所有顶点,即S=V为止。

  • 两个常见问题的解决:
1.如何存放最短路径?(最短路径中所有顶点都是最短路径)
采用一维数组path来保存:
![](https://img2020.cnblogs.com/blog/1776748/202005/1776748-20200516231441070-1145664932.png)

2.如何存放最短路径长度?
用一维数组dist[j]存储:dist[j]表示起点v到顶点j的最短路径长度并从源点开始根据最短顶点来调整。
  • eg:


  • Dijkstra算法的具体代码实现:(算法的时间复杂度为O(n2))

void Dijkstra(Graph G,int v)
{    
      int dist[MAXV],path[MAXV];
      int s[MAXV];
      int mindis,i,j,u;
      for (i=0;i<G.n;i++)
      {       
        dist[i]=G.edges[v][i];	//距离初始化
	s[i]=0;			//s[]置空
	if (G.edges[v][i]<INF)	//路径初始化
	       path[i]=v;  //顶点v到i有边时
	else
	      path[i]=-1;		
      }
      s[v]=1;	 		//源点v放入S中
      for (i=0;i<G.n;i++)	 	//循环n-1次
      {     
         mindis=INF;
	 for (j=0;j<G.n;j++)
	 {
             if (s[j]==0 && dist[j]<mindis) 
	     {        
                u=j;
		mindis=dist[j];
	     }
         }
	 s[u]=1;			//顶点u加入S中
	 for (j=0;j<G.n;j++)	//修改不在s中的顶点的距离
	 {
            if (s[j]==0)
	    {
                if (G.edges[u][j]<INF && dist[u]+G.edges[u][j]<dist[j])
	        {      
                  dist[j]=dist[u]+G.edges[u][j];
	   	   path[j]=u;
	        }
            }
         }
      }
    //输出最短路径
}
  • Dijkstra算法特点:
    1.不适用带负权值的带权图求单源最短路径。
    2.不适用求最长路径长度。
    3.最短路径长度是递增
    4.顶点u加入S后,不会再修改源点v到u的最短路径长度
    5.按Dijkstra算法,找第一个距离源点S最远的点A,这个距离在以后就不会改变。但A与S的最远距离一般不是直连。

Floyd算法(时间复杂性是O(n3)):

  • 算法思路:
有向图G=(V,E)采用邻接矩阵存储
二维数组A用于存放当前顶点之间的最短路径长度,分量A[i][j]表示当前顶点i到顶点j的最短路径长度。
递推产生一个矩阵序列A0,A1,…,Ak,…,An-1
     Ak+1[i][j]表示从顶点i到顶点j的路径上所经过的顶点编号k+1的最短路径长度。

  • Floyd算法的具体代码:
void Floyd(Graph G)		
{     
    int A[MAXVEX][MAXVEX];	//建立A数组
    int path[MAXVEX][MAXVEX];	//建立path数组
    int i, j, k;
   for (i=0;i<G.n;i++)   		
       for (j=0;j<G.n;j++) 
       {       
          A[i][j]=G.edges[i][j];
	 if (i!=j && G.edges[i][j]<INF)
	      path[i][j]=i; 	//i和j顶点之间有一条边时
       else			 
	      path[i][j]=-1;
       }

  for (k=0;k<G.n;k++)		//求Ak[i][j]
 {     
     for (i=0;i<G.n;i++)
       for (j=0;j<G.n;j++)
	    if (A[i][j]>A[i][k]+A[k][j])	//找到更短路径
	    {   
                A[i][j]=A[i][k]+A[k][j];	//修改路径长度
	       path[i][j]=path[k][j]; 	//修改最短路径为经过顶点k
          }
  }
}	

1.1.8 拓扑排序、关键路径

一.拓扑排序:

在一个有向图中找一个拓扑序列的过程称为拓扑排序。序列必须满足条件:
每个顶点出现且只出现一次。
若存在一条从顶点 A 到顶点 B 的路径,那么在序列中顶点 A 出现在顶点 B 的前面。
图中有回路,无法拓扑排序.拓扑排序可以用来检测图中是否有回路.

拓扑排序的过程:
1.从有向图中选取一个没有前驱的顶点,并输出之;
2.从有向图中删去此顶点以及所有以它为尾的弧;
3.重复上述两步,直至图空,或者图不空但找不到无前驱的顶点为止。


*拓扑排序:1->2->3->4->7->5->6

  • 具体代码:
typedef struct 	       	//表头结点类型
{     
      Vertex data;         	//顶点信息
      int count;           	//存放顶点入度
      ArcNode *firstarc;   	//指向第一条边
}VNode;

 void TopSort(AdjGraph *G)	//拓扑排序算法
{      
        int i,j;
        int St[MAXV],top=-1;	//栈St的指针为top
        ArcNode *p;
        for (i=0;i<G->n;i++)		//入度置初值0
	G->adjlist[i].count=0;
        for (i=0;i<G->n;i++)		//求所有顶点的入度
        {	
            p=G->adjlist[i].firstarc;
	    while (p!=NULL)
	   {        
                  G->adjlist[p->adjvex].count++;
	          p=p->nextarc;
	   }
        }
         
         for (i=0;i<G->n;i++)		//将入度为0的顶点进栈
	 if (G->adjlist[i].count==0)
	 {	
            top++;
	    St[top]=i;
	 }
         while (top>-1)			//栈不空循环
         {	  
            i=St[top];top--;			//出栈一个顶点i
	    printf("%d ",i);		//输出该顶点
	    p=G->adjlist[i].firstarc;		//找第一个邻接点
	    while (p!=NULL)		//将顶点i的出边邻接点的入度减1
	    {      
                 j=p->adjvex;
	         G->adjlist[j].count--;
	         if (G->adjlist[j].count==0)	//将入度为0的邻接点进栈
	         {      
                   top++;
		   St[top]=j;
	         }
	         p=p->nextarc;		//找下一个邻接点
	    }
        }
}

二.关键路径:

AOE-网:
用顶点表示事件,用有向边e表示活动,边的权c(e)表示活动持续时间。是带权的有向无环图
整个工程完成的时间为:从有向图的源点到汇点的最长路径。又叫关键路径(critical path).
“关键活动(key activity)”指的是:关键路径中的边.

AOE 网图关键词:

AOE网——带权的有向无环图
顶点--事件或状态
弧(有向边)--活动及发生的先后关系
权--活动持续的时间
起点--入度为0的顶点(只有一个)
终点--出度为0的顶点(只有一个)

eg:

拓扑序列为:
C4,C0,C3,C2,C1,C5

  • 求关键路径的过程:
    1.事件的最早开始和最迟开始时间
    事件v最早开始时间ve(v):v作为源点事件最早开始时间为0。
    2.事件v的最迟开始时间vl(v):定义在不影响整个工程进度的前提下,事件v必须发生的时间称为v的最迟开始时间

  • 活动:边 的最早开始时间和最迟开始时间:
    1.活动a(边)的最早开始时间e(a)指该活动起点x事件的最早开始时间,即:e(a)=ve(x)

2.活动a的最迟开始时间l(a)指该活动终点y事件的最迟开始时间与该活动所需时间之差,即:l(a)=vl(y)-c

  • 关键活动:d(a)=l(a)-e(a),若d(a)为0,则称活动a为关键活动。
    关键路径上的活动都是关键活动
    关键活动不存在富余时间,适当增加对关键活动的投资,减少关键活动的持续时间,缩短工程工期。

  • 求关键路径步骤:

1.对有向图拓扑排序
2.根据拓扑序列计算事件(顶点)的ve,vl数组
          ve(j) = Max{ve(i) + dut(<i,j>)}
           vl(i) = Min{vl(j) - dut(<i,j>)}
3.计算关键活动的e[],l[]。即边的最早、最迟时间
  e(i) = ve(j)
	l(i) = vl(k) - dut(<j, k>
4.找e=l边即为关键活动
5.关键活动连接起来就是关键路径

  • eg:

1.2.谈谈你对图的认识及学习体会

对图的认识:
图结构多对多结构其实是线性结构和非线性结构的大本营,线性结构一对一,树形结构一对多,图形结构多对多,图常见类型是有向图 无向图 稀疏图 稠密图 连通图 非连通图 完全图 非完全图等等。用邻接矩阵和邻接表进行存储。
学习体会:
感觉拓扑序列有点不太懂,有的结构是在下一个结构开讲之后才搞懂的,图的建立感觉最难,结构体定义代码量少了也不太容易记住。

2.阅读代码(来源:力扣(LeetCode))

2.1 题目及解题代码

  • 题目:

  • 解题代码:

/*
trustSize:表示对数
*trustColSize:表示每队是2个
*/

int findJudge(int N, int** trust, int trustSize, int* trustColSize) {

	int cnt[N + 2];
	memset(cnt, 0, sizeof(int) * (N + 2));

	// 如果被信任,则cnt++,如果信任了别人,则cnt置为-1,肯定不是法官,以后cnt不再增加
	for (int i = 0; i < trustSize; i++) {
		if (cnt[trust[i][1]] >= 0) {
			cnt[trust[i][1]]++;
		}
		cnt[trust[i][0]] = -1;
	}

	// 被信任数是N - 1的则是法官
	for (int i = 1; i <= N; i++) {
		if (cnt[i] == N - 1) {
			return i;
		}
	}
	return -1;
}




2.1.1 该题的设计思路

法官不相信任何人,说明法官不存在出度
所有人都信任法官,说明法官的入度为N-1
那么法官的出度加入度为N-1



时间复杂度:O(n)
空间复杂度:O(n)

2.1.2 该题的伪代码

int findJudge(int N, int** trust, int trustSize, int* trustColSize) {

	//定义出入度数组cnt

	// 如果被信任,则cnt++,如果信任了别人,则cnt置为-1,肯定不是法官,以后cnt不再增加
	for i = 0 to i = trustSize-1
		if 出度不为0即信任后继节点人
			后继者可能是法官,入度+1
		end if
		cnt[trust[i][0]] = -1;//前一个人相信别人,肯定不是法官
	end for

	// 被信任数是N - 1的则是法官
	for i = 1 to i = N
		if (入度为N - 1)
			return i;
		end if 
	}
	return -1;
}

2.1.3 运行结果

2.1.4分析该题目解题优势及难点。

优势:
根据图中节点的出度入度来判断节点人的身份,当所有人除了法官自己都相信法官时法官就可以判断出来了。
难点:
用图的方式来做这道题,其实我感觉方法比较多变叭,好像线性结构也可以做?

2.2 题目及解题代码

  • 题目:

  • 解题代码:

class Solution {
	public boolean isBipartite(int[][] graph) {
		int n = graph.length;
		int[] color = new int[n];
		Arrays.fill(color, -1);

		for (int start = 0; start < n; ++start) {
			if (color[start] == -1) {
				Stack<Integer> stack = new Stack();
				stack.push(start);
				color[start] = 0;

				while (!stack.empty()) {
					Integer node = stack.pop();
					for (int nei : graph[node]) {
						if (color[nei] == -1) {
							stack.push(nei);
							color[nei] = color[node] ^ 1;
						}
						else if (color[nei] == color[node]) {
							return false;
						}
					}
				}
			}
		}

		return true;
	}
}


2.2.1 该题的设计思路

使用数组(或者哈希表)记录每个节点的颜色: color[node]。颜色可以是 0, 1,或者未着色(-1 或者 null)。
搜索节点时,需要考虑图是非连通的情况。对每个未着色节点,从该节点开始深度优先搜索着色。每个邻接点都可以通过当前节点着相反的颜色。如果存在当前点和邻接点颜色相同,则着色失败。
使用栈完成深度优先搜索,栈类似于节点的 “todo list”,存储着下一个要访问节点的顺序。在 graph[node] 中,对每个未着色邻接点,着色该节点并将其放入到栈中。


时间复杂度:O(N+E)
空间复杂度:O(N)

2.2.2 该题的伪代码

class Solution {
	public boolean isBipartite(int[][] graph) {
		定义图节点数n
		颜色数组color[]
		Arrays.fill(color, -1);

		for (int start = 0; start < n; ++start) {
			if 邻接点未着色
				初始化栈stack
				未着色邻接点设为已着色并入栈

				while 栈不为空 
					取栈顶node
					for (int nei : graph[node]) //检查邻边
						if 当前节点未着色 
							nei节点入栈
							给color[nei]着相反色
						end if
						else if color[nei]和color[node]颜色相同;即弧两端节点颜色相同
							着色失败
						end if
					end for
				end while
			end if
		end for

		return true;
	}
}


2.2.3 运行结果

2.2.4分析该题目解题优势及难点

优势:
用红色蓝色标记模块,然后一条弧两边的节点不能是同一个颜色。
用颜色标记各节点之间的关系,到了统计合集的时候用颜色代表就可以了
难点:
需要注意图可能是不连通的:
1、有单独一个点的,没有任何连接,单独的点不算一个独立子集
2、有两个子图的,没有连接

2.3 题目及解题代码

  • 题目:

  • 解题代码:

int trap(vector<int>& height)
{
	if (height == null)
		return 0;
	int ans = 0;
	int size = height.size();
	vector<int> left_max(size), right_max(size);
	left_max[0] = height[0];
	for (int i = 1; i < size; i++) {
		left_max[i] = max(height[i], left_max[i - 1]);
	}
	right_max[size - 1] = height[size - 1];
	for (int i = size - 2; i >= 0; i--) {
		right_max[i] = max(height[i], right_max[i + 1]);
	}
	for (int i = 1; i < size - 1; i++) {
		ans += min(left_max[i], right_max[i]) - height[i];
	}
	return ans;
}

2.3.1 该题的设计思路

找到数组中从下标 i 到最左端最高的条形块高度 left_max.
找到数组中从下标 i 到最右端最高的条形块高度 right_max.
扫描数组 height 并更新答案:
累加 min(max_left[i], max_right[i])−height[i]到 ans 上

时间复杂度:O(n)。
空间复杂度:O(n)。

2.3.2 该题的伪代码

int trap(vector<int>& height)
{
	if height为 0(没有原始图)
		return 0;
	初始化可以接到的雨水量ans
	int size = height.size();
	定义左右两边最大高度数组left_max[]和right_max[]
	left_max[0] = height[0];//左边第一最大高度=当前最大高度
	for i = 1 to i = size-1 
		从当前元素向左扫描并更新:找到数组中从下标 i 到最左端最高的条形块高度left_max。
	end for
	right_max[size - 1] = height[size - 1];//右边第一最大高度=当前最大高度
	for i = size - 2 to i >= 0
		从当前元素向右扫描并更新:找到数组中从下标 i 到最左端最高的条形块高度right_max。
	end for
	for i = 1 to i = size - 2
		累加两边最大高度的较小值减去当前高度的值到ans上;
	end for
	return ans;
}

2.3.3 运行结果

2.3.4分析该题目解题优势及难点

优势:
提前存储可以接到雨水的最大值,不用每次都扫描获得最大值。

难点:
空间复杂度稍微大一点;
然后要注意雨水量是左边右边模块的更小值去减雨水的当前高度。