问题 A: 组合数
题目描述
输入
输出
样例输入 Copy
3 2
4 2
样例输出 Copy
3 2
6 4
因为求组合数的时候 long long 存不下,所以我们需要分解质因数再求解,它的就是把分子分母约去同时有的素因子以达到中间算阶乘的时候不会爆 long long
计算因子数用到了唯一分解定理
一个数 n 肯定能被分解成 n=p1^r1 * p2^r2 . . .*pn^rn
假设 p1p2,…pn 是它的素因子
假设 r1,r2,…rn 分别是 p1…pn 的幂次数
那么(1+r1)*(1+r2)….*(1+rn)就是他的因子数
/**
/*两个板子 快速幂+组合数
*/
#include<stdio.h>
#include <iostream>
#include <bits/stdc++.h>
#define maxn 200005
typedef long long ll;
using namespace std;
const ll mod=;
ll fac[maxn],inv[maxn];
ll pow_mod(ll a,ll n)
{
ll ret =;
while(n)
{
if(n&) ret=ret*a%mod;
a=a*a%mod;
n>>=;
}
return ret;
}
void init()
{
fac[]=;
for(int i=;i<maxn;i++)
{
fac[i]=fac[i-]*i%mod;
}
}
ll Cc(ll x, ll y)
{
return fac[x]*pow_mod(fac[y]*fac[x-y]%mod,mod-)%mod;
}
//ll pow(ll x,ll n)
//{
//ll temp(x),res(1);
//while(n)
//{
//if(n&1)
//{
//res *= temp;
//}
//temp *= temp;
//n>>=1;
//}
//return res;
//} long long C[][];
ll D(ll m,ll n)
{
for(int i=;i<;++i)
{
C[i][i]=;
C[i][]=;
}
for(int i=;i<;++i)
{
for(int j=i+;j<;++j)
{
C[j][i]=C[j-][i-]+C[j-][i];
}
}
// int m,n;
// while(cin>>m>>n)///C(m,n)
// {
// cout<<C[m][n]<<endl;
// }
return C[m][n];
}
#define N 440
int prime[N];
bool vis[N];
double fact(int n)//求阶乘
{
int i;
double sum;
sum=;
for(i=;i<=n;i++)
{
sum=sum*i;
}
return sum;
} int Prime()
{
int cnt = ;
for (int i = ; i <= N; ++i)
{
if (!vis[i])
{
prime[cnt++] = i;
}
for (int j = ; j < cnt&&i*prime[j] <= N; ++j)
{
vis[i*prime[j]] = ;
if (i%prime[j] == )break;
}
}
return cnt;
}
int num[];
int Fcnt;
void solve(int n,int y)
{
for (int i = ; i < Fcnt; ++i)
{
int c = , p = prime[i];
while (n / p )
{
c += n / p;
p *= prime[i];
}
num[i] = num[i] + y*c;
}
} //int main()
//{
// Fcnt=Prime();
// int n, m;
// while (scanf("%d%d", &n, &m) != EOF){
// memset(num, 0, sizeof(num));
// solve(n, 1);
// solve(m, -1);
// solve(n - m, -1);
// ll ans = 1;
// for (int i = 0; i < Fcnt; ++i)
// {
// ans *= (num[i] + 1);
// }
// printf("%lld\n", ans);
// }
//}
int main(){
ll p,q,k,a,b;
Fcnt=Prime();
int n,m;
double n1,m1,o1;
double fact(int n);
while(scanf("%lld%lld%lld%lld%lld",&p,&q,&n,&m,&b)!=EOF){
printf("%lld\n",p);
n1=fact(n);
m1=fact(m);
o1=fact(n-m);
int k = ;
memset(num, , sizeof(num));
solve(n, );
solve(m, -);
solve(n - m, -);
ll ans = ;
for(int i = ;i*i<=((int)(n1/(m1*o1)));i++){
if((int)(n1/(m1*o1))%i==)
k++;
}
printf("%.0f ",n1/(m1*o1),k); // ll CC = (long long)D(k,a);
ll aa = (long long)pow(q,m);
ll bb = (long long)pow(p,b);
ll sum = ((long long)(n1/(m1*o1)))*aa*bb;
// printf("%lld\n",CC);
printf("%lld\n",p);
printf("%lld\n",sum);
}
}
问题 A: 组合数的更多相关文章
- LCM性质 + 组合数 - HDU 5407 CRB and Candies
CRB and Candies Problem's Link Mean: 给定一个数n,求LCM(C(n,0),C(n,1),C(n,2)...C(n,n))的值,(n<=1e6). analy ...
- 计算一维组合数的java实现
背景很简单,就是从给定的m个不同的元素中选出n个,输出所有的组合情况! 例如:从1到m的自然数中,选择n(n<=m)个数,有多少种选择的组合,将其输出! 本方案的代码实现逻辑是比较成熟的方案: ...
- Noip2016提高组 组合数问题problem
Day2 T1 题目大意 告诉你组合数公式,其中n!=1*2*3*4*5*...*n:意思是从n个物体取出m个物体的方案数 现给定n.m.k,问在所有i(1<=i<=n),所有j(1< ...
- C++单元测试 之 gtest -- 组合数计算.
本文将介绍如何使用gtest进行单元测试. gtest是google单元测试框架.使用非常方便. 首先,下载gtest (有些google项目包含gtest,如 protobuf),复制目录即可使用. ...
- NOIP2011多项式系数[快速幂|组合数|逆元]
题目描述 给定一个多项式(by+ax)^k,请求出多项式展开后x^n*y^m 项的系数. 输入输出格式 输入格式: 输入文件名为factor.in. 共一行,包含5 个整数,分别为 a ,b ,k , ...
- AC日记——组合数问题 落谷 P2822 noip2016day2T1
题目描述 组合数表示的是从n个物品中选出m个物品的方案数.举个例子,从(1,2,3) 三个物品中选择两个物品可以有(1,2),(1,3),(2,3)这三种选择方法.根据组合数的定 义,我们可以给出计算 ...
- 【板子】gcd、exgcd、乘法逆元、快速幂、快速乘、筛素数、快速求逆元、组合数
1.gcd int gcd(int a,int b){ return b?gcd(b,a%b):a; } 2.扩展gcd )extend great common divisor ll exgcd(l ...
- 【BZOJ-4591】超能粒子炮·改 数论 + 组合数 + Lucas定理
4591: [Shoi2015]超能粒子炮·改 Time Limit: 10 Sec Memory Limit: 256 MBSubmit: 95 Solved: 33[Submit][Statu ...
- UOJ263 【NOIP2016】组合数问题
本文版权归ljh2000和博客园共有,欢迎转载,但须保留此声明,并给出原文链接,谢谢合作. 本文作者:ljh2000作者博客:http://www.cnblogs.com/ljh2000-jump/转 ...
- 组合数取模Lucas定理及快速幂取模
组合数取模就是求的值,根据,和的取值范围不同,采取的方法也不一样. 下面,我们来看常见的两种取值情况(m.n在64位整数型范围内) (1) , 此时较简单,在O(n2)可承受的情况下组合数的计算可以 ...
随机推荐
- java.lang.Void类源码解析_java - JAVA
文章来源:嗨学网 敏而好学论坛www.piaodoo.com 欢迎大家相互学习 在一次源码查看ThreadGroup的时候,看到一段代码,为以下: /* * @throws NullPointerEx ...
- 3D打印格式STL
STL格式及其转换 近期接触了3D打印的一些东西,也制作了一个vrml转stl的插件,对该领域多了一些认识. 目前尚没有打印机直接支持stl.obj等格式,在打印之前需使用厂家提供的软件将stl等格式 ...
- STM32开发板的TIM3开启和关闭
关闭定时器中断要考虑好多情况 1)关闭定时器时,定时器是否在处在工作状态 2)关闭定时器时,定时器是否正好进入中断,造成关闭程序出现断层,进而无法实现完整关闭程序,此时可以使用高一级别的外部中断强制进 ...
- Kohana Cache
The default cache group is loaded based on the Cache::$default setting. It is set to the file driver ...
- Http请求状态大全
一.HTTP状态码分类 HTTP状态码由三个十进制数字组成,第一个十进制数字定义了状态码的类型,后两个数字没有分类的作用.HTTP状态码共分为5种类型: 分类 分类描述 1** 信息 服务器收到请求, ...
- 官方文档翻译-Today
原文链接: Today 最后更新:2017-02-18 译文: 应用的今日视图扩展被称作小部件.小部件为用户提供快速访问重要信息的入口.例如,用户打开今日视图实时查看股票价格或天气情况,查看今天的日程 ...
- js刷新当前页面的5种方式
1. reload reload 方法,该方法强迫浏览器刷新当前页面.语法:location.reload([bForceGet]) 参数: bForceGet, 可选参数, 默认为 false, ...
- Charles抓取手机https请求
1.下载Charles工具,3.92破解版:http://pan.baidu.com/s/1cko2L4 密码:chmy 2.安装SSL证书,默认安装就可以 3.证书安装成功后,点击详细信息--> ...
- sqlserver高版本往低版本迁移
本例为sql2012 还原到sql2008. 要实现的功能是把sql2012的数据库备份到sql2008,数据库名字为Test,并且这两个数据库在不同的电脑中. 微软的软件设计方案基本上都是新版本兼容 ...
- @ResponseBody返回4种数据格式的数据
1.返回一个键值对或者集合 前端JS请求: //返回值为map的形式 $(".name").blur(function(){ $.ajax({ type:"Post&qu ...