我相信您可以将其重新表述为整数规划问题。
让:
x_ij = 1 if interval i's left endpoint is at j
0 otherwise
和:
y_ij = 1 if interval i's right endpoint is at j
0 otherwise
最后:
a_k = Total # of 1's in the interval [0, k]
b_k = Total # of 0's in the interval [0, k]
那么下面就相当于你的问题了:
maximize sum_ijk ( a_j * y_ij - a_k - x_ik # Ones inside
- b_j * y_ij - b_k - x_ik # Zeros inside
+ b_j * x_ij - b_k - y_(i-1)k # Zeros outside
- a_j * x_ij - a_k - y_(i-1)k # Ones outside
)
with respect to the constraints
sum_j x_ij = 1 for each i
sum_j y_ij = 1 for each i
0 <= x_ij <= 1 for each i, j
0 <= y_ij <= 1 for each i, j
sum_j * y_ij - j * x_ij > 0 for each i
sum_j * x_(i+1)j - j * y_ij > 0 for each i
由于每个x_ij 和y_ij 都是整数,这是一个整数规划问题。解除这个约束,你就会遇到一个线性规划问题,尽管这种情况下的结果很难解释。
为:
maximize sum_ijk ( a_j * y_ij - a_k - x_ik )
i 的总和是所有区间。每个术语 a_j * y_ij 仅针对 j 的一个值(该区间的右端点)为“on”。与a_j * x_ij 相同。那么差值就是一些a_k - a_r,也就是区间中1的总数。同样,其他三个术语计算正确和错误分类的出现次数。
对于约束:
sum_j x_ij = 1 for each i
sum_j y_ij = 1 for each i
说区间必须各有一个左端点和一个右端点,并且
sum j * y_ij - j * x_ij > 0 for each i
sum j * x_(i+1)j - j * y_ij > 0 for each i
表示左端点必须在右端点的左边,i+1st区间的右端点必须在ith区间的右端点的左边。