Noip2016day2 组合数问题problem

题目描述

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

说明

【样例1说明】

在所有可能的情况中，只有是2的倍数。

【子任务】

思路

这是一道数论题

首先要知道组合数的一般递推公式，它的递推公式和杨辉三角是一样的

c[i][j]=c[i-1][j-1]+c[i-1][j]

（解释：c[i][j]即为从i件物品中选j件的方案数。如果第i件物品不选，方案数就变为c[i-1][j]，如果选第i件物品，方案数就变为c[i-1][j-1]，总方案数就为两种情况的方案数之和）

为了不爆long long，每次求出c[i][j]后先模一下k

为了节约时间，进行二维求和，最后直接查找答案

#include<iostream>
using namespace std;
int n,m,t,k;
int c[][],s[][];
void get_c()
{
    for(int i=;i<=;i++)
    {
        for(int j=;j<=i;j++)
        {
            if(i==&&j==)c[i][j]=%k;
            else
            {
                c[i][j]=(c[i-][j]+c[i-][j-])%k;
            }
        }
    }
}
void get_s()
{
    if(c[][]==)s[][]=;
    for(int i=;i<=;i++)
    {
        for(int j=;j<=i;j++)
        {
            if(i==&&j==)continue;
            else
            {
                if((i==&&j)||(i==j)){
                    s[i][j]=s[i][j-];
                    if(c[i][j]==)s[i][j]++;
                }
                else if(i&&j==){
                    s[i][j]=s[i-][j];
                    if(c[i][j]==)s[i][j]++;
                }
                else if(i&&j)
                {
                    s[i][j]=s[i-][j]+s[i][j-]-s[i-][j-];
                    if(c[i][j]==)s[i][j]++;
                }
            }
        }
    }
}
int main()
{
    cin>>t>>k;
    get_c();
    get_s();
    while(t--)
    {
        cin>>n>>m;
        cout<<s[n][min(n,m)]<<endl;
    }
}