题目描述
组合数表示的是从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
数据范围
组合数表示的是从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.
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.