开头提示:本文章为作者原创,若有雷同,纯属巧合,谢谢!
注:组合数公式由于作者太懒而没有打出,因此带来的不便作者深表歉意!
【问题描述】
组合数表示的是从n个物品中选出m个物品的方案数。举个例子,从(1,2,3)三个物品中选择两个物品可以有(1,2),(1,3),(2,3)这三种选择方法。根据组合数的定义,我们可以给出计算组合数的一般公式:
(此处应有组合数公式)
其中n!=1*2*3*...*n。
小葱想知道如果给定n,m和 k,对于所有的0≤i≤n,0≤j≤min(i,m)有多少对(i,j)满足是k的倍数。
【输入格式】
从文件problem.in中读入数据。
第一行有两个整数t,k,其中t代表该测试点总共有多少组测试数据,k的意义见【问题描述】。
接下来t行每行两个整数n,m,其中n,m的意义见【问题描述】
【输出格式】
输出到文件problem.out中。
t行,每行一个整数代表所有的0≤i≤n,0≤j≤min(i,m)中有多少对(i,j)满足是k的倍数。
【输入输出样例1】
problem.in | problem.out |
1 2 3 3 | 1 |
【输入输出样例1说明】
在所有可能的情况中,只有(注:此处有组合数C(1)(2),1在上,2在下)=2是2的倍数。
【输入输出样例2】
problem.in | problem.out |
2 5 4 5 6 7 | 0 7 |
【子任务】
测试点 | n | m | k | t |
1 | ≤3 | ≤3 | =2 | =1 |
2 | =3 | ≤104 | ||
3 | ≤7 | ≤7 | =4 | =1 |
4 | =5 | ≤104 | ||
5 | ≤10 | ≤10 | =6 | =1 |
6 | =7 | ≤104 | ||
7 | ≤20 | ≤100 | =8 | =1 |
8 | =9 | ≤104 | ||
9 | ≤25 | ≤2000 | =10 | =1 |
10 | =11 | ≤104 | ||
11 | ≤60 | ≤20 | =12 | =1 |
12 | =13 | ≤104 | ||
13 | ≤100 | ≤25 | =14 | =1 |
14 | =15 | ≤104 | ||
15 | ≤60 | =16 | =1 | |
16 | =17 | ≤104 | ||
17 | ≤2000 | ≤100 | =18 | =1 |
18 | =19 | ≤104 | ||
19 | ≤2000 | =20 | =1 | |
20 | =21 | ≤104 |
这题还是比较简单的一道题,只要别上组合数公式的套就行了,貌似上了套只有20-30的分,所以这题可以用杨辉三角来解决,不过一定要先取余,不然很容易爆long long int,然后根据余数是否为0来判断是否是k的倍数,作者本人是用了另一个数组来记录答案,如果是k的倍数则满足该行该列的答案为上一行行该列的答案+1,否则则为上一行行该列的答案。例子如下:
杨辉三角(k=2时取余后) | 答案 |
1 1 1 1 0 1 1 1 1 1 1 0 0 0 1 1 1 0 0 1 1 | 0 0 0 0 1 0 0 1 0 0 0 2 1 1 0 0 2 2 2 0 0 |
如果要求到第3行第2列时的答案,则直接输出第3行前三列的和0+1+0=1即可(首行为第0行,首列为第0列)。
其实虽说是min(i,m),但是实际上用分列求和的方法并不用考虑i的值,因为在i<m时到i自动换到下一行,当i<m时三角并不会到m,当i>m时直接到m停就可以,所以直接<m就可以。
代码如下:
//noip2016提高组Day2 T1 (problem)
#include<iostream>
#include<fstream>
usingnamespace std;
ifstreamfin("problem.in");
ofstreamfout("problem.out");
structC//组合数的两个数据
{
int n;
int m;
};
intyh[2001][2001];
intans[2001][2001];
intmain()
{
memset(ans,0,sizeof(ans));//数组初值
int t,k,i=0,j,p;
fin>>t>>k;
C l[t];
for(i;i<t;i++)
fin>>l[i].n>>l[i].m;
for(i=0;i<2001;i++)//杨辉三角边处理
yh[i][0]=yh[i][i]=1%k;
for(i=2;i<2001;i++)//杨辉三角内部处理
for(j=1;j<i;j++)
yh[i][j]=(yh[i-1][j]%k+yh[i-1][j-1]%k)%k;
for(i=0;i<2001;i++)//判断是否为k的倍数
for(j=0;j<i;j++)
if(!yh[i][j])
ans[i][j]=ans[i-1][j]+1;//是,到当前行答案+1
else
ans[i][j]=ans[i-1][j];//不是,到当前行答案不变
for(i=0;i<t;i++)
{
p=0;
for(j=0;j<=l[i].m;j++) //到l[i].n行,各列求和
p+=ans[l[i].n][j];
fout<<p<<endl;
}
return 0;
}