第一步:標記化
處理運算式的第一步就是將其轉化為包含一個個獨立符號的列表。這一步很簡單,且不是本文的重點,因此在此處我省略了很多。
首先,我定義了一些標記(數字不在此中,它們是預設的標記)和一個標記類型:
token_map = {'+':'ADD', '-':'ADD', '*':'MUL', '/':'MUL', '(':'LPAR', ')':'RPAR'} Token = namedtuple('Token', ['name', 'value'])
下面就是我用來標記 `expr` 運算式的代碼:
split_expr = re.findall('[\d.]+|[%s]' % ''.join(token_map), expr)tokens = [Token(token_map.get(x, 'NUM'), x) for x in split_expr]
第一行是將運算式分割為基本標記的技巧,因此
'1.2 / ( 11+3)' --> ['1.2', '/', '(', '11', '+', '3', ')']
下一行命名標記,這樣分析器就能通過分類識別它們:
['1.2', '/', '(', '11', '+', '3', ')']->[Token(name='NUM', value='1.2'), Token(name='MUL', value='/'), Token(name='LPAR', value='('), Token(name='NUM', value='11'), Token(name='ADD', value='+'), Token(name='NUM', value='3'), Token(name='RPAR', value=')')]
任何不在 token_map 中的標記被假定為數字。我們的分詞器缺少稱為驗證的屬性,以防止非數字被接受,但幸運的是,運算器將在以後處理它。
就是這樣
第二步: 文法定義
我選擇的解析器實現自一個本地垂直解析器,其來源於LL解析器的一個簡單版本。它是一個最簡單的解析器實現,事實上,只有僅僅14行代碼。它是一種自上而下的解析器,這意味著解析器從最上層規則開始解析(like:expression),然後以遞迴方式嘗試按照其子規則方式解析,直至符合最下層的規則(like:number)。換句話解釋,當自底向上解析器(LR)逐步地收縮標記,使規則被包含在其它規則中,直到最後僅剩下一個規則,而自頂向下解析器(LL)逐步展開規則並進入到少數的抽象規則,直到它能夠完全符合輸入的標記。
在深入到實際的解析器實現之前,我們可對文法進行討論。在我之前發表的文章中,我使用過LR解析器,我可以像如下方式定義計算機文法(標記使用大寫字母表示):
add: add ADD mul | mul;mul: mul MUL atom | atom;atom: NUM | '(' add ')' | neg;neg: '-' atom;
(如果您還不理解上述文法,請閱讀我之前發表的文章)
現在我使用LL解析器,以如下方式定義計算機的文法:
rule_map = { 'add' : ['mul ADD add', 'mul'], 'mul' : ['atom MUL mul', 'atom'], 'atom': ['NUM', 'LPAR add RPAR', 'neg'], 'neg' : ['ADD atom'],}
大家可以看到,這裡有一個微妙的變化。有關"add and mul"的遞迴定義被反轉了。這是個非常重要的細節,我會向大家詳細說明這一點。
LR版本使用了左遞迴的模式。當LL解析器遇到遞迴的時候,它會嘗試去匹配規則。所以,當左遞迴發生是,解析器會進入無窮遞迴。甚至連聰明的LL解析器例如ANTLR也逃避不了這個問題,它會以友好的錯誤提示代替無窮的遞迴,而不像我們這個玩具解析器那樣。
左遞迴可以很容易的轉變為右遞迴,我就這麼做的。但是解析器並不是那麼簡單,它又會產生另一個問題:當左遞迴正確的解析 3-2-1 為(3-2)-1,而右遞迴卻錯誤的解析為3-(2-1)。我還沒想到一個簡單的解決辦法,所以為了讓事情簡單,我決定讓它繼續使用錯誤的解析格式,並在後面處理這個問題(請看步驟4)
第三步:解析為一個AST
演算法其實很簡單。我們會定義一個接收兩個參數的遞迴方法:第一個參數是我們要嘗試匹配的規則名稱,第二個參數是我們要保留的識別欄位表。我們從add(最上層規則)方法開始,其已包含完整的識別欄位表,遞迴調用已非常明確。方法將返回一個數組,其包含元素為:一個是當前匹配項,另一個是保留匹配的識別欄位表。我們將實現標識匹配功能,以使這段代碼可用(它們都是字串類型;一個是大寫格式,另一個是小寫格式)。
以下是解析器實現的代碼:
RuleMatch = namedtuple('RuleMatch', ['name', 'matched']) def match(rule_name, tokens): if tokens and rule_name == tokens[0].name: # 是否匹配標識? return RuleMatch(tokens[0], tokens[1:]) for expansion in rule_map.get(rule_name, ()): # 是否匹配規則? remaining_tokens = tokens matched_subrules = [] for subrule in expansion.split(): matched, remaining_tokens = match(subrule, remaining_tokens) if not matched: break # 運氣不好,跳出迴圈,處理下一個擴充定義! matched_subrules.append(matched) else: return RuleMatch(rule_name, matched_subrules), remaining_tokens return None, None # 無匹配結果
代碼4至5行說明:如果規則名稱(rule_name)確實是一個標識,並被包含在識別欄位表(tokens)中,同時檢查其是否匹配當前標識。如果是,運算式將返回匹配方法,識別欄位表任然進行使用。
代碼第6行說明:迭代將迴圈檢查是否匹配該規則名稱對應的子規則,通過遞迴實現每條子規則的匹配。如果規則名稱滿足匹配標識的條件,get()方法將返回一個空數組,同時代碼將返回空值(見16行)。
第9-15行,實現迭代當前的sub-rule,並嘗試順序地匹配他們。每次迭代都儘可能多的匹配標識。如果某一個標識無法匹配,我們就會放棄整個sub-rule。但是,如果所有的標識都匹配成功,我們就到達else語句,並返回rule_name的匹配值,還有剩下標識。
現在運行並看看1.2/(11+3)的結果。
>>> tokens = [Token(name='NUM', value='1.2'), Token(name='MUL', value='/'), Token(name='LPAR', value='('), Token (name='NUM', value='11'), Token(name='ADD', value='+'), Token(name='NUM', value='3'), Token(name='RPAR', value=')')] >>> match('add', tokens) (RuleMatch(name='add', matched=[RuleMatch(name='mul', matched=[RuleMatch(name='atom', matched=[Token(name='NUM', value='1.2')]), Token(name='MUL', value='/'), RuleMatch(name='mul', matched=[RuleMatch(name='atom', matched=[Token(name='LPAR', value='('), RuleMatch(name='add', matched=[RuleMatch(name='mul', matched=[RuleMatch(name='atom', matched=[Token(name='NUM', value='11')])]), Token(name='ADD', value='+'), RuleMatch(name='add', matched=[RuleMatch(name='mul', matched=[RuleMatch(name='atom', matched=[Token(name='NUM', value='3')])])])]), Token(name='RPAR', value=')')])])])]), [])
結果是一個tuple,當然我們並沒有看到有剩下的標識。匹配結果並不易於閱讀,所以讓我吧結果畫成一個圖:
add mul atom NUM '1.2' MUL '/' mul atom LPAR '(' add mul atom NUM '11' ADD '+' add mul atom NUM '3' RPAR ')'
這就是概念上的AST。通過你思維邏輯,或者在紙上描繪,想象解析器是如何運作的,這樣是個很好的鍛煉。我不敢說這樣是必須的,除非你想神交。你可以通過AST來協助你實現正確的演算法。
到目前為止,我們已經完成了可以處理二進位運算,一元運算,括弧和操作符優先權的解析器。
現在只剩下一個錯誤待解決,下面的步驟我們將解決這個錯誤。
第四步:後續處理
我的解析器並非在任何場合管用。最重要的一點是,它並不能處理左遞迴,迫使我把代碼寫成右遞迴方式。這樣導致,解析 8/4/2 這個運算式的時候,AST結果如下:
add mul atom NUM 8 MUL '/' mul atom NUM 4 MUL '/' mul atom NUM 2
如果我們嘗試通過AST計算結果,我們將會優先計算4/2,這當然是錯誤的。一些LL解析器選擇修正樹裡面的關聯性。這樣需要編寫多行代碼;)。這個不採納,我們需要使它扁平化。演算法很簡單:對於AST裡面的每個規則 1)需要修正 2)是一個二進位運算 (擁有sub-rules)3) 右邊的操作符同樣的規則:使後者扁平成前者。通過“扁平”,我意思是在其父節點的上下文中,通過節點的兒子代替這個節點。因為我們的穿越是DFS是後序的,意味著它從樹的邊緣開始,並一直到達樹根,效果將會累加。如下是代碼:
fix_assoc_rules = 'add', 'mul' def _recurse_tree(tree, func): return map(func, tree.matched) if tree.name in rule_map else tree[1] def flatten_right_associativity(tree): new = _recurse_tree(tree, flatten_right_associativity) if tree.name in fix_assoc_rules and len(new)==3 and new[2].name==tree.name: new[-1:] = new[-1].matched return RuleMatch(tree.name, new)
這段代碼可以讓任何結構的加法或乘法運算式變成一個平面列表(不會混淆)。括弧會破壞順序,當然,它們不會受到影響。
基於以上的這些,我可以把代碼重構成左關聯:
def build_left_associativity(tree): new_nodes = _recurse_tree(tree, build_left_associativity) if tree.name in fix_assoc_rules: while len(new_nodes)>3: new_nodes[:3] = [RuleMatch(tree.name, new_nodes[:3])] return RuleMatch(tree.name, new_nodes)
但是,我並不會這樣做。我需要更少的代碼,並且把計算代碼換成處理列表會比重構整棵樹需要更少的代碼。
第五步:運算器
對樹的運算非常簡單。只需用與後處理的代碼相似的方式對樹進行遍曆(即 DFS 後序),並按照其中的每條規則進行運算。對於運算器,因為我們使用了遞迴演算法,所以每條規則必須只包含數字和操作符。代碼如下:
bin_calc_map = {'*':mul, '/':div, '+':add, '-':sub}def calc_binary(x): while len(x) > 1: x[:3] = [ bin_calc_map[x[1]](x[0], x[2]) ] return x[0] calc_map = { 'NUM' : float, 'atom': lambda x: x[len(x)!=1], 'neg' : lambda (op,num): (num,-num)[op=='-'], 'mul' : calc_binary, 'add' : calc_binary,} def evaluate(tree): solutions = _recurse_tree(tree, evaluate) return calc_map.get(tree.name, lambda x:x)(solutions)
我使用 calc_binary 函數進行加法和減法運算(以及它們的同階運算)。它以左結合的方式計算資料行表中的這些運算,這使得我們的 LL文法不太容易擷取結果。
第六步:REPL
最樸實的REPL:
if __name__ == '__main__': while True: print( calc(raw_input('> ')) )
不要讓我解釋它 :)
附錄:將它們合并:一個70行的計算機
'''A Calculator Implemented With A Top-Down, Recursive-Descent Parser'''# Author: Erez Shinan, Dec 2012 import re, collectionsfrom operator import add,sub,mul,div Token = collections.namedtuple('Token', ['name', 'value'])RuleMatch = collections.namedtuple('RuleMatch', ['name', 'matched']) token_map = {'+':'ADD', '-':'ADD', '*':'MUL', '/':'MUL', '(':'LPAR', ')':'RPAR'}rule_map = { 'add' : ['mul ADD add', 'mul'], 'mul' : ['atom MUL mul', 'atom'], 'atom': ['NUM', 'LPAR add RPAR', 'neg'], 'neg' : ['ADD atom'],}fix_assoc_rules = 'add', 'mul' bin_calc_map = {'*':mul, '/':div, '+':add, '-':sub}def calc_binary(x): while len(x) > 1: x[:3] = [ bin_calc_map[x[1]](x[0], x[2]) ] return x[0] calc_map = { 'NUM' : float, 'atom': lambda x: x[len(x)!=1], 'neg' : lambda (op,num): (num,-num)[op=='-'], 'mul' : calc_binary, 'add' : calc_binary,} def match(rule_name, tokens): if tokens and rule_name == tokens[0].name: # Match a token? return tokens[0], tokens[1:] for expansion in rule_map.get(rule_name, ()): # Match a rule? remaining_tokens = tokens matched_subrules = [] for subrule in expansion.split(): matched, remaining_tokens = match(subrule, remaining_tokens) if not matched: break # no such luck. next expansion! matched_subrules.append(matched) else: return RuleMatch(rule_name, matched_subrules), remaining_tokens return None, None # match not found def _recurse_tree(tree, func): return map(func, tree.matched) if tree.name in rule_map else tree[1] def flatten_right_associativity(tree): new = _recurse_tree(tree, flatten_right_associativity) if tree.name in fix_assoc_rules and len(new)==3 and new[2].name==tree.name: new[-1:] = new[-1].matched return RuleMatch(tree.name, new) def evaluate(tree): solutions = _recurse_tree(tree, evaluate) return calc_map.get(tree.name, lambda x:x)(solutions) def calc(expr): split_expr = re.findall('[\d.]+|[%s]' % ''.join(token_map), expr) tokens = [Token(token_map.get(x, 'NUM'), x) for x in split_expr] tree = match('add', tokens)[0] tree = flatten_right_associativity( tree ) return evaluate(tree) if __name__ == '__main__': while True: print( calc(raw_input('> ')) )