生成满足约束的随机数的方法

第一种问题和方法

有三个变量$x_1，x_2，x_3$
生成满足约束的值：
$$x_1+x_2+x_3=n$$
$$a_1<x_1<b_1,a_2<x_2<b_2,a_3<x_3<b_3$$
思路：
1、$a_3<x_3<b_3$，所以$a_3<n-x_1-x_2<b_3$

2、由1得，$n-x_1-b_3<x_2<n-x_1-a_3$

3、由2得，$n-x_1-b_3>a_2,n-x_1-a_3<b_2$，所以，$x_1\in [n-(b_2+b_3),n-(a_2+a_3)]$

步骤：
1、生成$x_1\in [n-(b_2+b_3),n-(a_2+a_3)]$

2、生成$x_2\in [n-(x_1+b_3),n-(x_1+a_3)]$

3、生成$x_3=n-x_1-x_2$

# python 3.6
import random

# 递归生成随机数
def generate_randval(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
		generate_randval(x_low, x_up, x_sum, y)

if __name__=="__main__":
	x_low = [1, 2, 3]
	x_up = [3, 6, 9]
	x_sum =10  # 在6和18之间
	y = []
	generate_randval(x_low, x_up, x_sum, y)
	

第二种问题和方法

有三个变量$x_1，x_2，x_3$
生成满足约束的值：
$$x_1+x_2+x_3=n$$
$$a_1<x_1<b_1,a_2<x_2<b_2,a_3<x_3<b_3$$
$$a_1=a_2=a_3，b_1=b_2=b_3，c_1=c_2=c_3$$
思路:
1、$x_1=y_1/(y_1+y_2+y_3)\ast n$，所以$a_1<y_1/(y_1+y_2+y_3)\ast n<b_1$

2、若$a_1<y_1<b_1$，所以$0<1/(y_1+y_2+y_3)\ast n<1$

步骤：
1、生成$y_1\in [a_1,b_1]$,$y_2\in [a_2,b_2]$,$y_3\in [a_3,b_3]$

2、如果$y_1+y_2+y_3>n$, 生成$x_1=y_1/(y_1+y_2+y_3)\ast n$，$x_2=y_2/(y_1+y_2+y_3)\ast n$，$x_3=y_3/(y_1+y_2+y_3)\ast n$

3、如果$y_1+y_2+y_3<n$,生成$x_1=y_1\ast n/(y_1+y_2+y_3)$，$x_2=y_2\ast n/(y_1+y_2+y_3)$，$x_3=y_3\ast n/(y_1+y_2+y_3)$