ms_toollib 1.5.16

use ms_toollib::{
    cal_all_solution, cal_probability_cells_is_op, cal_probability_cells_not_mine,
    cal_probability_csp, cal_probability_enum, cal_probability_onboard, mark_board,
};

#[test]
fn board_1_works() {
    let board = vec![
        vec![10, 10, 10, 10, 10, 10, 10, 10],
        vec![10, 10, 10, 10, 10, 10, 10, 10],
        vec![10, 10, 10, 1, 10, 10, 10, 10],
        vec![10, 10, 10, 2, 10, 10, 10, 10],
        vec![10, 10, 10, 10, 10, 10, 10, 10],
        vec![10, 10, 10, 10, 10, 10, 8, 10],
        vec![10, 10, 10, 10, 10, 10, 10, 10],
        vec![10, 10, 10, 10, 10, 10, 10, 10],
    ];
    let ans = cal_probability_csp(&board, 10.0).unwrap();
    println!("{:?}", ans);
    assert_eq!(ans.2, [10, 10, 54]);
    let ans = cal_probability_enum(&board, 10.0).unwrap();
    println!("{:?}", ans);
    assert_eq!(ans.2, [10, 10, 54]);
}

#[test]
fn board_2_works() {
    let mut board = vec![
        vec![10, 10, 10, 10, 10, 10, 10, 10],
        vec![10, 10, 10, 10, 10, 10, 10, 10],
        vec![10, 10, 10, 10, 10, 10, 10, 10],
        vec![10, 10, 10, 10, 10, 10, 10, 10],
        vec![10, 10, 10, 10, 10, 10, 10, 10],
        vec![10, 10, 10, 10, 10, 10, 10, 10],
        vec![10, 10, 10, 10, 10, 10, 10, 10],
        vec![10, 10, 10, 10, 10, 10, 10, 10],
    ];
    let ans = cal_probability_csp(&board, 10.0).unwrap();
    println!("{:?}", ans);
    assert_eq!(ans.2, [0, 10, 64]);
    let ans = cal_probability_enum(&board, 10.0).unwrap();
    println!("{:?}", ans);
    assert_eq!(ans.2, [0, 10, 64]);
}

#[test]
fn board_3_works() {
    let mut board = vec![
        vec![10, 10, 10, 10, 10, 10, 10, 10],
        vec![10, 10, 10, 10, 10, 10, 10, 10],
        vec![10, 10, 10, 10, 10, 10, 10, 10],
        vec![10, 10, 10, 2, 10, 10, 10, 10],
        vec![10, 10, 10, 10, 10, 10, 10, 10],
        vec![10, 10, 10, 10, 10, 10, 10, 10],
        vec![10, 10, 10, 10, 10, 10, 10, 10],
        vec![10, 10, 10, 10, 10, 10, 10, 10],
    ];
    let ans = cal_probability_csp(&board, 10.0).unwrap();
    println!("{:?}", ans);
    assert_eq!(ans.2, [2, 10, 57]);
    let ans = cal_probability_csp(&board, 0.0).unwrap();
    println!("{:?}", ans);
    assert_eq!(ans.2, [2, 2, 57]);
    let ans = cal_probability_csp(&board, 60.0).unwrap();
    println!("{:?}", ans);
    assert_eq!(ans.2, [2, 57, 57]);
    let ans = cal_probability_enum(&board, 10.0).unwrap();
    println!("{:?}", ans);
    assert_eq!(ans.2, [2, 10, 57]);
}

#[test]
fn board_4_works() {
    let mut board = vec![
        vec![1, 10, 10, 10, 10, 10, 10, 10],
        vec![1, 10, 10, 10, 10, 10, 10, 10],
        vec![2, 10, 10, 10, 10, 10, 10, 10],
        vec![10, 10, 10, 10, 10, 10, 10, 10],
        vec![10, 10, 10, 10, 10, 10, 10, 10],
        vec![10, 10, 10, 10, 10, 10, 10, 10],
        vec![10, 10, 10, 10, 10, 10, 10, 10],
        vec![10, 10, 10, 10, 10, 10, 10, 10],
    ];
    let ans = cal_probability_csp(&board, 10.0).unwrap();
    println!("{:?}", ans);
    assert_eq!(ans.2, [2, 10, 59]);
    assert_eq!(ans.1, 0.14150943396226415);
    for (cell, p) in ans.0 {
        if cell == (0, 1) {
            assert_eq!(p, 0.07547169811320754)
        }
        if cell == (1, 1) {
            assert_eq!(p, 0.9245283018867925)
        }
        if cell == (2, 1) {
            assert_eq!(p, 0.0)
        }
        if cell == (3, 0) {
            assert_eq!(p, 0.5377358490566038)
        }
        if cell == (3, 1) {
            assert_eq!(p, 0.5377358490566038)
        }
    }
    let ans = cal_probability_enum(&board, 10.0).unwrap();
    println!("{:?}", ans);
    assert_eq!(ans.2, [2, 10, 59]);
    assert_eq!(ans.1, 0.14150943396226415);
    for (cell, p) in ans.0 {
        if cell == (0, 1) {
            assert_eq!(p, 0.07547169811320754)
        }
        if cell == (1, 1) {
            assert_eq!(p, 0.9245283018867925)
        }
        if cell == (2, 1) {
            assert_eq!(p, 0.0)
        }
        if cell == (3, 0) {
            assert_eq!(p, 0.5377358490566039)
        }
        if cell == (3, 1) {
            assert_eq!(p, 0.5377358490566039)
        }
    }
}

#[test]
fn board_5_works() {
    let mut board = vec![
        vec![10, 10, 10, 10, 10, 10, 10, 10],
        vec![10, 10, 1, 10, 10, 6, 10, 10],
        vec![10, 10, 2, 10, 10, 10, 1, 10],
        vec![10, 10, 3, 10, 10, 10, 10, 10],
        vec![10, 10, 10, 10, 10, 10, 10, 10],
        vec![10, 10, 10, 10, 5, 10, 10, 10],
        vec![10, 2, 10, 10, 10, 10, 1, 10],
        vec![10, 10, 10, 10, 10, 10, 10, 10],
    ];
    let _ = mark_board(&mut board, false);
    let ans = cal_probability_csp(&board, 29.0).unwrap();
    println!("{:?}", ans);
    assert_eq!(ans.2, [15, 29, 29]);
    let ans = cal_probability_enum(&board, 29.0).unwrap();
    println!("{:?}", ans);
    assert_eq!(ans.2, [15, 29, 29]);

    let ans = cal_probability_onboard(&board, 29.0);
    println!("{:?}", board);
    println!("{:?}", ans);
}

#[test]
fn board_6_works() {
    // 全10局面，传入比例0.2 → 钳位为12雷 → fallback
    let board = vec![
        vec![10, 10, 10, 10, 10, 10, 10, 10],
        vec![10, 10, 10, 10, 10, 10, 10, 10],
        vec![10, 10, 10, 10, 10, 10, 10, 10],
        vec![10, 10, 10, 10, 10, 10, 10, 10],
        vec![10, 10, 10, 10, 10, 10, 10, 10],
        vec![10, 10, 10, 10, 10, 10, 10, 10],
        vec![10, 10, 10, 10, 10, 10, 10, 10],
        vec![10, 10, 10, 10, 10, 10, 10, 10],
    ];
    let ans = cal_probability_csp(&board, 0.2).unwrap();
    println!("{:?}", ans);
    assert!(ans.0.is_empty());
    assert_eq!(ans.2, [0, 12, 64]);
    let ans = cal_probability_enum(&board, 0.2).unwrap();
    println!("{:?}", ans);
    assert_eq!(ans.2, [0, 12, 64]);
}

#[test]
fn board_7_works() {
    // 全10局面，雷数超上限 → 自动钳位到total_unopened
    let board = vec![
        vec![10, 10, 10, 10, 10, 10, 10, 10],
        vec![10, 10, 10, 10, 10, 10, 10, 10],
        vec![10, 10, 10, 10, 10, 10, 10, 10],
        vec![10, 10, 10, 10, 10, 10, 10, 10],
        vec![10, 10, 10, 10, 10, 10, 10, 10],
        vec![10, 10, 10, 10, 10, 10, 10, 10],
        vec![10, 10, 10, 10, 10, 10, 10, 10],
        vec![10, 10, 10, 10, 10, 10, 10, 10],
    ];
    let ans = cal_probability_csp(&board, 100.0).unwrap();
    println!("{:?}", ans);
    assert_eq!(ans.2, [0, 64, 64]);
    let ans = cal_probability_enum(&board, 100.0).unwrap();
    println!("{:?}", ans);
    assert_eq!(ans.2, [0, 64, 64]);
}

#[test]
fn board_8_works() {
    // 有约束局面，雷数不足下限 → 自动钳位到mine_min
    let board = vec![
        vec![10, 10, 10, 10, 10, 10, 10, 10],
        vec![10, 10, 10, 10, 10, 10, 10, 10],
        vec![10, 10, 10, 10, 10, 10, 10, 10],
        vec![10, 10, 10, 2, 10, 10, 10, 10],
        vec![10, 10, 10, 10, 10, 10, 10, 10],
        vec![10, 10, 10, 10, 10, 10, 10, 10],
        vec![10, 10, 10, 10, 10, 10, 10, 10],
        vec![10, 10, 10, 10, 10, 10, 10, 10],
    ];
    let ans = cal_probability_csp(&board, 1.0).unwrap();
    println!("{:?}", ans);
    assert_eq!(ans.2, [2, 2, 57]);
    let ans = cal_probability_csp(&board, 100.0).unwrap();
    println!("{:?}", ans);
    assert_eq!(ans.2, [2, 57, 57]);
    let ans = cal_probability_enum(&board, 1.0).unwrap();
    println!("{:?}", ans);
    assert_eq!(ans.2, [2, 2, 57]);
    let ans = cal_probability_enum(&board, 100.0).unwrap();
    println!("{:?}", ans);
    assert_eq!(ans.2, [2, 57, 57]);
}

#[test]
fn board_9_works() {
    // 有约束局面，传入比例 → 自动换算为雷数
    let board = vec![
        vec![10, 10, 10, 10, 10, 10, 10, 10],
        vec![10, 10, 10, 10, 10, 10, 10, 10],
        vec![10, 10, 10, 10, 10, 10, 10, 10],
        vec![10, 10, 10, 2, 10, 10, 10, 10],
        vec![10, 10, 10, 10, 10, 10, 10, 10],
        vec![10, 10, 10, 10, 10, 10, 10, 10],
        vec![10, 10, 10, 10, 10, 10, 10, 10],
        vec![10, 10, 10, 10, 10, 10, 10, 10],
    ];
    // 0.2 → 12雷
    let ans = cal_probability_csp(&board, 0.2).unwrap();
    println!("{:?}", ans);
    assert_eq!(ans.2, [2, 12, 57]);
    let ans = cal_probability_enum(&board, 0.2).unwrap();
    println!("{:?}", ans);
    assert_eq!(ans.2, [2, 12, 57]);
}

#[test]
fn board_10_works() {
    // 11标记（标雷）→ CSP当未知处理，不参与witness约束调整
    let mut board = vec![
        vec![10, 10, 10, 10, 10, 10, 10, 10],
        vec![10, 10, 10, 10, 10, 10, 10, 10],
        vec![10, 10, 10, 10, 10, 10, 10, 10],
        vec![10, 10, 10, 2, 10, 10, 10, 10],
        vec![10, 10, 10, 10, 10, 10, 10, 10],
        vec![10, 10, 10, 10, 10, 10, 10, 10],
        vec![10, 10, 10, 10, 10, 10, 10, 10],
        vec![10, 10, 10, 10, 10, 10, 10, 10],
    ];
    // 把一个邻居标雷
    board[3][2] = 11;
    let ans = cal_probability_csp(&board, 10.0).unwrap();
    println!("{:?}", ans);
    // 11视为已知雷→排除(3,2)，其余7格p=1/7
    assert!(!ans.0.iter().any(|&(pos, _)| pos == (3, 2)));
    assert_eq!(ans.0.len(), 7);
    for (_, p) in &ans.0 {
        assert!((p - 1.0 / 7.0).abs() < 1e-10);
    }
    assert_eq!(ans.2, [2, 10, 57]);
    let ans = cal_probability_enum(&board, 10.0).unwrap();
    println!("{:?}", ans);
    // ENUM: 11计入is_minenum，不在枚举结果中
    for (pos, p) in &ans.0 {
        assert_ne!(pos, &(3, 2));
    }
    assert_eq!(ans.0.len(), 7);
    assert_eq!(ans.2, [2, 10, 57]);
}

#[test]
fn board_11_works() {
    // 12标记（已知安全）→ CSP排除该格，不加入witness邻接
    let mut board = vec![
        vec![10, 10, 10, 10, 10, 10, 10, 10],
        vec![10, 10, 10, 10, 10, 10, 10, 10],
        vec![10, 10, 10, 10, 10, 10, 10, 10],
        vec![10, 10, 10, 2, 10, 10, 10, 10],
        vec![10, 10, 10, 10, 10, 10, 10, 10],
        vec![10, 10, 10, 10, 10, 10, 10, 10],
        vec![10, 10, 10, 10, 10, 10, 10, 10],
        vec![10, 10, 10, 10, 10, 10, 10, 10],
    ];
    // 把一个邻居标安全
    board[3][2] = 12;
    let ans = cal_probability_csp(&board, 10.0).unwrap();
    println!("{:?}", ans);
    // 12格不在结果中
    assert!(!ans.0.iter().any(|&(pos, _)| pos == (3, 2)));
    // 剩余7格，2雷 → 2/7
    assert_eq!(ans.0.len(), 7);
    for (_, p) in &ans.0 {
        assert!((p - 2.0 / 7.0).abs() < 1e-10);
    }
    assert_eq!(ans.2, [2, 10, 57]);
    let ans = cal_probability_enum(&board, 10.0).unwrap();
    println!("{:?}", ans);
    assert!(!ans.0.iter().any(|&(pos, _)| pos == (3, 2)));
    assert_eq!(ans.0.len(), 7);
    for (_, p) in &ans.0 {
        assert!((p - 2.0 / 7.0).abs() < 1e-10);
    }
    assert_eq!(ans.2, [2, 10, 57]);
}

#[test]
fn board_16_works() {
    let mut board = vec![
        vec![10, 10, 10, 10, 10, 10, 10, 10],
        vec![10, 10, 10, 10, 10, 10, 10, 10],
        vec![10, 10, 10, 10, 10, 10, 10, 10],
        vec![10, 10, 10, 2, 10, 10, 10, 10],
        vec![10, 10, 10, 10, 10, 10, 10, 10],
        vec![10, 10, 10, 10, 10, 10, 10, 10],
        vec![10, 10, 10, 10, 10, 10, 10, 10],
        vec![10, 10, 10, 10, 10, 10, 10, 10],
    ];
    let chording_cells = vec![(2, 2), (3, 2), (4, 2), (2, 4), (3, 4), (4, 4)];
    let p = cal_probability_cells_not_mine(&board, 10 as f64, &chording_cells);
    assert_eq!(p * 28.0, 1.0000000000000002);
}

#[test]
fn board_12_works() {
    // 全10+12标记 → fallback应排除12格
    let mut board = vec![
        vec![10, 10, 10, 10, 10, 10, 10, 10],
        vec![10, 10, 10, 10, 10, 10, 10, 10],
        vec![10, 10, 10, 10, 10, 10, 10, 10],
        vec![10, 10, 10, 10, 10, 10, 10, 10],
        vec![10, 10, 10, 10, 10, 10, 10, 10],
        vec![10, 10, 10, 10, 10, 10, 10, 10],
        vec![10, 10, 10, 10, 10, 10, 10, 10],
        vec![10, 10, 10, 10, 10, 10, 10, 10],
    ];
    board[0][0] = 12;
    board[7][7] = 12;
    let ans = cal_probability_csp(&board, 10.0).unwrap();
    println!("{:?}", ans);
    // 12格不在结果中
    assert!(!ans.0.iter().any(|&(pos, _)| pos == (0, 0)));
    assert!(!ans.0.iter().any(|&(pos, _)| pos == (7, 7)));
    // 剩余62格，10雷 → 10/62
    assert_eq!(ans.0.len(), 62);
    for (_, p) in &ans.0 {
        assert!((p - 10.0 / 62.0).abs() < 1e-10);
    }
    assert_eq!(ans.2, [0, 10, 62]);
    let ans = cal_probability_enum(&board, 10.0).unwrap();
    println!("{:?}", ans);
    // ENUM: 全10无witness → ans.0为空，概率通过p_unknow返回
    assert_eq!(ans.0.len(), 0);
    assert!((ans.1 - 10.0 / 62.0).abs() < 1e-10);
    assert_eq!(ans.2, [0, 10, 62]);
}

#[test]
fn board_13_works() {
    // 11标记（标雷）→ CSP当未知处理，不参与witness约束调整
    let mut board = vec![
        vec![10, 10, 10, 10, 10, 10, 10, 10],
        vec![10, 10, 10, 10, 10, 10, 10, 10],
        vec![10, 10, 10, 10, 10, 10, 10, 10],
        vec![10, 10, 10, 2, 10, 10, 10, 10],
        vec![10, 10, 10, 3, 10, 10, 10, 10],
        vec![10, 10, 10, 10, 10, 10, 10, 10],
        vec![10, 10, 10, 10, 10, 10, 10, 10],
        vec![10, 10, 10, 10, 10, 10, 10, 10],
    ];
    // 把一个邻居标雷
    board[3][2] = 11;
    board[0][0] = 11;
    board[3][1] = 11;
    board[6][6] = 12;
    board[6][6] = 12;
    board[7][0] = 12;
    let ans = cal_probability_csp(&board, 10.0).unwrap();
    println!("{:?}", ans);
    // 11视为已知雷→排除(3,2)，其余7格p=1/7
    let ans = cal_probability_enum(&board, 10.0).unwrap();
    println!("{:?}", ans);
}

// CSP 概率 > 1 的 bug 复现
#[test]
fn csp_prob_gt_1() {
    let board: Vec<Vec<i32>> = vec![
        vec![
            2, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10,
            10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10,
            10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 11, 1, 0, 1, 1,
        ],
        vec![
            10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10,
            10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10,
            10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 4, 3, 1, 3, 11,
        ],
        vec![
            10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10,
            10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10,
            10, 10, 10, 10, 10, 10, 10, 10, 3, 10, 10, 11, 4, 11, 4, 11,
        ],
        vec![
            10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10,
            10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10,
            10, 10, 10, 10, 10, 10, 5, 11, 2, 1, 3, 11, 4, 11, 3, 1,
        ],
        vec![
            10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10,
            10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10,
            10, 10, 10, 10, 10, 10, 10, 2, 1, 0, 1, 1, 2, 1, 1, 0,
        ],
        vec![
            10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10,
            10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10,
            10, 10, 10, 10, 10, 10, 10, 2, 1, 0, 0, 0, 1, 1, 1, 0,
        ],
        vec![
            10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10,
            10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10,
            10, 10, 10, 10, 10, 2, 3, 11, 2, 0, 0, 0, 1, 11, 2, 1,
        ],
        vec![
            10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10,
            10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10,
            10, 10, 10, 10, 10, 1, 2, 11, 4, 3, 2, 2, 2, 3, 10, 10,
        ],
        vec![
            10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10,
            10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10,
            10, 10, 10, 10, 10, 2, 2, 3, 11, 11, 11, 4, 11, 3, 10, 10,
        ],
        vec![
            10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10,
            10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10,
            10, 10, 10, 10, 10, 10, 10, 2, 10, 4, 10, 10, 10, 10, 10, 10,
        ],
        vec![
            10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10,
            10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10,
            10, 10, 10, 10, 10, 10, 3, 3, 2, 10, 10, 10, 10, 10, 10, 10,
        ],
        vec![
            10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10,
            10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10,
            10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10,
        ],
        vec![
            10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10,
            10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10,
            10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10,
        ],
        vec![
            10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10,
            10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10,
            10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10,
        ],
        vec![
            10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10,
            10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10,
            10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10,
        ],
        vec![
            10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10,
            10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10,
            10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10,
        ],
        vec![
            10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10,
            10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10,
            10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10,
        ],
        vec![
            10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10,
            10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10,
            10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10,
        ],
        vec![
            10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10,
            10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10,
            10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10,
        ],
        vec![
            10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10,
            10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10,
            10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10,
        ],
        vec![
            10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10,
            10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10,
            10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10,
        ],
        vec![
            10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10,
            10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10,
            10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10,
        ],
        vec![
            10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10,
            10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10,
            10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10,
        ],
        vec![
            10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10,
            10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10,
            10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10,
        ],
        vec![
            10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10,
            10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10,
            10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10,
        ],
        vec![
            10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10,
            10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10,
            10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10,
        ],
        vec![
            10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10,
            10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10,
            10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10,
        ],
        vec![
            10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10,
            10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10,
            10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10,
        ],
        vec![
            10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10,
            10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10,
            10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10,
        ],
        vec![
            10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10,
            10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10,
            10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10,
        ],
        vec![
            10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10,
            10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10,
            10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10,
        ],
        vec![
            10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10,
            10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10,
            10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10,
        ],
        vec![
            10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10,
            10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10,
            10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10,
        ],
        vec![
            10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10,
            10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10,
            10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10,
        ],
        vec![
            10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10,
            10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10,
            10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10,
        ],
    ];
    let ans = cal_probability_csp(&board, 540.0).unwrap();
    for &((r, c), p) in &ans.0 {
        assert!(
            p <= 1.0 + 1e-12,
            "CSP probability > 1 at ({},{}): {}",
            r,
            c,
            p
        );
    }
}