/*
 *
 * Copyright (c) 1994
 * Hewlett-Packard Company
 *
 * Permission to use, copy, modify, distribute and sell this software
 * and its documentation for any purpose is hereby granted without fee,
 * provided that the above copyright notice appear in all copies and
 * that both that copyright notice and this permission notice appear
 * in supporting documentation.  Hewlett-Packard Company makes no
 * representations about the suitability of this software for any
 * purpose.  It is provided "as is" without express or implied warranty.
 *
 */

/*
 * void randgen.seed(unsigned long j)  reseeds the generator with j
 *    
 * unsigned long ranggen(unsigned long limit)   
 *     returns a pseudo-random number n, 0 <= n < limit,
 *     uniformly distributed
 *
 * double randgen()
 *     returns a pseudo-random number d, 0.0 <= n < 1.0,
 *     uniformly distributed
 */

#include <stddef.h>
#include <climits>  // for ULONG_MAX

#define __SEED 161803398

class random_generator {
protected:
    unsigned long table[55];
    size_t index1;
    size_t index2;
public:
    unsigned long operator()(unsigned long limit) {
	index1 = (index1 + 1) % 55;
	index2 = (index2 + 1) % 55;
	table[index1] = table[index1] - table[index2];
	return table[index1] % limit;
    }
    inline double operator()() {
        return (double)operator()(ULONG_MAX) / ULONG_MAX;
    }
    void seed(unsigned long j);
    random_generator(unsigned long j) { seed(j); }
};

void random_generator::seed(unsigned long j) {
    unsigned long k = 1;
    table[54] = j;
    size_t i;
    for (i = 0; i < 54; i++) {
 	size_t ii = (21 * (i + 1) % 55) - 1;
        table[ii] = k;
        k = j - k;
        j = table[ii];
    }
    for (int loop = 0; loop < 4; loop++) {
        for (i = 0; i < 55; i++)
            table[i] = table[i] - table[(1 + i + 30) % 55];
    }
    index1 = 0;
    index2 = 31;
}

static random_generator randgen(__SEED);