/*
 * @(#) smt1_1.c - Random generator quality Chi^2 tester
 * (Task 1_1 at the Special Military Training course).
 * (c) 1996 Ivan Maidanski <ivmai@chat.ru> http://ivmai.chat.ru
 * Freeware program source. All rights reserved.
 **
 * Language: ANSI C
 * Tested with: Borland C++ v3.1
 * Last modified: 1996-12-02 22:25:00 GMT+03:00
 */

#include <stdio.h>

#define K 12

int m[1<<K];

int rnd_func(int x)
{
 return (x*5%(1<<K));
}

int main()
{
 int x,n,q,i;
 float chi,chi_alfa;
 printf("Input Chi^2(alfa,%d):\n",K-1);
 scanf("%f",&chi_alfa);
 printf("Input x0 [%d..%d],odd:\n",(1<<(K-2))-1,(1<<K)-1);
 scanf("%d",&x);
 printf("Input N (amount of generated numbers):\n");
 scanf("%d",&q);
 printf("Input n (amount of intervals):\n");
 scanf("%d",&n);
 for (i=0;i<n;i++)
  m[i]=0;
 for (i=0;i<q;i++)
  m[(long)(x=rnd_func(x))*n/(1<<K)]++;
 chi=0.0;
 for (i=0;i<n;i++)
  printf("%5d:%5d\n",i+1,m[i]);
 for (i=0;i<n;i++)
  chi+=(m[i]-(float)q/n)*(m[i]-(float)q/n)*n/q;
 printf("Chi^2: %f\n",chi);
 printf(chi<=chi_alfa? "Accepted\n" : "NOT ACCEPTED\n");
 return 0;
}