Альтернативный алгоритм генерации псевдослучайных чисел

namespace Random
{
    typedef long long rand_t;
    static rand_t srand = time(nullptr);
 
    rand_t random()
    {
        srand ^= srand << 21;
        srand ^= srand >> 35;
        srand ^= srand << 4;
        return srand;
    }
 
    rand_t random(rand_t min, rand_t max)
    {
        if (min > max) return max;
        max = max + 1 - min;
        rand_t out = random() % max;
        out = out < 0 ? -out : out;
        return out + min;
    }
};

srand(time(0));
//x-диапазон значений 
int a = rand() % (2 * x) - x;

#include <iostream>
#include <ctime>
#include <cmath>
 typedef long long rand_t;
 
    rand_t random()
    {
        static rand_t srand = time(NULL);
        srand ^= srand << 21;
        srand ^= srand >> 35;
        srand ^= srand << 4;
        return srand;
    }
 
    rand_t random(rand_t min, rand_t max)
    {
        if (min > max) return max;
        max = max + 1 - min;
        rand_t out = random() % max;
        out = out < 0 ? -out : out;
        return out + min;
    }
int main()
{
    const int N=100;
    int val[N+1]={0};
    for(int i=1;i<=2000000;i++)
    {
        val[random(0,100)]++;
    }
    double avg=0;
    for(int i=0;i<N;i++)
    {
        avg+=val[i];
        std::cout<<val[i]<<'\t';
        if( (i+1)%10==0)
            std::cout<<'\n';
    }
    avg/=N;
    double disp=0;
    for(int i=0;i<N;i++)
        disp+=(avg-val[i])*(avg-val[i]);
    disp/=N-1;
    double sigma=sqrt(disp);
    std::cout<<"\naverage="<<avg<<"\tsigma="<<sigma
              <<'\n'<<sigma/avg*100.0<<'%'<<std::endl;
    system("pause");
    return 0;
}