题目描述

解析

如果它不取等，那就和方程的解这道题一样了，但有了等号就很头疼
如何把等号去掉呢？
定义Bi=Ai+iB~i~=A~i~+iB i =A i +i那么我们就可以得到：0<B1<B2<...<Bm<=m+n0<B~1~<B~2~<...<B~m~<=m+n0<B 1 <B 2 <...<B m <=m+n,这样就可以转化为方程的解了
求C（n+m,n)即可
但是还有一个问题，就是它的数据规模太大，怎么办？
从洛谷学到了一种**分解质因数求组合数的方法，具体见下文代码吧

代码

#include<bits/stdc++.h>
using namespace std;
const int N=5e5+100;
const int mod=1000;
int t,k,x;
int n,m;
int ans[N],cnt;
void cheng(int x){for(int i=1;i<=cnt;i++) ans[i]*=x;int num=0;while(x){num++;x/=10;}cnt+=num+1;for(int i=1;i<=cnt;i++){ans[i+1]+=ans[i]/10;ans[i]%=10;}while(ans[cnt]==0) cnt--;cnt=min(cnt,102);
}
//int now[N];
//void chu(int x){//  int res=0;
//  for(int i=cnt;i>=1;i--){//      res*=10;res+=ans[i];
//      ans[i]=res/x;
//      res%=x;
//  }
//  while(ans[cnt]==0) cnt--;
//}
int p[N],v[N],c[N],tot;
int id[N];
void find_prime(){int top=n+m;for(int i=2;i<=top;i++){if(!v[i]){v[i]=i;p[++tot]=i;id[i]=tot;}for(int j=1;j<=tot;j++){int now=p[j];if(now>top/i||now>v[i]) continue;}}
}
void add_divide(int x){int top=floor(sqrt(x));for(int i=1;i<=tot;i++){int now=p[i];if(now>top) break;while(x%now==0){x/=now;c[i]++;}}if(x){c[id[x]]++;}
}
void minus_divide(int x){int top=floor(sqrt(x));for(int i=1;i<=tot;i++){int now=p[i];if(now>top) break;while(x%now==0){x/=now;c[i]--;}}if(x){c[id[x]]--;}
}
int main(){scanf("%d%d",&n,&m);find_prime();cnt=1;ans[1]=1;for(int i=2;i<=n+m;i++){add_divide(i);}for(int i=2;i<=n;i++){minus_divide(i);}for(int i=2;i<=m;i++){minus_divide(i);}for(int i=1;i<=tot;i++){for(int j=1;j<=c[i];j++){cheng(p[i]);}}for(int i=100;i>=1;i--) printf("%d",ans[i]);
}

YBTOJ：数列方案（组合数学）相关推荐

【BZOJ1485】[HNOI2009]有趣的数列（组合数学）
[BZOJ1485][HNOI2009]有趣的数列(组合数学) 题面 BZOJ 洛谷题解从小往大填数,要么填在最小的奇数位置,要么填在最小的偶数位置. 偶数位置填的数的个数不能超过奇数位置填的数的 ...
【XSY2808】董先生的休闲方案组合数学
题目描述有\(n\)个方案,编号为\(1\ldots n\). 最开始你不知道每个方案的编号. 你要按顺序提出这些方案. 每一个时刻你要做以下事情: 如果你阅读过下一个方案,就提出这个方案. 否则随 ...
ZZNU-OJ-2119 : 告辞，【卡特兰数列，组合数学】
2119 : 告辞时间限制:1 Sec 内存限制:256 MiB 提交:428 答案正确:102 提交状态编辑讨论区题目描述整个世界都在散发着恋爱的恶臭,只有spring依旧保持着单身贵族 ...
组合数学 —— 卡特兰数列（Catalan）
[概述] 卡特兰数列是组合数学中一个常出现在各种计数问题中出现的数列,其前几项为 : 1, 1, 2, 5, 14, 42, 132, 429, 1430, 4862, 16796, 58786, 2 ...
【Java系列】斐波拉契数列
介绍斐波那契数列(Fibonacci sequence),又称黄金分割数列,因数学家莱昂纳多·斐波那契(Leonardo Fibonacci)以兔子繁殖为例子而引入,故又称为"兔子数列&q ...
暑期训练日志----2018.7.31
训练第二天. 早上主持开会,讨论了半个多小时,发现很多人达到了传说的境界--会用就行,在我讲线性同余方程求法的时候,很多人认为那就是个模版,通解特解套用就行,但是为什么这么样套,以及定理2是怎么应用的 ...
以LOL为起点，谈谈游戏数值设计核心思路
文档在今年三月份我动笔写了一小部分,但当时思路凌乱,行文梗阻,于是丢在一边构思了半年,现在又回过头来慢慢写,希望能写好写完吧,初衷是希望即时萌新也能看懂,但是好像并不能行--本来几个数据就能弄好的东西 ...
游戏数值——LOL篇以LOL为起点-说游戏数值设计核心思路
附文文档在今年三月份我动笔写了一小部分,但当时思路凌乱,行文梗阻,于是丢在一边构思了半年,现在又回过头来慢慢写,希望能写好写完吧,初衷是希望即时萌新也能看懂,但是好像并不能行--本来几 ...
【组合数学】生成函数 ( 正整数拆分 | 重复有序拆分 | 不重复有序拆分 | 重复有序拆分方案数证明 )
文章目录一.重复有序拆分二.不重复有序拆分 1.无序拆分基本模型 2.全排列三.重复有序拆分方案数证明参考博客 : 按照顺序看 [组合数学]生成函数简要介绍 ( 生成函数定义 | 牛顿二项式 ...

YBTOJ：数列方案（组合数学）

文章目录

题目描述

解析

代码

YBTOJ：数列方案（组合数学）相关推荐

最新文章

热门文章