#include <stdio.h>
#include <stdlib.h>

#define N 6

double gosa(int a[], double b[])
{
  int i;
  double res;
  res = 0;
  for (i=0; i<N; i++) { 
    res = res + fabs(((double) a[i])-b[i]);
  }
  return (res);
}

int main(void)
{
    int x,i,j,k, total_life, denomi;
    int occur[N];
    int res_to_test[N]={2, 5, 2, 8, 17, 6};
    double prob_ratio[N] = {1, 2, 1, 3, 6, 3};
    double prob_list[N];
    double prob_dist[N] = {0,0,0,0,0,0};
    double ideal[N];
    int lessn, moren;

    printf("input seed>");
    scanf("%d",&x);
    srand(x);

    total_life = 0;
    for (i=0; i<N; i++) total_life += res_to_test[i];

    denomi = 0;
    for (i=0; i<N; i++) denomi += prob_ratio[i];

    prob_list[0] = prob_ratio[0]/denomi;
    ideal[0] = prob_list[0] * total_life;
    for (i=1; i<N; i++) {
      prob_list[i] = prob_list[i-1] + (prob_ratio[i]/denomi);
      ideal[i] = (prob_ratio[i]/denomi) *  total_life;
    }

    /*
    for (i=0; i<N; i++) {
      printf("prob_list[%d]=%6f\n",i,prob_list[i]);
    }
    */

    lessn = 0;
    moren = 0;
    for (j=0; j<1000; j++) {
      for (k=0; k< N; k++) occur[k]=0;
      
      for (i=0; i<total_life; i++) {
	  double ran;
	  ran = ((double) rand())/RAND_MAX;
	  for (k=0; k<N; k++) {
	    if (ran < prob_list[k]) {
	      occur[k]++;
	      break;
	    }
	  }
      }

      /*
      printf("occur=");
      for (k=0; k<N ; k++) {
	printf("=%3d, ", occur[k]);
      } 
      printf("\n");
      */

      if (gosa(res_to_test, ideal) <= gosa(occur, ideal)) { 
	moren ++; 
      } else {
	lessn ++;
      }
    }
    printf("(lessn,moren)=%8d, %8d\n", lessn, moren) ;
}