//  粒子の凝集シミュレーション（アスキーアート）
//  粒子は領域の上からランダムにやって来て下、左、右に移動する
//  粒子は領域の底に到達したら吸着する
//  粒子は他の粒子に接すると吸着する
#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <time.h>

#define  PMAX  300                                      // 粒子の数
#define  XMAX  80                                       // 水平方向MAX
#define  YMAX  20                                       // 垂直方向MAX

void  main( )
{
    unsigned long  seed;                                // 乱数の種
    char  buf[YMAX][XMAX];                              // 1画面分のマス
    char  pc='*';                                       // 粒子の文字
    int   n,x,y;                                        // 粒子位置
    int   ir;                                           // 吸着
    int   i;

    printf("粒子の凝集シミュレーション（ＡＡ）\n");

    memset(&buf[0][0],' ',XMAX*YMAX);                   // 配列を空白で初期化
    time(&seed);                                        // 時刻取得(秒)
    srand(seed);                                        // 乱数の種をまく

    for( i=0; i<PMAX; i++ )                             // 粒子発生
      {
        x=rand( )%XMAX; y=0;                            // 最初の位置

        if ( buf[y][x]==pc ) continue;                  // すでに粒子あり

        while( 1 )
          {
            ir=1;                                       // 吸着(仮)

            if ( y==YMAX-1 ) break;                     // 領域の底

            if ( buf[y+1][x]==pc ) break;               // 下に粒子あり

            if ( x<XMAX-1 )
              { if ( buf[y+1][x+1]==pc ) break; }       // 右下に粒子あり

            if ( x<XMAX-1 )
              { if ( buf[y][x+1]==pc ) break; }         // 右に粒子あり

            if ( x<XMAX-1 && y>0 )
              { if ( buf[y-1][x+1]==pc ) break; }       // 右上に粒子あり

            if ( y>0 )
              { if ( buf[y-1][x]==pc ) break; }         // 上に粒子あり

            if ( x>0 && y>0 )
              { if ( buf[y-1][x-1]==pc ) break; }       // 左上に粒子あり

            if ( x>0 )
              { if ( buf[y][x-1]==pc ) break; }         // 左に粒子あり

            if ( x>0 )                                  // 左下に粒子あり
              { if ( buf[y+1][x-1]==pc ) break; }

            ir=0;                                       // 吸着なし

            //--移動--
            n=rand( )%3;                                // 3方向

            switch( n )
              {
                case 0: y++; break;                     // 下
                case 1: x++; break;                     // 右
                case 2: x--; break;                     // 左
              }

            if ( y<0 || y>=YMAX ) break;                // 領域の外
            if ( x<0 || x>=XMAX ) break;                // 領域の外
          }

        if ( ir==1 ) buf[y][x]=pc;                      // 粒子吸着
      }

    //--結果表示--
    for( y=0; y<YMAX; y++ )
      { printf("%.*s\n",XMAX,&buf[y][0]); }             // 1行分表示
}