题意: 求组合数C(p, q) / C(r, s)结果保留5为小数. 分析: 先用筛法求出10000以内的质数,然后计算每个素数对应的指数,最后再根据指数计算答案. #include <cstdio> #include <cmath> #include <cstring> ; int pri[maxn], cnt, e[maxn]; //e记录每个素数的质数 ]; void add_interger(int n, int d) //乘以n的d次幂 { ; i <…

白书P171 对m,n!分解,质因子指数取min #include<iostream> #include<algorithm> #include<cstdio> #include<cstring> #include<cstdlib> #include<string> #include<sstream> #include<vector> #define rep(i,j,k) for(register int i=…

经过紫书的分析,已经将问题转化为求组合数C(n-1, 0)~C(n-1, n-1)中能够被m整除的个数,并输出编号(这n个数的编号从1开始) 首先将m分解质因数,然后记录下每个质因子对应的指数. 由组合恒等式,我们可以递推C(n, k)的质因数的个数. 一个没什么用的小优化:因为杨辉三角每一行都是对称的,所以我们可以求出前一半答案,然后根据对称性求出后一半的答案. 需要注意的是,如果答案中有类似C(10, 5)的数,就不要再对称了,会重复的. 这个优化貌似也只能优化0.05s左右. #inclu…

#include<iostream> #include<stdio.h> #include<algorithm> #include<string.h> #include<math.h> #define ll long long using namespace std; int main(){ int n; ; ){ k++; scanf("%d",&n); )break; ); ; ll ans=; ;i<=m;…

题意:给出G和L,求最小的a使得gcd(a,b)=G,lcm(a,b)=L 显然a>=G,所以a取G,b要满足质因子质数为L的同次数,b取L //此处应有代码…

UVA - 10375 Choose and divide Choose and divide Time Limit: 1000MS   Memory Limit: 65536K Total Submissions: 4053   Accepted: 1318 Description The binomial coefficient C(m,n) is defined as m! C(m,n) = -------- n!(m-n)! Given four natural numbers p, q…

题意:求C(p,q)/C(r,s),4个数均小于10000,答案不大于10^8 思路:根据答案的范围猜测,不需要使用高精度.根据唯一分解定理,每一个数都可以分解成若干素数相乘.先求出10000以内的所有素数,用a数组表示唯一分解式中个素数的指数,求出每个分子部分的素因子,并且相应的素数的指数加一.分母则减一.最后求解唯一分解式的值. #include<stdio.h> #include<string.h> #include<math.h> ; int pr[N],p[N…

UVa 10375 Choose and divide 题目: http://acm.hust.edu.cn/vjudge/problem/viewProblem.action?id=19601 思路: maxn=10000 如果计算maxn!再保存的话显然装不下. 但答案由阶乘的积或商组成,所以可以用唯一分解定理求解.大题思路就是把目前答案的质因子的指数用数组e保存,乘除都对e操作. 需要注意的是筛法求素数优化后的写法. 代码: #include<iostream> #include<…

uva10375 Choose and Divide(唯一分解定理) 题意: 已知C(m,n)=m! / (n!*(m-n!)),输入整数p,q,r,s(p>=q,r>=s,p,q,r,s<=10000), 计算C(p,q)/C(r,s).输出保证不超过10^8,保留5位小数. 分析: 本题时间上基本上没有太大的限制,可以暴力求解C(m,n); #include<cstdio> #include<cstring> #include<iostream>…

首先说一下什么是唯一分解定理 唯一分解定理:任何一个大于1的自然数N,如果N不是质数,那么N可以分解成有限个素数的乘积:例:N=(p1^a1)*(p2^a2)*(p3^a3)......其中p1<p2<p3...... 所以当我们求两个很大的数相除时  唯一分解定理是一个不错的选择,不会爆范围 下面具体说一下怎么求唯一分解定理: 首先我们需要知道所有的素数:  用埃式算法打表求得: void is_prime() { cnt=; ;i<=maxn;i++) { if(!vis[i]) {…