P2822 组合数问题

题目描述

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

​​  C​n​m​​=​m!(n−m)!​​n!​​
其中n! = 1 × 2 × · · · × n//不清楚不要怪我，暂时找不到图。。。

小葱想知道如果给定n,m和k，对于所有的0 <= i <= n，0 <= j <= min(i,m)有多少对 (i,j)满足C_i^jC
​i
​j
​​ 是k的倍数。

输入输出格式

输入格式：
第一行有两个整数t,k，其中t代表该测试点总共有多少组测试数据，k的意义见 【问题描述】。

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

输出格式：
t行，每行一个整数代表答案。

输入输出样例

输入样例#1：
1 2
3 3
输出样例#1：
1

输入样例#2：
2 5
4 5
6 7
输出样例#2：
0
7

这其实是一道纯粹的数论题，然而洛谷上的通过率还不到10%。。。

代码如下（pascal）

const maxn=2001;
var
 f,flag:array[0..maxn,0..maxn]of longint;
 i,j,k,t,n,m:longint;
begin
  readln(t,k);
  for i:=1 to maxn-1 do begin
   f[i,1]:=i mod k; f[i,i]:=1; end;
  for i:=3 to maxn-1 do
   for j:=2 to i-1 do
    f[i,j]:=(f[i-1,j]+f[i-1,j-1]) mod k;//对杨辉三角的维护，防止太大导致爆0~
  for i:=1 to maxn-1 do
   for j:=1  to  i do
     if f[i,j]mod k=0 then flag[i,j]:=1;//也可以不要这个mod，因为上面已经mod了，能被k整除的对其进行flag标记
     fillchar(f,sizeof(f),0);
  for i:=1 to maxn-1 do
    for j:=1 to maxn-1 do
    f[i,j]:=f[i-1,j]+f[i,j-1]-f[i-1,j-1]+flag[i,j];//从f[1,1]开始找能被k整除的数的个数
  for i:=1 to t do begin
   readln(n,m);
   writeln(f[n,m]);
   end;
end.

其实最主要的地方就是找到杨辉三角后进行对数组的维护~只有做了这题你才能发现数学的美妙╮(╯_╰)╭，唉，oi选手就是这么命苦啊！不仅要把编程学好，还要把其他科学好。。。（语数英……）蒟蒻表示十分无奈。