Yer another goto 10 labirint generator

rand()%2

Версия функций rand() и srand() в библиотеке С для Linux использует тот же генератор чисел, что и в функциях random() и srandom(), так что младшие биты в числе случайны настолько, насколько и старшие. В то же время в старой реализации rand() младшие биты являются гораздо менее случайными, чем старшие 

В книге Numerical Recipes in C: The Art of Scientific Computing (William H. Press, Brian P. Flannery, Saul A. Teukolsky, William T. Vetterling; New York: Cambridge University Press, 1992 (2nd ed., p. 277)) даны следующие комментарии:

    "Если Вы желаете получить случайное число в промежутке от 1 до 10, Вы всегда должны использовать старшие биты, например:

        j=1+(int) (10.0*rand()/(RAND_MAX+1.0));

    не следует выполнять такое действие:

        j=1+(rand() % 10);

    (т.е. использовать младшие биты)." 

double y = (double) rand();
y = y / (double) pow(10, 1+floor(log10(y)));

double y1 = (double) rand();
double y2 = (double) rand();
int l1 = 1.0+(int) floor(log10(y1));
int l2 = 1.0+(int) floor(log10(y2));
y1 = y1 / (double) pow(10, (double)l1);
y1 += y2 / (double) pow(10, (double)(l1+l2));

yes 'c=(╱ ╲);printf ${c[RANDOM%2]}'|bash

cat 2.c 
#include <stdio.h>		// printf
#include <stdlib.h>		// random number functions
#include <sys/time.h>		// gettimeofday
#include <limits.h>		// LONG_MAX
#include <math.h>		// floor
#include <unistd.h>		// read
#include <sys/types.h>		// open
#include <sys/stat.h>		// open
#include <fcntl.h>		// open

double dtime(){
	double t;
	struct timeval tv;
	gettimeofday(&tv, NULL);
	t = tv.tv_sec + ((double)tv.tv_usec)/1e6;
	return t;
}

void urandom_ini(){
	double tt = dtime() * 1e6;
	double mx = (double)LONG_MAX;
	srand48((long)(tt - mx * floor(tt/mx)));
}

void dev_random_ini(){
	long r_ini;
	int fd = open("/dev/random", O_RDONLY);
	if(-1 == fd){fprintf(stderr,"Error open /dev/random!\n"); urandom_ini(); return;}
	if(sizeof(long) != read(fd, &r_ini, sizeof(long))){
		fprintf(stderr,"Error read /dev/random!\n"); urandom_ini();
		close(fd); return;
	}
	close(fd);
	srand48(r_ini);
}

int main(int argc, char **argv){
	int i, hist[10];
	long part = LONG_MAX / 10, R;
	dev_random_ini();
	for(i = 0; i < 10; i++) hist[i] = 0;
	for(i = 0; i < 1000000; i++){
		R = lrand48();
		hist[(int)(R / part)]++;
	}
	for(i = 0; i < 10; i++)
		printf("%ld < X < %ld: %d\n", part*i, part*(i+1), hist[i]);
	return 0;
}

gcc -lm 2.c && ./a.out 
0 < X < 922337203685477580: 1000000
922337203685477580 < X < 1844674407370955160: 0
1844674407370955160 < X < 2767011611056432740: 0
2767011611056432740 < X < 3689348814741910320: 0
3689348814741910320 < X < 4611686018427387900: 0
4611686018427387900 < X < 5534023222112865480: 0
5534023222112865480 < X < 6456360425798343060: 0
6456360425798343060 < X < 7378697629483820640: 0
7378697629483820640 < X < 8301034833169298220: 0
8301034833169298220 < X < 9223372036854775800: 0

cat 3.c 
#include <stdio.h>		// printf
#include <stdlib.h>		// random number functions
#include <sys/time.h>		// gettimeofday
#include <limits.h>		// LONG_MAX
#include <math.h>		// floor
#include <unistd.h>		// read
#include <sys/types.h>		// open
#include <sys/stat.h>		// open
#include <fcntl.h>		// open

double dtime(){
	double t;
	struct timeval tv;
	gettimeofday(&tv, NULL);
	t = tv.tv_sec + ((double)tv.tv_usec)/1e6;
	return t;
}

void urandom_ini(){
	double tt = dtime() * 1e6;
	double mx = (double)UINT_MAX;
	srand((unsigned int)(tt - mx * floor(tt/mx)));
}

void dev_random_ini(){
	unsigned int r_ini;
	int fd = open("/dev/random", O_RDONLY);
	if(-1 == fd){fprintf(stderr,"Error open /dev/random!\n"); urandom_ini(); return;}
	if(sizeof(int) != read(fd, &r_ini, sizeof(int))){
		fprintf(stderr,"Error read /dev/random!\n"); urandom_ini();
		close(fd); return;
	}
	close(fd);
	srand(r_ini);
}

typedef union {
	long L;
	double D;
	unsigned int I[2];
} LD;

long rnd(){
	LD r;
	r.I[0] = (rand() << 1 | rand());
	r.I[1] = rand();
	return r.L;
}

int main(int argc, char **argv){
	int i, hist[10], s = 0;
	long part = LONG_MAX / 10, R;
	dev_random_ini();
	for(i = 0; i < 10; i++) hist[i] = 0;
	for(i = 0; i < 1000000; i++){
		R = rnd();
		hist[(int)(R / part)]++;
	}
	for(i = 0; i < 10; i++){
		printf("%ld < X < %ld: %d\n", part*i, part*(i+1), hist[i]);
		s += hist[i];
	}
	printf("Total: %d numbers\n", s);
	return 0;
}

gcc -lm 3.c && ./a.out 
0 < X < 922337203685477580: 99989
922337203685477580 < X < 1844674407370955160: 100490
1844674407370955160 < X < 2767011611056432740: 100048
2767011611056432740 < X < 3689348814741910320: 99839
3689348814741910320 < X < 4611686018427387900: 99679
4611686018427387900 < X < 5534023222112865480: 99907
5534023222112865480 < X < 6456360425798343060: 100391
6456360425798343060 < X < 7378697629483820640: 99900
7378697629483820640 < X < 8301034833169298220: 99698
8301034833169298220 < X < 9223372036854775800: 100059
Total: 1000000 numbers

Make distribution of 10000000 random numbers
0 < X < 922337203685477580: 999239
922337203685477580 < X < 1844674407370955160: 1000186
1844674407370955160 < X < 2767011611056432740: 1000296
2767011611056432740 < X < 3689348814741910320: 1001116
3689348814741910320 < X < 4611686018427387900: 1000714
4611686018427387900 < X < 5534023222112865480: 998833
5534023222112865480 < X < 6456360425798343060: 999500
6456360425798343060 < X < 7378697629483820640: 1000769
7378697629483820640 < X < 8301034833169298220: 1000420
8301034833169298220 < X < 9223372036854775800: 998927
Total: 10000000 numbers

LCG(2^48, 25214903917, 11) | 4.1 | 0.65 |  4 |  21 | -- *"LCG(2^48, 25214903917, 11) is the drand48" 
                                                         from the Unix standard library.
LCG(2^31–1, 16807, 0)      | 3.8 | 3.6  |  3 |  42 | -- *rand() в BSD-шном libc
LCG(2^32, 69069, 1)        | 3.2 | 0.67 | 11 | 106 | -- *часто используется из-за простоты запоминания
Unix-random-32             | 4.7 | 1.6  |  5 | 101 | -- *это все rand() в glibc. "The Unix-random’s are RNGs of 
Unix-random-64             | 4.7 | 1.5  |  4 |  57 | --  different sizes derived from BSD Unix and available 
Unix-random-128            | 4.7 | 1.5  |  2 |  13 | 19  as function random() on several Unix-Linux platforms. 
Unix-random-256            | 4.7 | 1.5  |  1 |   8 | 11  They are LFib(232, 7, 3, +), LFib(232, 15, 1, +), 
                                                         LFib(232, 31, 3, +) and LFib(232, 63, 1, +), with 
                                                         the least significant bit of each random number dropped."
ran3 / LFib(10^9,55,24,−)  | 2.2 | 0.9  |    |  11 | 17 *вариации используются много где, напр., Mono / .NET
Java.util.Random           | 6.3 | 0.76 |  1 |   9 | 21
MT19937                    | 4.3 | 1.6  |    |   2 |  2 *очень много где, включая c++11, Питон, Ruby и т.д.
LFib(2^48, 607, 273, +)    | 2.4 | 1.4  |    |   2 |  2 *есть в бусте