Kota's Blog

コンピューターの掛け算と割り算

2021-03-29 Computer Science

この記事ではコンピュータの掛け算と割り算に使われるシフト演算の説明と計算をします。筆者が情報数学の基礎を勉強していて「すげー」と思った部分を書いているので間違っている点があればコメントで指摘していただけると喜びます。

掛け算と割り算

よくある誤解として、コンピューターの掛け算は足し算をループして実装していると思われがちですが、そんなことはありません。
値は全て2進数で扱われるので足し算であれば0010+0010=01000010+0010=0100みたいな計算になるのですが、例えば4×24\times2の計算を

0100+0100=1000=8

みたいに一つずつ足しながらやるかと言えば、こんなクソコーダーみたいな実装はしないのです。
よく見ると、4の2進数01000100と8の2進数10001000だと、1の位置が1ビット左に移動してますね。このように、nnビット分移動した時に乗数が2n2^nになることをシフトといい、これを使った計算をシフト演算というのです。コンピューターはこのシフト演算を使って掛け算と割り算を計算しています。
もう一つの例として8×48\times4の計算をクソコーダー風にすると

1000+1000+1000+1000=00100000=32

となります。これも、8の2進数を0000100000001000と考えると2ビット分移動しており、22=4=乗数2^2=4=乗数が成り立っています。

シフト演算とは?

シフト演算には論理シフトと算術シフトという2種類の演算が存在します。

論理シフト

論理シフトは上の演算でも利用していた通り、レジスタ上のビットを全て左、右へ移動する演算です。移動して空いた端のビットには0が入ります。
上で紹介したシフト演算はこの論理シフトの左シフトといい、左へnnビットシフトすると2n2^n倍されるという演算です。これを使って掛け算が行われています。割り算には右シフトといい、右へnnビットシフトすると2n2^{-n}倍されるという演算が利用されています。

算術シフト

算術シフトは論理シフトと違い、符号を考慮します。レジスタにおいて、左端のビットは1であれば負の数、0であれば正の数といった感じで符号を表します。その符号ビットだけは残したまま他の全てのビットを論理シフトと同じようにずらす演算が算術演算です。

論理シフトしてみよう

ここまで簡単にシフト演算を説明してきましたが、実際に論理シフトを使った計算をしてみましょう。算術シフトに関しては論理シフトと計算方法は近いので割愛します。

問題1

32ビットのレジスタに16進数A6E2が入っている。これを3ビット右に論理シフトした値を16進数で求めよ

問題のレジスタを2進数で表すと0000 0000 0000 0000 1010(←A) 0110(←6) 1110(←E) 0010(←2)になります。ここから3ビット分右に論理シフト、つまり3個分全ての値を右にずらしてみましょう。シフト後の値は0000 0000 0000 0000 0001 0100 1101 1100、これを16進数に直すと、

0001=23×0+22×0+21×0+20×1=10001=2^3\times0+2^2\times0+2^1\times0+2^0\times1=1
0100=23×0+22×1+21×0+20×0=40100=2^3\times0+2^2\times1+2^1\times0+2^0\times0=4
1101=23×1+22×1+21×0+20×1=13=D1101=2^3\times1+2^2\times1+2^1\times0+2^0\times1=13=D
1100=23×1+22×1+21×0+20×0=12=C1100=2^3\times1+2^2\times1+2^1\times0+2^0\times0=12=C

よって答えは14DCになります。


問題2

レジスタに格納されている正の整数xを10倍にする処理を以下から選べ。この時桁が溢れることはない。
1. xを2ビット左にシフトした値にxを加算し、更に1ビット左にシフトする
2. xを2ビット左にシフトした値にxを加算し、更に2ビット左にシフトする
3. xを3ビット左にシフトした値と、xを2ビット左にシフトした値を加算する
4. xを3ビット左にシフトした値にxを加算し、更に1ビット左にシフトする

桁が溢れるというのはシフト演算の過程で1のビットが端から漏れるということです。
先述のシフト演算に関する前提知識があれば、式を作って計算してみれば簡単に答えがわかります。

選択肢1

xを2ビット左にシフトするということはつまり値は222^2倍されることになります。そしてその値にxを加算し(+x+x)、更に1ビット左にシフトする、つまり212^1倍した結果を調べれば10倍かどうか分かるというわけです。

x=(x×22+x)×21=5x×2=10xx'=(x\times2^2+x)\times2^1=5x\times2=10x

ご覧の通り10倍になっているので選択肢1が正解です。試験などで解いた時はこれで次の問題に移って問題ないのですが一応他の選択肢も計算しておきましょう。

選択肢2

x=(x×22+x)×22=5x×4=20xx'=(x\times2^2+x)\times2^2=5x\times4=20x

選択肢3

x=x×23+x×22=8x+4x=12xx'=x\times2^3+x\times2^2=8x+4x=12x

選択肢4

x=(x×23+x)×21=9x×2=18xx'=(x\times2^3+x)\times2^1=9x\times2=18x


問題3

以下の16ビットの固定小数点レジスタの内容を3ビット論理左シフトしたものをa、2ビット論理右シフトしたものをbとした時、aがbの何倍になるのか求めよ。
0000 0100 1011 1000

3ビット論理左シフトしたaa0010 0101 1100 0000、2ビット論理右シフトしたbb0000 0001 0010 1110となります。この問題では桁溢れが起きないため、単純にaaは元のレジスタの内容から232^3倍値が増え、bb222^{-2}倍値が増えた物であると考えられるのです。この時比べたいのはaabbですから、元のレジスタの内容をxxとすると
ab=23xx22=23x×22x=25=32\frac{a}{b}=\frac{2^3x}{\frac{x}{2^2}}=\frac{2^3x\times2^2}{x}=2^5=32
よってaabbの32倍という結果になる。


おまけ-足し算と引き算

コンピューターの足し算は論理演算子を組み合わせて実現します。具体的には以下のような論理回路で計算を行います。
image.png
足し算であれば何不自由なく計算できるコンピューターですが、引き算は少し複雑です。繰り下がりに該当する論理回路が存在しないため足し算のような回路を組むことができないのです。ではどう引き算を実現しているのかというと、補数という値を利用します。
補数とは「任意の数nnに足せば繰り上がりする数」のことです。10進数で言えばn+補数=10n+補数=10の数値ということで、1の補数は9、7の補数は3といった具合に補数が決定します。

補数を使って引き算する

足し算を使って引き算、と言われると違和感しかありませんが補数を使った計算は至って簡単です。例として969-6の計算をしてみたいと思います。
969-6の計算において、補数を求めるnnは6になります(減数)。つまり補数は4になるのですが、ここで被減数(引かれる数)と補数を足してみましょう。

9+4=13

この13の1の位に注目すると、3…なんと元々の引き算と同じ答えが出ているのです。被減数と補数を足して最上位の桁を取り除く、この方法でコンピューターは引き算をしています。感動

あれ、でも補数ってどうやって導くの?

賢い方は気づいたはず、上記の方法で引き算をしようとしても、そもそも補数を求めるのに引き算使ってるから意味ないじゃん、と。なんとなくブートストラップ問題的な、鶏卵な問題になりそうですが、実は2進数を機械的に扱うことで解決してしまうんです。

まず補数を算出する対象xxを4とします。この4を2進数で表すと0100となります。**この2進数を反転させたビットに1を足すと、補数が求まってしますのです!**計算して確かめてみましょう。

1011+1=11001011+1=1100

1001は10進数に直すと12になり、本来4の補数は6であるべきなので間違いのように思えますが、この6は2の補数と言い、2進数で計算する上でxxと足しても繰り上がりをしない最大の数なのです。つまり今までこのおまけ章で扱ってきた補数は10進数用、10の補数だったというわけです。2の補数を10進数で使っても上手くいくはずがないので、被減数を8(=1000)として2進数で4の補数を足してみましょう。

1000+1100=101001000+1100=10100
最上位の桁を無視すると0100=40100=4、10進数で計算した84=48-4=4と一致し、足し算だけで引き算を実現することができました。

終わりに

僕はC言語やJavaを書いたことがないので実際にシフト演算をコードで実装したことはありませんが、いつの日か低レイヤーの鬼になってレジスタの実装とかしてみたいですね。

参考文献

算術シフト - Wikipedia
論理シフト - Wikipedia
2の補数
シフト演算