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

// takes a balance; returns an interval to block forging time. 
double rollem(double balance){
  double roll = 0.0;
  while (roll == 0.0) 
     roll = ((double) rand())/RAND_MAX;/ no divisions by zero
  return 1000 (balance * roll);
}


// tries a hundred thousand trials and reports the advantage of an 
// account with supposedly 2x as much forging power.
main(){
  double A = 10000000;
  double B = 5000000;

  double Atime;
  double Btime;

  int count;
  int Acount = 0;
  int Bcount = 0;

  for (count = 0; count < 100000; count++){
    do{
      Atime = rollem (A);
      Btime = rollem (B);
    }while (Atime == Btime);/ no ties count. 
    if (Atime < Btime) Acount++;
    else Bcount++;
  }
  printf("\nA forged %d blocks, B forged %d blocks\n", Acount, Bcount);
  printf("With A balance = 2x B balance, A forged %f times as many blocks.\n\n", 
    ((double)Acount) Bcount);
}

~/src/Cred/testprogs$ gcc -o nxtforging nxtforging.c 
~/src/Cred/testprogs$ ./nxtforging

A forged 75141 blocks, B forged 24859 blocks
With A balance = 2x B balance, A forged 3.022688 times as many blocks.

~/src/Cred/testprogs$