除法表达式(数论初步) By ACReaper
来源:互联网 发布:培育发展新动力优化 编辑:程序博客网 时间:2024/05/16 05:12
这道题目,读了许久,都没明白过来为什么判断表达式的值能否为整数,只要判断E = x1/(x2/x3/....xk) = x1 * x3 * ...*xk / x2是否为整数。后来才明白,其实是这样的,无论你怎么加括号下x2都必然在分母位置,要证明这个可以用数学归纳法去证明,这里不再证明,但你可以试试几个,稍微验证一下。而且其它的数都可以放在分子位置因为我们都可以这样加括号x1/(x2/x3.....xk)。也就是说,如果一个表达式(加括号后)不是x2的倍数,则,这个(加完括号的)表达式值肯定不是整数,而如果一个表达式的值是x2的倍数则,我们一定能把这个表达式写成成x1 * x3 *.....xk/ x2的形式。换句话说,我们只要判断x1 * x3 *.....xk/ x2这个值就可以知道它是不是能写成值为整数的表达式。可能你会疑惑,请注意所有加括号后值能为整数的,必定也能改写为x1 * x3 *.....xk/ x2,这是个充分必要条件!
so为我们可以直接用约分方法来写这道题目
#include <stdio.h>#define MAXN 10001int A[MAXN];int n;int gcd(int a, int b){return b == 0? a : gcd(b, a % b);}int judge(int A[]){A[2] /= gcd(A[2],A[1]);for(int i = 3; i <= n; i++){A[2] /= gcd(A[i],A[2]);}return A[2] == 1;}int main(){while(scanf("%d",&n) != EOF){for(int i = 1; i<= n; i++){scanf("%d",A + i);}if(judge(A)){printf("Yes\n");}else{printf("No\n");}}return 0;}
这里用到了欧几里德算法,别说不知道啊~~~~~~。
还有一个很重要的定理gcd(a,b) * lcm(a,b) = a * b;so我们可以通过gcd求lcm,lcm(a,b) = a / gcd(a,b) * b,这样写是为了防止乘法溢出!
- 除法表达式(数论初步) By ACReaper
- 幂取模(数论初步) By ACReaper
- 扩展欧几里德算法(数论初步) By ACReaper
- 模线性方程(数论初步) By ACReaper
- 无平方因子的数(数论初步) By ACReaper
- 简单素数筛选法介绍(数论初步) By ACReaper
- 直线上的点(数论初步) By ACReaper
- 同余与模运算(数论初步) By ACReaper
- 大整数取模(数论初步) By ACReaper
- 离散概率初步 By ACReaper
- 表达式树(图论) By ACReaper
- 数论之除法表达式
- 【数论】除法表达式
- 约数的个数(数论中的计数问题) By ACReaper
- 编码与解码(数论中的计数问题) By ACReaper
- 小于n且与n互素的个数(数论中的计数问题) By ACReaper
- 帕斯卡恒等式 By ACReaper
- 范德蒙很等式 By ACReaper
- 云计算应用之Dropcam HD
- (个人笔记) java.util.logging.Logger 写日志
- mysql授权+IP登录
- boost中的range与slice函数简介
- MapReduce 程序模板 (采用 新/旧API)
- 除法表达式(数论初步) By ACReaper
- Debian 5.0.5的介绍与详细安装图解
- MSSQL如何删除字段的所有约束和索引
- windows下查看linux文件的工具--Explore2fs
- vs2010下cppunit 配置和实例
- poj2362 - Square
- VMware中安装CentOS6.4 Minimal并启用共享文件夹
- myeclipse里web project和web service project的区别
- Oracle IMP/EXP 导入导出