如何随机化相似二进制向量并确保无重复元素组合？

问题分析 #

你需要处理两个等长、0和1数量相同的二进制向量A、B，要求：

1. 二者独立随机化；
2. 对应位置的元素组合(A(i), B(i))无重复（即四种组合(0,0)、(0,1)、(1,0)、(1,1)各出现一次，反映到A+B上就是包含一个0、两个1、一个2）；
3. 不能出现A与B完全相同的情况。

你的原代码存在两个核心问题：

• 循环条件仅检查(0,0)的数量，未覆盖(1,1)的数量限制，也没校验A、B是否完全相同；
• 每次循环同时打乱A和B，满足条件的概率低，容易陷入死循环。

高效解决方案：约束性构造B #

不需要依赖随机循环碰运气，先固定A的随机化结果，再根据A的结构直接构造符合要求的B，逻辑如下：

1. 随机化A；
2. 提取A中0和1的索引；
3. 从A的0索引中选1个、1索引中选1个作为B的0的位置，保证B的0和1数量与原向量一致；
4. 剩余位置设为1，这样自然满足所有组合各出现一次的要求。

代码实现

% 初始化原始向量
A = [1 1 0 0];
B_template = A;

% 第一步：随机化A
A = A(randperm(length(A)));

% 提取A中0和1的索引
zero_indices = find(A == 0);
one_indices = find(A == 1);

% 从A的0、1索引中各选1个，作为B的0的位置
b_zero_pos = [zero_indices(randi(length(zero_indices))), one_indices(randi(length(one_indices)))];

% 构造符合要求的B
B = zeros(size(A));
B(b_zero_pos) = 0;
B(setdiff(1:length(A), b_zero_pos)) = 1;

% 验证结果（可选）
disp('随机化后的A：'); disp(A);
disp('构造后的B：'); disp(B);
disp('A+B的结果：'); disp(A+B);

代码说明 #

• 该方法直接通过逻辑构造满足条件的B，完全避免死循环，效率远高于随机循环；
• 最终的A+B必然包含一个0（对应(0,0)）、两个1（对应(0,1)和(1,0)）、一个2（对应(1,1)），严格符合你对组合无重复的要求；
• 由于B的0位置同时覆盖了A的0和1索引，A与B不可能完全相同。

原代码的修正思路（如果坚持用循环） #

若一定要保留循环逻辑，需修正循环条件，同时优化随机策略（比如只打乱B，固定A），示例代码如下：

A = [1 1 0 0];
B = A;

% 先随机化A
A = A(randperm(length(A)));

% 循环直到满足条件
while true
    % 仅随机化B
    B = B(randperm(length(B)));
    % 检查条件：(0,0)和(1,1)的数量都为1，且A与B不完全相同
    count_00 = sum(A == 0 & B == 0);
    count_11 = sum(A == 1 & B == 1);
    if count_00 == 1 && count_11 == 1 && ~isequal(A, B)
        break;
    end
end

disp('修正循环后的A：'); disp(A);
disp('修正循环后的B：'); disp(B);
disp('A+B的结果：'); disp(A+B);

不过这种方法仍存在随机概率问题，向量越长，满足条件的概率越低，不推荐使用。

内容的提问来源于stack exchange，提问作者Nic