从random随机数看验证码重复数字

本文分享自天翼云开发者社区《从random随机数看验证码重复数字》.作者：王****淋

## 推导：

1. 有6个不同数字

- (0个重复数字: 1+1+1+1+1+1):

$$p(6) = C_{10}^6 * 6! \div 10^6$$

2. 有5个不同数字

- (某个数字重复了2次: 1+1+1+1+2):

$$p(5) = (\frac{C_{10}^5 * C_5^1}{2!*1!*1!*1!*1!} ) * 6! \div 10^6$$

$$ = (\frac{C_{10}^5 * C_5^1}{2!} ) * 6! \div 10^6$$

3. 有4个不同数字

- (某个数字重复了3次: 1+1+1+3，或某2个数字各重复2次: 1+1+2+2)。(省略$1!$):

$$p(4) = (\frac{C_{10}^4 * C_4^1}{3!} + \frac{C_{10}^4 * C_4^2}{2!*2!}) * 6! \div 10^6$$

4. 有3个不同数字

- (某个数字重复了4次: 1+1+4，或某个数字重复3次+某数字重复2次: 1+3+2, 或3个数字哥重复2次: 2+2+2):

$$p(3) = (\frac{C_{10}^3 * C_3^1}{4!} + \frac{C_{10}^3 * C_3^1 * C_2^1}{3!*2!} +\frac{C_{10}^3 * C_3^3}{2!*2!*2!}) * 6! \div 10^6$$

5. 有2个不同数字

- (某个数字重复了5次: 1+5，或某个数字重复4次+某数字重复2次: 4+2, 或2个数字哥重复3次: 3+3): (省略$1!$):

$$p(2) = (\frac{C_{10}^2 * C_2^1}{5!} + \frac{C_{10}^2 * C_2^1}{4!*2!} + \frac{C_{10}^2 * C_2^2}{3!*3!}) * 6! \div 10^6$$

6. 有1个不同数字(6个相同的重复数字):

$$p(1) = \frac{C_{10}^1}{6!} * 6! \div 10^6$$

## 蒙特卡洛

```py

import random

# 6位验证码, 相同数字的概率: 蒙特卡洛

def mtkr(count:int=10, is_show:bool=True, debug:bool = False):

ans = [0 for i in range(6)]

for i in range(count):

array_i = [random.randint(0, 9) for i in range(6)]

n = len(list(set(array_i)))

if debug:

print(n, ": ", array_i)

ans[6-n] += 1

# end_for

if is_show:

print('-' * 30)

print('蒙特卡洛: count:', count, "; \n模拟结果: ")

for i in range(6):

print("unique", 6-i, ": ", ans[i]/count)

return [i/count for i in ans]

ans = mtkr(100000, is_show=True, debug=False)

print(ans)

res_plot(ans)

# 结果：

#模拟结果:

# unique 6 : 0.15148

# unique 5 : 0.45396

# unique 4 : 0.32748

# unique 3 : 0.0643

# unique 2 : 0.00278 u

# nique 1 : 0.0

```