1999 年度「計算機基礎論３Ｂ」 1999-11-05

今回もたくさんありますが、 順番に取り組んでください。 でも課題１の締め切りは 1999 年 11 月 11 日（木） 16 時 30 分ですからお忘れなく。

HTML についてもう少し

課題１に関連して、 「ここが一行あいてしまうんですが」 「ここにもう少し空白をいれたいんですが」 といった質問をいただきました。

• <br> とすればそこで強制改行する。
• <pre> と </pre> ではさまれた部分は 「整形済みテキスト」ということで打った通り表示される。 ただし、「<」「>」「&」は“解釈”されるので完全に文字どおりにはならないが。
• 半角スペース「 」は「空白」扱いなので何個打ってもすきまは増えないが、 全角スペース「　」は「白い文字」扱いなので打てば打っただけ空白ができる。

ところで、Netscape のほかには lynx が 実習用 WS にインストールされているので、 興味のある人は試してみてください。 これはもっとも“地味”なブラウザのひとつだと思われます。

lynx の使い方

lynx は、上でも少し説明しましたが、 テキストベースの地味なブラウザです。 この WS では使いにくいようですが、 いちおう使い方を説明しておきます。 時間のある人だけ試してください。 パソコン通信から大学のマシンにログインしたときなどには使える唯一のブラウザとなるでしょう。

※ ……って言ったってもうパソコン通信なんてはやらないよねぇ。

シェルツールまたはコマンドツール （スクロールをオフにしておくほうがいいみたい）の中に 「lynx ~iwase/public_html/99comp/index.html」 のように、 見たいファイル名などを引数につけて打ちこんで起動します。 使いにくいといったのはリンクの張ってあるところ （Netscape なら色を変えて下線付きで表示されるところ） といまカーソルがあたっているところ （Netscape ならマウスカーソルが当っているところ） が同じく白黒反転文字で表示されるからです。 よって、 「とにかく起動してからリンクをたどろう」 と思わないで、 いきなり見たいファイルを引数にとって起動するのがよいでしょう。 自分の書いたホームページを lynx で見ると、 ちょっと違って見えておもしろい（？）かも知れません。 スペースバーで次の画面、 「↑」で前の画面、 Ctrl+R で再読み込みです。 終了は「Q」。

Netscape の画面の印刷

「ファイル」「印刷(P)...」と選び、 「印刷コマンド」の中に「lpr -PNPR215」とします。 用紙サイズは「A4」だと思います。

はじめに

今回は、Ｃ言語の続きとして、くり返しと配列について学びます。 きょうのプログラムもサブディレクトリ ~iwase/prog に置いてありますので、 必要なかたはコピーして利用してください。

くり返しを含むＣプログラム

いままでのプログラムは関数電卓でも代用できるものでした。 それに対し「これこれを 1000 回くり返せ」のようにくり返しを含む計算は、 コンピュータを使うほうがずっと楽です。

Ｃ言語でくり返しを実行させる方法には for と while と do-while の三つがあります。 ここでは for と while だけを説明し、do-while はとりあげません。 まずは抽象的に説明しますが、 そのあとで実例が出てきますから、 合わせて理解してください。

    for (文１; 式２; 文３) {
        ...
    }

• まず文１（初期設定）を実行。
• 式２（継続条件）が真ならば中カッコの中を実行する； 文３（再初期化）を実行し式２のチェックにもどる。
• 式２が偽ならば，中カッコはとばして次へ進む。

もう一つの構文

    while (条件) {
        ...
    }

• 条件が真なら中カッコの中を実行し、また条件のチェックへもどる。
• 条件が偽なら中カッコはとばして次へ進む。

== sum.c ======================================================================
#include <stdio.h>  /* ユーザが数を入力すると，1 からその数までの和を出力 */

main() {
    int n, sum, limit;

    printf("いくつまで足しますか？\n");
    scanf("%d", &limit);

    sum = 0;

    for (n = 1; n <= limit; n++) {
        sum = sum + n;
    }
    printf("1 から %d まで足すと %d です.\n", limit, sum);
}
===============================================================================

ここまでで自作プログラムを試してみたい人は、 1 から始めて 1, 2, 4, 8, 16, ... と倍々になってゆく数列を表示するプログラムを書いてみましょう。 2 の何乗で答えがおかしくなるでしょうか？

== arale.c ====================================================================
#include <stdio.h>  /* アラレ数 */

main() {
    int n;

    printf("自然数を入れてください.\n");
    scanf("%d", &n);
    while (n != 1) {
        if (n % 2 == 0) {
            n = n / 2;          /* n が偶数のとき */
        } else {
            n = 3 * n + 1;      /* n が奇数のとき */
        }
        printf("%d ", n);
    }
}
===============================================================================

• ある項が偶数のときはその項を 2 で割ったものを次の項とする
• ある項が奇数のときはその項の 3 倍に 1 を加えたものを次の項とする

このプログラムは、ユーザが入力した数から始めて 1 になるまで、 この数列を出力します。

== exp.c =====================================================================
#include <stdio.h>  /* 自然対数の底 e の計算 */

main() {
    int n;
    double sum, kaizyo;

    sum = 1;
    kaizyo = 1;
    n = 1;

    for (n = 1; n < 24; n++) {
        kaizyo = kaizyo * n;
        sum = sum + 1 / kaizyo;
        printf("第 %d 項までの和は %.25f です.\n", n, sum);
    }
}
===============================================================================

自然対数の底 e を e = 1/0! + 1/1! + 1/2! + 1/3! + ... として計算しています。 Ｃ言語には階乗を表わす演算子はないので自力で計算しなければなりません。 ここでは n! を計算した結果を捨てないで次の (n+1)! の計算に使っています。

なお、「kaizyo = kaizyo * n」のところでは double 型と int 型を掛けています。 こういった場合、結果は double 型になります。

== tansindo.c =================================================================
#include <stdio.h>  /* 単振動 */

#define T 0.01          /* 時刻のきぎみ */

main() {
    double t, x, v;     /* 時刻，変位，速度 */

    x = 1; v = 0;       /* 初期条件 */

    for (t = 0; t < 10; t = t + T) {
        printf("時刻 %5.2f, 変位 %9.6f, 速度 %9.6f\n", t, x, v);
        x = x + v * T;  /* T だけ後の位置を計算 */
        v = v - x * T;  /* T だけ後の速度を計算 */
    }
}
===============================================================================

近似計算で微分方程式を解くプログラムの例です。 ここでは x を変位、v を速度、a を加速度とするとき微分方程式 「a = -x」を初期値 x = 1, v = 0 のもとで解いてみました。

時刻は T = 0.01 きぎみにし、 その間は等速運動をするとして T だけ後の位置を計算しているのが最後から二つめの文、 その間は等加速運動をするとして T だけ後の速度を計算しているのが最後の文です。

「#define T 0.01」というのは 「以後 T が出てきたら 0.01 で置きかえろ」ということでした。 「#define」を使わなくても T の代わりに 0.01 と書けば同じことですが、 時刻のきざみを変えてこのプログラムを試してみたいとき 「#define」を使っていれば「#define T 0.01」の「0.01」 だけを書きかえればよいのに対し、 「#define」を使っていなければ三箇所に出てくる T の値をすべて書きかえなければなりません。

このプログラムは非常にいい加減な近似をしています。 もっと精度のよいやり方は別の講義などにゆずります。

理論的にこの微分方程式を解けば x = cos(t) になります。 それとの差を並べて出力してみるのも面白いかも知れません。

※ ここで「出力リダイレクト」の説明をしておこう。 「a.out」の代わりに「a.out >> filename」とすると、 filename という名前のファイルの内容の最後に a.out の出力結果を追加する。 （ファイルがなければ新規に作られる。） これを利用すると、 長い出力を行なうプログラムの結果をあとからエディタでゆっくり見ることができる。 a.out でなく、ほかのコマンドでも同様。

「かつ」や「または」は次のように書きます。 これらは for や while の条件部にも使えます。

    if ((x > 0) && (x < 8)) {       /* 「かつ」 */
        ........
    }

    if ((x < 0) || (x > 2)) {       /* 「または」 */
        ........
    }

配列

「配列」というのは，変数の列です。

== sieve.c ====================================================================
#include <stdio.h>  /* エラトステネスのふるいで 10000 以下の素数を全て求める */

#define UNERASED 0              /* ←英語として正しいかどうか */
#define ERASED 1                /*   責任はもてません :-)     */

int a[10000+1];                         /* 配列 a[0], ... , a[10000] */

main() {
    int n, i;

    for (n = 2; n <= 10000; n++) {      /* 配列の初期化 */
        a[n] = UNERASED;
    }

    for (n = 2; n <= 100; n++) {        /* ふるいにかける */
        for (i = 2; n * i <= 10000; i++) {
            a[n*i] = ERASED;
        }
    }

    for (n = 2; n <= 10000; n++) {      /* 結果の出力 */
        if (a[n] == UNERASED) {
            printf("%d は素数.\n", n);
        } else {
            printf("%d は合成数.\n", n);
        }
    }
}
===============================================================================

エラトステネスのふるい （Eratosthenes's sieve）を使って 10000 以下の素数を全て求めるプログラムです。

数学科のみなさんには説明するまでもないでしょうが、 紙と鉛筆でやる場合はこうやります。 1 から 10000 までの自然数を書いておき、 まず 2 の倍数を斜線で消します。 次に 3 の倍数を消します。 次に 4 の倍数ですが、 これはすでに 2 の倍数として消されているので消さなくてよい……。 こうやっていって残ったものが素数でした。

「int a[10000+1];」は 10001 個の int 型変数の配列 a[0], ... , a[10000] を宣言する構文です。 Ｃ言語の配列は必ず添字 0 から始まるので、 自然数 n に a[n] を対応させるとしたら a[10000] まで用意するためにはサイズ 10001 の配列が必要です。 「int a[10001];」と書いてもいいのですが、 「a[0] があるから +1 が必要だよ」 ということを強調するためわざわざこう書きました。

そして、「斜線で消す」という代わりに、 a[5] に 0 が代入されていたら 5 は消されていない、 1 が代入されていたら消されている、としました。

「#define」で始まる二行ですが、 0 と 1 のどちらが消されている印でどちらが消されていない印だったかわからなくなるといけないので、 こうやって名前をつけました。 「#define」にはこういう用法もあります。 この UNERASED と ERASED は変数名ではなく定数名なので大文字で書いています。

三つの for 文が何をやっているかは、 簡単にコメントしておいたので読めばわかると思います。 二つ目の for 文はその中にさらに for 文を含んでいます。 なお、そこでは合成数の倍数も消していますので、 このプログラムは最も効率のよいものではありません。

ここまでに習ったことを組み合わせると、いろいろなプログラムが書けます。

いつだったかは、 課題として「第一項、第二項が 1 である Fibonacci 数列の最初の十数項について、 その値、およびその項と前の項との比の値を出力するプログラムを書け」 というのを出しました。 みなさんもやってみてはいかがでしょう？