满足条件的随机数的产生及python程序

有三个变量满足

$p_1+p_2+p_3=n$

$a_1<p_1<b_1$

$a_2<p_2<b_2$

$a_3<p_3<b_3$

则思路：

• step 1：

$a_3<p_3<b_3$，

根据 $p_1+p_2+p_3=n$，

$a_3<n-p_1-p_2<b_3$

$n-(b_3+p_1)<p_2<n-(a_3+p_1)$

• step 2:

由于$a_2<p_2<b_2$

综上$a_2\le n-(a_3+p_1)$

而且$b_2\ge n-(b_3+p_1)$

可得

$n-(a_2+a_3)\le p_1\le n-(b_2+b_3)$

结合$a_1<p_1<b_1$

可先随机取满足条件的$p_1$的值

• step 3:

然后再递推回去得到满足条件的$p_2$的值，进而可以根据$p_1+p_2+p_3=1$确定$p_3$的值

---

Python实现：
#%%

import random

def g(x_low, x_up, x_sum, y):
if len(x_low) == 1:
y.append(x_sum)
print(x_sum)
else:
a = max(x_sum-sum(x_up[1:len(x_up)]), x_low[0])
b = min(x_sum-sum(x_low[1:len(x_low)]), x_up[0])
temp = random.uniform(a,b)
y.append(temp)
print(temp)
x_low = x_low[1:len(x_low)]
x_up = x_up[1:len(x_up)]
x_sum = x_sum - temp
g(x_low, x_up, x_sum, y)
#return y

x_low = [0, 0, 0] #即$a_1,a_2,a_3=0$
x_up = [1, 1, 1] #即$b_1,b_2,b_3=1$
x_sum =1 #即$p_1+p_2+p_3=1$
y = []
Y = []
i=0

g(x_low, x_up, x_sum, y)