Jeans & Development コンピューターのことなどを綴ったメモ (旧：目から鱗 w/SQLite)

ちょっと訳有って、円周率（π、パイ、pi）を、コンピューターに計算させてみた。とりあえず、例によって、PHPで。

円周率の計算は、手作りコンピューターでやろうとしていることなので、アルゴリズムとしては最も簡単なものを利用した。最初に作成したスクリプトは、以下のとおり。

for ($max=1;$max<=8192;$max=$max*2) {
    $t=time();
    echo "$max: ".calc_pi($max)."  (".(time()-$t).")\n";
}

function calc_pi($max) {
    $total=0;
    for ($x=0;$x<$max;$x++) {
        for ($y=0;$y<$max;$y++) {
            if ($x*$x + $y*$y < $max * $max) $total++;
        }
    }
    $pi=4 * ( $total/($max*$max) );
    return $pi;
}

正方形の内部にまんべんなく点を打ったとき、その点から一つの頂点までの距離が正方形の辺の長さより短いかどうかをピタゴラスの定理を用いて判断することで、円の内部かどうかを調べ、円の面積から円周率を計算している。

実行結果は、以下のようになった。

1: 4  (0)
2: 4  (0)
4: 3.75  (0)
8: 3.5  (0)
16: 3.34375  (0)
32: 3.25390625  (0)
64: 3.19921875  (0)
128: 3.170166015625  (0)
256: 3.15673828125  (0)
512: 3.1490936279297  (0)
1024: 3.1453971862793  (1)
2048: 3.1435289382935  (2)
4096: 3.1425535678864  (10)
8192: 3.1420764327049  (40)

<円周率の計算 のつづき>

Windows NTのcmd.exeでは、ファイルをcmdのウインドウにdrag&dropすることで、ファイルのフルパスを入力することができた。これは、XP以降ではできなくなったようで（もしかしたら、Vista以降かな？）、不便を感じたので、SendToからファイルパスをクリップボードにコピーするスクリプトを書いてみた。

<VBScriptで、クリップボードに文字列をコピーする のつづき>

本来は日記ブログに書くべき内容なのだが、少し思うところあってこちらに書いてみる。

先日、NHKでニュースを見ていたときのこと。高齢者ドライバーの話になっていた。Webでの当該ニュースへのリンクを見つけられないのだが、大元はたぶんこれ。実際にニュースとして聞いたものを要約すると次のような感じであった。

愛媛県では高齢者率が高く、２５％の人が高齢者である。高齢者ドライバーによる交通事故も年々増加し、去年は事故全体の１６％が高齢者ドライバーによるものであった。増加する高齢者ドライバーの事故に関して、愛媛県でどのような対応をしているのか、追った。（以下略）

さっと聞き流すともっともな内容であるが、良く考えるとこれは少し変な話である。全体で２５％の方が高齢者であるのに対し、事故を起こした高齢者ドライバーの率が１６％であるというのは、高齢者ドライバーが優秀であることを示しているのではないか？

愛媛県の人口１４万人のうち、３万五千人が高齢者、高齢者以外は１０万五千人である。仮に交通事故の総件数が１００だと仮定すると、高齢者ドライバーによる事故が１６件、高齢者でないドライバーによる事故が８４件。これらの値から、高齢者ではおよそ２２００人に一人が事故を起こし、高齢者以外ではおよそ１３００人に一人が事故を起こしていることになる。このデータを見る限り、高齢者ドライバーは優秀ドライバーなのだ。

<変な話 のつづき>