数据流分析共性问题
之前分别介绍了可用表达式分析 和活跃变量分析 两类经典的数据流分析实例,不难发现,它们的分析过程间存在相似性,都是基于图的表示,并通过建立数据流方程组求解的分析框架,这样统一的分析框架,为求解数据流方程组提供高效的通用算法提供了可能。
不同分析实例的相似性
我们可以比较可用表达式分析与活跃变量分析的数据流方程组形式化定义:
AE l = { ∅ if l = init ( c ) ⋂ { φ l ′ ( AE l ′ ) ∣ ( l ′ , l ) ∈ flow ( c ) } otherwise \text{AE}_l =
\begin{cases}
\emptyset & \text{if } \textcolor{red}{l = \text{init}(c)} \\
\textcolor{red}{\bigcap} \{ \varphi_{l'}(\text{AE}_{l'}) \mid (l', l) \in \text{flow}(c) \} & \text{otherwise}
\end{cases} AE l = { ∅ ⋂ { φ l ′ ( AE l ′ ) ∣ ( l ′ , l ) ∈ flow ( c )} if l = init ( c ) otherwise
LV l = { Var c if l ∈ final ( c ) ⋃ { φ l ′ ( LV l ′ ) ∣ ( l , l ′ ) ∈ flow ( c ) } otherwise \text{LV}_l =
\begin{cases}
\text{Var}_c & \text{if } \textcolor{red}{l \in \text{final}(c)} \\
\textcolor{red}{\bigcup} \{ \varphi_{l'}(\text{LV}_{l'}) \mid (l, l') \in \text{flow}(c) \} & \text{otherwise}
\end{cases} LV l = { Var c ⋃ { φ l ′ ( LV l ′ ) ∣ ( l , l ′ ) ∈ flow ( c )} if l ∈ final ( c ) otherwise
它们在表示上存在相似性,实际上我们可以定义一个统一的形式表达来定义这样的分析过程:
对于 c ∈ Cmd c \in \text{Cmd} c ∈ Cmd 且 l ∈ Lab c l \in \text{Lab}_c l ∈ Lab c ,分析信息(A I \mathrm{AI} AI )由下面形式的方程描述:
AI l = { ι if l ∈ E ⨆ { φ l ′ ( AI l ′ ) ∣ ( l ′ , l ) ∈ F } otherwise \text{AI}_l = \begin{cases}
\iota & \text{if } l \in E \\
\bigsqcup \{ \varphi_{l'}(\text{AI}_{l'}) \mid (l', l) \in F \} & \text{otherwise}
\end{cases} AI l = { ι ⨆ { φ l ′ ( AI l ′ ) ∣ ( l ′ , l ) ∈ F } if l ∈ E otherwise
其中:
组合操作子 ⨆ \bigsqcup ⨆ :对应公式中的组合操作,有路径量化方式确定;
块B ′ B' B ′ 的迁移函数φ l ′ \varphi_{l'} φ l ′ :用于处理分析信息的迁移;
流关系:flow ( c ) \text{flow}(c) flow ( c ) 或 flow R ( c ) \text{flow}^\text{R}(c) flow R ( c ) (即 f l o w \mathrm{flow} flow 关系的逆),对应公式中的 ( l ′ , l ) ∈ F (l', l) \in F ( l ′ , l ) ∈ F ;
端结点的标签集合E E E :{ init ( c ) } \{ \text{init}(c) \} { init ( c )} 或 final ( c ) \text{final}(c) final ( c ) ,用于界定公式中 l ∈ E l \in E l ∈ E 的情况。
当 F = f l o w ( c ) F = \mathrm{flow(c)} F = flow ( c ) 时属于前向分析的信息流;当 F = f l o w R ( c ) F = \mathrm{flow^R(c)} F = flo w R ( c ) 时属于后向分析的信息流。而组合操作子 ⊔ \sqcup ⊔ 则根据路径上的量化方式,如可用表达式分析的 M u s t \mathrm{Must} Must 类型使用 ⨆ = ⋂ \bigsqcup = \bigcap ⨆ = ⋂ ,而如活跃变量分析的 M a y \mathrm{May} May 类型则使用 ⨆ = ⋃ \bigsqcup = \bigcup ⨆ = ⋃ 。
分析框架
不同的数据流分析实例在相同的形式化表达下构建数据流方程组,在相同的分析框架下,通过不动点迭代来求解数据流方程组。
把方程组的解刻画为一个转换子的不动点;
引入偏序关系来比较分析结果;
建立最小上界或最大下界作为组合操作子;
保证迁移函数的单调性;
确保不动点迭代的终止性(满足递增链条件);
通过 w o r k l i s t \mathrm{worklist} worklist 算法来对不动点迭代进行优化。
理论基础
序理论(O r d e r T h e o r y \mathrm{Order\ Theory} Order Theory )
偏序(P a r t i a l O r d e r \mathrm{Partial\ Order} Partial Order )
IMPORTANT 偏序关系是定义在集合 D D D 上的一种二元关系 ⊑ \sqsubseteq ⊑ ,需要满足下面三个条件:
自反性:∀ a ∈ D \forall a \in D ∀ a ∈ D 有 a ⊑ a a \sqsubseteq a a ⊑ a ;
传递性:若 a ⊑ b a \sqsubseteq b a ⊑ b 且 b ⊑ c b \sqsubseteq c b ⊑ c ,则 a ⊑ c a \sqsubseteq c a ⊑ c ;
反对称性:若 a ⊑ b a \sqsubseteq b a ⊑ b 且 b ⊑ a b \sqsubseteq a b ⊑ a ,则 a = b a = b a = b .
若 ⊑ \sqsubseteq ⊑ 构成集合 D D D 上的偏序关系,则称 ( D , ⊑ ) (D, \sqsubseteq) ( D , ⊑ ) 为偏序集(P o s e t \mathrm{Poset} Poset )。
TIP 设 S S S 是一个非空集合,二元关系 R R R 是 A A A 中元素的有序对(二元序偶)集合,即 R ⊆ S × S R \subseteq S\times S R ⊆ S × S (S × S S \times S S × S 表示 S S S 中所有的有序对 < a , b > <a, b> < a , b > 的集合);若 < a , b > ∈ R <a, b> \in R < a , b >∈ R ,则称”a a a 与 b b b 具有关系 R R R “,记作 a R b aRb a R b 。
不要混淆偏序关系与偏序集。两者都是集合,但偏序关系是”元素间的关系集合”,而偏序集则是”带关系的元素集合载体”,即前者描述一种关系,后者描述一个带有某种关系的集合。
WARNING 需要注意,偏序集中并非所有元素都能够比较关系,但是只要能够比较,就满足自反、传递以及反对称的规则。
示例
( Z , ≤ ) (\mathbb{Z}, \le) ( Z , ≤ )
( P ( x ) , ⊆ ) (\mathscr{P}(x), \subseteq) ( P ( x ) , ⊆ ) ,其中 P \mathscr{P} P 表示幂集
( Z 2 , ⊑ ) (\mathbb{Z}^2, \sqsubseteq) ( Z 2 , ⊑ ) ,其中,( a , b ) ⊑ ( a ′ , b ′ ) ⇔ a ≥ a ′ ∧ b ≥ b ′ (a, b) \sqsubseteq (a', b') \Leftrightarrow a \ge a' \land b \ge b' ( a , b ) ⊑ ( a ′ , b ′ ) ⇔ a ≥ a ′ ∧ b ≥ b ′
前面提到,偏序集中并非所有元素都能够比较关系,如 ( P ( { 1 , 2 } ) , ⊆ ) (\mathscr{P}(\{1, 2\}), \subseteq) ( P ({ 1 , 2 }) , ⊆ ) 中的元素 { 1 } \{1\} { 1 } 和 { 2 } \{2\} { 2 } 之间就不能进行比较。示例3表示了区间包含关系。
格(l a t t i c e \mathrm{lattice} lattice )
IMPORTANT 设偏序集 ( D , ⊑ ) (D, \sqsubseteq) ( D , ⊑ ) ,若集合 D D D 中任意两个元素 a , b a, b a , b 都有最小上界(即上确界,L U B \mathrm{LUB} LUB )和最大下界(即下确界,G L B \mathrm{GLB} GLB ),则称之为格,该格记作 ( D , ⊑ , ⊔ , ⊓ ) (D, \sqsubseteq, \sqcup, \sqcap) ( D , ⊑ , ⊔ , ⊓ ) 。其中:
最小上界:记为 a ⊔ b a\sqcup b a ⊔ b ∀ a , b ∈ D \forall a, b \in D ∀ a , b ∈ D ,存在唯一元素 c ∈ D c \in D c ∈ D ,满足:
a ⊑ c a\sqsubseteq c a ⊑ c 且 b ⊑ c b \sqsubseteq c b ⊑ c (即 c c c 是 a a a 和 b b b 的上界);
若存在 d ∈ D d \in D d ∈ D 使得 a ⊑ d a\sqsubseteq d a ⊑ d 且 b ⊑ d b \sqsubseteq d b ⊑ d ,则 c ⊑ d c \sqsubseteq d c ⊑ d (即 c c c 是所有上界中最”小”的那个).
最大下界:记为 a ⊓ b a\sqcap b a ⊓ b ∀ a , b ∈ D \forall a, b \in D ∀ a , b ∈ D ,存在唯一元素 c ∈ D c \in D c ∈ D ,满足:
c ⊑ a c\sqsubseteq a c ⊑ a 且 c ⊑ b c \sqsubseteq b c ⊑ b (即 c c c 是 a a a 和 b b b 的下界);
若存在 d ∈ D d \in D d ∈ D 使得 d ⊑ a d\sqsubseteq a d ⊑ a 且 d ⊑ b d \sqsubseteq b d ⊑ b ,则 d ⊑ c d \sqsubseteq c d ⊑ c (即 c c c 是所有上界中最”大”的那个).
示例
整数 ( Z , ≤ , m a x , m i n ) (\mathbb{Z}, \le, max, min) ( Z , ≤ , ma x , min )
整数区间 ( { [ a , b ] ∣ a , b ∈ Z , a ≤ b } ∪ 0 , ⊆ , ⊔ , ⊓ ) (\{[a, b] | a, b \in \mathbb{Z}, a \le b\} \cup {0}, \subseteq, \sqcup, \sqcap) ({[ a , b ] ∣ a , b ∈ Z , a ≤ b } ∪ 0 , ⊆ , ⊔ , ⊓ ) ,其中:
[ a , b ] ⊔ [ a ′ , b ′ ] = def [ min ( a , a ′ ) , max ( b , b ′ ) ] [a, b] \sqcup [a', b'] \stackrel{\text{def}}{=} [\min(a, a'), \max(b, b')] [ a , b ] ⊔ [ a ′ , b ′ ] = def [ min ( a , a ′ ) , max ( b , b ′ )]
[ a , b ] ⊓ [ a ′ , b ′ ] = def [ max ( a , a ′ ) , min ( b , b ′ ) ] , or ∅ if max ( a , a ′ ) > min ( b , b ′ ) [a, b] \sqcap [a', b'] \stackrel{\text{def}}{=} [\max(a, a'), \min(b, b')], \text{ or } \emptyset \text{ if } \max(a, a') > \min(b, b') [ a , b ] ⊓ [ a ′ , b ′ ] = def [ max ( a , a ′ ) , min ( b , b ′ )] , or ∅ if max ( a , a ′ ) > min ( b , b ′ )
如,整数区间 ( { [ a , b ] ∣ a , b ∈ { 0 , 1 , 2 } , a ≤ b } ∪ 0 , ⊆ , ⊔ , ⊓ ) (\{[a, b] | a, b \in \{0, 1, 2\}, a \le b\} \cup {0}, \subseteq, \sqcup, \sqcap) ({[ a , b ] ∣ a , b ∈ { 0 , 1 , 2 } , a ≤ b } ∪ 0 , ⊆ , ⊔ , ⊓ ) ,其中每条边表示偏序关系。
完全格(C o m p l e t e L a t t i c e \mathrm{Complete\ Lattice} Complete Lattice )
IMPORTANT 设偏序集 ( D , ⊑ ) (D, \sqsubseteq) ( D , ⊑ ) ,若集合 D D D 的任意子集 X X X (无论有限还是无限)都存在最小上界 ⊔ X \sqcup X ⊔ X 和最大下界 ⊓ X \sqcap X ⊓ X ,特别地,对集合 D D D 总存在最小元 ⊥ = ⊓ ∅ \bot = \sqcap \empty ⊥ = ⊓ ∅ 和最大元 ⊤ = ⊔ D \top = \sqcup D ⊤ = ⊔ D ,此时称之为完全格,记作 ( D , ⊑ , ⊔ , ⊓ , ⊥ , ⊤ ) (D, \sqsubseteq, \sqcup, \sqcap, \bot, \top) ( D , ⊑ , ⊔ , ⊓ , ⊥ , ⊤ ) 。
实际上,定义中任意子集 X X X 都存在最小上界和最大上界是比格的定义中任意两个元素都存在最小上界和最大下界更严格的条件,因此完全格首先是格。此外对于集合 D D D 而言,其最小元 ⊥ = ⊓ D = ⊓ ∅ \bot = \sqcap D = \sqcap \empty ⊥ = ⊓ D = ⊓ ∅ ,其最大元 ⊤ = ⊔ D \top = \sqcup D ⊤ = ⊔ D 。
TIP 对于 ∅ \empty ∅ 而言,⊥ = ⊤ = ⊔ ∅ = ⊓ ∅ \bot = \top = \sqcup \empty = \sqcap \empty ⊥ = ⊤ = ⊔ ∅ = ⊓ ∅ .
示例
幂集 ( P ( S ) , ⊆ , ∪ , ∩ , ∅ , S ) (\mathscr{P}(S), \subseteq, \cup, \cap, \empty, S) ( P ( S ) , ⊆ , ∪ , ∩ , ∅ , S )
如 ( P { 0 , 1 , 2 } , ⊆ , ∪ , ∩ , ∅ , { 0 , 1 , 2 } ) (\mathscr{P}\{0, 1, 2\}, \subseteq, \cup, \cap, \empty, \{0, 1, 2\}) ( P { 0 , 1 , 2 } , ⊆ , ∪ , ∩ , ∅ , { 0 , 1 , 2 }) :
示例
其中,对于子集 { { 1 } , { 2 } , { 0 , 1 } } \{\{1\}, \{2\}, \{0, 1\}\} {{ 1 } , { 2 } , { 0 , 1 }} ,其最小上界为 { 0 , 1 , 2 } \{0, 1, 2\} { 0 , 1 , 2 } ,最大下界为 ∅ \empty ∅ 。
带无穷界的整数区间
( I , ⊆ , ⊔ , ⊓ , ∅ , [ − ∞ , + ∞ ] ) (\mathbb{I}, \subseteq, \sqcup, \sqcap, \emptyset, [-\infty, +\infty]) ( I , ⊆ , ⊔ , ⊓ , ∅ , [ − ∞ , + ∞ ])
I = def { [ a , b ] ∣ a ∈ Z ∪ { − ∞ } , b ∈ Z ∪ { + ∞ } , a ≤ b } ⊔ i ∈ I [ a i , b i ] = def [ min i ∈ I a i , max i ∈ I b i ] ⊓ i ∈ I [ a i , b i ] = def [ max i ∈ I a i , min i ∈ I b i ] , or ∅ if max > min \begin{aligned}
&\mathbb{I} \stackrel{\text{def}}{=} \{ [a, b] \mid a \in \mathbb{Z} \cup \{-\infty\}, b \in \mathbb{Z} \cup \{+\infty\}, a \leq b \} \\
&\sqcup_{i \in I} [a_i, b_i] \stackrel{\text{def}}{=} [\min_{i \in I} a_i, \max_{i \in I} b_i] \\
&\sqcap_{i \in I} [a_i, b_i] \stackrel{\text{def}}{=} [\max_{i \in I} a_i, \min_{i \in I} b_i] \text{, or } \emptyset \text{ if } \max > \min
\end{aligned} I = def {[ a , b ] ∣ a ∈ Z ∪ { − ∞ } , b ∈ Z ∪ { + ∞ } , a ≤ b } ⊔ i ∈ I [ a i , b i ] = def [ i ∈ I min a i , i ∈ I max b i ] ⊓ i ∈ I [ a i , b i ] = def [ i ∈ I max a i , i ∈ I min b i ] , or ∅ if max > min
链(C h a i n \mathrm{Chain} Chain )
IMPORTANT 设偏序集 ( D , ⊑ ) (D, \sqsubseteq) ( D , ⊑ ) ,取集合 D D D 的非空子集 S S S ,若 S S S 中的任意两个元素关于 ⊑ \sqsubseteq ⊑ 皆可比,则称该子集 S S S 为 ( D , ⊑ ) (D, \sqsubseteq) ( D , ⊑ ) 的一条链。
例如,∀ S ⊆ N \forall S \subseteq \mathbb{N} ∀ S ⊆ N 是 ( N , ⊆ ) (\mathbb{N}, \subseteq) ( N , ⊆ ) 上的一条链(高度可能无穷)。
递增链
IMPORTANT 设序列 ( d i ) i ∈ N (d_i)_{i \in \mathbb{N}} ( d i ) i ∈ N 是链,如果对于 ∀ i ∈ N \forall i \in \mathbb{N} ∀ i ∈ N 都有:d i ⊑ d i + 1 d_i \sqsubseteq d_{i+1} d i ⊑ d i + 1 ,则称该链为递增链。
递增链条件(A s c e n d i n g C h a i n C o n d i t i o n , A C C \mathrm{Ascending\ Chain\ Condition, ACC} Ascending Chain Condition , ACC )
IMPORTANT 设偏序集 ( D , ⊑ ) (D, \sqsubseteq) ( D , ⊑ ) ,如果对其中的任意递增链 d 0 ⊑ d 1 ⊑ . . . d_0 \sqsubseteq d_1 \sqsubseteq ... d 0 ⊑ d 1 ⊑ ... 最终都会稳定(即 ∃ n ∈ N \exist n \in \mathbb{N} ∃ n ∈ N 使得对所有 m > n m > n m > n ,有 a m = a n a_m = a_n a m = a n ),则称 D D D 满足递增链条件。
TIP 这里只要求 ( D , ⊑ ) (D, \sqsubseteq) ( D , ⊑ ) 满足偏序集,因此该偏序集中可以没有任何递增链,甚至可以不存在任何链。
有穷高度的链肯定满足递增链条件A C C ACC A C C ,但是满足递增链条件的链不一定是有穷高度(比如有可能存在不稳定的递减链)。
TIP 链 C C C 的高度定义为:链 C C C 中存在的最长严格递增链,则该最长严格递增链的长度称为链 C C C 的高度,记作 h t ( C ) \mathrm{ht}(C) ht ( C ) 。
其中,严格递增链是在递增链定义的基础上要求 d i ⊏ d i + 1 d_i \sqsubset d_{i+1} d i ⊏ d i + 1 ,即要求 d i ≠ d i + 1 d_i \not = d_{i+1} d i = d i + 1 。
因此,对于递增链在趋于稳定后,后半段是稳定值 a m a_m a m 的无限重复,并未构成严格递增链关系,此时该链的高度也并非无穷。
IMPORTANT 完全格且满足递增链条件是不动点迭代终止的充分条件!
完全偏序(C o m p l e t e P a r t i a l O r d e r , C P O \mathrm{Complete\ Partial\ Order, CPO} Complete Partial Order , CPO )
IMPORTANT 设偏序集 ( D , ⊑ ) (D, \sqsubseteq) ( D , ⊑ ) ,若:
D D D 存在最小元,记作 ⊥ D \bot_D ⊥ D ;
对 D D D 中每条链 S S S ,其最小上界 ⊔ S \sqcup S ⊔ S 存在。
则称 ( D , ⊒ ) (D, \sqsupseteq) ( D , ⊒ ) 为完全偏序(C P O \mathrm{CPO} CPO ),记作 ( D , ⊑ , ⊔ , ⊥ D ) (D, \sqsubseteq, \sqcup, \bot_D) ( D , ⊑ , ⊔ , ⊥ D )
示例
( N , ⊑ ) (\mathbb{N}, \sqsubseteq) ( N , ⊑ )
不是 C P O \mathrm{CPO} CPO ,N \mathrm{N} N 本身构成链,其最小上界为 ∞ \infin ∞ ,但是 ∞ ∉ N \infin \not \in \mathrm{N} ∞ ∈ N ;
( N ∪ { ∞ } , ⊑ ) (\mathbb{N} \cup \{\infin\}, \sqsubseteq) ( N ∪ { ∞ } , ⊑ )
是 C P O \mathrm{CPO} CPO 。
不动点理论
函数
单调性
IMPORTANT 设偏序集 ( D , ⊑ ) (D, \sqsubseteq) ( D , ⊑ ) 和 ( D ′ , ⊑ ′ ) (D', \sqsubseteq') ( D ′ , ⊑ ′ ) 上的全函数 f : D → D ′ f: D \rightarrow D' f : D → D ′ ,若 ∀ x , y ∈ D \forall x, y \in D ∀ x , y ∈ D ,且只要 x ⊑ y x \sqsubseteq y x ⊑ y 则有 f ( x ) ⊑ ′ f ( y ) f(x) \sqsubseteq' f(y) f ( x ) ⊑ ′ f ( y ) ,则该全函数 f f f 是单调的。
这里的偏序集 ( D , ⊑ ) (D, \sqsubseteq) ( D , ⊑ ) 和 ( D ′ , ⊑ ′ ) (D', \sqsubseteq') ( D ′ , ⊑ ′ ) 可以理解成定义域和值域。需要注意的是,f ( x ) ⊑ ′ f ( y ) f(x) \sqsubseteq' f(y) f ( x ) ⊑ ′ f ( y ) 中 ⊑ \sqsubseteq ⊑ 需要映射到 ⊑ ′ \sqsubseteq' ⊑ ′ 。若 D D D 中存在链 s s s ,单调性确保了该链经函数 f f f 映射到 D ′ D' D ′ 后构成新链 s ′ s' s ′ (根据单调性定义很容易证明)。
上确界连续(⊔ − \sqcup- ⊔ − 连续)
IMPORTANT 定义一 :
设 ( D , ⊑ ) (D, \sqsubseteq) ( D , ⊑ ) 和 ( D ′ , ⊑ ′ ) (D', \sqsubseteq') ( D ′ , ⊑ ′ ) 是完全偏序(C P O \mathrm{CPO} CPO ),存在全函数 f : D → D ′ f: D \rightarrow D' f : D → D ′ ,若对 D D D 中的每条链 s s s ,⊔ f ( s ) \sqcup f(s) ⊔ f ( s ) 存在且有 f ( ⊔ s ) = ⊔ f ( s ) f(\sqcup s) = \sqcup f(s) f ( ⊔ s ) = ⊔ f ( s ) ,则称 f f f 是 ⊔ − \sqcup- ⊔ − 连续(上确界连续)的
上面的定义中其实隐含了函数单调性:∀ x , y ∈ D \forall x, y \in D ∀ x , y ∈ D 且 x ⊑ y x \sqsubseteq y x ⊑ y ,则 x , y x,y x , y 在同一条链 { x , y } \{x, y\} { x , y } 上,而根据定义 f ( y ) = f ( ⊔ { x , y } ) = ⊔ f ( { x , y } ) = ⊔ { f ( x ) , f ( y ) } f(y) = f(\sqcup \{x, y\}) = \sqcup f(\{x, y\}) = \sqcup \{f(x), f(y)\} f ( y ) = f ( ⊔ { x , y }) = ⊔ f ({ x , y }) = ⊔ { f ( x ) , f ( y )} ,而 f ( x ) ⊑ ′ ⊔ { f ( x ) , f ( y ) } f(x) \sqsubseteq' \sqcup \{f(x), f(y)\} f ( x ) ⊑ ′ ⊔ { f ( x ) , f ( y )} ,从而 f ( x ) ⊑ ′ f ( y ) f(x) \sqsubseteq' f(y) f ( x ) ⊑ ′ f ( y ) ,故得证函数 f f f 是单调的。
函数 f f f 单调性保证了 D D D 中的链 s s s 经函数 f f f 映射到 ( D ′ , ⊑ ′ ) (D', \sqsubseteq') ( D ′ , ⊑ ′ ) 后仍然构成链 s ′ s' s ′ 。于是上面的定义可以理解为:对 D D D 中的任意一条链 s s s ,f f f 作用在链 s s s 的最小上界(上确界)上的结果,等于 f f f 作用在链 s s s 中的每个元素构成的新链 s ′ = f ( s ) s' = f(s) s ′ = f ( s ) 的最小上界(上确界),也即 f ( ⊔ s ) = ⊔ f ( s ) f(\sqcup s) = \sqcup f(s) f ( ⊔ s ) = ⊔ f ( s ) 。
IMPORTANT 定义二 :
函数 f f f 是 ⊔ − \sqcup- ⊔ − 连续的,当且仅当 f f f 是单调的且对 D D D 中的每条链 s s s ,有 f ( ⊔ s ) ⊑ ′ ⊔ f ( s ) f(\sqcup s) \sqsubseteq' \sqcup f(s) f ( ⊔ s ) ⊑ ′ ⊔ f ( s ) 。
这里和定义一相对比,条件给出函数 f f f 是单调的,并且 f ( ⊔ s ) ⊑ ′ ⊔ f ( s ) f(\sqcup s) \sqsubseteq' \sqcup f(s) f ( ⊔ s ) ⊑ ′ ⊔ f ( s ) ,因此实际上我们只需要基于函数 f f f 的单调性证明 ⊔ f ( s ) ⊑ ′ f ( ⊔ s ) \sqcup f(s) \sqsubseteq' f(\sqcup s) ⊔ f ( s ) ⊑ ′ f ( ⊔ s ) 即可得到定义一中的条件,也就能够证明函数 f f f 是 ⊔ − \sqcup- ⊔ − 连续。
证明如下:∀ x ∈ s \forall x \in s ∀ x ∈ s ,有 x ⊑ ⊔ s x \sqsubseteq \sqcup s x ⊑ ⊔ s ,于是 f ( x ) ⊑ ′ f ( ⊔ s ) f(x) \sqsubseteq' f(\sqcup s) f ( x ) ⊑ ′ f ( ⊔ s ) ,即对链 s s s 中任意元素 x x x 映射到 D ′ D' D ′ 中都有上界 f ( ⊔ s ) f(\sqcup s) f ( ⊔ s ) ,于是 ⊔ f ( s ) = ⊔ { f ( x 1 ) , . . . , f ( x k ) } ⊑ ′ f ( ⊔ s ) \sqcup f(s) = \sqcup \{f(x_1), ..., f(x_k)\} \sqsubseteq' f(\sqcup s) ⊔ f ( s ) = ⊔ { f ( x 1 ) , ... , f ( x k )} ⊑ ′ f ( ⊔ s ) 。根据条件 f ( ⊔ s ) ⊑ ′ ⊔ f ( s ) f(\sqcup s) \sqsubseteq' \sqcup f(s) f ( ⊔ s ) ⊑ ′ ⊔ f ( s ) ,于是有
f ( ⊔ s ) = ⊔ f ( s ) f(\sqcup s) = \sqcup f(s) f ( ⊔ s ) = ⊔ f ( s )
再根据定义一,可知函数 f f f 是 ⊔ − \sqcup- ⊔ − 连续的。
TIP 若 D D D 只包含有穷链,则 f f f 是 ⊔ − \sqcup- ⊔ − 连续的当且仅当 f f f 是单调的。
提示一下,既然是有穷链,则 ⊔ s ∈ s \sqcup s \in s ⊔ s ∈ s ,而根据单调性即可证明 ⊔ f ( s ) = f ( ⊔ s ) \sqcup f(s) = f(\sqcup s) ⊔ f ( s ) = f ( ⊔ s ) 。
示例
考虑从 ( N ∪ { ∞ } , ≤ ) (\mathbb{N}\cup \{\infin\}, \le) ( N ∪ { ∞ } , ≤ ) 到 ( { F , T } , ⊑ ) (\{F, T\}, \sqsubseteq) ({ F , T } , ⊑ ) 的函数 f ( x ) f(x) f ( x ) ,其中 F ⊑ T F \sqsubseteq T F ⊑ T ,判断函数 f ( x ) f(x) f ( x ) 单调性与 ⊔ − \sqcup- ⊔ − 连续性:
f ( x ) : i f x ≥ 2 t h e n T e l s e F f(x): \mathrm{if\ x\ge 2\ then\ T\ else\ F} f ( x ) : if x ≥ 2 then T else F
f ( x ) : i f x = = ∞ t h e n T e l s e F f(x): \mathrm{if\ x == \infin\ then\ T\ else\ F} f ( x ) : if x == ∞ then T else F
显然,两者中的 f ( x ) f(x) f ( x ) 均单调。但是 2.中 f ( x ) f(x) f ( x ) 不是 ⊔ − \sqcup- ⊔ − 连续的,理由如下:
根据定义,N ∪ { ∞ } \mathbb{N}\cup \{\infin\} N ∪ { ∞ } 中任意链 s s s 均有,⊔ f ( s ) \sqcup f(s) ⊔ f ( s ) 存在且有 f ( ⊔ s ) = ⊔ f ( s ) f(\sqcup s) = \sqcup f(s) f ( ⊔ s ) = ⊔ f ( s ) ,而对于2.中的 f ( x ) f(x) f ( x ) ,不妨取 s = N s = \mathbb{N} s = N ,由于 ∞ ∉ N \infin \not \in \mathbb{N} ∞ ∈ N ,于是 ⊔ f ( s ) = F ⊑ T = f ( ⊔ s ) \sqcup f(s) = F \sqsubseteq T = f(\sqcup s) ⊔ f ( s ) = F ⊑ T = f ( ⊔ s ) ,故得证。
不动点理论
基础概念
IMPORTANT 设 f f f 是偏序集( D , ⊑ ) (D, \sqsubseteq) ( D , ⊑ ) 上的函数,则 D D D 中元素 x x x 称为函数 f f f 的:
不动点(f i x p o i n t \mathrm{fixpoint} fixpoint ):若 f ( x ) = x \color{red}{\mathrm{f(x) = x}} f ( x ) = x ;
前不动点(p r e − f i x p o i n t \mathrm{pre-fixpoint} pre − fixpoint ):若 x ⊑ f ( x ) \color{red}{\mathrm{x \sqsubseteq f(x)}} x ⊑ f ( x ) ;
后不动点(p o s t − f i x p o i n t \mathrm{post-fixpoint} post − fixpoint ):若 f ( x ) ⊑ x \color{red}{\mathrm{f(x) \sqsubseteq x}} f ( x ) ⊑ x .
图示
图中的 p r e pre p r e 和 p o s t post p os t 都是相对于不动点 x 0 x_0 x 0 ,实际上 p o s t post p os t 所指部分相对于 x 1 x_1 x 1 也是前不动点。
IMPORTANT 不动点集合 f p f = { x ∈ D : f ( x ) = x } \mathrm{fp}\ f = \{x \in D: f(x) = x\} fp f = { x ∈ D : f ( x ) = x } ;
x ∈ f p f x \in \mathrm{fp}\ f x ∈ fp f ,若对任意 y ∈ f p f y \in \mathrm{fp}\ f y ∈ fp f 皆有 x ⊑ y x \sqsubseteq y x ⊑ y ,则称 x x x 为函数 f f f 的最小不动点,记作 l f p f \mathrm{lfp}\ f lfp f ;
x ∈ f p f x \in \mathrm{fp}\ f x ∈ fp f ,若对任意 y ∈ f p f y \in \mathrm{fp}\ f y ∈ fp f 皆有 y ⊑ x y \sqsubseteq x y ⊑ x ,则称 x x x 为函数 f f f 的最大不动点,记作 g f p f \mathrm{gfp}\ f gfp f ;
记 l f p d f \mathrm{lfp}_d\ f lfp d f 为关于偏序 ⊑ \sqsubseteq ⊑ 比 d d d 大的最小不动点;
记 g f p d f \mathrm{gfp}_d\ f gfp d f 为关于偏序 ⊑ \sqsubseteq ⊑ 比 d d d 小的最大不动点;
应用示例
以不动点理论在联立递归方程求解中的应用为例:
问题背景:联立递归方程
考虑由两个递归方程组成的方程组:
{ x 1 = f ( x 1 , x 2 ) x 2 = g ( x 1 , x 2 ) \begin{cases}
&x_1 = f(x_1, x_2) \\
&x_2 = g(x_1, x_2) \\
\end{cases} { x 1 = f ( x 1 , x 2 ) x 2 = g ( x 1 , x 2 )
其中,x 1 , x 2 x_1, x_2 x 1 , x 2 是格 x x x 中的元素;f , g f, g f , g 是定义在该格上的函数。
转化为不动点问题
为求解上述方程组,我们可以将其抽象为一个函数的不动点问题:
定义函数 F : X × X → X × X F: X \times X \rightarrow X \times X F : X × X → X × X ,其中 X × X X \times X X × X 表示原格的 “乘积格”,保持偏序结构,函数 F F F 形式为:
F ( [ x 1 x 2 ] ) = [ f ( x 1 , x 2 ) g ( x 1 , x 2 ) ] F(
\begin{bmatrix}
x_1 \\
x_2 \\
\end{bmatrix}
) =
\begin{bmatrix}
f(x_1, x_2) \\
g(x_1, x_2) \\
\end{bmatrix} F ( [ x 1 x 2 ] ) = [ f ( x 1 , x 2 ) g ( x 1 , x 2 ) ]
此时,方程组的解 ( x 1 , x 2 ) (x_1, x_2) ( x 1 , x 2 ) 满足:
F ( [ x 1 x 2 ] ) = [ x 1 x 2 ] F(
\begin{bmatrix}
x_1 \\
x_2 \\
\end{bmatrix}
) =
\begin{bmatrix}
x_1 \\
x_2 \\
\end{bmatrix} F ( [ x 1 x 2 ] ) = [ x 1 x 2 ]
即 ( x 1 , x 2 ) (x_1, x_2) ( x 1 , x 2 ) 是函数 F F F 的不动点,于是原方程组最小解为 l f p F \mathrm{lfp}\ F lfp F 。
T a r s k i \mathrm{Tarski} Tarski 不动点定理
IMPORTANT 完全格( D , ⊑ , ⊔ , ⊓ , ⊥ , ⊤ ) (D, \sqsubseteq, \sqcup, \sqcap, \bot, \top) ( D , ⊑ , ⊔ , ⊓ , ⊥ , ⊤ ) 上的单调函数 f : D → D f: D \rightarrow D f : D → D 的不动点集合 f p \mathrm{fp} fp 构成一个非空完全格,且:
l f p f = ⊓ p o s t f p f = ⊓ { x ∈ D : f ( x ) ⊑ x } g f p f = ⊔ p r e f p f = ⊔ { x ∈ D : f ( x ) ⊒ x } \begin{align*}
& \mathrm{lfp}\ f = \sqcap \mathrm{postfp}\ f = \sqcap \{x \in D: f(x) \sqsubseteq x\} \\
& \mathrm{gfp}\ f = \sqcup \mathrm{prefp}\ f = \sqcup \{x \in D: f(x) \sqsupseteq x\} \\
\end{align*} lfp f = ⊓ postfp f = ⊓ { x ∈ D : f ( x ) ⊑ x } gfp f = ⊔ prefp f = ⊔ { x ∈ D : f ( x ) ⊒ x } 其中,由不动点构成的非空完全格为 ( fp f , ⊑ , λ S . lfp ⊔ S , λ S . gfp ⊓ S , lfp f , gfp f ) (\text{fp } f, \sqsubseteq, \lambda S.\text{lfp}_{\sqcup S}, \lambda S.\text{gfp}_{\sqcap S}, \text{lfp } f, \text{gfp } f) ( fp f , ⊑ , λ S . lfp ⊔ S , λ S . gfp ⊓ S , lfp f , gfp f ) ,其中 λ S \lambda S λ S 构造一个以 “集合 S S S ” 为输入、以 “对应的不动点” 为输出的函数,于是 λ S . lfp ⊔ S \lambda S.\text{lfp}_{\sqcup S} λ S . lfp ⊔ S 表示对应的最小不动点,λ S . gfp ⊓ S \lambda S.\text{gfp}_{\sqcap S} λ S . gfp ⊓ S 表示对应的最大不动点。
TIP T a r s k i \mathrm{Tarski} Tarski 不动点定理证明了不动点存在性,但并未给出求解不动点的方法。
K l e e n e \mathrm{Kleene} Kleene 不动点定理
IMPORTANT 设 ( D , ⊑ ) (D, \sqsubseteq) ( D , ⊑ ) 为完全偏序 C P O CPO C P O ,若函数 f : D → D f: D \rightarrow D f : D → D 是 ⊔ − \sqcup- ⊔ − 连续的,则有:
l f p f = ⊔ { f i ( ⊥ ) ∣ i ∈ N } \mathrm{lfp}\ f = \sqcup \{f^i(\bot)| i \in \mathbb{N}\} lfp f = ⊔ { f i ( ⊥ ) ∣ i ∈ N } 其中,i ∈ N i \in \mathbb{N} i ∈ N 表示迭代次数。
图示
TIP K l e e n e \mathrm{Kleene} Kleene 不动点定理给出了实际求解不动点的方法:从 C P O CPO C P O 的 ⊥ \bot ⊥ 开始,反复应用函数 f f f 进行迭代来构造序列:⊥ , f ( ⊥ ) , f ( f ( ⊥ ) ) , . . . , f n ( ⊥ ) \bot, f(\bot), f(f(\bot)), ..., f^n(\bot) ⊥ , f ( ⊥ ) , f ( f ( ⊥ )) , ... , f n ( ⊥ ) ,该序列的”最小上确界”就是函数 f f f 的最小不动点。
满足递增链条件的完全偏序 C P O CPO C P O 上的不动点迭代
K l e e n e \mathrm{Kleene} Kleene 不动点定理并未要求迭代过程必须在有限步内终止,其核心实际上是通过无穷迭代的极限来定义不动点,而这种极限在多数情况下无法通过实际计算抵达。于是需要更加严格条件来保证不动点迭代的终止性。
IMPORTANT 若 C P O ( D , ⊑ , ⊔ , ⊥ ) CPO(D, \sqsubseteq, \sqcup, \bot) C P O ( D , ⊑ , ⊔ , ⊥ ) 上没有无穷递增链,f f f 是 D D D 上的单调函数,那么 K l e e n e \mathrm{Kleene} Kleene 迭代序列 a i + 1 = f ( a i ) a_{i+1} = f(a_i) a i + 1 = f ( a i ) 将会在有穷时间内收敛于 l f p ⊥ f \mathrm{lfp}_\bot \ f lfp ⊥ f 。
TIP 有穷高度的链肯定满足递增链条件A C C ACC A C C ,但是满足递增链条件的链不一定是有穷高度(比如有可能存在不稳定的递减链)。
而对 C P O CPO C P O 而言存在最小元 ⊥ \bot ⊥ ,于是满足递增链条件的 C P O CPO C P O 一定没有无穷递增链。
在一些问题上,如:
有穷格 :即格中元素个数是有穷的,如符号格;
满足递增链条件的无穷格 :即格中元素个数是无穷的,但是格的高度是有穷的,如常量格.
上述定理提供了迭代方法来精确地计算最小不动点,并且保证了迭代计算的终止性。
具体应用
数据流系统
之前我们学习的可用表达式分析 和活跃变量分析 也属于数据流系统,这里给出其形式化定义:
一个数据流系统 S = ( L a b , E , F , ( D , ⊑ ) , ι , φ ) \mathcal{S} = (Lab, E, F, (D, \sqsubseteq), \iota, \varphi) S = ( L ab , E , F , ( D , ⊑ ) , ι , φ ) 由以下部分组成:
有限的(程序)标签集 L a b Lab L ab (此处为 L a b c Lab_c L a b c );
极值标签集 E ⊆ L a b E \subseteq Lab E ⊆ L ab (此处为 { init ( c ) } \{\text{init}(c)\} { init ( c )} 或 final ( c ) \text{final}(c) final ( c ) );
流关系 F ⊆ L a b × L a b F \subseteq Lab \times Lab F ⊆ L ab × L ab (此处为 flow ( c ) \text{flow}(c) flow ( c ) 或 flow R ( c ) \text{flow}^R(c) flow R ( c ) );
满足递增链条件(A C C ACC A C C )的完全格 ( D , ⊑ ) (D, \sqsubseteq) ( D , ⊑ ) (带有上确界(LUB)算子 ⨆ \bigsqcup ⨆ 和最小元 ⊥ \bot ⊥ );
极值标签对应的极值值 ι ∈ D \iota \in D ι ∈ D ;
一组单调传递函数 { φ l ∣ l ∈ L a b } \{\varphi_l \mid l \in Lab\} { φ l ∣ l ∈ L ab } ,类型为 φ l : D → D \varphi_l: D \rightarrow D φ l : D → D 。
于是我们将可用表达式分析与活跃变量分析映射到上面的形式化定义中:
问题 可用表达式(Available Expressions) 活跃变量(Live Variables) E E E { init ( c ) } \{\text{init}(c)\} { init ( c )} final ( c ) \text{final}(c) final ( c ) F F F flow ( c ) \text{flow}(c) flow ( c ) flow R ( c ) \text{flow}^R(c) flow R ( c ) D D D 2 C E x p c 2^{CExp_c} 2 C E x p c 2 V a r c 2^{Var_c} 2 V a r c ⊑ \color{red}{\sqsubseteq} ⊑ ⊇ \color{red}{\supseteq} ⊇ ⊆ \color{red}{\subseteq} ⊆ ⨆ \bigsqcup ⨆ ⋂ \bigcap ⋂ ⋃ \bigcup ⋃ ⊥ \color{red}{\bot} ⊥ C E x p c CExp_c C E x p c ∅ \emptyset ∅ ι \iota ι ∅ \emptyset ∅ V a r c Var_c V a r c φ l \varphi_l φ l φ l ( d ) = ( d ∖ kill ( B l ) ) ∪ gen ( B l ) \varphi_l(d) = (d \setminus \text{kill}(B^l)) \cup \text{gen}(B^l) φ l ( d ) = ( d ∖ kill ( B l )) ∪ gen ( B l ) -
注意表中标红的部分。对于可用表达式分析,偏序关系 ⊑ \sqsubseteq ⊑ 定义为 ⊇ \color{red}{\supseteq} ⊇ (即集合包含关系的逆),这里可以这样理解:不动点迭代是 “沿递增链向上” 的过程(即从 ⊥ \bot ⊥ 出发,每一步 d i ⊑ d i + 1 d_i \sqsubseteq d_{i+1} d i ⊑ d i + 1 逐步逼近,直到稳定),最终达到的稳定的结果即为最小不动点,因此最终稳定的结果应该是一个精确的结果,反映到可用表达式集合中表现为:集合越小,精度越高。所以 “向上迭代” 需要对应集合从大到小收缩,即从最粗略的 ⊥ \bot ⊥ 逐步收缩到最精确的最小不动点。于是 “迭代向上”就必须要把偏序 ⊒ \sqsupseteq ⊒ 定义为集合的 ⊇ \supseteq ⊇ 。
数据流方程组与不动点
数据流方程组(d a t a f l o w e q u a t i o n s y s t e m \mathrm{dataflow\ equation\ system} dataflow equation system )
给定数据流系统 S = ( L a b , E , F , ( D , ⊑ ) , ι , φ ) \mathcal{S} = (Lab, E, F, (D, \sqsubseteq), \iota, \varphi) S = ( L ab , E , F , ( D , ⊑ ) , ι , φ ) ,其中 L a b = { 1 , … , n } Lab = \{1, \dots, n\} L ab = { 1 , … , n } (不失一般性)
S \mathcal{S} S 确定方程组(其中 l ∈ L a b l \in Lab l ∈ L ab ):
Al l = { ι i f l ∈ E ⨆ { φ l ′ ( Al l ′ ) ∣ ( l ′ , l ) ∈ F } o t h e r w i s e \text{Al}_l = \begin{cases}
\iota & if\ l \in E \\
\bigsqcup \{\varphi_{l'}(\text{Al}_{l'}) \mid (l', l) \in F\} & otherwise
\end{cases} Al l = { ι ⨆ { φ l ′ ( Al l ′ ) ∣ ( l ′ , l ) ∈ F } i f l ∈ E o t h er w i se
( d 1 , … , d n ) ∈ D n (d_1, \dots, d_n) \in D^n ( d 1 , … , d n ) ∈ D n 被称为解 ,若:
d l = { ι i f l ∈ E ⨆ { φ l ′ ( d l ′ ) ∣ ( l ′ , l ) ∈ F } o t h e r w i s e d_l = \begin{cases}
\iota & if\ l \in E \\
\bigsqcup \{\varphi_{l'}(d_{l'}) \mid (l', l) \in F\} & otherwise
\end{cases} d l = { ι ⨆ { φ l ′ ( d l ′ ) ∣ ( l ′ , l ) ∈ F } i f l ∈ E o t h er w i se
S \mathcal{S} S 确定变换:
Φ S : D n → D n : ( d 1 , … , d n ) → ( d 1 ′ , … , d n ′ ) \Phi_{\mathcal{S}}: D^n \rightarrow D^n: (d_1, \dots, d_n) \rightarrow (d_1', \dots, d_n') Φ S : D n → D n : ( d 1 , … , d n ) → ( d 1 ′ , … , d n ′ )
其中
d l ′ : = { ι i f l ∈ E ⨆ { φ l ′ ( d l ′ ) ∣ ( l ′ , l ) ∈ F } o t h e r w i s e d_l' := \begin{cases}
\iota & if\ l \in E \\
\bigsqcup \{\varphi_{l'}(d_{l'}) \mid (l', l) \in F\} & otherwise
\end{cases} d l ′ := { ι ⨆ { φ l ′ ( d l ′ ) ∣ ( l ′ , l ) ∈ F } i f l ∈ E o t h er w i se
( d 1 , … , d n ) ∈ D n (d_1, \dots, d_n) \in D^n ( d 1 , … , d n ) ∈ D n 求解该方程组 当且仅当它是 Φ S \Phi_{\mathcal{S}} Φ S 的不动点 。
通过不动点迭代来求解数据流问题
要点分析
( D , ⊑ ) (D, \sqsubseteq) ( D , ⊑ ) 是完全格 ,确保 Φ S \Phi_S Φ S 良定义。
由于 ( D , ⊑ ) (D, \sqsubseteq) ( D , ⊑ ) 是满足 A C C ACC A C C 的完全格,因此 ( D n , ⊑ n ) (D^n, \sqsubseteq^n) ( D n , ⊑ n ) 也是完全格(其中 ( d 1 , … , d n ) ⊑ n ( d 1 ′ , … , d n ′ ) (d_1, \dots, d_n) \sqsubseteq^n (d_1', \dots, d_n') ( d 1 , … , d n ) ⊑ n ( d 1 ′ , … , d n ′ ) 当且仅当对所有 1 ≤ i ≤ n 1 \leq i \leq n 1 ≤ i ≤ n ,有 d i ⊑ d i ′ d_i \sqsubseteq d_i' d i ⊑ d i ′ )。
传递函数 φ l \varphi_l φ l 在 ( D , ⊑ ) (D, \sqsubseteq) ( D , ⊑ ) 中的单调性,蕴含 Φ S \Phi_S Φ S 在 ( D n , ⊑ n ) (D^n, \sqsubseteq^n) ( D n , ⊑ n ) 中的单调性(因为上确界算子 ⨆ \bigsqcup ⨆ 也具有单调性)。
因此,(最小)不动点可通过迭代有效计算:
fix ( Φ S ) = ⨆ { Φ S k ( ⊥ D n ) ∣ k ∈ N } \text{fix}(\Phi_S) = \bigsqcup \{\Phi_S^k(\bot_{D^n}) \mid k \in \mathbb{N}\} fix ( Φ S ) = ⨆ { Φ S k ( ⊥ D n ) ∣ k ∈ N }
其中 ⊥ D n = ( ⊥ D , … , ⊥ D ) ⏟ n 次 \bot_{D^n} = \underbrace{(\bot_D, \dots, \bot_D)}_{n \text{ 次}} ⊥ D n = n 次 ( ⊥ D , … , ⊥ D ) 。
若 ( D , ⊑ ) (D, \sqsubseteq) ( D , ⊑ ) 的高度为 m m m
⟹ ( D n , ⊑ n ) \implies (D^n, \sqsubseteq^n) ⟹ ( D n , ⊑ n ) 的高度为 m ⋅ n m \cdot n m ⋅ n
⟹ \implies ⟹ 不动点迭代最多需要 m ⋅ n m \cdot n m ⋅ n 步。
Comments