NOIP2016组合数问题 杨辉三角形+组合数

于 2017-10-18 21:44:53 发布
阅读量826 收藏
点赞数
CC 4.0 BY-SA版权
分类专栏： ------数论-------
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/LWD_D/article/details/78278326
本文介绍了一种利用组合数学原理进行快速查询的方法。通过预处理C(n,m)的值，并结合二维前缀和技巧，实现了O(1)的时间复杂度查询特定条件下的组合数。文章提供了完整的C++代码实现。

摘要生成于 C知道 ，由 DeepSeek-R1 满血版支持， 前往体验 >

【还好不是在考场上遇到这道题】
我会说在这之前我都不记得有这么一个公式？总之蒟蒻的数论真的太差了……

C(n,m)=C(n-1,m)+C(n-1,m-1)

知道这个公式后一切都变的简单了。先用O(n*m)预处理出c[i][j]（在代码里写的是f[i][j]) 的值，然后用二维前缀和预处理出答案，在查询时就可以O(1)回答啦！

#include<cstdio>
#include<cstring> 
using namespace std;
int T,k,f[2000+5][2000+5],sum[2005][2005];
int main()
{
    memset(f,-1,sizeof(f));
    scanf("%d%d",&T,&k);
    for(int i=1;i<=2000;i++){
        f[i][1]=i%k;f[i][i]=1;
    }
    for(int i=1;i<=2000;i++)
    for(int j=1;j<=i;j++){
        if(f[i][j]!=-1) continue;
        f[i][j]=(f[i-1][j]+f[i-1][j-1])%k;
    }
    for(int i=1;i<=2000;i++)
    for(int j=1;j<=2000;j++)
        sum[i][j]=sum[i][j-1]+sum[i-1][j]-sum[i-1][j-1]+(f[i][j]==0);

    while(T--){
        int n,m; 
        scanf("%d%d",&n,&m);
        printf("%d\n",sum[n][m]);
    }
    return 0;
}