|
|
/*========================================================*/ ＜＜＜あゆしゃのＣ言語プログラミング＞＞＞ /*========================================================*/ 第４６１回 平方根の少数 発行 ２００４年６月１４日（月曜日） 発行数 約２９００ {magclick} /*========================================================*/ はじめに （ 決り文句 ） /*========================================================*/ ・このメールマガジンは、主にまぐまぐさんから発行しています。 ・ジャンルは、マルチメディアのプログラム、Ｃ言語です。 ・横６０文字で作成し、インデントは大抵半角スペース４つです。 ・ここで扱うプログラムは、C言語と半光年以内のものです。 ・登録解除は、メルマガのホームページでお願いします。 ・過去ログはバックナンバー（下欄参照）を活用して下さい。 ・内容は私が感じたもので、最新の技術も、へたれもあります。 ・わかりやすさを優先させる為、たまに嘘があるかもしれません。 ・セキュリティ突破のため、暗号化された単語があります。 /*========================================================*/ ご挨拶 /*========================================================*/ こんにちは。あゆしゃです。 先日、株の売買ゲームに登録しました。 野村のバーチャル株式投資倶楽部 http://www2.nomura.co.jp/vstock/VirtualServlet? 仮想のお金（初期値：１００万円）と、実際の株式の状態を元に 資金を増やすゲームです。 参加費は無料で、扱われる株式の値段は本物ですので、株の お勉強に最適です。 ★Ayusyaの成績 木曜日 999,865円です。（49,945番／90,783人中） 金曜日 1,002,444円です。（21,299番／91,327人中） 土曜日 1,001,669円です。（21,913番／91,794人中） 毎日、株価を見るのが楽しみになりました。 結構、簡単に儲かるものですね、金額は１００円単位で少ないの ですが。 2386サイトサポート・インスティテュートには、偶然にもぼろ儲 けさせていただきました。 １日で＋４００円。２株買って売ったので、 ＋８００円−７６円（手数料）＝７３４円の儲けでした。 本番では１００株単位らしいので、３８万円の元手で１日に ★４万円の儲け。。ぐは。。！ {magclick} /*========================================================*/ 今回のお題 ＜＜ 平方根の少数 ＞＞ /*========================================================*/ 平方根シリーズの冒頭にて、 「平方根とは、桁数を叩き割った数字です」 という説明をしました。 覚えていますか？ ３ヶ月前ですが。 ★私は忘れていました。 しかしふと、少数を計算したいなぁと思い、いろいろと考えて いたところ、 この原則が、重要でした。 つまり、少数の桁数も叩き割られる、ということです。 /*========================================================*/ 具体的に見てみましょうか。 ルート12345 は 111.10805 です。 １６進数でいうと、ルート0x3039 は 0x6F です、となります。 で、は。 固定小数点数として、下位８ビットが少数であると考えると、 0x30.39＝48.22265625 であり、 ルート48.22265625 は、6.94425 となります。 しかし、0x0.6F は 0.43359375 です。 求めるべき 6.94425 とは異なります。 ★ここで、小数部の桁数を叩き割るのです。 ８ビットだった小数部を４ビットにします。 0x06.F は、6.9375 です。 求めるべき 6.94425 と似ています！ ビットの桁落ちがあるので、多少誤差が出ますが、平方根が 求められていますよね？ /*========================================================*/ あゆしゃ法に限らず、すべての整数版平方根アルゴリズムは、 ★固定少数点数で小数を計算できるのです。 r = ( UINT )( 123456.789 * 256 ); // 小数部８ビット ret = ays_sqrt( r ); ret <<= 4; // 小数部が４ビットになっているので８に戻す という感じです。 /*========================================================*/ 計算はガウスさん、１６進数変換はWindowsの電卓を使用しまし た。 ガウスさんの答えで小数桁が多いものは、５桁だけ拾いました。 Windowsの電卓で 0x0.6F を１０進数に変換するには、 １．6F を１０進数の 111 にする ２．２５６で割る とすると求められます。 小数部が８ビットなので２５６で割るのです。 {magclick} /*========================================================*/ さいごに /*========================================================*/ 新C言語使いにおくるチョー基本講座 第７回。 数字の先頭に０ｘをつけると、１６進数になります。 だから、０ｘ１０と書くと、１０進数では１６になります。 ちなみに、数字の先頭を０にすると、８進数になります。 だから、０１０と書くと、１０進数では８になります。 ★８進数はバグの元です。使わないほうがいいです。 ★０４０６１６など、１０進数のつもりで０を使うのは危険です！ ８進数はともかく、１６進数はよく使用します。 これぐらいは暗記しておいたほうが便利です。 １００は０ｘ６４ ０ｘ１０は１６ ０ｘ１００は２５６ ０ｘ４００は１０２４ ０ｘ１０００は４０９６ ０ｘ１００００は６５５３６ これぐらいかなぁ。 ２の乗数も、１乗（１）から１６乗（６５５３６）まで覚えて おいたほうがいいです。 私は１乗から１６乗までは覚えていますが、１７乗以上は覚えて いません。 また、覚えているといっても、１３乗は？といきなり聞かれると 指で数えなければ答えられません。 それでも生きていけます。 {magclick} /*========================================================*/ 次回予告 /*========================================================*/ 次回は６月１８日（水曜日）に、第４６２回をお送りします。 お題は「少数の読み込み」 いっそのこと、LC で固定小数点数を扱うようにしましょう。 そのためには少数の文字列の読み込みと、文字列への変換が 必要です。 お楽しみに！ /*========================================================*/ 最後の決り文句 /*========================================================*/ このメールマガジンは、まぐまぐさんから発行しています。 このメールマガジンを解除したい場合は、まぐまぐさんをご利用 ください。このメルマガのまぐまぐアイディーは最後にあります。 このメールマガジンには広告が挿入されていますか？ このメールマガジンの内容に文面の引用はありませんか？ めーらっくすの場合はめーらっくすの利用方に従ってください。 このメールマガジンの内容の、転用、流用、宣伝、リンク、 宝くじが当たりますようにぃ！！ なんて大歓迎です。 {magclick} /*========================================================*/ /*========================================================*/ 発行者 あゆしゃ ホームページ::あゆしゃの世界 http://ayusya.hp.infoseek.co.jp/ ご意見・ご感想・ご質問メール mailto:ayusya@flamenco.plala.or.jp まぐまぐ::アイディー 0000020674 まぐまぐ::登録と解除 http://www.mag2.com/m/0000020674.htm まぐまぐ::バックナンバー http://jazz.tegami.com/backnumber/frame.cgi?id=0000020674 めーらっくす::アイディー MM3E1AEE285AB4F めーらっくす::登録と解除 http://www.mailux.com/mm_dsp.php?mm_id=MM3E1AEE285AB4F めーらっくす::バックナンバー★最近のものならこちらが便利★ http://www.mailux.com/mm_bno_list.php?mm_id=MM3E1AEE285AB4F |