Luogu P2822 组合数问题

https://www.luogu.org/problemnew/show/P2822

题目描述

组合数 C_n^mCnm​ 表示的是从 nn 个物品中选出 mm 个物品的方案数。举个例子，从 (1,2,3)(1,2,3) 三个物品中选择两个物品可以有 (1,2),(1,3),(2,3)(1,2),(1,3),(2,3) 这三种选择方法。根据组合数的定义，我们可以给出计算组合数 C_n^mCnm​ 的一般公式：

C_n^m=\frac{n!}{m!(n-m)!}Cnm​=m!(n−m)!n!​

其中 n!=1\times2\times\cdots\times nn!=1×2×⋯×n；特别地，定义 0!=10!=1。

小葱想知道如果给定 n,mn,m 和 kk，对于所有的 0\leq i\leq n,0\leq j\leq \min \left ( i, m \right )0≤i≤n,0≤j≤min(i,m) 有多少对 (i,j)(i,j) 满足 C_i^jCij​ 是 kk 的倍数。

输入输出格式

输入格式：

 

第一行有两个整数 t,kt,k，其中 tt 代表该测试点总共有多少组测试数据，kk 的意义见问题描述。

接下来 tt 行每行两个整数 n,mn,m，其中 n,mn,m 的意义见问题描述。

 

输出格式：

 

共 tt 行，每行一个整数代表所有的 0\leq i\leq n,0\leq j\leq \min \left ( i, m \right )0≤i≤n,0≤j≤min(i,m) 中有多少对 (i,j)(i,j) 满足 C_i^jCij​ 是 kk 的倍数。

 

输入输出样例

输入样例#1： 复制

1 2
3 3

输出样例#1： 复制

1

输入样例#2： 复制

2 5
4 5
6 7

输出样例#2： 复制

0
7

说明

【样例1说明】

在所有可能的情况中，只有C_2^1 = 2C21​=2是2的倍数。

【子任务】

 

我竟然只有25分

初始化出错了：（，我怎么这么菜

#include<cstdio>
#include<iostream>
using namespace std;
int read()
{
    int ret=0; char ch=getchar();
    while(ch<'0'||ch>'9') ch=getchar();
    while(ch>='0'&&ch<='9')
        ret=(ret<<1)+(ret<<3)+ch-'0',ch=getchar();
    return ret;
}
 
const int N=2005;
int t,p,n,m,c[N][N],a[N][N];
 
int main()
{
    //freopen("problem.in","r",stdin);
    //freopen("problem.out","w",stdout);
    t=read(),p=read();
    for(int i=1;i<=2000;i++) c[0][i]=c[i][i]=1;
    for(int j=2;j<=2000;j++)
        for(int i=1;i<j;i++)
            c[i][j]=(c[i-1][j-1]+c[i][j-1])%p;
    for(int i=1;i<=2000;i++)
        for(int j=1;j<=2000;j++)
            a[i][j]=a[i-1][j]+a[i][j-1]-a[i-1][j-1]+(i<j&&c[i][j]==0);
    
	while(t--)
        n=read(),m=read(),
        printf("%d\n",a[m][n]);
    return 0;       
}