唯一分解定理 (Choose and Divide Uva10375)
来源:互联网 发布:广州乐乎城市青年社区 编辑:程序博客网 时间:2024/05/22 18:56
The binomial coefficient C(m; n) is defined as
C(m; n) = m!
(m - n)! n!
Given four natural numbers p, q, r, and s, compute the the result of dividing C(p; q) by C(r; s).
Input
Input consists of a sequence of lines. Each line contains four non-negative integer numbers giving values
for p, q, r, and s, respectively, separated by a single space. All the numbers will be smaller than 10,000
with p ≥ q and r ≥ s.
Output
For each line of input, print a single line containing a real number with 5 digits of precision in the fraction, giving the number as described above. You may assume the result is not greater than 100,000,000.
Sample Input
10 5 14 9
93 45 84 59
145 95 143 92
995 487 996 488
2000 1000 1999 999
9998 4999 9996 4998
Sample Output
0.12587
505606.46055
1.28223
0.48996
2.00000
3.99960
C(m; n) = m!
(m - n)! n!
Given four natural numbers p, q, r, and s, compute the the result of dividing C(p; q) by C(r; s).
Input
Input consists of a sequence of lines. Each line contains four non-negative integer numbers giving values
for p, q, r, and s, respectively, separated by a single space. All the numbers will be smaller than 10,000
with p ≥ q and r ≥ s.
Output
For each line of input, print a single line containing a real number with 5 digits of precision in the fraction, giving the number as described above. You may assume the result is not greater than 100,000,000.
Sample Input
10 5 14 9
93 45 84 59
145 95 143 92
995 487 996 488
2000 1000 1999 999
9998 4999 9996 4998
Sample Output
0.12587
505606.46055
1.28223
0.48996
2.00000
3.99960
简单翻译一下,告诉你p,q,r,s求C(p,q)/C(r,s)
简单的暴力是不可能过得,所以我们要用数学的方法,我们运用唯一分解法,将质因数分解之后存在数组中,这样乘除就是这个因数的个数的加减,最后将所有因数乘起来就可以了~
#include<bits/stdc++.h>using namespace std;vector<int> prime;int book[20000],e[20000];void add_integer(int n,int d){for(int i=0;i<prime.size();i++){while(n%prime[i]==0){e[i]+=d;n/=prime[i];}if(n==1) break;}}void add(int n,int d){for(int i=1;i<=n;i++)add_integer(i,d);}void init(){book[1]=1;for(int i=2;i<=101;i++)if(!book[i])for(int j=i*i;j<=10000;j+=i)book[j]=1;for(int i=1;i<=10000;i++)if(!book[i])prime.push_back(i);}double Pow(int a,int b){int ok=0;if(b==0) return 1;if(b<0) ok=1,b=-b;double ans=1;for(int i=1;i<=b;i++)ans*=double(a);if(ok) return double(1/ans);return ans;}int main(){init();int p,q,r,s;while(cin>>p>>q>>r>>s){memset(e,0,sizeof(e));add(p,1);add(q,-1);add(p-q,-1);add(r,-1);add(s,1);//作为除数的话和前面的正好相反 add(r-s,1);double ans=1;for(int i=0;i<prime.size();i++)ans*=Pow(prime[i],e[i]);printf("%.5lf\n",ans);}return 0;}
阅读全文
0 0
- uva10375 Choose and Divide(唯一分解定理)
- UVA10375 Choose and Divide (唯一分解定理)
- 唯一分解定理 (Choose and Divide Uva10375)
- uva10375(唯一分解定理) Choose and divide
- Uva10375 Choose and divide(数论、唯一分解定理、快速幂)
- 【数论】Choose and Divide, UVa10375 【组合数学】【唯一分解定理】【精度】
- uva10375 Choose and divide 筛素数,唯一分解原理
- UVA 10375 Choose and divide(唯一分解定理)
- UVA 10375 Choose and divide(唯一分解定理)
- UVA 10375 Choose and divide (唯一分解定理)
- UVA10375 Choose and divide 质因数分解
- Choose and divide唯一分解定理/分解/直接 UVA
- uva10375(唯一分解定理)
- UVa10375(唯一分解定理)
- uva10375(唯一分解定理)
- uvaoj 10375 Choose and divide 唯一分解定理
- UVA 10375 - Choose and divide【唯一分解定理】
- Uva 10375 Choose and divide 素数唯一分解定理
- 扩展欧几里得模板(洛谷1082 同余方程NOIP 2012 提高组 第二天 第一题)
- python流程控制语句
- C/C++中struct与class的区别
- 快速幂取模+数学知识(Colossal Fibonacci Numbers! uva 11582)
- 扩展欧几里得应用解不定方程(uva12169,Disgruntled Judge,NWERC 2008)
- 唯一分解定理 (Choose and Divide Uva10375)
- 唯一分解定理(Minimum Sum LCM,Uva10791)
- Java中不容易被注意到的地方
- 欧拉函数phi值的计算模板
- 欧拉函数phi的应用(洛谷2158 [SDOI2008]仪仗队)
- 数学方法(GCD XOR Uva12716)
- 字符串转换成整数+不用加减乘除做加法
- 组合数杨辉三角(Irrelevant Elements uva1635)
- 堆or优先队列(洛谷1631 序列合并)