上述问题中，我对01变量y_kr再进行补充，这个变量表示的是如果工件k在r位置上加工则等于1，否则等于0，对于这个变量也有了一些约束，如何把约束加到上述代码中呢，请重新编写一下，谢谢

好的，以下是重新编写后的示例代码：

import numpy as np

# 定义字典S
S = {
    ('job1', 'position1'): 0.2,
    ('job1', 'position2'): 0.5,
    # ...
}

# 定义参数
F = 0.5
a = -0.5

# 定义计算函数
def calculate_p_kr(p_k, r):
    return p_k * (F + (1 - F) * r ** a)

# 定义约束
constraints = {
    ('job1', 'position1'): [2, 3, 4],
    ('job1', 'position2'): [1, 3, 5],
    # ...
}

# 生成列表
u_jk_list = []
for j in range(1, 11):
    for k in range(1, 11):
        s_jk = S.get(('job{}'.format(j), 'position{}'.format(k)), 0)
        p_kr = calculate_p_kr(0.5, k)
        y_kr = 1 if k in constraints.get(('job{}'.format(j), 'position{}'.format(k)), []) else 0
        u_jk = s_jk + p_kr * y_kr + (1 - y_kr) * 0.1
        u_jk_list.append(u_jk)

# 生成矩阵
u_jk_matrix = np.array(u_jk_list).reshape(10, 10)
print(u_jk_matrix)

在这个示例代码中，我们添加了一个约束字典`constraints`，其中每个键值对表示了一个位置上允许加工的工件编号列表。在计算y_kr的时候，我们根据工件编号是否在对应位置上的约束列表中来设置y_kr的值。如果工件编号在约束列表中，则y_kr为1，否则为0。最后，我们使用numpy库将这个列表转换为一个10*10的矩阵，并且将这个矩阵打印出来。