CJOJ 2255 【NOIP2016】组合数问题 / Luogu 2822 组合数问题 (递推)

Description

组合数$$C^m_n$$表示的是从n个物品中选出m个物品的方案数。举个例子,从(1,2,3) 三个物品中选择两个物品可以有(1,2),(1,3),(2,3)这三种选择方法。根据组合数的定 义,我们可以给出计算组合数的一般公式:

\[C^m_n=\frac{n!}{m!(n-m)!}
\]

其中n! = 1 × 2 × · · · × n

小葱想知道如果给定n,m和k,对于所有的0 <= i <= n,0 <= j <= min(i,m)有多少对 (i,j)满足是k的倍数。

Input

第一行有两个整数t,k,其中t代表该测试点总共有多少组测试数据。

接下来t行每行两个整数n,m。

Output

t行,每行一个整数代表答案。

Sample Input

1 2

3 3

Sample Output

1

Http

Luogu:https://www.luogu.org/problem/show?pid=2822

CJOJ:http://oj.changjun.com.cn/problem/detail/pid/2255

Source

递推

解决思路

不要被题目给出的那个式子给迷惑了!要是真按照那么算,肯定会爆掉!

这道题我们可以离线地来做。我们知道,计算组合数除了上面题目中给出的通项式之外,还有一个递推式:

\[C^m_n=C^{m-1}_{n}+C^{m-1}_{n-1}
\]

那么我们可以开一个二维数组存下所有的组合数(题目范围只有2000嘛)。至于是否整除的问题,我们可以在每一步计算的时候%k,这样既保证了不会爆掉long long 又能求出是否整除(若最后算出来是0就是可以整除啦)。

我们再把能整除的置为1,不能整除的置为0,求一下二维前缀和,然后就可以O(1)地回答询问啦。

代码

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
using namespace std; const int maxN=2000;
const int inf=2147483647; int n,m,mod;
int C[maxN+10][maxN+10];
int If[maxN+10][maxN+10]; int main()
{
int TT;
cin>>TT>>mod;
memset(C,0,sizeof(C));
for (int i=0;i<=maxN;i++)//递推算出所有的组合数并取膜
{
C[i][0]=1;
for (int j=1;j<=i;j++)
C[i][j]=(C[i-1][j]+C[i-1][j-1])%mod;
}
for (int i=0;i<=maxN;i++)
for (int j=0;j<=i;j++)//构造二维前缀和初值
If[i][j]=(bool)(C[i][j])^1;
/*for (int i=0;i<=10;i++)
{
for (int j=0;j<=10;j++)
cout<<If[i][j]<<' ';
cout<<endl;
}
cout<<endl;*/
for (int i=1;i<=maxN;i++)//计算二维前缀和
for (int j=1;j<=maxN;j++)
If[i][j]=If[i][j]+If[i-1][j]+If[i][j-1]-If[i-1][j-1];
/*for (int i=0;i<=10;i++)
{
for (int j=0;j<=10;j++)
cout<<If[i][j]<<' ';
cout<<endl;
}
cout<<endl;*/
for (int ti=1;ti<=TT;ti++)//回答询问
{
cin>>n>>m;
cout<<If[n][m]<<endl;
}
return 0;
}

最新文章

  1. vsftpd安装配置 530 Permission denied.错误
  2. wampserver的安装以及使用
  3. 基础知识:Q&amp;A
  4. java 性能优化:35 个小细节,让你提升 java 代码的运行效率
  5. 和以往印象不同的Java
  6. Centos7安装完毕后重启提示Initial setup of CentOS Linux 7 (core)的解决方法
  7. 理解Java接口
  8. jquery左右链接类似frameset的插件
  9. cwRsync window下的跨服务器的文件同步
  10. Swift游戏实战-跑酷熊猫 12 与平台的碰撞
  11. Android进阶——声波振幅显示
  12. STL 常用的一些容器总结
  13. Centos下删除文件名乱码文件
  14. zend studio里面这块注释是用什么快捷键按出来的?
  15. Cannot find a valid baseurl for repo: base
  16. CentOS openssh升级到openssh可用
  17. csrf jsonp
  18. Learning to Rank(转)
  19. Vue中的事件与常见的问题处理
  20. 当ORACLE归档日志满后如何正确删除归档日志

热门文章

  1. mysqldump命令详解
  2. python3-如何正常使用HTMLTestRunner.py,生成自动化测试报告
  3. 简单vue项目脚手架(vue+webpack2.0+vuex+vue-router)
  4. (转) Unicode(UTF-8, UTF-16)令人混淆的概念
  5. Hibernate入门(一)
  6. mac 安装mysqldb组件包及mac中安装mysql-python遇到的问题
  7. SQL SERVER 使用BULK Insert将txt文件中的数据批量插入表中(1)
  8. 关于HTML学习整理(一)
  9. Xmanager连接到RedHat Enterprise Linux 6.8
  10. salesforce零基础学习(七十三)ProcessInstanceWorkItem/ProcessInstanceStep/ProcessInstanceHistory浅谈