random — Generate pseudo-random numbers (2024)

Source code: Lib/random.py

This module implements pseudo-random number generators for variousdistributions.

For integers, there is uniform selection from a range. For sequences, there isuniform selection of a random element, a function to generate a randompermutation of a list in-place, and a function for random sampling withoutreplacement.

On the real line, there are functions to compute uniform, normal (Gaussian),lognormal, negative exponential, gamma, and beta distributions. For generatingdistributions of angles, the von Mises distribution is available.

Almost all module functions depend on the basic function random(), whichgenerates a random float uniformly in the half-open range 0.0 <= X < 1.0.Python uses the Mersenne Twister as the core generator. It produces 53-bit precisionfloats and has a period of 2**19937-1. The underlying implementation in C isboth fast and threadsafe. The Mersenne Twister is one of the most extensivelytested random number generators in existence. However, being completelydeterministic, it is not suitable for all purposes, and is completely unsuitablefor cryptographic purposes.

The functions supplied by this module are actually bound methods of a hiddeninstance of the random.Random class. You can instantiate your owninstances of Random to get generators that don’t share state.

Class Random can also be subclassed if you want to use a differentbasic generator of your own devising: see the documentation on that class formore details.

The random module also provides the SystemRandom class whichuses the system function os.urandom() to generate random numbersfrom sources provided by the operating system.


The pseudo-random generators of this module should not be used forsecurity purposes. For security or cryptographic uses, see thesecrets module.

See also

M. Matsumoto and T. Nishimura, “Mersenne Twister: A 623-dimensionallyequidistributed uniform pseudorandom number generator”, ACM Transactions onModeling and Computer Simulation Vol. 8, No. 1, January pp.3–30 1998.

Complementary-Multiply-with-Carry recipe for a compatible alternativerandom number generator with a long period and comparatively simple updateoperations.

Bookkeeping functions

random.seed(a=None, version=2)

Initialize the random number generator.

If a is omitted or None, the current system time is used. Ifrandomness sources are provided by the operating system, they are usedinstead of the system time (see the os.urandom() function for detailson availability).

If a is an int, it is used directly.

With version 2 (the default), a str, bytes, or bytearrayobject gets converted to an int and all of its bits are used.

With version 1 (provided for reproducing random sequences from older versionsof Python), the algorithm for str and bytes generates anarrower range of seeds.

Changed in version 3.2: Moved to the version 2 scheme which uses all of the bits in a string seed.

Changed in version 3.11: The seed must be one of the following types:None, int, float, str,bytes, or bytearray.


Return an object capturing the current internal state of the generator. Thisobject can be passed to setstate() to restore the state.


state should have been obtained from a previous call to getstate(), andsetstate() restores the internal state of the generator to what it was atthe time getstate() was called.

Functions for bytes


Generate n random bytes.

This method should not be used for generating security tokens.Use secrets.token_bytes() instead.

Added in version 3.9.

Functions for integers

random.randrange(start, stop[, step])

Return a randomly selected element from range(start, stop, step).

This is roughly equivalent to choice(range(start, stop, step)) butsupports arbitrarily large ranges and is optimized for common cases.

The positional argument pattern matches the range() function.

Keyword arguments should not be used because they can be interpretedin unexpected ways. For example randrange(start=100) is interpretedas randrange(0, 100, 1).

Changed in version 3.2: randrange() is more sophisticated about producing equally distributedvalues. Formerly it used a style like int(random()*n) which could produceslightly uneven distributions.

Changed in version 3.12: Automatic conversion of non-integer types is no longer supported.Calls such as randrange(10.0) and randrange(Fraction(10, 1))now raise a TypeError.

random.randint(a, b)

Return a random integer N such that a <= N <= b. Alias forrandrange(a, b+1).


Returns a non-negative Python integer with k random bits. This methodis supplied with the Mersenne Twister generator and some other generatorsmay also provide it as an optional part of the API. When available,getrandbits() enables randrange() to handle arbitrarily largeranges.

Changed in version 3.9: This method now accepts zero for k.

Functions for sequences


Return a random element from the non-empty sequence seq. If seq is empty,raises IndexError.

random.choices(population, weights=None, *, cum_weights=None, k=1)

Return a k sized list of elements chosen from the population with replacement.If the population is empty, raises IndexError.

If a weights sequence is specified, selections are made according to therelative weights. Alternatively, if a cum_weights sequence is given, theselections are made according to the cumulative weights (perhaps computedusing itertools.accumulate()). For example, the relative weights[10, 5, 30, 5] are equivalent to the cumulative weights[10, 15, 45, 50]. Internally, the relative weights are converted tocumulative weights before making selections, so supplying the cumulativeweights saves work.

If neither weights nor cum_weights are specified, selections are madewith equal probability. If a weights sequence is supplied, it must bethe same length as the population sequence. It is a TypeErrorto specify both weights and cum_weights.

The weights or cum_weights can use any numeric type that interoperateswith the float values returned by random() (that includesintegers, floats, and fractions but excludes decimals). Weights are assumedto be non-negative and finite. A ValueError is raised if allweights are zero.

For a given seed, the choices() function with equal weightingtypically produces a different sequence than repeated calls tochoice(). The algorithm used by choices() uses floatingpoint arithmetic for internal consistency and speed. The algorithm usedby choice() defaults to integer arithmetic with repeated selectionsto avoid small biases from round-off error.

Added in version 3.6.

Changed in version 3.9: Raises a ValueError if all weights are zero.


Shuffle the sequence x in place.

To shuffle an immutable sequence and return a new shuffled list, usesample(x, k=len(x)) instead.

Note that even for small len(x), the total number of permutations of xcan quickly grow larger than the period of most random number generators.This implies that most permutations of a long sequence can never begenerated. For example, a sequence of length 2080 is the largest thatcan fit within the period of the Mersenne Twister random number generator.

Changed in version 3.11: Removed the optional parameter random.

random.sample(population, k, *, counts=None)

Return a k length list of unique elements chosen from the populationsequence. Used for random sampling without replacement.

Returns a new list containing elements from the population while leaving theoriginal population unchanged. The resulting list is in selection order so thatall sub-slices will also be valid random samples. This allows raffle winners(the sample) to be partitioned into grand prize and second place winners (thesubslices).

Members of the population need not be hashable or unique. If the populationcontains repeats, then each occurrence is a possible selection in the sample.

Repeated elements can be specified one at a time or with the optionalkeyword-only counts parameter. For example, sample(['red', 'blue'],counts=[4, 2], k=5) is equivalent to sample(['red', 'red', 'red', 'red','blue', 'blue'], k=5).

To choose a sample from a range of integers, use a range() object as anargument. This is especially fast and space efficient for sampling from a largepopulation: sample(range(10000000), k=60).

If the sample size is larger than the population size, a ValueErroris raised.

Changed in version 3.9: Added the counts parameter.

Changed in version 3.11: The population must be a sequence. Automatic conversion of setsto lists is no longer supported.

Discrete distributions

The following function generates a discrete distribution.

random.binomialvariate(n=1, p=0.5)

Binomial distribution.Return the number of successes for n independent trials with theprobability of success in each trial being p:

Mathematically equivalent to:

sum(random() < p for i in range(n))

The number of trials n should be a non-negative integer.The probability of success p should be between 0.0 <= p <= 1.0.The result is an integer in the range 0 <= X <= n.

Added in version 3.12.

Real-valued distributions

The following functions generate specific real-valued distributions. Functionparameters are named after the corresponding variables in the distribution’sequation, as used in common mathematical practice; most of these equations canbe found in any statistics text.


Return the next random floating point number in the range 0.0 <= X < 1.0

random.uniform(a, b)

Return a random floating point number N such that a <= N <= b fora <= b and b <= N <= a for b < a.

The end-point value b may or may not be included in the rangedepending on floating-point rounding in the expressiona + (b-a) * random().

random.triangular(low, high, mode)

Return a random floating point number N such that low <= N <= high andwith the specified mode between those bounds. The low and high boundsdefault to zero and one. The mode argument defaults to the midpointbetween the bounds, giving a symmetric distribution.

random.betavariate(alpha, beta)

Beta distribution. Conditions on the parameters are alpha > 0 andbeta > 0. Returned values range between 0 and 1.


Exponential distribution. lambd is 1.0 divided by the desiredmean. It should be nonzero. (The parameter would be called“lambda”, but that is a reserved word in Python.) Returned valuesrange from 0 to positive infinity if lambd is positive, and fromnegative infinity to 0 if lambd is negative.

Changed in version 3.12: Added the default value for lambd.

random.gammavariate(alpha, beta)

Gamma distribution. (Not the gamma function!) The shape andscale parameters, alpha and beta, must have positive values.(Calling conventions vary and some sources define ‘beta’as the inverse of the scale).

The probability distribution function is:

 x ** (alpha - 1) * math.exp(-x / beta)pdf(x) = -------------------------------------- math.gamma(alpha) * beta ** alpha
random.gauss(mu=0.0, sigma=1.0)

Normal distribution, also called the Gaussian distribution.mu is the mean,and sigma is the standard deviation. This is slightly faster thanthe normalvariate() function defined below.

Multithreading note: When two threads call this functionsimultaneously, it is possible that they will receive thesame return value. This can be avoided in three ways.1) Have each thread use a different instance of the randomnumber generator. 2) Put locks around all calls. 3) Use theslower, but thread-safe normalvariate() function instead.

Changed in version 3.11: mu and sigma now have default arguments.

random.lognormvariate(mu, sigma)

Log normal distribution. If you take the natural logarithm of thisdistribution, you’ll get a normal distribution with mean mu and standarddeviation sigma. mu can have any value, and sigma must be greater thanzero.

random.normalvariate(mu=0.0, sigma=1.0)

Normal distribution. mu is the mean, and sigma is the standard deviation.

Changed in version 3.11: mu and sigma now have default arguments.

random.vonmisesvariate(mu, kappa)

mu is the mean angle, expressed in radians between 0 and 2*pi, and kappais the concentration parameter, which must be greater than or equal to zero. Ifkappa is equal to zero, this distribution reduces to a uniform random angleover the range 0 to 2*pi.


Pareto distribution. alpha is the shape parameter.

random.weibullvariate(alpha, beta)

Weibull distribution. alpha is the scale parameter and beta is the shapeparameter.

Alternative Generator

class random.Random([seed])

Class that implements the default pseudo-random number generator used by therandom module.

Changed in version 3.11: Formerly the seed could be any hashable object. Now it is limited to:None, int, float, str,bytes, or bytearray.

Subclasses of Random should override the following methods if theywish to make use of a different basic generator:

seed(a=None, version=2)

Override this method in subclasses to customise the seed()behaviour of Random instances.


Override this method in subclasses to customise the getstate()behaviour of Random instances.


Override this method in subclasses to customise the setstate()behaviour of Random instances.


Override this method in subclasses to customise the random()behaviour of Random instances.

Optionally, a custom generator subclass can also supply the following method:


Override this method in subclasses to customise thegetrandbits() behaviour of Random instances.

class random.SystemRandom([seed])

Class that uses the os.urandom() function for generating random numbersfrom sources provided by the operating system. Not available on all systems.Does not rely on software state, and sequences are not reproducible. Accordingly,the seed() method has no effect and is ignored.The getstate() and setstate() methods raiseNotImplementedError if called.

Notes on Reproducibility

Sometimes it is useful to be able to reproduce the sequences given by apseudo-random number generator. By reusing a seed value, the same sequence should bereproducible from run to run as long as multiple threads are not running.

Most of the random module’s algorithms and seeding functions are subject tochange across Python versions, but two aspects are guaranteed not to change:

  • If a new seeding method is added, then a backward compatible seeder will beoffered.

  • The generator’s random() method will continue to produce the samesequence when the compatible seeder is given the same seed.


Basic examples:

>>> random() # Random float: 0.0 <= x < 1.00.37444887175646646>>> uniform(2.5, 10.0) # Random float: 2.5 <= x <= 10.03.1800146073117523>>> expovariate(1 / 5) # Interval between arrivals averaging 5 seconds5.148957571865031>>> randrange(10) # Integer from 0 to 9 inclusive7>>> randrange(0, 101, 2) # Even integer from 0 to 100 inclusive26>>> choice(['win', 'lose', 'draw']) # Single random element from a sequence'draw'>>> deck = 'ace two three four'.split()>>> shuffle(deck) # Shuffle a list>>> deck['four', 'two', 'ace', 'three']>>> sample([10, 20, 30, 40, 50], k=4) # Four samples without replacement[40, 10, 50, 30]


>>> # Six roulette wheel spins (weighted sampling with replacement)>>> choices(['red', 'black', 'green'], [18, 18, 2], k=6)['red', 'green', 'black', 'black', 'red', 'black']>>> # Deal 20 cards without replacement from a deck>>> # of 52 playing cards, and determine the proportion of cards>>> # with a ten-value: ten, jack, queen, or king.>>> deal = sample(['tens', 'low cards'], counts=[16, 36], k=20)>>> deal.count('tens') / 200.15>>> # Estimate the probability of getting 5 or more heads from 7 spins>>> # of a biased coin that settles on heads 60% of the time.>>> sum(binomialvariate(n=7, p=0.6) >= 5 for i in range(10_000)) / 10_0000.4169>>> # Probability of the median of 5 samples being in middle two quartiles>>> def trial():...  return 2_500 <= sorted(choices(range(10_000), k=5))[2] < 7_500...>>> sum(trial() for i in range(10_000)) / 10_0000.7958

Example of statistical bootstrapping using resamplingwith replacement to estimate a confidence interval for the mean of a sample:

# https://www.thoughtco.com/example-of-bootstrapping-3126155from statistics import fmean as meanfrom random import choicesdata = [41, 50, 29, 37, 81, 30, 73, 63, 20, 35, 68, 22, 60, 31, 95]means = sorted(mean(choices(data, k=len(data))) for i in range(100))print(f'The sample mean of {mean(data):.1f} has a 90% confidence ' f'interval from {means[5]:.1f} to {means[94]:.1f}')

Example of a resampling permutation testto determine the statistical significance or p-value of an observed differencebetween the effects of a drug versus a placebo:

# Example from "Statistics is Easy" by Dennis Shasha and Manda Wilsonfrom statistics import fmean as meanfrom random import shuffledrug = [54, 73, 53, 70, 73, 68, 52, 65, 65]placebo = [54, 51, 58, 44, 55, 52, 42, 47, 58, 46]observed_diff = mean(drug) - mean(placebo)n = 10_000count = 0combined = drug + placebofor i in range(n): shuffle(combined) new_diff = mean(combined[:len(drug)]) - mean(combined[len(drug):]) count += (new_diff >= observed_diff)print(f'{n} label reshufflings produced only {count} instances with a difference')print(f'at least as extreme as the observed difference of {observed_diff:.1f}.')print(f'The one-sided p-value of {count / n:.4f} leads us to reject the null')print(f'hypothesis that there is no difference between the drug and the placebo.')

Simulation of arrival times and service deliveries for a multiserver queue:

from heapq import heapify, heapreplacefrom random import expovariate, gaussfrom statistics import mean, quantilesaverage_arrival_interval = 5.6average_service_time = 15.0stdev_service_time = 3.5num_servers = 3waits = []arrival_time = 0.0servers = [0.0] * num_servers # time when each server becomes availableheapify(servers)for i in range(1_000_000): arrival_time += expovariate(1.0 / average_arrival_interval) next_server_available = servers[0] wait = max(0.0, next_server_available - arrival_time) waits.append(wait) service_duration = max(0.0, gauss(average_service_time, stdev_service_time)) service_completed = arrival_time + wait + service_duration heapreplace(servers, service_completed)print(f'Mean wait: {mean(waits):.1f} Max wait: {max(waits):.1f}')print('Quartiles:', [round(q, 1) for q in quantiles(waits)])

See also

Statistics for Hackersa video tutorial byJake Vanderplason statistical analysis using just a few fundamental conceptsincluding simulation, sampling, shuffling, and cross-validation.

Economics Simulationa simulation of a marketplace byPeter Norvig that shows effectiveuse of many of the tools and distributions provided by this module(gauss, uniform, sample, betavariate, choice, triangular, and randrange).

A Concrete Introduction to Probability (using Python)a tutorial by Peter Norvig coveringthe basics of probability theory, how to write simulations, andhow to perform data analysis using Python.


These recipes show how to efficiently make random selectionsfrom the combinatoric iterators in the itertools module:

def random_product(*args, repeat=1): "Random selection from itertools.product(*args, **kwds)" pools = [tuple(pool) for pool in args] * repeat return tuple(map(random.choice, pools))def random_permutation(iterable, r=None): "Random selection from itertools.permutations(iterable, r)" pool = tuple(iterable) r = len(pool) if r is None else r return tuple(random.sample(pool, r))def random_combination(iterable, r): "Random selection from itertools.combinations(iterable, r)" pool = tuple(iterable) n = len(pool) indices = sorted(random.sample(range(n), r)) return tuple(pool[i] for i in indices)def random_combination_with_replacement(iterable, r): "Choose r elements with replacement. Order the result to match the iterable." # Result will be in set(itertools.combinations_with_replacement(iterable, r)). pool = tuple(iterable) n = len(pool) indices = sorted(random.choices(range(n), k=r)) return tuple(pool[i] for i in indices)

The default random() returns multiples of 2⁻⁵³ in the range0.0 ≤ x < 1.0. All such numbers are evenly spaced and are exactlyrepresentable as Python floats. However, many other representablefloats in that interval are not possible selections. For example,0.05954861408025609 isn’t an integer multiple of 2⁻⁵³.

The following recipe takes a different approach. All floats in theinterval are possible selections. The mantissa comes from a uniformdistribution of integers in the range 2⁵² ≤ mantissa < 2⁵³. Theexponent comes from a geometric distribution where exponents smallerthan -53 occur half as often as the next larger exponent.

from random import Randomfrom math import ldexpclass FullRandom(Random): def random(self): mantissa = 0x10_0000_0000_0000 | self.getrandbits(52) exponent = -53 x = 0 while not x: x = self.getrandbits(32) exponent += x.bit_length() - 32 return ldexp(mantissa, exponent)

All real valued distributionsin the class will use the new method:

>>> fr = FullRandom()>>> fr.random()0.05954861408025609>>> fr.expovariate(0.25)8.87925541791544

The recipe is conceptually equivalent to an algorithm that chooses fromall the multiples of 2⁻¹⁰⁷⁴ in the range 0.0 ≤ x < 1.0. All suchnumbers are evenly spaced, but most have to be rounded down to thenearest representable Python float. (The value 2⁻¹⁰⁷⁴ is the smallestpositive unnormalized float and is equal to math.ulp(0.0).)

See also

Generating Pseudo-random Floating-Point Values apaper by Allen B. Downey describing ways to generate morefine-grained floats than normally generated by random().

random — Generate pseudo-random numbers (2024)


How to generate random numbers in pseudocode? ›

Example Algorithm for Pseudo-Random Number Generator
  1. Accept some initial input number, that is a seed or key.
  2. Apply that seed in a sequence of mathematical operations to generate the result. ...
  3. Use that resulting random number as the seed for the next iteration.
  4. Repeat the process to emulate randomness.
Oct 26, 2020

Do random number generators produce pseudo-random numbers? ›

Computerized random number generators produce what are called pseudo-random numbers, in that the numbers are produced by a deterministic process that gives exactly the same sequence of numbers every time.

How to make a pseudo-random number? ›

The algorithm is as follows: take any number, square it, remove the middle digits of the resulting number as the "random number", then use that number as the seed for the next iteration.

What is the difference between random and pseudorandom numbers? ›

Pseudorandom Number Generation. Software-generated random numbers only are pseudorandom. They are not truly random because the computer uses an algorithm based on a distribution, and are not secure because they rely on deterministic, predictable algorithms.

Which command is used to generate pseudo-random numbers? ›

If u is a uniform random number on (0,1) , then using X = F - 1 ( U ) generates a random number X from a continuous distribution with specified cdf F . rng('default') % For reproducibility mu = 1; X = expinv(rand(1e4,1),mu);

How to write a pseudo code? ›

How to Write Pseudocode. Always capitalize the initial word (often one of the main six constructs). Make only one statement per line. Indent to show hierarchy, improve readability and show nested constructs.

What are the two functions to generate the pseudo random number? ›

The Excel RAND and RANDBETWEEN functions generate pseudo-random numbers from the Uniform distribution, aka rectangular distribution, where there is equal probability for all values that a random variable can take on.

Which method produces a pseudorandom number? ›

Pseudo Random Number Generator(PRNG) refers to an algorithm that uses mathematical formulas to produce sequences of random numbers. PRNGs generate a sequence of numbers approximating the properties of random numbers. A PRNG starts from an arbitrary starting state using a seed state.

How does RNG actually work? ›

Random number generators are typically software, pseudo random number generators. Their outputs are not truly random numbers. Instead they rely on algorithms to mimic the selection of a value to approximate true randomness.

What is a pseudorandom generator? ›

A pseudo-random number generator (PRNG) is a program written for, and used in, probability and statistics applications when large quantities of random digits are needed. Most of these programs produce endless strings of single-digit numbers, usually in base 10, known as the decimal system.

Why is it called pseudo-random numbers? ›

Pseudo means false, in the sense that the number are not really random! They are generated according to a deterministic algorithm whose aim is to imitate as closely as possible what randomness would look like.

What is a pseudo-random code called? ›

Pseudo-Random Codes

The first pseudo-random code is called the C/A (Coarse Acquisition) code. It modulates the L1 carrier. It repeats every 1023 bits and modulates at a 1MHz rate.

What is a disadvantage of a pseudorandom number generator? ›

The Bad Guys

Since PRNGs are deterministic, it is comparable to knowing the poker hands before they are dealt. Due to the vast cryptographic applications of CSPRNGs, many different parties may be working against their security.

Is there a truly random number generator? ›

RANDOM.ORG offers true random numbers to anyone on the Internet. The randomness comes from atmospheric noise, which for many purposes is better than the pseudo-random number algorithms typically used in computer programs.

How do computers generate pseudorandom numbers? ›

Computers generate random numbers in two main categories: true random numbers and pseudo-random numbers. True random numbers are generated by measuring physical phenomena like atmospheric noise or radioactive decay, while algorithms and a seed value generate pseudo-random numbers.

How do you generate a random number in code? ›

The rand function in C

In C rand() returns a value between 0 and a large integer (called by the name "RAND_MAX" found in stdlib. h). To convert from the C format to the Matlab format, in C we would say: "rand() / (float)RAND_MAX". To create a boolean value (true/false) for a flip of a coin!

How do you make a random number generator in Makecode? ›

Select Set "Variable" to 0

Click category Math and select the pick random 0 to 10 block. This will create your random number generator. Type 1 to set the random number generator to choose between two different numbers (1 & 0).

What is the formula to create random number? ›

The formula is =RAND()*(b−a)+a. If you want a number between 20 and 60 — where a is the lower bound value and b is the upper bound value — the formula would be =RAND()*(60−20)+20. Note that this formula will also not return a value equal to that of the upper bound value.

What are the two functions to generate the pseudo-random number? ›

The Excel RAND and RANDBETWEEN functions generate pseudo-random numbers from the Uniform distribution, aka rectangular distribution, where there is equal probability for all values that a random variable can take on.

Top Articles
Latest Posts
Article information

Author: Fr. Dewey Fisher

Last Updated:

Views: 6469

Rating: 4.1 / 5 (62 voted)

Reviews: 93% of readers found this page helpful

Author information

Name: Fr. Dewey Fisher

Birthday: 1993-03-26

Address: 917 Hyun Views, Rogahnmouth, KY 91013-8827

Phone: +5938540192553

Job: Administration Developer

Hobby: Embroidery, Horseback riding, Juggling, Urban exploration, Skiing, Cycling, Handball

Introduction: My name is Fr. Dewey Fisher, I am a powerful, open, faithful, combative, spotless, faithful, fair person who loves writing and wants to share my knowledge and understanding with you.