|
|
/*========================================================*/ ＜＜＜あゆしゃのＣ言語プログラミング＞＞＞ /*========================================================*/ 第４３６回 大型計算機 割り算 発行 ２００４年３月３１日（水曜日） 発行数 約２９００ {magclick} /*========================================================*/ はじめに （ 決り文句 ） /*========================================================*/ ・このメールマガジンはまぐまぐさんから発行しています。 ・ジャンルは、マルチメディアのプログラム、Ｃ言語です。 ・このメールマガジンは、横６０文字で作成しています。 また、インデントはすべて半角スペース４つで構成しています。 ・ここで扱うプログラムは、C言語と半光年以内のものです。 ・登録解除は、まぐまぐさんのホームページでお願いします。 ・まぐまぐさんのバックナンバー（下欄参照）を活用して下さい。 ・ここは私の復習の場です。内容は約１ヶ月内外に私が勉強した 内容になっています。最新の技術があれば、へたれもあります。 ・わかりやすさを優先させる為、たまに嘘があるかもしれません。 ・セキュリティ突破のため、暗号化された単語があります。 /*========================================================*/ ご挨拶 /*========================================================*/ こんにちは。あゆしゃです。 少し時間が空いてしまいました。 別に、ドdラrクエをやっていたわけでは、まったくそうです。 山登りにてこずりまして、ようやくグランバニアです。 {magclick} /*========================================================*/ 今回のお題 ＜＜ 大型計算機 割り算 ＞＞ /*========================================================*/ のまえに、前回、全然間違いでした。 CXInt& mul( CXInt& a ); for( i = 0; i < SIZE * 8; i++ ) { if( *a.m_buff & 1 ) add( a ); shl(); } return *this; こんな感じでしょうか、と紹介しましたが、これでは全然駄目 でした。 どうも、何かのために盲目であったようです。 紹介した、UINTの掛け算は、こんな感じでした。 UINT ret = 0; int i = 32; while( --i >= 0 ) { if( b & 1 ) ret += a; b >>= 1, a <<= 1; } return ret; ループカウンタを使っているところが手抜きですが、 見比べてみると、シフトの数が足りません。変数の数が足りませ ん。何か、すごく変です！ 次のは、動作確認済みの掛け算です。 UINT mul_uint4( UINT a, UINT b ) { UINT ret = 0; do { if( b & 1 ) ret += a; a <<= 1; } while( b >>= 1 ); return ret; } ループカウンタを使う必要のない形にして、簡略化しました。 これにあわせてCXIntの掛け算を作りましょう。 CXInt& mul( CXInt& b ); CXInt a( *this ); clear(); do { if( *b.m_buff & 1 ) add( a ); a.shl(); b.shr(); } while( b.noZero() ); return *this; こんな感じでしょうね。 /*========================================================*/ 前回までで、なんとか掛け算関数まで作りました。 class CXInt; enum { SIZE = 1024; } UCHAR m_buff[ SIZE ]; /*========================================================*/ CXInt& add( CXInt& a ); CXInt& sub( CXInt& a ); CXInt& mul( CXInt& b ); /*========================================================*/ 今回は、割り算です。 以前紹介したUINT版の割り算は、動作確認してみると、いろいろ と悪い点がありました。 動作確認した結果、以下のほうが良いです。 UINT div_uint4( UINT a, UINT b ) { UINT ret = 0; UINT mask = 1; while( a > b && ( I4 )b >= 0 ) { b <<= 1; if( ( I4 )( mask <<= 1 ) < 0 ) break; } while( a ) { if( a >= b ) ret += mask, a -= b; b >>= 1; if( ! ( mask >>= 1 ) ) break; } return ret; } 基本的な考え方は変わっていませんが、０を上手く演算できない などの「バグ」を修正してあるほか、ループカウンタを使わない 工夫を施してあります。 これを踏まえて、CXInt版です。 CXInt& div( CXInt& b ); CXInt a( *this ); clear(); CXInt mask; *mask.m_buff = 1; while( a.cmp( b ) > 0 && b.noNeg() ) { b.shl(); mask.shl(); if( mask.isNeg() ) break; } while( a.noZero() ) { if( a.cmp( b ) >= 0 ) { add( mask ); a.sub( b ); } b.shr(); mask.shr(); if( mask.isZero() ) break; } return *this; ふぅ、こんな感じでしょうか。 {magclick} /*========================================================*/ さいごに /*========================================================*/ {magclick} /*========================================================*/ 次回予告 /*========================================================*/ 次回は４月５日（月曜日）に、第４３７回を送ります。 お題は「大型計算機 シフトなど」 平方根の前に、掛け算割り算で必要となる、雑多な関数です。 お楽しみに！ /*========================================================*/ 最後の決り文句 /*========================================================*/ このメールマガジンは、まぐまぐさんから発行しています。 このメールマガジンを解除したい場合は、まぐまぐさんをご利用 ください。このメルマガのまぐまぐアイディーは最後にあります。 このメールマガジンには広告が挿入されていますか？ このメールマガジンの内容に文面の引用はありませんか？ めーらっくすの場合はめーらっくすの利用方に従ってください。 このメールマガジンの内容の、転用、流用、宣伝、リンク、 パパ父かわい〜 なんて大歓迎です。 {magclick} /*========================================================*/ /*========================================================*/ 発行者 あゆしゃ まぐまぐアイディー 0000020674 まぐまぐバックナンバー http://jazz.tegami.com/backnumber/frame.cgi?id=0000020674 あゆしゃの世界 http://ayusya.hp.infoseek.co.jp/ 登録と解除 http://www.mag2.com/m/0000020674.htm ご意見・ご感想・ご質問メール mailto:ayusya@flamenco.plala.or.jp めーらっくす ＜＜過去ログがタイトル別になっています＞＞ http://www.mailux.com/mm_dsp.php?mm_id=MM3E1AEE285AB4F |