软件工程

本类阅读TOP10

·PHP4 + MYSQL + APACHE 在 WIN 系统下的安装、配置
·Linux 入门常用命令(1)
·Linux 入门常用命令(2)
·使用 DCPROMO/FORCEREMOVAL 命令强制将 Active Directory 域控制器降级
·DirectShow学习(八): CBaseRender类及相应Pin类的源代码分析
·基于ICE方式SIP信令穿透Symmetric NAT技术研究
·Windows 2003网络负载均衡的实现
·一网打尽Win十四种系统故障解决方法
·数百种 Windows 软件的免费替代品列表
·收藏---行百里半九十

分类导航
VC语言Delphi
VB语言ASP
PerlJava
Script数据库
其他语言游戏开发
文件格式网站制作
软件工程.NET开发
看似简单的求R^n...

作者:未知 来源:月光软件站 加入时间:2005-2-28 月光软件站

求出R^n 其中R为实数 0.0 < R < 99.999  n是一个整数 n<=25
大家很容易想到,这还不容易么?可是如果要求最精确的结果呢?就是说小数点之后有多少输出多少...并且计算时间不多于1s...是的,这是一道ACM的试题
我暂时没有考虑精确度的问题,感觉陷阱在1sec上面,于是没有敲出来代码就开始想,可能这是个坏习惯...如果仅仅是r*r*r.....这样岂不是一道白吃的问题?于是想了好久找到如下解决R^n的节省资源的优化办法...

1划分法

a^n = (a^(n/2))^2 n是奇数
a^n = a*a^(n-1)   n是偶数
例子:
普通算法:
   2*2*2*2*2*2*2*2 ~~ 8 次乘法  时间复杂度~~ O(n).
划分法:
   2^8 == (a^4)^2 == ((a^2)^2)^2 ~~ 3 次乘法 时间复杂度~~ O(log n).

O(log n) is much faster than O(n).  :)

long square(long n) { return n*n; }

long fastexp(long base,long power) {
  if (power == 0)
    return 1;
  else if (power%2 == 0)
    return square(fastexp(base,power/2));
  else
    return base * (fastexp(base,power-1));
}

2公式法

a^n = exp(log(a)*n)    pow(a,n)
code:

#include <stdio.h>
#include <math.h>

void main() {
  printf("%lf\n",exp(log(8.0)*1/3.0));
  printf("%lf\n",pow(8.0,1/3.0));
}

但是似乎本题的目的也不在于此 要计算出那么多的位数恐怕不是普通的数据类型能够办到的 在学数据结构最基础的数组的时候 曾经有一个例子 用数组求大数的幂 看来此处可能用得到
那之前的努力白费了么?当然没有...起码我们知道了以后有这种计算的时候如何节省内存和cpu了。。。
待续...
ps:如果有人知道更简便的方法或者用普通数据类型就能解决的麻烦指点我一下 我很笨地:)
参考资料:http://www.comp.nus.edu.sg/~stevenha/index.html




相关文章

相关软件