noip2016D2T1 组合数问题

最新推荐文章于 2018-08-20 11:03:00 发布
最新推荐文章于 2018-08-20 11:03:00 发布
阅读量206 收藏
点赞数
分类专栏： 数论
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/xh_0319/article/details/77926283
版权

题目描述

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

C_n^m=\frac{n!}{m!(n - m)!}C​n​m​​=​m!(n−m)!​​n!​​

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

小葱想知道如果给定n,m和k，对于所有的0 <= i <= n，0 <= j <= min(i,m)有多少对 (i,j)满足C_i^jC​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说明】

在所有可能的情况中，只有C_2^1 = 2C​2​1​​=2是2的倍数。

【子任务】


TLE了两个点

没想好前缀和怎么写【叹气

90分~~

#include<bits/stdc++.h>
using namespace std;
template <typename T> void read(T &x){
	x=0;int f=1;char ch=getchar();
	for(;!isdigit(ch);ch=getchar())if(ch=='-')f=-1;
	for(;isdigit(ch);ch=getchar())x=(x<<1)+(x<<3)+ch-'0';
	x*=f;
}

typedef long long LL;

int T,K;
int n,m;

LL f[2010][2010];

int main(){
	read(T),read(K);
	memset(f,0,sizeof(f));
	f[1][1]=1;
	for(int i=1;i<=2001;++i){
		f[i][1]=i%K;
		f[i][i]=1;
	}
	for(int i=2;i<=2001;++i)
		for(int j=2;j<=i-1;++j){
			f[i][j]=f[i-1][j-1]+f[i-1][j];
			if(f[i][j]>=K) f[i][j]%=K;
		}
	/*for(int i=1;i<=6;++i){
		for(int j=1;j<=7;++j)cout<<"dp["<<i<<"]["<<j<<"]="<<dp[i][j]<<" "; 
		cout<<endl;
	}*/
	for(int t=1;t<=T;++t){
		read(n),read(m);
		LL ans=0;
		/*for(int i=1;i<=n;++i){
			for(int j=1;j<=m;++j) cout<<"f["<<i<<"]["<<j<<"]="<<f[i][j]<<" "; 
			cout<<endl;
		}*/
		for(int i=1;i<=n;++i){
			for(int j=1;j<=min(i,m);++j)
				if(f[i][j]==0){
					//cout<<i<<" "<<j<<endl;
					ans++;
				}
		}
		cout<<ans<<endl;
	}
 	return 0;
}

加了一个比较弱的前缀和 95分


#include<bits/stdc++.h>
using namespace std;
template <typename T> void read(T &x){
	x=0;int f=1;char ch=getchar();
	for(;!isdigit(ch);ch=getchar())if(ch=='-')f=-1;
	for(;isdigit(ch);ch=getchar())x=(x<<1)+(x<<3)+ch-'0';
	x*=f;
}

typedef long long LL;

int T,K;
int n,m;

LL f[2010][2010];
LL dp[2010][2010];

int main(){
	read(T),read(K);
	memset(f,0,sizeof(f));
	memset(dp,0,sizeof(dp));
	f[1][1]=1;
	for(int i=1;i<=2001;++i){
		f[i][1]=i%K;
		f[i][i]=1;
	}
	for(int i=2;i<=2001;++i)
		for(int j=2;j<=i-1;++j){
			f[i][j]=f[i-1][j-1]+f[i-1][j];
			if(f[i][j]>=K) f[i][j]%=K;
		}
	for(int i=1;i<=2001;++i){
		if(f[i][1]==0) dp[i][1]=1;
		for(int j=2;j<=i;++j)
			if(f[i][j]==0) dp[i][j]=dp[i][j-1]+1;
			else dp[i][j]=dp[i][j-1];
		for(int k=i+1;k<=2001;++k)
			dp[i][k]=dp[i][k-1];
	}
	/*for(int i=1;i<=6;++i){
		for(int j=1;j<=7;++j)cout<<"dp["<<i<<"]["<<j<<"]="<<dp[i][j]<<" "; 
		cout<<endl;
	}*/
	for(int t=1;t<=T;++t){
		read(n),read(m);
		LL ans=0;
		/*for(int i=1;i<=n;++i){
			for(int j=1;j<=m;++j) cout<<"f["<<i<<"]["<<j<<"]="<<f[i][j]<<" "; 
			cout<<endl;
		}*/
		/*for(int i=1;i<=n;++i){
			for(int j=1;j<=min(i,m);++j)
				if(f[i][j]==0){
					//cout<<i<<" "<<j<<endl;
					ans++;
				}
		}
		cout<<ans<<endl;*/ 
		for(int i=1;i<=n;++i) ans+=dp[i][m];
		cout<<ans<<endl;
	}
 	return 0;
}


二维前缀和（标算

#include<bits/stdc++.h>
using namespace std;
template <typename T> void read(T &x){
	x=0;int f=1;char ch=getchar();
	for(;!isdigit(ch);ch=getchar())if(ch=='-')f=-1;
	for(;isdigit(ch);ch=getchar())x=(x<<1)+(x<<3)+ch-'0';
	x*=f;
}

typedef long long LL;

int T,K;
int n,m;

LL f[2010][2010];
LL dp[2010][2010];

int main(){
	read(T),read(K);
	memset(f,0,sizeof(f));
	memset(dp,0,sizeof(dp));
	f[1][1]=1;
	for(int i=1;i<=2001;++i){
		f[i][1]=i%K;
		f[i][i]=1;
	}
	for(int i=2;i<=2001;++i)
		for(int j=2;j<=i-1;++j){
			f[i][j]=f[i-1][j-1]+f[i-1][j];
			if(f[i][j]>=K) f[i][j]%=K;
		}
	/*for(int i=1;i<=2001;++i){
		if(f[i][1]==0) dp[i][1]=1;
		for(int j=2;j<=i;++j)
			if(f[i][j]==0) dp[i][j]=dp[i][j-1]+1;
			else dp[i][j]=dp[i][j-1];
		for(int k=i+1;k<=2001;++k)
			dp[i][k]=dp[i][k-1];
	}*/
	for(int i=2;i<=2001;++i){
		for(int j=1;j<=i;++j){
			if(f[i][j]==0) dp[i][j]=dp[i-1][j]+dp[i][j-1]-dp[i-1][j-1]+1;
			else dp[i][j]=dp[i-1][j]+dp[i][j-1]-dp[i-1][j-1];
		}
		for(int k=i+1;k<=2001;++k)
			dp[i][k]=dp[i][k-1];
	} 
	/*for(int i=1;i<=6;++i){
		for(int j=1;j<=7;++j)cout<<"dp["<<i<<"]["<<j<<"]="<<dp[i][j]<<" "; 
		cout<<endl;
	}*/
	for(int t=1;t<=T;++t){
		read(n),read(m);
		LL ans=0;
		/*for(int i=1;i<=n;++i){
			for(int j=1;j<=m;++j) cout<<"f["<<i<<"]["<<j<<"]="<<f[i][j]<<" "; 
			cout<<endl;
		}*/
		/*for(int i=1;i<=n;++i){
			for(int j=1;j<=min(i,m);++j)
				if(f[i][j]==0){
					//cout<<i<<" "<<j<<endl;
					ans++;
				}
		}
		cout<<ans<<endl;*/ 
		/*for(int i=1;i<=n;++i) ans+=dp[i][m];
		cout<<ans<<endl;*/
		cout<<dp[n][m]<<endl;
	}
 	return 0;
}





酒酿小河桃
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
NOIP2016D2T1-组合数问题
LiuKairui的博客
07-19 354
问题描述 组合数 CmnCnmC_n^m 表示的是从 n 个物品中选出 m 个物品的方案数。举个例子，从 (1,2,3) 三个物品中选择两个物品可以有 (1,2),(1,3),(2,3) 这三种选择方法。根据组合数的定义，我们可以给出计算组合数 CmnCnmC_n^m的一般公式： Cmn=n!m!(n−m)!Cnm=n!m!(n−m)!C_n^m=\frac{n!}{m!(n-m)!} 其中...
NOIP2016D2T1组合数问题
luoyuef的博客
10-31 355
题目复制变成乱码0.0 这道题刚开始看上去是很懵很懵的，仔细读题，还是很懵。给的组合数公式基本上没什么用（应该都会算吧。。？） 然后自己手写了一遍C（1，1） 到C（5，5）发现组合数递推式f【i，j】=f【i-1】【j】+f【i-1】【j-1】； 然后我想这题不是水题就秒了嘛？开个long long算一下有几个能被整除，询问用二维前缀和维护，然而并不是这样。。开了long long还是爆了0.0
[NOIP2016提高组]D2T1：组合数问题
TYB的博客
08-05 850
暴力
【NOIP2016】D2 T1 组合数问题
风凌残度的博客
12-08 451
数论 杨辉三角形 组合
[NOIP]2016 D2T1 组合数问题
Le Petit Prince
08-25 452
组合数问题具体描述： 组合数Cmn表示从n个物品中选m个物品的方案数，根据组合数的定义，我们可以知道计算组合数的一般公式： Cmn=n！/（m！*（n-m）！） 现在，小葱想知道如果给的n，m，k，  对于0ji是k的倍数 输入描述 Input Description 第一行有两个整数t, k，其中t代表该测试点总共有多少组测试数据，k的意义见【问题描述】。 
解题报告——noip 2016提高组
11-14
解题报告——noip 2016提高组,初赛详细解析。
NOIP2016普及组与提高组试题：NOIP2016复赛普及组试题&NOIP2016复赛提高组试题
11-20
《NOIP2016普及组与提高组试题解析》 全国青少年信息学奥林匹克联赛（National Olympiad in Informatics in Provinces，简称NOIP）是中国计算机学会主办的一项旨在选拔优秀中学生进入更高层次信息学竞赛的活动。...
NOIP2016 提高组 复赛 换教室
11-26
NOIP2016提高组复赛换教室
NOIP2016复赛数据
10-24
【NOIP2016复赛数据】是中国信息学奥赛（NOIP）2016年全国决赛阶段的比赛数据，对于参赛者来说，这些数据是理解和分析比赛问题、进行算法设计与编程实践的重要资源。NOIP，全称为全国青少年信息学奥林匹克联赛，是由...
NOIP2016提高组试题（C++）
09-16
NOIP2016提高组试题（C++）,初赛原题，很好的资源，可以用于教学
NOIP2016 D2T1 组合数问题
baben4194的博客
08-20 123
洛谷P2822 数学真重要啊…… 其实解这一题的关键就是组合恒等式：C(n,m)=C(n-1,m)+C(n-1,m-1)，然后再知道组合数的矩阵（杨辉三角）和题中n,m的关系就很容易解决了（然而做这题之前我并不知道组合恒等式于是杯具了） 由上文提到的几何恒等式，我们可以将组合数打成如下矩阵（C(i,j)行数代表i，列数代表j（均从0开始）） 1 1 1 1 2 1...
洛谷 P2822 [NOIP2016 D2T1] 组合数问题 [90分伪题解]
依然一笑作春温。
11-24 1673
唯一分解定理
[NOIP2016D2T1]组合数问题
ymzqwq的博客
11-06 388
传送门 一道水题debug一上午 我一定要发上来纪念一下 //话说最近在做历年NOIP，T1和T2有一道必WA。。。。弱成这样我大概没救了QAQ 其实就是C[n][m]=C[n−1][m−1]+C[n−1][m]C[n][m]=C[n−1][m−1]+C[n−1][m]C[n][m]=C[n-1][m-1]+C[n-1][m] 然后把前缀和预处理出来，就可以O（1）O（1）O（1）出答案啦...
【noip2016】d2解题报告
weixin_30266829的博客
11-10 62
t1：组合数问题 之前有ac代码来着……现在找不到了，回头再补个 t2：蚯蚓，三个队列维护，前两天刚做 之前博客：http://www.cnblogs.com/Loi-dfkdsmbd/p/7799604.html t3：愤怒的小鸟 丢个dalao博客：http://blog.csdn.net/qq_36312502/article/details/78482945 转载于:https:...
NOIP2016D2T2-蚯蚓
LiuKairui的博客
07-17 394
问题描述 本题中，我们将用符号⌊c⌋ 表示对 c 向下取整，例如：⌊3.0⌋=⌊3.1⌋=⌊3.9⌋=3 。 蛐蛐国最近蚯蚓成灾了！隔壁跳蚤国的跳蚤也拿蚯蚓们没办法，蛐蛐国王只好去请神刀手来帮他们消灭蚯蚓。 蛐蛐国里现在共有 n 只蚯蚓（ n 为正整数）。每只蚯蚓拥有长度，我们设第 i 只蚯蚓的长度为 ai (i=1,2,…,n )，并保证所有的长度都是非负整数（即：可能存在长度为 0 的蚯...
[NOIP2016D2T2]借教室
OIER成长之路
07-20 1507
二分或线段树 NOIP2016D2T2
NOIP2016D2T3-愤怒的小鸟
LiuKairui的博客
07-21 320
在外面，不方便写博客，代码刚过了，存一下，回去慢慢补… 小鸟 #include &amp;lt;iostream&amp;gt; #include &amp;lt;cstdio&amp;gt; #include &amp;lt;cstring&amp;gt; #include &amp;lt;algorithm&amp;gt; #include &amp;lt;cmath&amp;gt; using namespace st
BZOJ4870: [Shoi2017]组合数问题
neither_nor
04-24 1880
人话题意：从nk个物品里选模k余r个物品，问方案数模P 大傻逼题，绝对傻逼题，但是考场上没几个人想出来…… 太傻逼了 f[i][j]表示前i个物品选模p余j个，然后矩乘即可 #include #include #include #include #include #include #include #include #include #include #include #include #
【组合数学+转化问题】BZOJ4005[JLOI2015]骗我呢
MEET YOU FIRST TIME
02-01 1278
一道不错的组合数学 转化问题的思路很优秀
[BZOJ4870][Shoi2017]组合数问题（矩阵乘法）
女装的你，如此好看！
06-02 292
Source https://www.lydsy.com/JudgeOnline/problem.php?id=4870 Solution 首先，组合数有一个重要的性质： Cmn=Cmn−1+Cm−1n−1Cnm=Cn−1m+Cn−1m−1C_n^m=C_{n-1}^m+C_{n-1}^{m-1} 设： f(m,r)=∑i=0∞Cik+rmf(m,r)=∑i=0∞Cmik+rf(m...

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值