用邻接矩阵或关联矩阵表示图,称为图的代数表示。用矩阵表示图,主要有两个优点:

  • (1) 能够把图输入到计算机中
  • (2) 可以用代数方法研究图论

Spectral Graph Theory 是研究图和代数的一个领域

# 邻接矩阵

Adjaceny matrix is n×nn\times n symmetric matrix

# 定义

GGnn 阶图,V={v1,z2,,vn}V=\{v_1,z_2,\cdots,v_n\} ,邻接矩阵 A(G)=(aij)A(G)= (a_{ij}) ,其中:

aij={l,vivj间边数0,vivj不邻接ai上的环数,i=ja_{ij}=\begin{cases} l,\ \ \ v_i\text{与}v_j\text{间边数}\\ 0,\ \ \ v_i\text{与}v_j\text{不邻接}\\ a_i\text{上的环数},\ \ \ i=j \end{cases}

# 代数性质

非负性 —— 由邻接矩阵定义知 aija_{ij} 是非负整数,即邻接矩阵是非负整数矩阵

(无向图) 对称性 —— 显然 aij=ajia_{ij}=a_{ji} ,所以邻接矩阵是对称矩阵。

同一图的不同形式的邻接矩阵是相似矩阵。这是因为,同图的两种不同形式矩阵可以通过交换行和相应的列变成致。

如果 G 为简单图,则 A(G)A(G)布尔矩阵;行和 (列和) 等于对应顶点的度数;矩阵元素总和为图的总度数,也就是 GG 的边数的 2 倍

# 连通的充要条件

GG 连通的充要条件是:A(G)A(G) 不能与如下矩阵相似:

(A1100A22)\begin{pmatrix} A_{11} & 0\\ 0 & A_{22} \end{pmatrix} \quad

# 性质 (定理)

Ak(G)=(aij(k))A^{k}(G)=(a_{ij}^{(k)}) ,则 aij(k)a_{ij}^{(k)} 表示顶点 viv_ivjv_j 的途径长度为 kk 的途径数目;


证明:对 kk 作数学归纳法证明

k=1k=1 时,由邻接矩阵定义;

设结论对 k1k-1 时成立,考虑 kk 的情况:

Ak=Ak1A=(ai1(k1)aj1+ai2(k1)aj2++ain(k1)ajn)n×n=(aij(k))A^k=A^{k-1}A=(a_{i1}^{(k-1)}a_{j1}+a_{i2}^{(k-1)}a_{j2}+\cdots+a_{in}^{(k-1)}a_{jn})_{n\times n}=(a_{ij}^{(k)})

同时,考虑顶点 viv_ivjv_j 的途径长度为 kk 的途径,设 vmv_mviv_ivjv_j 的途径中点,且该点与 vjv_j 邻接。则 viv_ivjv_j 的经过 vmv_m 且途径长度为 kk 的途径数目为:

aim(k1)amja_{im}^{(k-1)}a_{mj}

所以 viv_ivjv_j 的途径长度为 kk 的途径数目为:

ai1(k1)aj1+ai2(k1)aj2++ain(k1)ajn=aij(k)a_{i1}^{(k-1)}a_{j1}+a_{i2}^{(k-1)}a_{j2}+\cdots+a_{in}^{(k-1)}a_{jn}=a_{ij}^{(k)}


推论:

GG 是连通的,则对于 iji\neq jviv_ivjv_j 的距离是使 AnA^{n}aij(k)0a_{ij}^{(k)}\neq0 的最小整数

这是显然的

# 关联矩阵

Incidence matrix is n×mn\times m matrix of a (n,m)(n,m) graph

# 定义

GG(n,m)(n,m) 图。定义 GG 的关联矩阵: M(G)=(aij)n×mM(G)=(a_{ij})_{n\times m} ,其中:aija_{ij}viv_ieje_j 关联的次数 (取值为 0,1, 或 2 (环)).

  • 关联矩阵的元素为 0,1 或 2 ;
  • 关联矩阵的每列和为 2 ; 每行的和为对应顶点度数;

# 连通的充要条件

无环图 GG 连通的充分必要条件𝑅(𝑀)=n1𝑅(𝑀)=n-1


证明充分性:

GG 不连通,假设 GG 有两个连通分支 G1G_1G2G_2 , 又设 G1G_1G2G_2 的关联矩阵分别为 𝑀1𝑀_1𝑀2𝑀_2 , 则 GG 的关联矩阵可以写为:

M(G)=(M1M2)M(G)= \begin{pmatrix} M_1 & \\ & M_2 \end{pmatrix}

于是 𝑅(𝑀)=V(G1)1+V(G2)1=n2𝑅(𝑀)=V(G_1)-1+V(G_2)-1=n-2 ,矛盾!所以 GG 一定连通。


证明必要性:

令:

M=(m1m2mn)M= \begin{pmatrix} m_1\\ m_2\\ \vdots\\ m_n \end{pmatrix}

由于 𝑀𝑀 中每列恰有两个 “1” 元素,所有行向量的和为 0 (模 2 加法运算,即为偶数),所以有 R(M)n1R(M)\leq n-1

另一方面,在 MM 中任意去掉一行 mkm_k ,由于 GG连通的,因此,mkm_k 中存在元素 “1” ,这样,MM 中去掉行 mkm_k ,后的行按模 2 相加不等于零,即它们是线性无关的。所以有 R(M)n1R(M)\geq n-1

因此,R(M)=n1R(M)=n-1

# 基本关联矩阵

GG 的关联矩阵中删掉任意一行后得到的矩阵可以完全决定 GG ,该矩阵称为 GG基本关联矩阵。删掉的行对应的顶点称为该基本关联矩阵的参考点

图的关联矩阵及其性质是网络图论的基础,在电路分析中有重要应用

图的关联矩阵比邻接矩阵大得多,不便于计算机存储。但二者都有各自的应用特点

# 定义 2

Incidence matrix CC is m×nm\times n matrix of a (n,m)(n,m) graph

用来定义拉普拉斯矩阵 L=CTCL=C^{T}C

  • 每行代表一条边,每列代表一个顶点
  • 每行和为 0 ,有一个 1 和 -1 ,分别表示边的起点和终点
  • 其余元素全部为 0

C=(e1Te2TemT)C= \begin{pmatrix} e_1^{T}\\ e_2^{T}\\ \vdots\\ e_m^{T} \end{pmatrix}

# 度矩阵

Degree matrix is n×nn\times n diagonal matrix

# 定义

度矩阵是对角阵,对角上的元素为各个顶点的度。

⚠️有向图中,可能会有不同的定义。一般度数为各个顶点入度与出度之和,例如拉普拉斯矩阵引入的度矩阵。

# 拉普拉斯矩阵

Laplacian matrix (LL) is n×nn\times n symmetric matrix

拉普拉斯矩阵是对称矩阵,保留了图的 “无向” 的属性,丢失了 “有向” 的属性

GG 是顶点集合为 V(G)={v1,v2,,vn}V(G)=\{v_1,v_2,\cdots,v_n\} 的图, DD 为图的度矩阵A=aijA=a_{ij}GG邻接矩阵CC 为图的关联矩阵

定义:

L=DA=CTCL=D-A=C^{T}C

拉普拉斯矩阵 L=(lij)L=(l_{ij})nn 阶方阵,其中:

lij={d(vi),i=jaij,ijl_{ij}=\begin{cases} d(v_i), & i=j\\ -a_{ij}, & i\neq j \end{cases}

  • 所有特征值为非负实数: 0λ0λ1λn10\leq\lambda_0\leq\lambda_1\leq\cdots\leq\lambda_{n-1}
  • 特征向量为实 (正交) 向量
  • 每一行、列之和都等于 0
  • xTLx0x^{T}Lx\geq0 是半正定二次型,或者说拉普拉斯矩阵 LL 是半正定矩阵

[Fiedler, 1973] GGkk 个连通分支当且仅当:

λ0=λ1==λk1=0\lambda_0=\lambda_1=\cdots=\lambda_{k-1}=0

GG 显然至少有一个连通分支,所以 λ0=0\lambda_0=0

GG 的割边 (cut edges) 数量等于:

14xTLx=14(i,j)E(xixj)2\frac{1}{4}x^{T}Lx=\frac{1}{4}\sum_{(i,j)\in E}(x_i-x_j)^2

(等式右侧推导在下面)

对任意对称矩阵 MM 定义 λ2\lambda_2

λ2=minxxTMxxTx\lambda_2=\min_{x}\frac{x^{T}Mx}{x^{T}x}

考虑拉普拉斯矩阵 LL 对应的二次型 xTLxx^{T}Lx

xTLx=i,j=1nLijxixj=i,j=1n(DijAij)xixj=i=1nDiixi2(i,j)E2xixj=(i,j)E(xi2+xj22xixj)=(i,j)E(xixj)2\begin{aligned} x^{T}Lx &=\sum_{i,j=1}^{n}L_{ij}x_ix_j\\ &=\sum_{i,j=1}^{n}(D_{ij}-A_{ij})x_ix_j\\ &=\sum_{i=1}^{n}D_{ii}x_i^2-\sum_{(i,j)\in E}2x_ix_j\\ &=\sum_{(i,j)\in E}(x_i^2+x_j^2-2x_ix_j)\\ &=\sum_{(i,j)\in E}(x_i-x_j)^2 \end{aligned}

此时:

λ2=minAll labelings of so thatxi=0(i,j)E(xixj)2ixi2\lambda_2=\min_{\text{All labelings of so that }\sum x_i=0}\frac{\sum_{(i,j)\in E}(x_i-x_j)^2}{\sum_ix_i^2}

# 子图空间

# 回路空间

https://lfool.github.io/LFool-Notes/other/ 图论.html# 回路系统简介

# 断集空间