/* life.tc - 3/13/2010 - lrb

int dim,pop;char cell(50*50),dell(50*50)

blinker [
 cell((dim/2-1)*dim+dim/2-1)='*'
 cell((dim/2-1)*dim+dim/2)='*'
 cell((dim/2-1)*dim+dim/2+1)='*'
 ]

glider [
 cell((dim/2-1)*dim+dim/2)='*'
 cell(dim/2*dim+dim/2+1)='*'
 cell((dim/2+1)*dim+dim/2-1)='*'
 cell((dim/2+1)*dim+dim/2)='*'
 cell((dim/2+1)*dim+dim/2+1)='*'
 ]

longline int l [
 int i,j;j=(dim-l)/2;for (i=j;i<j+l;++i) cell((dim/2)*dim+i)='*'
 ]

manual [
 int r,c
 pl "";pl "";ps "0 to end";pl "";pl ""
 ps "row (0:49) ";r=gn
 while (r) [ps "column (0:49) ";c=gn;cell(r*dim+c)='*';ps "row (0:49) ";r=gn]
 ]

life [
 char b,bb(1),lifefile(30),rcell(1);int k,p,slast,r,c,gen,l,lpop,periodic
 rcell(0)=' ';rcell(1)='*';dim=50;p=50
 pl "";ps "life.tc - 3/13/2010";pl ""
 pl "";ps "please visit http://primepuzzle.com/tc/life.html";pl ""
 while (b!='q') [
  periodic=pop=0;posc 43,55;ps "gens desired (2:999) ";gen=gn
  ceol;posc 44,1;ps "b/g/l/f/m/a/s/r/w/p/q/? ";bb(0)=b=getchar
  if !((b=='r')+(b=='w')+(b=='p')+(b=='s')) memset cell,dim*dim,' '
  if (b=='b') [blinker;b='p']
  if (b=='g') [glider;b='p']
  if (b=='l') [pl "";ps " line length (3:46) ";l=gn;longline l;b='p']
  if (b=='m') [manual;b='p']
  if (b=='r') [pl "";ps "filename ";gs lifefile;readfile(lifefile,cell,cell+dim*dim-1,1)]
  if (b=='w') [pl "";ps "filename ";gs lifefile;writefile(lifefile,cell,cell+dim*dim-1,1)]
  if (b=='f') [
   for (c=0;c<dim;++c) cell(c)=cell((dim-1)*dim+c)='*'
   for (r=0;r<dim;++r) cell(dim*r)=cell(dim*(r+1)-1)='*'
   b='p'
   ]
  if (b=='a') [
   pl "";ps "Enter a seed number ";slast=last=gn
   if !(last%2) ++last /* make number odd
   seed=last;for (c=0;c<dim*dim;++c) cell(c)=rcell(random (0,1))
   b='p'
   ]
  if (b=='s') [pl "";ps "Enter speed number 1:100 ";p=gn;b='p']
  if (b=='?') [
   pl "";pl "";ps "(b)link (g)lide (l)ine (f)ill (m)anual"
   pl "r(a)ndom (s)peed (r)ead (w)rite (p)lay (q)uit"
   ]
  if (b=='p') [
   off;hilo '1';cls;display;posc 47,55
   ps "gen ";pn 0;ps " population";pn pop
   posc 48,55;ps "speed";pn p;if slast [ps " seed";pn slast]
   pause p;lpop=pop
   for (k=1;k<=gen-1;++k) [
    core;pause p;strcpy cell,dell;display;posc 47,55
    ps "gen";pn k;ps " population";pn pop;ceol
    posc 48,55;ps "speed";pn p;if slast [ps " seed";pn slast]
    if !pop [posc 47,35;color '4';beep 880,5;ps "poof!";break]
    if (pop==lpop) [
     ++periodic;if (periodic==13) [
      color '4';posc 41,55;ps "periodic or static!";break
      ]
     ]
     else [lpop=pop;periodic=1]
    ]
   on;hilo '0';color '7';pl ""
   ]
   if !index("bglfmasrwpq?",12,bb,1) [display;posc 47,55;ps "huh?";ceol]
 ]]

posc int r,c [MC '\'-'A',1;MC '\'-1,1;pos r;MC ';',1;pos c;MC 'H',1]

pos int rc [
 if (rc>9) [MC rc/10+48,1;rc=rc%10]
 MC rc+48,1
 ]

ceol [MC '\'-1-'@',1;MC '\'-1,1;MC 'K',1]

display [
 int i,j;pop=0
 for (i=0;i<dim;++i) [
  for (j=0;j<dim;++j) [
   posc i,j;putchar cell(dim*i+j)
   if (cell(dim*i+j)=='*') ++pop
 ]]]

core [
 int r,c,count
 for (r=0;r<dim;++r) [
  for (c=0;c<dim;++c) [
   count=neighbor(r,c)
   dell(dim*r+c)=cell(dim*r+c) /* default is no change
   if (cell(dim*r+c)=='*') [
    if ((count>3)+(count<2)) dell(dim*r+c)=' '
    ]
   else [
    if (count==3) dell(dim*r+c)='*'
 ]]]]

neighbor int row,col [
 int count
 if !((row-1<0)+(col-1<0)) count=count+(cell(dim*(row-1)+col-1)=='*')
 if !((row-1<0)) count=count+(cell(dim*(row-1)+col)=='*')
 if !((row-1<0)+(col+1>(dim-1))) count=count+(cell(dim*(row-1)+col+1)=='*')
 if !((col-1<0)) count=count+(cell(dim*row+col-1)=='*')
 if !((col+1>(dim-1))) count=count+(cell(dim*row+col+1)=='*')
 if !((row+1>(dim-1))+(col-1<0)) count=count+(cell(dim*(row+1)+col-1)=='*')
 if !((row+1>(dim-1))) count=count+(cell(dim*(row+1)+col)=='*')
 if !((row+1>(dim-1))+(col+1>(dim-1))) count=count+(cell(dim*(row+1)+col+1)=='*')
 return count 
 ]

pause int u [
 int p,q,r
 for (p=0;p<u;++p) [
  for (q=0;q<u;++q) [
   for (r=0;r<u;++r) [;]
 ]]]