calcRPN関数で電卓を作ろう

calcRPN
1。逆ポーランド式の解釈はどうやる?
[b]このワークシートは[url=https://www.geogebra.org/m/twxxx3yq]Math by Code[/url]の一部です。[br][/b][br](1*2)+(4-3)という中置式から木を作り、逆ポーランド記法ReversePolandNotation, RPNで出力すると、[br][b]1, 2, *, 4 ,3,- になることがわかった。[br]カンマが邪魔なので省くと、1 2 * 4 3 -だ。[br]では、逆ポーランド記法を通常の式に戻すことをやってみよう。[br][br][/b][br][color=#9900ff][b][u][size=150]質問:演算子が後置の逆ポーランド記法RPN「4 3 2 + 1[b][u][size=150] ー[/size][/u][/b] *」を普通の中置式にするとどうなるでしょうか。[br][/size][/u][/b][/color][br]4 3 2 + 1 - *[br]演算子だけとりだすと、+, - , *だ。[br]だから+をして、ーをして、*をしている。[br]では何と何に作用させているんだろう。[br][br]4 [u]3 2[u]+[/u][/u] 1 ー *[br]まず、+の直前の2数3と2を取り出して、(3+2)にして戻す。[br]4 [u](3+2) 1 ー [/u]*[br]次は、ーの直前の2数(3+2)と1 を取り出して、(3+2)ー1にして戻そう。[br][u]4 (3+2) [u]ー[/u]1[u]*[/u] [/u][br]最後に*の直前の2数4 と (3+2) ー1 を取り出して、4 *( (3+2) ー1 )にして戻そう。[br]だから、4 * ( (3+2) ー1 )となるね。[br]この手順は目で見ながら実行したのだけれど、空白を区切りとして、数値を演算子を読み取ることで[br]ルールから今の作業を実行したい。[br][br][color=#9900ff][b][u][size=150]質問:後置式を中置式に変換して返すプログラムはどうやればできるだろうか。[br][/size][/u][/b][/color][br]さて、さて、データ構造といえばスタックかキューかハッシュの出番かな?[br]左から順に、つまり入れた順に使うわけではないから、スタックでしょうか。[br]stackにSさんと名前をつけよう。[br][br]RPN=「4 [u]3 2[u]+[/u][/u] 1 ー *」を空白区切りでSさんにわたそう。[br]Sさんは数値ならどんどん上に積んでいく。[br]Sに左から4、3、2の順を渡すと、上か順に[b][u]x、y、残りをysと名付けるなら、[/u][/b]x=2, y=3, ys=4となる。[br][color=#0000ff]もとの式の左右の順番とx,yの順番が逆になってしまうので計算順に注意しよう。[br][/color]次は(+)が出たから、上からx=2,y=3を取り出して、y+x=3+2=5をSさんに戻す。[br]残った4の上に5が積まれるので、上から順にx=5,y=4となる。[br]次の1はSさんに渡す。Sさんは上から順にx=1,y=5,ys=4となる。[br]次は(ー)が出たから、上からx=1,y=5を取り出して、y-x=5-1=4を戻す。[br]Sさんは上から順にx=4,y=4となる。[br]次も(*)が出たから、上からx=4, y=4を取り出し、y*x=4*4=16を戻す。[br]x=16が入っている。[br]RPNが空になったとき、Sさんのスタックが1つの数値だけになりそれが答えになるね。[br]
2。Haskellで電卓を作る。
複数の演算子があると理解しにくいので、まずは、1つだけで作ってみよう。[br][br]RPN=文字列を受取り、和を返す関数calcDiff( RPN)を作ってみよう。[br][br]文字列RPNをリストLに分解するwords関数を使う。[br]Lを左から順に読み取って処理は、左畳み込みfoldl関数にfoldStack関数にLと空のリストを渡す。[br]最後に畳み込んだあとのリストの先頭をhead関数で読み込む。[br]この3つの関数をパイプラインでつないだ、合成関数がcalcDiffです。[br]ただし、foldStack関数だけは未定義なので、where句で末尾にケース分けして反応を定義します。[br]foldStack関数が使うスタック(だたのリストです)Sは最初は空です。S=[ ][br]・Lからの読み取りが数値numなら, Sの中身xsの先頭にnumを追加する。Sの先頭から順にx,y,ysとなる。[br]・Lからの読み取りが”-”なら,Sの先頭からx,yなので、y-xを計算してSの先頭に追加する。[br]最後にSの先頭(head)を読み取る。[br][br]これをHaskellでかくとこうなります。[br][color=#0000ff][b]calcDiff = head . foldl foldStack [] . words [br] where [br] foldStack (x:y:ys) "-" =(y-x):ys[br] foldStack xs num = read num :xs[br][/b][/color]これを使うには、[br]コマンドラインでカレントディレクトリをプログラムファイルのある場所に移動し、[br]ghciと打ちHaskellを立ち上げてから、プログラムファイルをロード(:l ファイル名)します。[br][b]ghci>calcDiff " 10 - 2 "[br]8[br][b]ghci>[/b]calcDiff " 2 - 10 "[br]-8[br][/b]となるので、x,y の逆順で(y-x)を計算して残りysに積んでいるのが効いてますね。[br][br][b][color=#9900ff][u][size=150]質問:RPN文字列を読み取り、たし算以外の四則計算などもできる関数calcRPNはどうやれば作れますか。[br][/size][/u][/color][/b][br]簡単です。読み取った演算記号が-以外の場合の式をwhere句に、好きなだけ追加するだけです。[br]たとえば、sum や ln など、haskellが読める関数がそのままかけます。[br][br][color=#0000ff][b][size=150]calcRPN = head . foldl foldStack [] . words [br] where [br] foldStack (x:y:ys) "+" =(y+x):ys[br] foldStack (x:y:ys) "-" =(y-x):ys[br] foldStack (x:y:ys) "*" =(y*x):ys[br] foldStack (x:y:ys) "/" =(y/x):ys[br] foldStack (x:y:ys) "^" =(y**x):ys[br] foldStack (x:xs) "ln" = log x:xs[br] foldStack xs "sum" =[sum xs][br] foldStack xs num = read num :xs[br][/size][/b][/color][br][b]ghci>[/b]calcRPN "4 3 2 + 1 - *"[br]16 [br][b]ghci>[/b]calcRPN "2 10 ^"[br]1024.0[br][b]ghci>[/b]calcRPN "2.71828 ln"[br]0.999999327347282[br][b]ghci>[/b]calcRPN "10 20 20 sum"[br]60.0[br][br]やはり、haskellのような言語は、関数を組み合わせる、加工型、関数型プログラムには向いてますね。[br]すばらしい。
3.juliaで電卓を作る。
[color=#9900ff][b][size=150]質問:関数名を同じにしてjuliaでも電卓を作るにはどうしたらよいでしょうか。[br][/size][/b][/color][br][b]juliaにもfoldl,foldr関数がある[/b]のですが、簡単な演算のくり返しでは作動しますが、[br]一般の関数は難しいようです。[br]しかたないので、whileで代用しましょう。[br][br][b]words[/b]関数と[b]foldStack[/b]関数を作っておき、単体で作動テストをしてから、[br]それらを組み込んだ[b]calcRPN[/b]関数を作ります。[br][br][color=#0000ff]function [b]words[/b](x) # 1スペースで区切る[br] return split(x," ")[br]end[br]function [b]foldStack[/b](lst)[br] ope = ["*","+","-","/"][br] opeS = Set(ope)[br] Stack =[][br] item = popfirst!(lst)[br] # y x ope を抜き出して Stack の上からx ,y の順につまれた数で[br] # z= y ope x を計算します。[br] while item ∉ opeS[br] pushfirst!(Stack,item)[br] item = popfirst!(lst)[br] end[br] if item in ope[br] x = parse.(Int,popfirst!(Stack))[br] y = parse.(Int,popfirst!(Stack))[br] z = 0[br] if item == ope[1][br] z = y * x[br] elseif item == ope[2][br] z = y + x[br] elseif item == ope[3][br] z = y - x[br] elseif item == ope[4][br] z = y / x[br] end[br] end[br] #取り出した3要素y,x,opeを1要素zにしてリストlstに戻し、[br] #計算しなかったStackの残りも、リストlstにもどします。[br] pushfirst!(Stack,z)[br] while length(Stack)>0[br] p = string(popfirst!(Stack))[br] pushfirst!(lst,p)[br] end[br] return lst[br]end[br][br][/color][color=#0000ff]function [b]calcRPN[/b](RPN)[br] res = RPN |> words[br] while length(res)>2 [br] res = foldStack(res)[br] end[br] return parse.(Int,res)[br]end[br][br]"4 3 2 + 1 - *" |> calcRPN[br][/color]#==========================[br][OUT] [br]1-element Vector{Int64}:[br] 16[br][br]RPNの数式ならば加減乗除が自由にできます。[br]二項演算の種類を増やしたり、単項演算,多項演算もできるようにすることも可能ですが、[br]あちこち行数が増えてしまいます。[br]まあ、juliaは演算「//」ですぐに分数計算に切り替わるので、分数電卓を作りたかったら向いているかもしれません。[br][br][color=#9900ff][b][u][size=150]質問:このcalcRPN関数のもとになったfoldStack関数を改造して、もっと計算できる演算の種類をふやしてみましょう。[br]また、geogebraでcalcRPN電卓を作るにはどうしたらよいでしょうか。[br][br][/size][/u][/b][/color]geogebraで作るにはRPN式をスペースで区切りリストlstにします。[br]演算個数opedを読み取ります。使われている演算記号リストopeinを作ります。[br]最初の演算記号の位置を読み取ったら、opePos1に格納します。[br][br]opePos1までの3要素をrpn1というリストにします。[br]rpn1を中置順に直して、textToNumberで計算して、[br]もとのリストから3個削除して、計算結果v1を挿入します。[br]すると、lstからopein(1)だけを計算したリストlst1が作られます。[br][br]もし、演算個数opedが1より大きかったら、同様にしてlst1のrpn2を計算結果v2に置き換えたlst2を[br]作ります。opedが1より大きくなかったら、v2はv1のままでlst2はlst1のままにします。[br][br]もし、演算個数opedが2より大きかったら、同様にしてrpn3を計算して、結果answを得ます。大きくなかったら、answはv2のままです。
calcRPN

Information: calcRPN関数で電卓を作ろう