随机整数生成答案

【问题标题】：Random integer generating随机整数生成
【发布时间】：2018-10-09 20:28:56
【问题描述】：

我遇到了一个奇怪的问题。也许有人可以指导我查阅相关文献。

所以，在 Python 中，我创建了这个方法，它附加随机整数来设置，直到出现重复值。当生成的整数对于特定的集合不是唯一的时，方法会刹车：

import random

def count_no_repeat(i,j):
    random_set = set()
    while True:
        new_number = random.randint(i,j)
        if new_number in random_set:
            break
        random_set.add(new_number)
    return len(random_set) + 1

然后，我重复这个方法一千次来计算：需要多少步骤才能产生非之前的价值

stats = []
for _ in range(1000):
    stats.append(count_no_repeat(1,n))

n - 整数生成器有上限。

得到了这样的结果：对于 n = 100：

对于 n = 1000：

对于 n = 10000：

对于 n = 100000：

所以，对于这个实验的中位数：

增长相对缓慢；
停留在图上的位置（对于 10000 次实验也是如此）；

谁能提供帮助，并说，为什么会这样？谢谢！

【问题讨论】：

好像是模拟生日攻击。
感谢嵌入 img 的
您的问题是为什么会发生这种情况？
您能否澄清您需要帮助的问题或主题？
这可能在Math 上问得更好，尽管你也可能会在那里得到一个“数学”的答案。

标签： python random statistics distribution

【解决方案1】：

您正在计算广义生日问题的 PDF。基本上都在https://en.wikipedia.org/wiki/Birthday_problem。唯一的问题是 Wiki 页面正在讨论问题的 CDF（请参见此处的第一张图），您正在采样 PDF，p(n, k) - p(n, k-1) 的值。这是您的采样图（蓝色）与 PDF（橙色），如果您需要代码，请告诉我

更新

无论如何，最好将代码放在这里，以免丢失。所有阶乘均作为 Gamma 函数计算，pbar/p 的表达式通过对数完成，避免溢出，因此调用 Gamma 函数的对数 lgamma。

import matplotlib.pyplot as plt
import numpy as np
import math
import random

def pbar(k, n): # as in wiki article, computed via log/exp
    l = 0
    try:
        l = math.lgamma(n + 1)  - math.lgamma(n - k + 1) - k*math.log(n)
    except ValueError:
        l = -50
    return math.exp(l)

def p(k, n):
    return 1.0 - pbar(k, n)

def count_no_repeat(i, j): # original sampling code
    random_set = set()
    while True:
        new_number = random.randint(i,j)
        if new_number in random_set:
            break
        random_set.add(new_number)
    return len(random_set) + 1

# 100 of numbers, 1mln of samples
n = 100
N = 1000000
stats = np.zeros(n+2, dtype = np.float32)
meds  = []

for _ in range(0, N):
    q = count_no_repeat(1, n)
    stats[q] += 1
    meds.append(q)

print(np.median(meds))

stats /= float(N)
x = np.linspace(0, n+1, n+2)

# computing PDF
z = []
for k in x:
    if k == 0:
        z.append(0)
    else:
        z.append(p(k, n) - p(k-1, n))

plt.plot(x, stats, 'o')
plt.plot(x, z)
plt.show()

【讨论】：