NOIP2016Day2T1组合数问题(problem)

开头提示：本文章为作者原创，若有雷同，纯属巧合，谢谢！

注：组合数公式由于作者太懒而没有打出，因此带来的不便作者深表歉意！

【问题描述】

组合数表示的是从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;

}