C++ 随机生成 RxC 列联表及源码实现

3.3 本项目采用的核心思想

👉 逐单元格随机分配 + 剩余量约束

核心原则：

1. 总剩余数量逐步减少
2. 每个单元格分配一个不超过剩余量的随机值
3. 最后一个单元格自动确定

这种方法具备以下优势：

• 实现简单
• 不会出现非法解
• 非常适合教学
• 易于扩展为'固定行和 / 列和'版本

实现思路

4.1 算法总体流程

输入 R, C, N
初始化 R×C 矩阵为 0
remaining = N
for i in [0, R-1]:
    for j in [0, C-1]:
        if 是最后一个单元格:
            table[i][j] = remaining
        else:
            在 [0, remaining] 中随机取值
            table[i][j] = value
            remaining -= value
返回 table

4.2 为什么算法一定正确？

1. 非负性：随机值下界为 0
2. 整数性：整型随机数
3. 总和守恒：每一步都维护 remaining
4. 终止必然性：最后一个格子确定值

4.3 随机性的说明

该方法并不保证'所有可能表等概率'，但在以下场景完全足够：

• 随机测试
• Monte Carlo 仿真
• 教学演示
• 压力测试

📌 若需要严格等概率采样，则必须引入更复杂的组合计数或 MCMC 方法，属于进阶内容。

完整实现代码

头文件 (contingency_table.h)

/******************************************************
 * File: contingency_table.h
 ******************************************************/
#ifndef CONTINGENCY_TABLE_H
#define CONTINGENCY_TABLE_H

#include <vector>

/*
 * 随机生成一个 R×C 列联表
 * 所有元素为非负整数
 * 所有元素之和等于 N
 */
std::vector<std::vector<int>> generateContingencyTable(int R, int C, int N);

#endif

实现文件 (contingency_table.cpp)

/******************************************************
 * File: contingency_table.cpp
 ******************************************************/
#include "contingency_table.h"
#include <random>
#include <stdexcept>

std::vector<std::vector<int>> generateContingencyTable(int R, int C, int N) {
    if (R <= 0 || C <= 0 || N < 0) throw std::invalid_argument("Invalid parameters");
    
    std::vector<std::vector<int>> table(R, std::vector<int>(C, 0));
    std::random_device rd;
    std::mt19937 gen(rd());
    int remaining = N;
    
    for (int i = 0; i < R; ++i) {
        for (int j = 0; j < C; ++j) {
            /* 最后一个单元格，剩余量全部填入 */
            if (i == R - 1 && j == C - 1) {
                table[i][j] = remaining;
            } else {
                std::uniform_int_distribution<> dist(0, remaining);
                int value = dist(gen);
                table[i][j] = value;
                remaining -= value;
            }
        }
    }
    return table;
}

主程序 (main.cpp)

/******************************************************
 * File: main.cpp
 ******************************************************/
#include <iostream>
#include "contingency_table.h"

int main() {
    int R = 3;
    int C = 4;
    int N = 20;
    
    auto table = generateContingencyTable(R, C, N);
    
    std::cout << "Generated contingency table:\n";
    for (const auto& row : table) {
        for (int v : row) std::cout << v << " ";
        std::cout << "\n";
    }
    
    return 0;
}

代码解读

6.1 generateContingencyTable

• 核心生成函数
• 逐单元格随机分配整数
• 通过 remaining 变量保证总和约束
• 最后一个单元格自动补齐

6.2 随机数生成部分

• 使用 std::mt19937 作为高质量伪随机数引擎
• uniform_int_distribution 确保整数均匀分布
• 每次运行生成不同表（可扩展为固定种子）

总结

通过本项目，我们完成了：

1. 列联表的数学建模理解
2. 随机生成问题的工程化拆解
3. 一个稳定、可复用的 C++ 实现
4. 为卡方检验、Monte Carlo 模拟打下基础

该实现具备：

• 清晰逻辑
• 教学友好
• 易于扩展
• 工程可用

常见问题

Q1：生成结果是'均匀分布'的吗？

A：不是严格组合意义下的均匀，但足够随机用于仿真和测试。

Q2：如何固定行和或列和？

A：可先对每一行（或列）分配剩余量，属于自然扩展。

Q3：是否可以生成稀疏表？

A：可以通过限制随机上界或引入零概率控制。

扩展方向

1. 固定行和 / 列和的列联表生成
2. 基于多项分布的概率生成
3. 等概率采样（Stars and Bars）
4. Fisher 精确检验专用生成器
5. 并行 Monte Carlo 仿真