Документация Python неофициальный перевод

random.md

485 строк · 42.1 КБ · обычная страница · сырой текст · скачать

1> **Источник:** https://python-all.ru/3.12/library/random.html2>3> «Документация Python на русском» – неофициальный перевод официальной документации Python: версии от 2.6 до 3.16, полнотекстовый поиск, английский оригинал рядом с переводом. Эта Markdown-версия страницы предназначена для работы с LLM: вставьте её в ChatGPT, Claude или Cursor.45---67# `random` – Генерация псевдослучайных чисел89**Исходный код:** [Lib/random.py](https://python-all.ru/src/3.12/Lib/random.py)1011---1213Этот модуль реализует генераторы псевдослучайных чисел для различных распределений.1415Для целых чисел есть равномерный выбор из диапазона. Для последовательностей есть равномерный выбор случайного элемента, функция для генерации случайной перестановки списка на месте и функция для случайной выборки без возвращения.1617На вещественной прямой есть функции для вычисления равномерного, нормального (гауссовского), логнормального, отрицательного экспоненциального, гамма и бета распределений. Для генерации распределений углов доступно распределение фон Мизеса.1819Почти все функции модуля зависят от базовой функции [`random()`](https://python-all.ru/3.12/library/random.html#random.random), которая генерирует случайное число с плавающей запятой, равномерно распределённое в полуоткрытом интервале `0.0 <= X < 1.0`. Python использует Вихрь Мерсенна в качестве основного генератора. Он производит числа с плавающей запятой точностью 53 бита и имеет период 2\*\*19937-1. Базовая реализация на C является быстрой и потокобезопасной. Вихрь Мерсенна – один из наиболее тщательно протестированных генераторов случайных чисел из существующих. Однако, будучи полностью детерминированным, он не подходит для всех целей и совершенно непригоден для криптографических целей.2021Функции, предоставляемые этим модулем, на самом деле являются привязанными методами скрытого экземпляра класса [`random.Random`](https://python-all.ru/3.12/library/random.html#random.Random). Можно создать собственные экземпляры [`Random`](https://python-all.ru/3.12/library/random.html#random.Random), чтобы получить генераторы, не разделяющие состояние.2223Класс [`Random`](https://python-all.ru/3.12/library/random.html#random.Random) также можно подклассифицировать, если требуется использовать другой базовый генератор собственного изобретения: см. документацию по этому классу для получения дополнительных сведений.2425Модуль [`random`](https://python-all.ru/3.12/library/random.html#module-random) также предоставляет класс [`SystemRandom`](https://python-all.ru/3.12/library/random.html#random.SystemRandom), который использует системную функцию [`os.urandom()`](https://python-all.ru/3.12/library/os.html#os.urandom) для генерации случайных чисел из источников, предоставляемых операционной системой.2627> **Предупреждение**28>29> Псевдослучайные генераторы этого модуля не следует использовать для целей безопасности. Для безопасности или криптографического использования см. модуль [`secrets`](https://python-all.ru/3.12/library/secrets.html#module-secrets).3031> **См. также**32>33> М. Мацумото и Т. Нисимура, «Вихрь Мерсенна: 623-мерно равномерно распределенный генератор псевдослучайных чисел», ACM Transactions on Modeling and Computer Simulation том 8, № 1, январь, стр. 3–30, 1998.34>35> [Complementary-Multiply-with-Carry recipe](https://python-all.ru/3.12/library/random.html) для совместимого альтернативного генератора случайных чисел с длинным периодом и сравнительно простыми операциями обновления.3637## Учётные функции3839#### `random.seed(a=None, version=2)`4041Инициализирует генератор случайных чисел.4243Если *a* опущен или равен `None`, используется текущее системное время. Если источники случайности предоставляются операционной системой, они используются вместо системного времени (см. функцию [`os.urandom()`](https://python-all.ru/3.12/library/os.html#os.urandom) для получения подробной информации о доступности).4445Если *a* – целое число, оно используется напрямую.4647В версии 2 (по умолчанию) объект [`str`](https://python-all.ru/3.12/library/stdtypes.html#str), [`bytes`](https://python-all.ru/3.12/library/stdtypes.html#bytes) или [`bytearray`](https://python-all.ru/3.12/library/stdtypes.html#bytearray) преобразуется в [`int`](https://python-all.ru/3.12/library/functions.html#int) и используются все его биты.4849В версии 1 (предназначенной для воспроизведения случайных последовательностей из старых версий Python) алгоритм для [`str`](https://python-all.ru/3.12/library/stdtypes.html#str) и [`bytes`](https://python-all.ru/3.12/library/stdtypes.html#bytes) генерирует более узкий диапазон начальных значений (семян).5051Изменено в версии 3.2: Переход на схему версии 2, которая использует все биты строкового начального значения.5253Изменено в версии 3.11: Параметр *seed* должен быть одного из следующих типов: `None`, [`int`](https://python-all.ru/3.12/library/functions.html#int), [`float`](https://python-all.ru/3.12/library/functions.html#float), [`str`](https://python-all.ru/3.12/library/stdtypes.html#str), [`bytes`](https://python-all.ru/3.12/library/stdtypes.html#bytes) или [`bytearray`](https://python-all.ru/3.12/library/stdtypes.html#bytearray).5455#### `random.getstate()`5657Возвращает объект, фиксирующий текущее внутреннее состояние генератора. Этот объект можно передать в [`setstate()`](https://python-all.ru/3.12/library/random.html#random.setstate) для восстановления состояния.5859#### `random.setstate(state)`6061*state* должен быть получен из предыдущего вызова [`getstate()`](https://python-all.ru/3.12/library/random.html#random.getstate), а [`setstate()`](https://python-all.ru/3.12/library/random.html#random.setstate) восстанавливает внутреннее состояние генератора в то, которое было на момент вызова [`getstate()`](https://python-all.ru/3.12/library/random.html#random.getstate).6263## Функции для байтов6465#### `random.randbytes(n)`6667Генерирует *n* случайных байтов.6869Этот метод не следует использовать для генерации токенов безопасности. Вместо этого используйте [`secrets.token_bytes()`](https://python-all.ru/3.12/library/secrets.html#secrets.token_bytes).7071Добавлено в версии 3.9.7273## Функции для целых чисел7475#### `random.randrange(stop)`7677#### `random.randrange(start, stop[, step])`7879Возвращает случайно выбранный элемент из `range(start, stop, step)`.8081Это примерно эквивалентно `choice(range(start, stop, step))`, но поддерживает произвольно большие диапазоны и оптимизировано для распространённых случаев.8283Шаблон позиционных аргументов соответствует функции [`range()`](https://python-all.ru/3.12/library/stdtypes.html#range).8485Именованные аргументы не следует использовать, так как они могут быть интерпретированы неожиданным образом. Например, `randrange(start=100)` интерпретируется как `randrange(0, 100, 1)`.8687Изменено в версии 3.2: [`randrange()`](https://python-all.ru/3.12/library/random.html#random.randrange) стала более совершенной в генерации равномерно распределённых значений. Ранее использовался подход, подобный `int(random()*n)`, который мог давать несколько неравномерные распределения.8889Изменено в версии 3.12: Автоматическое преобразование нецелых типов больше не поддерживается. Вызовы, такие как `randrange(10.0)` и `randrange(Fraction(10, 1))`, теперь вызывают [`TypeError`](https://python-all.ru/3.12/library/exceptions.html#TypeError).9091#### `random.randint(a, b)`9293Возвращает случайное целое *N* такое, что `a <= N <= b`. Псевдоним для `randrange(a, b+1)`.9495#### `random.getrandbits(k)`9697Возвращает неотрицательное целое число Python с *k* случайными битами. Этот метод предоставляется генератором Mersenne Twister, и некоторые другие генераторы также могут реализовывать его как опциональную часть API. Если он доступен, [`getrandbits()`](https://python-all.ru/3.12/library/random.html#random.getrandbits) позволяет [`randrange()`](https://python-all.ru/3.12/library/random.html#random.randrange) работать с произвольно большими диапазонами.9899Изменено в версии 3.9: Этот метод теперь принимает ноль для *k*.100101## Функции для последовательностей102103#### `random.choice(seq)`104105Возвращает случайный элемент из непустой последовательности *seq*. Если *seq* пуста, возбуждает [`IndexError`](https://python-all.ru/3.12/library/exceptions.html#IndexError).106107#### `random.choices(population, weights=None, *, cum_weights=None, k=1)`108109Возвращает список размера *k* из элементов, выбранных из *population* с возвращением. Если *population* пуста, возбуждает [`IndexError`](https://python-all.ru/3.12/library/exceptions.html#IndexError).110111Если задана последовательность *weights*, выборки производятся в соответствии с относительными весами. Альтернативно, если задана последовательность *cum\_weights*, выборки производятся по накопленным весам (возможно, вычисленным с помощью [`itertools.accumulate()`](https://python-all.ru/3.12/library/itertools.html#itertools.accumulate)). Например, относительные веса `[10, 5, 30, 5]` эквивалентны накопленным весам `[10, 15, 45, 50]`. Внутренне относительные веса преобразуются в накопленные перед выполнением выборок, поэтому передача накопленных весов экономит работу.112113Если не указаны ни *weights*, ни *cum\_weights*, выборки производятся с равной вероятностью. Если передана последовательность весов, её длина должна равняться длине последовательности *population*. Указание и *weights*, и *cum\_weights* приводит к [`TypeError`](https://python-all.ru/3.12/library/exceptions.html#TypeError).114115*weights* или *cum\_weights* могут использовать любой числовой тип, совместимый со значениями [`float`](https://python-all.ru/3.12/library/functions.html#float), возвращаемыми [`random()`](https://python-all.ru/3.12/library/random.html#module-random) (сюда входят целые числа, числа с плавающей запятой и дроби, но не десятичные числа). Предполагается, что веса неотрицательны и конечны. Если все веса равны нулю, возбуждается [`ValueError`](https://python-all.ru/3.12/library/exceptions.html#ValueError).116117Для заданного начального значения функция [`choices()`](https://python-all.ru/3.12/library/random.html#random.choices) с равными весами обычно выдаёт другую последовательность, чем многократные вызовы [`choice()`](https://python-all.ru/3.12/library/random.html#random.choice). Алгоритм, используемый в [`choices()`](https://python-all.ru/3.12/library/random.html#random.choices), применяет арифметику с плавающей запятой для внутренней согласованности и скорости. Алгоритм, используемый в [`choice()`](https://python-all.ru/3.12/library/random.html#random.choice), по умолчанию использует целочисленную арифметику с повторными выборками, чтобы избежать малых смещений из-за ошибок округления.118119Добавлено в версии 3.6.120121Изменено в версии 3.9: Возбуждает [`ValueError`](https://python-all.ru/3.12/library/exceptions.html#ValueError), если все веса равны нулю.122123#### `random.shuffle(x)`124125Перемешивает последовательность *x* на месте.126127Чтобы перемешать неизменяемую последовательность и вернуть новый перемешанный список, используйте `sample(x, k=len(x))`.128129Обратите внимание, что даже для небольших `len(x)` общее число перестановок *x* может быстро превысить период большинства генераторов случайных чисел. Это означает, что большинство перестановок длинной последовательности никогда не могут быть сгенерированы. Например, последовательность длины 2080 – это самая длинная, которая умещается в период генератора случайных чисел Mersenne Twister.130131Изменено в версии 3.11: Удалён необязательный параметр *random*.132133#### `random.sample(population, k, *, counts=None)`134135Возвращает список длины *k* из уникальных элементов, выбранных из последовательности population. Используется для случайной выборки без возвращения.136137Возвращает новый список, содержащий элементы из population, не изменяя исходную population. Результирующий список упорядочен по порядку выбора, так что все подсписки также будут корректными случайными выборками. Это позволяет разделить победителей розыгрыша (выборку) на обладателей главного и второго призов (подсписки).138139Элементы population не обязаны быть [хэшируемыми](https://python-all.ru/3.12/glossary.html#term-hashable) или уникальными. Если population содержит повторения, то каждое вхождение является возможным элементом выборки.140141Повторяющиеся элементы можно указывать по одному или с помощью необязательного параметра *counts*, передаваемого только по ключу. Например, `sample(['red', 'blue'], counts=[4, 2], k=5)` эквивалентно `sample(['red', 'red', 'red', 'red', 'blue', 'blue'], k=5)`.142143Чтобы выбрать выборку из диапазона целых чисел, используйте объект [`range()`](https://python-all.ru/3.12/library/stdtypes.html#range) в качестве аргумента. Это особенно быстро и эффективно по памяти при выборке из большой популяции: `sample(range(10000000), k=60)`.144145Если размер выборки превышает размер популяции, возбуждается исключение [`ValueError`](https://python-all.ru/3.12/library/exceptions.html#ValueError).146147Изменено в версии 3.9: Добавлен параметр *counts*.148149Изменено в версии 3.11: Параметр *population* должен быть последовательностью. Автоматическое преобразование множеств в списки больше не поддерживается.150151## Дискретные распределения152153Следующая функция генерирует дискретное распределение.154155#### `random.binomialvariate(n=1, p=0.5)`156157[Биномиальное распределение](https://python-all.ru/3.12/library/random.html). Возвращает количество успехов для *n* независимых испытаний, где вероятность успеха в каждом испытании равна *p*:158159Математически эквивалентно:160161```python162sum(random() < p for i in range(n))163```164165Количество испытаний *n* должно быть неотрицательным целым числом. Вероятность успеха *p* должна быть в диапазоне `0.0 <= p <= 1.0`. Результат – целое число в диапазоне `0 <= X <= n`.166167Добавлено в версии 3.12.168169## Вещественные распределения170171Следующие функции генерируют конкретные вещественные распределения. Имена параметров функций соответствуют переменным в уравнении распределения, как принято в математической практике; большинство этих уравнений можно найти в любом учебнике по статистике.172173#### `random.random()`174175Возвращает следующее случайное число с плавающей запятой в диапазоне `0.0 <= X < 1.0`176177#### `random.uniform(a, b)`178179Возвращает случайное число с плавающей запятой *N* такое, что `a <= N <= b` при `a <= b` и `b <= N <= a` при `b < a`.180181Значение конечной точки `b` может как входить, так и не входить в диапазон в зависимости от округления чисел с плавающей запятой в выражении `a + (b-a) * random()`.182183#### `random.triangular(low, high, mode)`184185Возвращает случайное число с плавающей запятой *N* такое, что `low <= N <= high` и с указанной *модой* между этими границами. Границы *low* и *high* по умолчанию равны нулю и единице. Аргумент *mode* по умолчанию равен середине между границами, что даёт симметричное распределение.186187#### `random.betavariate(alpha, beta)`188189Бета-распределение. Условия на параметры: `alpha > 0` и `beta > 0`. Возвращаемые значения находятся в диапазоне от 0 до 1.190191#### `random.expovariate(lambd=1.0)`192193Экспоненциальное распределение. *lambd* – это 1.0, делённая на желаемое среднее. Оно должно быть ненулевым. (Параметр мог бы называться «lambda», но это зарезервированное слово в Python.) Возвращаемые значения находятся в диапазоне от 0 до положительной бесконечности, если *lambd* положительно, и от отрицательной бесконечности до 0, если *lambd* отрицательно.194195Изменено в версии 3.12: Добавлено значение по умолчанию для `lambd`.196197#### `random.gammavariate(alpha, beta)`198199Гамма-распределение. (*Не* гамма-функция!) Параметры формы и масштаба, *alpha* и *beta*, должны иметь положительные значения. (Соглашения о вызовах различаются, и некоторые источники определяют 'beta' как обратную величину масштаба.)200201Функция плотности распределения вероятности:202203```python204          x ** (alpha - 1) * math.exp(-x / beta)205pdf(x) =  --------------------------------------206            math.gamma(alpha) * beta ** alpha207```208209#### `random.gauss(mu=0.0, sigma=1.0)`210211Нормальное распределение, также называемое распределением Гаусса. *mu* – это среднее, а *sigma* – это стандартное отклонение. Эта функция немного быстрее, чем функция [`normalvariate()`](https://python-all.ru/3.12/library/random.html#random.normalvariate), определённая ниже.212213Примечание о многопоточности: когда два потока одновременно вызывают эту функцию, они могут получить одно и то же возвращаемое значение. Этого можно избежать тремя способами: 1) Пусть каждый поток использует свой экземпляр генератора случайных чисел. 2) Установите блокировки вокруг всех вызовов. 3) Используйте более медленную, но потокобезопасную функцию [`normalvariate()`](https://python-all.ru/3.12/library/random.html#random.normalvariate).214215Изменено в версии 3.11: *mu* и *sigma* теперь имеют аргументы по умолчанию.216217#### `random.lognormvariate(mu, sigma)`218219Логнормальное распределение. Если взять натуральный логарифм этого распределения, получится нормальное распределение со средним *mu* и стандартным отклонением *sigma*. *mu* может принимать любое значение, а *sigma* должна быть больше нуля.220221#### `random.normalvariate(mu=0.0, sigma=1.0)`222223Нормальное распределение. *mu* – это среднее, а *sigma* – стандартное отклонение.224225Изменено в версии 3.11: *mu* и *sigma* теперь имеют аргументы по умолчанию.226227#### `random.vonmisesvariate(mu, kappa)`228229*mu* – это средний угол, выраженный в радианах от 0 до 2\**pi*, а *kappa* – параметр концентрации, который должен быть больше или равен нулю. Если *kappa* равен нулю, это распределение сводится к равномерному случайному углу в диапазоне от 0 до 2\**pi*.230231#### `random.paretovariate(alpha)`232233Распределение Парето. *alpha* – параметр формы.234235#### `random.weibullvariate(alpha, beta)`236237Распределение Вейбулла. *alpha* – параметр масштаба, а *beta* – параметр формы.238239## Альтернативный генератор240241#### `class random.Random([seed])`242243Класс, реализующий генератор псевдослучайных чисел по умолчанию, используемый модулем [`random`](https://python-all.ru/3.12/library/random.html#module-random).244245Изменено в версии 3.11: Ранее *seed* мог быть любым хешируемым объектом. Теперь он ограничен: `None`, [`int`](https://python-all.ru/3.12/library/functions.html#int), [`float`](https://python-all.ru/3.12/library/functions.html#float), [`str`](https://python-all.ru/3.12/library/stdtypes.html#str), [`bytes`](https://python-all.ru/3.12/library/stdtypes.html#bytes) или [`bytearray`](https://python-all.ru/3.12/library/stdtypes.html#bytearray).246247Подклассы `Random` должны переопределить следующие методы, если они хотят использовать другой базовый генератор:248249#### `seed(a=None, version=2)`250251Переопределите этот метод в подклассах, чтобы настроить [`seed()`](https://python-all.ru/3.12/library/random.html#random.seed) поведение экземпляров `Random`.252253#### `getstate()`254255Переопределите этот метод в подклассах, чтобы настроить [`getstate()`](https://python-all.ru/3.12/library/random.html#random.getstate) поведение экземпляров `Random`.256257#### `setstate(state)`258259Переопределите этот метод в подклассах, чтобы настроить [`setstate()`](https://python-all.ru/3.12/library/random.html#random.setstate) поведение экземпляров `Random`.260261#### `random()`262263Переопределите этот метод в подклассах, чтобы настроить [`random()`](https://python-all.ru/3.12/library/random.html#random.random) поведение экземпляров `Random`.264265При желании подкласс пользовательского генератора может также предоставить следующий метод:266267#### `getrandbits(k)`268269Переопределите этот метод в подклассах, чтобы настроить [`getrandbits()`](https://python-all.ru/3.12/library/random.html#random.getrandbits) поведение экземпляров `Random`.270271#### `class random.SystemRandom([seed])`272273Класс, использующий функцию [`os.urandom()`](https://python-all.ru/3.12/library/os.html#os.urandom) для генерации случайных чисел из источников, предоставляемых операционной системой. Доступен не на всех системах. Не зависит от состояния программного обеспечения, и последовательности невоспроизводимы. Соответственно, метод [`seed()`](https://python-all.ru/3.12/library/random.html#random.seed) не действует и игнорируется. Методы [`getstate()`](https://python-all.ru/3.12/library/random.html#random.getstate) и [`setstate()`](https://python-all.ru/3.12/library/random.html#random.setstate) вызывают [`NotImplementedError`](https://python-all.ru/3.12/library/exceptions.html#NotImplementedError) при вызове.274275## Примечания о воспроизводимости276277Иногда бывает полезно воспроизводить последовательности, выдаваемые генератором псевдослучайных чисел. При повторном использовании одного и того же начального значения (seed) одна и та же последовательность должна воспроизводиться от запуска к запуску, если не запущено несколько потоков.278279Большинство алгоритмов и функций инициализации генератора модуля random могут меняться в разных версиях Python, но два аспекта гарантированно остаются неизменными:280281- Если будет добавлен новый метод инициализации, будет предложен обратно совместимый инициализатор.282- Метод [`random()`](https://python-all.ru/3.12/library/random.html#random.Random.random) генератора продолжит выдавать ту же последовательность, если совместимому инициализатору передать то же начальное значение.283284## Примеры285286Простые примеры:287288```python289>>> random()                          # Случайное число с плавающей точкой: 0.0 <= x < 1.02900.37444887175646646291292>>> uniform(2.5, 10.0)                # Случайное число с плавающей точкой: 2.5 <= x <= 10.02933.1800146073117523294295>>> expovariate(1 / 5)                # Интервал между поступлениями в среднем 5 секунд2965.148957571865031297298>>> randrange(10)                     # Целое число от 0 до 9 включительно2997300301>>> randrange(0, 101, 2)              # Чётное целое число от 0 до 100 включительно30226303304>>> choice(['win', 'lose', 'draw'])   # Один случайный элемент из последовательности305'draw'306307>>> deck = 'ace two three four'.split()308>>> shuffle(deck)                     # Перемешать список309>>> deck310['four', 'two', 'ace', 'three']311312>>> sample([10, 20, 30, 40, 50], k=4) # Четыре выборки без возвращения313[40, 10, 50, 30]314```315316Симуляции:317318```python319>>> # Шесть вращений колеса рулетки (взвешенная выборка с возвращением)320>>> choices(['red', 'black', 'green'], [18, 18, 2], k=6)321['red', 'green', 'black', 'black', 'red', 'black']322323>>> # Раздать 20 карт без возвращения из колоды324>>> # из 52 игральных карт и определить долю карт325>>> # со значением десять: десятка, валет, дама или король.326>>> deal = sample(['tens', 'low cards'], counts=[16, 36], k=20)327>>> deal.count('tens') / 203280.15329330>>> # Оценить вероятность выпадения 5 или более орлов за 7 подбрасываний331>>> # несимметричной монеты, которая выпадает орлом в 60% случаев.332>>> sum(binomialvariate(n=7, p=0.6) >= 5 for i in range(10_000)) / 10_0003330.4169334335>>> # Вероятность того, что медиана 5 выборок окажется в двух средних квартилях336>>> def trial():337...     return 2_500 <= sorted(choices(range(10_000), k=5))[2] < 7_500338...339>>> sum(trial() for i in range(10_000)) / 10_0003400.7958341```342343Пример [статистического бутстрепа](https://python-all.ru/3.12/library/random.html) с использованием повторной выборки с возвращением для оценки доверительного интервала для среднего значения выборки:344345```python346# https://www.thoughtco.com/example-of-bootstrapping-3126155347from statistics import fmean as mean348from random import choices349350data = [41, 50, 29, 37, 81, 30, 73, 63, 20, 35, 68, 22, 60, 31, 95]351means = sorted(mean(choices(data, k=len(data))) for i in range(100))352print(f'The sample mean of {mean(data):.1f} has a 90% confidence '353      f'interval from {means[5]:.1f} to {means[94]:.1f}')354```355356Пример [пермутационного теста с повторной выборкой](https://python-all.ru/3.12/library/random.html) для определения статистической значимости или [p-значения](https://python-all.ru/3.12/library/random.html) наблюдаемого различия между эффектами препарата и плацебо:357358```python359# Пример из книги "Statistics is Easy" Денниса Шаши и Манды Уилсон360from statistics import fmean as mean361from random import shuffle362363drug = [54, 73, 53, 70, 73, 68, 52, 65, 65]364placebo = [54, 51, 58, 44, 55, 52, 42, 47, 58, 46]365observed_diff = mean(drug) - mean(placebo)366367n = 10_000368count = 0369combined = drug + placebo370for i in range(n):371    shuffle(combined)372    new_diff = mean(combined[:len(drug)]) - mean(combined[len(drug):])373    count += (new_diff >= observed_diff)374375print(f'{n} label reshufflings produced only {count} instances with a difference')376print(f'at least as extreme as the observed difference of {observed_diff:.1f}.')377print(f'The one-sided p-value of {count / n:.4f} leads us to reject the null')378print(f'hypothesis that there is no difference between the drug and the placebo.')379```380381Моделирование времени поступления и обслуживания запросов для многоканальной очереди:382383```python384from heapq import heapify, heapreplace385from random import expovariate, gauss386from statistics import mean, quantiles387388average_arrival_interval = 5.6389average_service_time = 15.0390stdev_service_time = 3.5391num_servers = 3392393waits = []394arrival_time = 0.0395servers = [0.0] * num_servers  # время, когда каждый сервер становится доступным396heapify(servers)397for i in range(1_000_000):398    arrival_time += expovariate(1.0 / average_arrival_interval)399    next_server_available = servers[0]400    wait = max(0.0, next_server_available - arrival_time)401    waits.append(wait)402    service_duration = max(0.0, gauss(average_service_time, stdev_service_time))403    service_completed = arrival_time + wait + service_duration404    heapreplace(servers, service_completed)405406print(f'Mean wait: {mean(waits):.1f}   Max wait: {max(waits):.1f}')407print('Quartiles:', [round(q, 1) for q in quantiles(waits)])408```409410> **См. также**411>412> [Statistics for Hackers](https://python-all.ru/3.12/library/random.html) – видеоурок от [Jake Vanderplas](https://python-all.ru/3.12/library/random.html) по статистическому анализу, использующему всего несколько фундаментальных понятий, включая моделирование, выборку, перемешивание и кросс-валидацию.413>414> [Economics Simulation](https://python-all.ru/3.12/library/random.html) – симуляция рынка от [Peter Norvig](https://python-all.ru/3.12/library/random.html), демонстрирующая эффективное использование многих инструментов и распределений, предоставляемых этим модулем (gauss, uniform, sample, betavariate, choice, triangular и randrange).415>416> [A Concrete Introduction to Probability (using Python)](https://python-all.ru/3.12/library/random.html) – учебное пособие от [Peter Norvig](https://python-all.ru/3.12/library/random.html), охватывающее основы теории вероятностей, написание симуляций и выполнение анализа данных с помощью Python.417418## Рецепты419420Эти рецепты показывают, как эффективно делать случайный выбор из комбинаторных итераторов модуля [`itertools`](https://python-all.ru/3.12/library/itertools.html#module-itertools):421422```python423def random_product(*args, repeat=1):424    "Random selection from itertools.product(*args, **kwds)"425    pools = [tuple(pool) for pool in args] * repeat426    return tuple(map(random.choice, pools))427428def random_permutation(iterable, r=None):429    "Random selection from itertools.permutations(iterable, r)"430    pool = tuple(iterable)431    r = len(pool) if r is None else r432    return tuple(random.sample(pool, r))433434def random_combination(iterable, r):435    "Random selection from itertools.combinations(iterable, r)"436    pool = tuple(iterable)437    n = len(pool)438    indices = sorted(random.sample(range(n), r))439    return tuple(pool[i] for i in indices)440441def random_combination_with_replacement(iterable, r):442    "Choose r elements with replacement.  Order the result to match the iterable."443    # Результат будет в set(itertools.combinations_with_replacement(iterable, r)).444    pool = tuple(iterable)445    n = len(pool)446    indices = sorted(random.choices(range(n), k=r))447    return tuple(pool[i] for i in indices)448```449450По умолчанию [`random()`](https://python-all.ru/3.12/library/random.html#random.random) возвращает числа, кратные 2⁻⁵³, в диапазоне *0.0 ≤ x \< 1.0*. Все такие числа равномерно распределены и точно представимы как числа с плавающей запятой Python. Однако многие другие представимые числа с плавающей запятой в этом интервале не могут быть выбраны. Например, `0.05954861408025609` не является целым кратным 2⁻⁵³.451452Следующий рецепт использует другой подход. Все числа с плавающей запятой в интервале могут быть выбраны. Мантисса берётся из равномерного распределения целых чисел в диапазоне *2⁵² ≤ mantissa \< 2⁵³*. Показатель степени берётся из геометрического распределения, где показатели меньше *-53* встречаются вдвое реже, чем следующий больший показатель.453454```python455from random import Random456from math import ldexp457458class FullRandom(Random):459460    def random(self):461        mantissa = 0x10_0000_0000_0000 | self.getrandbits(52)462        exponent = -53463        x = 0464        while not x:465            x = self.getrandbits(32)466            exponent += x.bit_length() - 32467        return ldexp(mantissa, exponent)468```469470Все [вещественнозначные распределения](https://python-all.ru/3.12/library/random.html#real-valued-distributions) в классе будут использовать новый метод:471472```python473>>> fr = FullRandom()474>>> fr.random()4750.05954861408025609476>>> fr.expovariate(0.25)4778.87925541791544478```479480Этот рецепт концептуально эквивалентен алгоритму, который выбирает из всех чисел, кратных 2⁻¹⁰⁷⁴, в диапазоне *0.0 ≤ x \< 1.0*. Все такие числа равномерно распределены, но большинство приходится округлять вниз до ближайшего представимого числа с плавающей запятой Python. (Значение 2⁻¹⁰⁷⁴ – это наименьшее положительное ненормализованное число с плавающей запятой, равное `math.ulp(0.0)`.)481482> **См. также**483>484> [Generating Pseudo-random Floating-Point Values](https://python-all.ru/3.12/library/random.html) – статья Аллена Б. Дауни, описывающая способы генерации более мелкозернистых чисел с плавающей запятой, чем обычно генерирует [`random()`](https://python-all.ru/3.12/library/random.html#random.random).485