[C++]Генерация случайных чисел типа double во всем диапазоне

[Peek-a-boo]

Нужна помощь в написании программы, генерирующей псевдослучайные числа типа double на всем диапазоне, используя функцию rand().
Я пытался реализовать это, заполняя отдельно каждый из восьми байт случайными числами:

#include "stdio.h"
#include "stdlib.h"
int main()
{
	double x;
        int seed;
	int *px;
	px=(int*)&x;
        scanf ("%i",seed);
	srand(seed);
	*px=rand()%256+(rand()%256)*(1<<8)+(rand()%256)*(1<<16)+(rand()%256)*(1<<24)+(rand()%256)*(1<<32)+(rand()%256)*(1<<40)+(rand()%256)*(1<<48)+(rand()%256)*(1<<56);
	printf ("%lf", x);printf("\n");
	return 0;
}

Независимо от значения seed всегда выводится чилсо вида:
-92559598516065116000000000000000000 000000000000000000000000000000000.0 00000
Также нужно сделать проверку на получение Not a Number'а - насчет этого у меня вообще никаких идей нет.