weak pseuodrandom number generator

A pseudorandom number generator recently proposed here, namely: int rand1(int seedval) { return (seed * 183041 % 183319 + 1); } needs some cleaning up. It should be something like: unsigned long rand1(unsigned long n) { return ( ( ( n * 183041L) % 183319 ) + 1 ); } where n is initally set to some seed value. However, this is particularly weak, and quickly degenerates into a cycle, usually of length 208, as the following program will confirm: #include #include unsigned long a[15000]; unsigned long rand1( unsigned long n ); void main(int argc, char **argv) { unsigned int i=0, j; if ( argc < 2 ) return; a[0] = atol(argv[1]); while ( ++i < 15000 ) { a[i] = rand1(a[i-1]); for ( j=0; j