|
|
/*========================================================*/ ＜＜＜あゆしゃのＣ言語プログラミング＞＞＞ /*========================================================*/ 第５７９回 加算表作成処理３ 発行 ２００５年６月２８日（月曜日） 発行数 約２６００ {magclick} /*========================================================*/ はじめに （ 決り文句 ） /*========================================================*/ ・このメールマガジンは、主にまぐまぐさんから発行しています。 ・ジャンルは、マルチメディアのプログラム、Ｃ言語です。 ・横６０文字で作成し、インデントは大抵半角スペース４つです。 ・ここで扱うプログラムは、C言語と半光年以内のものです。 ・登録解除は、メルマガのホームページでお願いします。 ・過去ログはバックナンバー（下欄参照）を活用して下さい。 ・内容は私が感じたもので、最新の技術も、へたれもあります。 ・わかりやすさを優先させる為、たまに嘘があるかもしれません。 ・セキュリティ突破のため、暗号化された単語があります。 /*========================================================*/ ご挨拶 /*========================================================*/ こんにちは。あゆしゃです。 先日、大型計算機をバージョンアップしました。 宣伝文句としては、「計算精度の向上」をうたっていますが、 じつは、 ・・・ ・・・・・・大嘘です。 確かに、「結果として」精度は上がりましたが・・・ それはともかく、今回は相互紹介があります。 ●「とおってもわかりやすい経済用語」8910部 ────────────────────────────────────── ●「とおってもわかりやすい経済用語」 「なあんだこんなことだったのか！！ 」新聞・ニュースが良く分かります。 よく読むが、よく聞くが何となく意味の分からない言葉を分かりやすく説明します。 今まで新聞はテレビ欄しか見なかったというあなたも経済欄がこんなに簡単に！！ 【登録】→ http://www.mag2.com/m/0000113177.htm ────────────────────────────────────── 以前紹介したものと同じです。 売掛債権（うりかけさいけん）が読めないというレベルの あゆしゃにとって、まだまだ敷居の高いメルマガではありますが、 最近、攻殻S2がアニメ放送されていることもあって、 そういった技術（？）分野に少し興味を抱いている今日この頃 です。 {magclick} /*========================================================*/ 今回のお題 ＜＜ 加算表作成処理３ ＞＞ /*========================================================*/ さて現在は少し話題を変えて、素因数分解のプログラムについて お話します。 ６００桁の数字の素因数分解に成功すると、２０００万円の 懸賞金がもらえます。 http://www.rsasecurity.com/rsalabs/node.asp?id=2093 素因数分解とは、そんなに難しいのでしょうか？ /*========================================================*/ 前回のプログラムで、７まで指定した場合の結果です。 const int c357[] = { 2,2,4,2,4,2,4,6,2,6,4,2,4,6,6,2,6,4,2,6,4,6,8,4,2,4,2,4,8,6, 4,6,2,4,6,2,6,6,4,2,4,6,2,6,4,2,4,2,10,2,10,2,4,2,4,6,2,6,4, 2,4,6,6,2,6,4,2,6,4,6,8,4,2, 0}; ３と５と７の最小公倍数の３倍は、３１５です。 それ以後の加算表の内容を○×素数表で数えてみると、 ４，２，４，８，６，４，６，２，４，６，２、 と続いています。 この数列が登場するのは１列目の右から５つ目以降です。 つまり、 const int c357[] = { 2,2,4,2,4,2,4,6,2,6,4,2,4,6,6,2,6,4,2,6,4,6,8,4,2, // この部分がループしている 4,2,4,8,6,4,6,2,4,6,2,6,6,4,2,4,6,2,6,4,2,4,2,10,2,10, 2,4,2,4,6,2,6,4,2,4,6,6,2,6,4,2,6,4,6,8,4,2, 0}; ということになります。 /*========================================================*/ このループの先頭の位置は、いったい何なのでしょうか？ 数えてみると、ループしている先頭の４の数字によって求め られる素数は、１０３から１０７になるところ、ということが わかりました。 つまり、１０５という、３と５と７の最小公倍数の位置です。 ３と５の場合はどうでしょうか？ const int c35[] = { 2,2,4,2,4,2,4,6,2,6,4,2, 0}; これは、「４，２，４，６」以降がループしている、ということ でした。 「４」がきたときの素数を数えて調べると、１３から１７へ 変わるときです。 つまり、３と５の最小公倍数の１５をまたいでいるところです。 /*========================================================*/ ★法則性がつかめました！ この「ループが始まる場所」を示す定数も出力するように しましょう。 int last = 3; a = 5; int cnt = 0; int col = 0; do { if( data[ a ] == 0 ) { fprintf( fp, "%d,", a - last ); col += a - last < 10 ? 2 : 3; if( col >= 58 ) { fprintf( fp, EN ); col = 0; } last = a; if( last < minimax ) cnt++; } a += 2; } while( a < minimax * 3 ); // 最後のデータを出力 fprintf( fp, EN "0};" EN ); // ループの先頭位置を出力 fprintf( fp, "const int* %s_head = %s + %d;" EN, name, name, cnt ); cntに、ループの先頭位置のインデックスが入ります。 長くなった場合の自動改行もいれました。 これによる出力結果です。 const int c35[] = { 2,2,4,2,4,2,4,6,2,6,4,2, 0}; const int* c35_head = c35 + 4; const int c357[] = { 2,2,4,2,4,2,4,6,2,6,4,2,4,6,6,2,6,4,2,6,4,6,8,4,2,4,2,4,8, 6,4,6,2,4,6,2,6,6,4,2,4,6,2,6,4,2,4,2,10,2,10,2,4,2,4,6,2, 6,4,2,4,6,6,2,6,4,2,6,4,6,8,4,2, 0}; const int* c357_head = c357 + 25; うん、場所はあっていますね。 コンパイルも通ります。 {magclick} /*========================================================*/ さいごに /*========================================================*/ ちなみに頭文字のｃに、特に意味はありません。 {magclick} /*========================================================*/ 次回予告 /*========================================================*/ 次回は６月３０日（水曜日）に、第５８０回をお送りします。 お題は「大型計算機 1.0.0.7」 大型計算機で変更した内容について紹介します。 お楽しみに！ /*========================================================*/ 最後の決り文句 /*========================================================*/ このメールマガジンは、まぐまぐさんから発行しています。 このメールマガジンを解除したい場合は、まぐまぐさんをご利用 ください。このメルマガのまぐまぐアイディーは最後にあります。 このメールマガジンには広告が挿入されていますか？ このメールマガジンの内容に文面の引用はありませんか？ めーらっくすの場合はめーらっくすの利用方に従ってください。 このメールマガジンの内容の、転用、流用、宣伝、リンク、 限界マックス！ なんて大歓迎です。 {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 |