Standard rand()
function gives numbers not big enough for me: I need unsigned long long
ones. How do we get really big random numbers? I tried modifying a simple hash function but it's too big, takes too long to run and never produces numbers which are less than 1e5!!
6 Answers
You can easily do this with std::uniform_int_distribution<unsigned long long>
.
Simple example code (taken from here, modified to use unsigned long long
):
#include <random>
#include <iostream>
int main()
{
std::random_device rd;
std::mt19937 gen(rd());
std::uniform_int_distribution<unsigned long long> dis(lowerBorder, upperBorder);
for (int n=0; n<10; ++n)
std::cout << dis(gen) << ' ';
std::cout << '\n';
}
Note that the seeding of the mersenne twister as done here for demo purposes is not perfect, for example see here.
-
1
-
-
@Alexandru No, that will not improve randomness. You can try if it improves speed on your environment, that might potentially work.– Baum mit Augen ♦Dec 20, 2017 at 12:41
Here's a portable C99 solution that returns a random 64-bit number:
unsigned long long llrand() {
unsigned long long r = 0;
for (int i = 0; i < 5; ++i) {
r = (r << 15) | (rand() & 0x7FFF);
}
return r & 0xFFFFFFFFFFFFFFFFULL;
}
Explanation: rand()
returns integers in the range 0 to RAND_MAX
and RAND_MAX
is only guaranteed to be at least 32,767 (15 random bits). long long
is guaranteed to have 64 bits but may be larger.
-
4Be really weary of this. It may seem reasonable but the numbers it produces may not (actually, probably will not) pass many statistical tests of randomness. Jan 23, 2015 at 18:11
-
1@AlterMann As I said,
long long
might be larger than 64 bits. If it is 64 bits, the compiler should make it a noop. Jan 23, 2015 at 18:17 -
2@LưuVĩnhPhúc MSVC (which has a reasonable market share) has
#define RAND_MAX 0x7fff
Jan 23, 2015 at 18:28 -
1@NikBougalis can you provide a citation or an example of a test that this would fail on (that rand() by itself would pass)– FoonMar 30, 2015 at 18:58
-
2@NikBougalis Still not seeing issue; assuming rand() is uniformly distributed, not seeing any number that can't be generated or any that would be generated more or less often. With that said, ran this against dieharder 3.31.1; there were two "good" tests that came back weak [not failed] (as opposed to just using the output of rand(), which had many failures, presumably due to rand() not return 32 bits of randomness at a time, or with uchar x = rand() % 256 and outputting 8 bits at a time; only 1 suspect test came back weak (rest passed)– FoonMar 31, 2015 at 1:49
If you want just to produce unsigned long long from value returned by rand() and do not care about the characteristics of the result consider the following function that must be compiler version and platform independent (because no "magic numbers" are used):
// this header has RAND_MAX value
#include <stdlib.h>
// and this header has ULLONG_MAX
#include <limits.h>
unsigned long long ullrand()
// Produces pseudo-random numbers from 0 to ULLONG_MAX
// by filling all bits of unsigned long long integer number
// with bits of several "small" integer numbers generated by rand()
{
unsigned long long myrndnum = 0; // at the beginning just zero
unsigned long long counter = ULLONG_MAX; // at the beginning we have all bits set as 1
// ... and while at least one bit is still set to 1
while(counter > 0) {
myrndnum = (myrndnum * (RAND_MAX + 1)) + rand(); // fill some bits from rand()
counter /= (RAND_MAX + 1); // decrease number of 1-bits in counter
}
// Return the result
return myrndnum;
}
But if you want some sequence of random numbers with certain predetermined characteristics, you should look in some specific guides or math books. E.g. https://www.gnu.org/software/gsl/manual/html_node/Random-number-generator-algorithms.html
You didn't ask for a specific OS and the answers here are really good, but on Linux (and probably on other OSes too) you can also read from a random device.
Example:
#include <stdio.h>
#include <assert.h>
#define RANDDEV "/dev/urandom"
unsigned long long bigrand(void) {
FILE *rdp;
unsigned long long num;
rdp = fopen(RANDDEV, "rb");
assert(rdp);
assert(fread(&num, sizeof(num), 1, rdp) == 1);
fclose(rdp);
return num;
}
Written on mobile, may have bugs. :P
-
1Maybe you should turn off buffering, or use read(), to avoid wasting random bytes? Jan 24, 2015 at 22:49
You can also use the boost library (taken from link):
#include <ctime> // std::time
#include <boost/random/mersenne_twister.hpp>
#include <boost/random/linear_congruential.hpp>
#include <boost/random/uniform_real.hpp>
#include <boost/random/variate_generator.hpp>
#include <boost/generator_iterator.hpp>
int main()
{
long long my_min = 1;
long long my_max = 1e5;
boost::mt19937 generator(static_cast<unsigned int>(std::time(0)));
boost::variate_generator<boost::mt19937&, boost::uniform_real<> >
die_gen(generator, boost::uniform_real<> (my_min, my_max));
boost::generator_iterator<boost::variate_generator<boost::mt19937&, boost::uniform_real<> > > die(&die_gen);
std::cout<<"Generated random numbers: \n";
for (int i=0; i <10 ; i++)
{
std::cout<< static_cast<long long>(*die++) << std::endl;
}
return 0;
}
try this:
long N=1000000;
long randNumber;
for(long i=0;i<N;i++)
randNumber=i+rand()
-
1While this code may answer the question, it would be better to include some context, explaining how it works and when to use it. Code-only answers are not useful in the long run. Mar 28, 2016 at 18:38
-
1This is a terrible method to generate random numbers. It's slow, badly distributed, and causes undefined behavior due to signed overflow.– interjayFeb 13, 2017 at 16:39
int
-wide number of bits at a time into theunsigned long long
using bit masks.rand()
and using some of the better RNGs c++ provides en.cppreference.com/w/cpp/numeric/random