【プロジェクトオイラー】数学の難問をプログラムの力で解く!!Part1(オイラープロジェクトの紹介)
こんにちは。じばるどーねのとぼです。
皆さん、数学は好きですか?
プログラミングとかやったりしますか?
わたし、数学とかプログラムで難しい問題を解く行為が(双方全くつよくないしむしろよわいのだけれど。精進しろ)好きで、競技プログラミングとかもちょこちょこやってるんですよ。
とはいっても今やっているのは日本でとても有名なAtCoderだけなんですが。世界的に有名なCodeforcesやCodeJamなどにも今後挑戦してみたいなと思ってはいます。英語と夜に弱いので、日本時刻で真夜中に行われる英語のコンテストであるCodeforcesなどはきついですが...。
ここでいう「競技プログラミング」ってのは、プログラミングの問題を制限時間の中でいかに早く、多く解けるかを競う、比較的短時間のコンテストです。AtCoderなどの競技プログラミングでは、各問題にはいくつかのデータセットがあって、各データセットでは入力データが与えられるのでそれに対応した正しい出力データを出したときのみ正解(AC)となります(部分点がある問題もあります)。また、各問題に対して正しい出力を出す場合でも、各問題に与えられた実行時間の制限時間(1秒とか3秒とか)を超えてしまうとTLE(時間制限エラー)となってしまいます。
一方で、1つの問題に対して1つの解答があり、かつ各問題に実行時間制限がないようなプログラミングの問題を提供するサイトも存在します。その中の1つに、「Project Euler」(プロジェクト オイラー) というものがあります。かの偉大な数学者オイラーの名を冠している通り、整数問題や幾何問題などの数学の難問をプログラムの力を使って解こう!というものです。
Project Eulerに挑戦してみよう!!
TLで微妙にこのプロジェクトオイラーが流行っていて、数学とプログラミングが好きな私にはうってつけの内容だったので早速やってみることにしました!!
では早速、Google先生に[Project Euler]で検索をかけてみましょう。
...。
_人人人人人人人_
> 全 文 英 語 <
 ̄Y^Y^Y^Y^Y^Y^Y^ ̄
...はい。このように、少しスタートの威圧感が強いですが、まずは登録をしましょう!
右上の「Register」とかかれたボタンを押して、Username(ニックネーム)やPasswordを決めて登録をしましょう!
登録が終わったら「Archives」(アーカイブ)と書かれたボタンをクリックすれば早速問題一覧を見ることができます!!
...こんな感じで問題一覧が表示されています。現在なんと問題数は600を超えています!気長に楽しみましょう。
ちなみに各問題の「Difficulty」ってとこにあるオレンジ色のメーターが、その名の通り各問題の難易度の指標になっています。この難易度メーターは5%~100%まであり、メーターが低いほど簡単なものになっています(私は現状40%以上の問題を解けていません...)。見ての通り最初の方の問題は簡単なものが多いので、最初は順番に解いていきましょう。
まあ実際に、最初の問題を見てみましょう。
和訳すると、「1000未満の数で3または5の倍数である数をすべて足すといくつになるか求めなさい」ってことですね。
(英語でbelow~やunder~は全て~未満となることに注意してください。以下であれば、"5 or under(5以下)"のように"OO or under/below"のような形で表されます。)
数学の問題として取り組んだらちょっと計算がめんどくさいですが...、プログラムの問題としてはよくある「Fizz Buzz問題」と同程度で、簡単に解けるものですね。
問題が解けたら、下の回答フォームからAnswerを入力しましょう。
(Confirmation Codeを入力するのを忘れずに!たまに1と7とか見分け付きづらいときがありますが、そういう時は画像をクリックすると新しい画像に更新されます!)
入力した答えが正解であれば、
不正解であれば、
といったように、全身白タイツの男が回答の正誤に一喜一憂してくれます。かわいらしいですね。
「英語アレルギーで、問題が解けない!!」
という方も、もしかするといらっしゃるかもしれません。でも大丈夫です。プロジェクトオイラーには有志の方々によって作られた日本語版問題Wikiがあるのです!
それがこちら↓↓
ただし、後半の問題にはまだ翻訳されていないものもあるので、英語ガチ勢の方は是非翻訳完成目指してください(投げ)。
なお、私は、「むしろ原文で問題解いた方が英語の勉強にもなる!!」と思いながら、2問に1問くらいのペースで和訳をガン見しています(。
だって細かい条件が英語だと理解しにくいんだもん...(甘え)。
ちなみに、問題に正解すると自分が何番目に正解したかわかるだけでなく、その問題に関する議論を他の正解者とすることができるようになります(議論スレッドに入れるようになる)。このスレッドでは、最初の100個のコメントと最新の100個のコメントが表示されます。自分のソースコードと他の人のソースコードを比較したり評価しあったりするのも成長につながると思います。(なお、議論は英語で())
ちなみに、私のProblem1に用いたソースコード(C++)のメイン関数はこんな感じです↓↓
int sum = 0;
for (int i = 1; i < 1000; i++)
if (i % 3 == 0 || i % 5 == 0) sum += i;
cout << sum << endl;
return 0;
}
うん。おいしい。
まあ、何の変哲もない感じですね。Problem1なので、オイラープロジェクトの数ある問題の中でも一番簡単なものになっています。
もう少し難易度の高い問題(Difficultyはやはり5%だが)を見てみましょう。
Problem14。これは、有名な数学の未解決問題の1つ「コラッツ予想」に関する実験をする問題です。数学好きとしては題材がコラッツ予想ってだけでワクワクしますね。
コラッツ予想は、任意の正の整数は
・偶数なら2で割る
・奇数なら3倍して1を足す
という操作を繰り返すと最終的に1になる、というものです。
「ほんとかなぁ」(ゴロリ風に)
と思うわけですが、今までのところ反例となる整数は見つかっていないため正しいのではないかと考えられている命題です。
で、今回は正の整数を100万未満に限定したとき、1に至るもまでの操作の回数(ステップ)が一番長くなるのはどの整数だろうか、という問題になります。
いやー、楽しいですね。
興味のある方は是非自分で実装して解いてみてください。
数行先でネタバレします(ソースコード載せます)
こんな感じで実装しました。
//EP14.コラッツ問題
int main() {
int num[1000001] = { 0 };
num[1] = 0, num[2] = 1;
int id = 1;
for (int i = 3; i < 1000000; ++i) {
if (i % 2 == 0) num[i] = 1 + num[i / 2];
else {
ll t = i;
int s = 0;
while (t >= i) {
if (t % 2) t = 3 * t + 1;
else t /= 2;
s++;
}
num[i] = s + num[t];
}
if (ans < num[i]) {
ans = num[i];
id = i;
}
}
cout << id << endl;
return 0;
}
さっきよりは長いコードになりました(それはそう)。
配列numがその名の通り各整数を表しています。さすがに1~100万の整数すべてに対してはじめから操作をシミュレーションしていったら時間がかかりすぎてしまうので、偶奇に場合分けして漸化式に落とし込んでいます。
num[1]=0
iを2以上の整数として
num[i]= 1 + num[i/2] (iが偶数の時)
s + num[t] (iが奇数の時)
(ここでsはtがi未満の値になるまで2で割ったり3倍+1したりする操作を繰り返した回数、tはその結果最後の操作で遷移したi未満の整数)
と漸化式を立てました。奇数状態が続いたら必ずしも高速に動くとは保証できない気もしますが...
まあ、細かいことは気にしない
おしまい。
Part2へ続く!!