NOIP2016 提高组Day2 T1 组合数问题

题目描述
组合数表示的是从n个物品中选出m个物品的方案数。举个例子，从（1,2,3) 三个物品中选择两个物品可以有（1,2),(1,3),(2,3)这三种选择方法。根据组合数的定 义，我们可以给出计算组合数的一般公式： 
 
其中n! = 1×2×···×n 
小葱想知道如果给定n,m和k，对于所有的0<=i<= n，0<=j<= min(i,m)有多少对 (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 
数据范围

 

分析

这道题其实是杨辉三角问题，因为最后的答案只要对k取余，所以可以先算出每个C[n,m]对k取余的值，用flag[n,m]储存它能否整除k，最后的答案统计f[1..n,1..m]即可。

CODE

const
  maxn=5000;
var
  f,flag:array[0..maxn,0..maxn]of longint;
  i,j,t,k,n,m:longint;
begin
  readln(t,k);
  for i:=1 to maxn do 
  begin 
    a[i,1]:=i mod k;
    a[i,i]:=1; 
  end;
  for i:=3 to maxn do
    for j:=2 to i-1 do
      a[i,j]:=(a[i-1,j]+a[i-1,j-1]) mod k;
  for i:=1 to maxn do
    for j:=1 to i do
      if a[i,j]=0 then
        flag[i,j]:=1;
        
    fillchar(a,sizeof(a),0);
    
  for i:=1 to maxn do
    for j:=1 to maxn do
      a[i,j]:=a[i-1,j]+a[i,j-1]-a[i-1,j-1]+flag[i,j];
  for i:=1 to t do
  begin
    read(n,m);
    writeln(a[n,m]);
  end;
end.