UVA 1649 Binomial coefficients
https://vjudge.net/problem/UVA-1649
题意:
输入m,求所有的C(n,k)=m
m<=1e15
如果枚举n,那么C(n,k)先递增后递减
如果枚举k,那么C(n,k)单调递增
所以可以枚举k,二分n,直至C(n,k)=m
k枚举到什么时候?
根据公式 C(n,k)=C(n,n-k)
所以只管那个小的k
k<n-k 即 k<n/2,
也就是对于每个n,k只算到n/2
所以 当C(k*2,k)>m 时,停止枚举
然后对于这个k,二分n
边界:l=k*2,r=m
#include<queue>
#include<cstdio>
#include<iostream>
using namespace std;
#define pairr pair<long long,long long>
#define make(a,b) make_pair(a,b)
priority_queue<pairr,vector<pairr>,greater<pairr> >q;
long long n;
long long C(long long m,long long k)
{
long long x=;
for(int i=;i<=k;i++)
{
if(x/i>n/(m-i+)) return n+;
x*=m-i+;
x/=i;
}
return x;
}
void solve()
{
long long l,r,mid,tmp;
for(int k=;C(k<<,k)<=n;k++)
{
l=k<<;
r=n;
while(l<=r)
{
mid=l+r>>;
tmp=C(mid,k);
if(tmp<n) l=mid+;
else if(tmp==n)
{
q.push(make(mid,k));
if(mid!=k<<) q.push(make(mid,mid-k));
break;
}
else r=mid-;
}
}
}
int main()
{
int T;
scanf("%d",&T);
long long siz;
while(T--)
{
scanf("%lld",&n);
solve();
siz=q.size();
printf("%d\n",siz);
while(siz--)
{
printf("(%lld,%lld)",q.top().first,q.top().second);
if(siz) printf(" ");
else printf("\n");
q.pop();
}
}
}
UVA 1649 Binomial coefficients的更多相关文章
- UVa 1649 Binomial coefficients 数学
		题意: \(C(n, k) = m(2 \leq m \leq 10^{15})\),给出\(m\)求所有可能的\(n\)和\(k\). 分析: 设\(minK = min(k, n - k)\),容 ... 
- UVA - 1649 Binomial coefficients (组合数+二分)
		题意:求使得C(n,k)=m的所有的n,k 根据杨辉三角可以看出,当k固定时,C(n,k)是相对于n递增的:当n固定且k<=n/2时,C(n,k)是相对于k递增的,因此可以枚举其中的一个,然后二 ... 
- 51nod1245 Binomial Coefficients Revenge
		题目来源: HackerRank 基准时间限制:2 秒 空间限制:131072 KB 分值: 640 C(M,N) = M! / N! / (M - N)! (组合数).给出M和质数p,求C(M,0 ... 
- 51nod 1245 Binomial Coefficients Revenge
		Description C(M,N) = M! / N! / (M - N)! (组合数).给出M和质数p,求C(M,0), C(M,1)......C(M,M)这M + 1个数中,有多少数不是p的倍 ... 
- 紫书 习题 10-21 UVa 1649 (组合数)
		C(n, k) = m, 固定k,枚举k 这里用到了组合数的一个性质 当k固定的时候,C(2 * k, k) 最小 C(m, k)最大(对于这道题而言是这样,因为大于m 就最终答案不可能为m了) 所以 ... 
- 【51nod 1245】Binomial Coefficients Revenge
		题目大意 C(M,N) = M! / N! / (M - N)! (组合数).给出M和质数p,求C(M,0), C(M,1)......C(M,M)这M + 1个数中,有多少数不是p的倍数,有多少是p ... 
- Some series and integrals involving the Riemann zeta function binomial coefficients and the harmonic numbers
		链接:http://pan.baidu.com/s/1eSNkz4Y 
- 99 Lisp Problems 列表处理(P1~P28)
		L-99: Ninety-Nine Lisp Problems 列表处理类问题的解答,用Scheme实现,首先定义几个在后续解题中用到的公共过程: ; common procedure (define ... 
- UVA10375 Choose and divide 质因数分解
		质因数分解: Choose and divide Time Limit: 3000MS Memory Limit: Unknown 64bit IO Format: %lld & %l ... 
随机推荐
- 自测之Lesson13:共享内存
			题目:创建一个64K的共享内存. 实现代码: #include <stdio.h> #include <sys/ipc.h> #include <sys/shm.h> ... 
- Java 集合框架之Collection
			此图是 java 中 Collection 相关的接口与类的关系的类图.其中,类只是集合框架的一部分,比较常用的一部分. 第一次画类图,着实很费劲,不过收获也不小. 下面是相关接口和类的解释说明.文字 ... 
- 由作业题引发对C++引用的一些思考
			首先分析一段代码: #include <bits/c++config.h> #include <ostream> #include <iostream> #incl ... 
- 第二章 shell的语法
			变量:字符串.数字.环境和参数 获取变量内容可以在变量前使用$字符,使用echo指令可以将变量内容输出到终端. wuchao@wuchao-Lenovo:~$ var=hello wuchao@wuc ... 
- 解决打包遇到的_mssql问题
			明明上一篇还说着打包好了,这一个又出现了问题,真是让人揪心呀!错误如下: Traceback (most recent call last): File "macc.py", li ... 
- vi/sed等遵循的搜索正则语法
			转自:http://blog.csdn.net/lanxinju/article/details/5731843 一.查找 查找命令 /pattern<Enter> :向下查找patter ... 
- 第三部分shell编程3(shell脚本2)
			7. if 判断一些特殊用法 if [ -z $a ] 这个表示当变量a的值为空时会怎么样if grep -q '123' 1.txt; then 表示如果1.txt中含有'123'的行时会怎么样if ... 
- (转)Linux常用性能检测命令
			一.uptime Uptime命令的显示结果包括服务器已经运行了多长时间,有多少登陆用户和对服务器性能的总体评估(load average).load average值分别记录了上个1分钟,5 ... 
- 用php实现一个双向队列 如何实现?
			PHP面试题作业 class DuiLie { private $array = array();//声明空数组 public function setFirst($item){ return arr ... 
- SQL SERVER技术内幕之10 事务并发
			1.事务 1.1事务的定义 事务是作为单个工作单元而执行的一系列操作.定义事务边界有显式和隐式两种.显式事务的定义以BEGIN TRAN作为开始,以COMMIT TRAN提交事务,以ROLLBACK ... 
