1. 参考论文
Reweighted Random Walks for Graph Matching【1】
The pagerank citation ranking: Bringing order to the web【2】
A spectral technique for correspondence problems using pairwise constraints【3】
2. 论文概述
论文提出了全新的思路——马尔科夫链去解决图匹配问题，提出了“基本版” Random Walk 并从随机游走角度恰好解释了【3】的谱方法SM的理论。随后根据【2】的思路改进成 RRM，通过实验证明其效果还是蛮不错的。
3. 论文介绍
3.1 传统问题描述
传统的图匹配问题，优化的目标函数为：
通过限制条件可以看出这是 IQP (integer quadratic program) 问题，其中 w 是两个图的相似度矩阵，x 是排列矩阵的列向量化（vec）
3.2 Markov角度问题描述
这是论文最核心的地方，即如何将传统的图匹配问转化为随机过程Markov的问题。我们先看平时场见的相似度矩阵 w 的形式：
观察上图(g)，可以发现不管行、列，其坐标每一项代表含义都是两图中两个点的假设对应关系，如 1a 。我们最终的目的也是找到类似如此的对应关系，只不过这里借助的是相似度矩阵。
那么，作者就是从相似度矩阵出发，将其改造成随机游走模型，如图：
两个图P、Q的点的对应关系就当作Grw的结点，而相似度矩阵的对角线上的值（即一阶点的对应关系）当作Grw的 attribute，相似度矩阵非对角线上的 attribute 当作Grw的两个顶点（如13--22）构成的边的值。这样一来，我们将相似度矩阵就变换成了Grw的形式，同时，整个图匹配问题，变成了在改图上的随机游走过程，最终的目标成了在Grw上选点的过程—— To select the nodes in Grw（ selecting reliable nodes in the Grw）
而如何选点，论文借鉴了【2】.
3.3 转移矩阵
Markov 过程需要一个转移矩阵（概率矩阵），那么对应到这里，怎么构建呢？首先可以明确的是，肯定是从 w 入手，将其改造成转移矩阵——每行和为1，每一项 ≥ 0（满足）。
对于前者条件，一种常见的思路是：
即求 w 每一行的和，然后改行分别处以该数——但是其弊处明显，文中提到这种做法会放大不匹配点的概率，即相似度矩阵本来是相似度值的大小，通过除法，会削弱相似度大的值，反之放大原本小的值的匹配概率（如 1 0 2；3 0 6；---->均为 1/3 0 2/3）
作者采用将 w 矩阵每一项均除以 max(sum(和))，即以此做法保护原来相似度值较大的点。思路是这样，然后为了满足Markov的sum（行）=1的特点，需要改造转移矩阵：
其中左式 P 就是转移矩阵，其中dmax = max(sum(行和))，d(向量)是每行和，式子第一项就是为了满足转移矩阵行和为1，同时底部第二行是为了保证P是个方阵，原本是 xP，既然 P 是多了一项，所以改造x，引入（x xabs）P，即如上图右边式子，上标n代表转移次数
而引入 xabs，其实就是对Grw再此进行改造，增加一个虚拟结点（absorb 吸收点），因为 P 的转移概率是（0000....000 1）, 即恰好代表 nm*nm+1 这个结点不会转移到其他点，只能自转移，概率为1。
如下图，就是引入该 abs 的图
3.4 为什么求解Markov链的稳态就是结果
论文最终就是通过 P 转移矩阵，求得 x 的稳态（稳态Markov链），我想因为 x 的初始分布（应该）可以平均分配（即随机游走开始（1/mn*mn 1/mn*mn.....）），而Grw的edge代表了相似度，进一步该相似度被转化成了P（即概率表示），那么最终得到的 x 的值（再还原回矩阵）代表了两两结点如何匹配的score，然后使用匈牙利算法即可选出该x还原回来的矩阵的最佳匹配score
至于为什么稳态（转移矩阵n次幂）一定存在：
