Noip 2016 组合数问题 - 矩阵前缀和

最新推荐文章于 2025-10-25 23:50:14 发布

原创最新推荐文章于 2025-10-25 23:50:14 发布 · 543 阅读

0 ·

CC 4.0 BY-SA版权

NOIP 同时被 3 个专栏收录

169 篇文章

订阅专栏

模拟

22 篇文章

订阅专栏

数学

18 篇文章

订阅专栏

本文介绍了一种利用杨辉三角特性结合矩阵前缀和的方法来解决特定数学问题的算法实现。通过预计算杨辉三角并进行模运算，避免了直接计算组合数时可能遇到的大数问题。此外，通过构建辅助数组并使用矩阵前缀和技术快速定位符合条件的元素。

C n m = C n - 1 m - 1 + C n m - 1

$C_{m}^{n} = C_{m-1}^{n-1}+C_{m-1}^{n}$
杨辉三角，只不过需要注意

C1m=mCm1=m $C_m^1=m$
矩阵前缀和即可
注意%完后那些符合要求的数值上为0，但是全局变量初始化全为0，为了避免麻烦用另一个数组把那些算出来为0的全标记为1，那些不在计算范围之内的还有不合要求的都设为0，然后求这个数组的矩阵前缀和，注意数组不要访问越界。。。
一开始我循环都访问到2005了答案就错了

#include <cstdio>
#include <iostream>
#include <algorithm>
using namespace std;
#define debug(x) cout<<#x<<"="<<x<<endl;
int t,n,m,k,ans,san[2005][2005],f[2005][2005];
int main() {
    cin >> t >> k;
    san[1][0] = 1;
    san[1][1] = 1;
    for(int i=2; i<=2002; i++){
        san[i][0] = 1;
        for(int j=1; j <= i; j++) {
            san[i][j] = (san[i-1][j] + san[i-1][j-1]) % k;
            if(!san[i][j]) f[i][j] = 1;
        }
    }
    for(int i=1; i<=2002; i++) {
        for(int j=1; j<=2002; j++) {
            f[i][j] = f[i][j] + f[i-1][j] + f[i][j-1] - f[i-1][j-1];
        }
    }
    for(int o=1; o<=t; o++) {
        cin >> n >> m;
        cout << f[n][m] << endl;
    }
    return 0;
}