Exponentiation(高精度幂计算)(标签:链表,字符串,快速幂计算)

题目描述

对数值很大、精度很高的数进行高精度计算是一类十分常见的问题。比如,对国债进行计算就是属于这类问题。
现在要你解决的问题是:对一个实数R( 0.0 < R < 99.999 ),要求写程序精确计算 R 的 n 次方(Rn),其中n 是整数并且 0 < n <= 25。

Input

T输入包括多组 R 和 n。 R 的值占第 1 到第 6 列,n 的值占第 8 和第 9 列。

Output

对于每组输入,要求输出一行,该行包含精确的 R 的 n 次方。输出需要去掉前导的 0 后不要的 0 。如果输出是整数,不要输出小数点。

Sample Input

95.123 12
0.4321 20
5.1234 15
6.7592 9
98.999 10
1.0100 12

Sample Output

548815620517731830194541.899025343415715973535967221869852721
.00000005148554641076956121994511276767154838481760200726351203835429763013462401
43992025569.928573701266488041146654993318703707511666295476720493953024
29448126.764121021618164430206909037173276672
90429072743629540498.107596019456651774561044010001
1.126825030131969720661201

解题思路

对于此类高精度计算,一定不能用常规的+、-、*、/运算符运算,因其会导致溢出问题,所以考虑将原数据的每一位分开,并且顺序放入线性表中,模拟竖式的计算过程进行计算。在选择线性表的种类时,根据实际需要选择数组或链表。在C++ STL标准模板库中,以数组实现的线性表的容器是vector,其容量可以随数据大小变化,即在开始时会根据初始数据的大小开辟一个比数据的长度稍大的空间,比如初始数组中只有一个数,但是容器会先开辟8个数据的容量。若后续对容器进行数据量的增加,会填满这8个空位,那么容器会自动判断是否有空位,容量不够的话,会重新开辟一块更大容量的空间,并且将原数组中的内容拷贝到新的空间。这就导致,在此问题中若用数组实现线性表时,会浪费更多的程序运行时间,所以此处选择链表形式的线性表,在STL中,使用线性表需要添加头文件:

#include <list>

将题目中输入的s转存到list当中,在转存时,为了方便操作,将高位与低位倒置放在链表中(因为模拟竖式运算时,需要从低位向高位计算)。并且同时记录原始数据中小数点的位置。

list<int> origin, res(, ); //origin为转换为链表之后存储底数的变量
int dot_pos = , res_dot = ; //dot_pos为原底数小数点位置, res_dot记录结果的小数点位置
for (int i = s.size() - ; i >= ; --i) {
if (s[i] == '.') dot_pos = s.size() - i - ; //记录是否在原底数中有小数点的存在,记录逆序小数点位置
else origin.push_back(s[i]-''); //转化char为int类型
}

接下来运用快速幂计算的算法思路来求解s的n次幂。普通的幂计算是用递归(递推)方式求解,计算n次幂需要n次递归计算,时间复杂度为O(n)。而快速幂算法的简单思路为:若计算s^n,将n表示为二进制形式,如5可表示为101,即1+4,则,s^5=s^(1+4)=s^1*s^4=1*s^1*s^4。根据以上所述,依次计算s^1, s^2, s^4, ...,将s^1与s^4累乘到1上,即可得到结果。时间复杂度为O(logn)。详细代码如下:

/*使用快速幂运算求高精度幂*/
while (n > ) {
if (n & ) {
multiple(res, origin);
res_dot += dot_pos; //更新结果的小数点位置
}
multiple(origin, origin);
if (dot_pos > ) dot_pos <<= ;
n >>= ;
}

其中,multiple计算两个存储在链表中的两个数的乘积,结果通过第一个参数的引用返回。multiple函数的定义如下:

/*计算存储在两个链表中的两个数的乘积,结果通过Num1的引用返回*/
void multiple(list<int>& num1, list<int> num2) { //模拟乘法竖式进行计算,规定num2为乘数,num1为被乘数
list<int> temp, origin_num1(num1); //temp存储乘数与被乘数其中一位的临时结果,origin_num1为num1的拷贝,以防后续将num1清空导致原数据丢失
int key_index = ; //记录两位数相乘时被乘数的左移位数
num1.erase(num1.begin(), num1.end());
for (auto key : origin_num1) {
temp.erase(temp.begin(), temp.end());
for (int i = ; i < key_index; ++i) temp.push_back();
multiple_one(num2, key, temp);
add(num1, temp);
++key_index;
}
}

注意,因为传入的两个数num1与num2可能是同一个变量,所以此处实现时,为了函数体中能够保留num2原参数值,没有传入引用。传入引用时,C++编译器直接对实参进行操作,类似指针,而不加&符号时,C++在初始化函数后,拷贝一份实参的实体,并对这个新的实体进行操作,关于C++引用的详细内容,可参考这里

上述代码中的multiple_one与add函数分别计算一位数与多位数相乘以及两个多位数相加,实现思想与multiple函数相同,即模拟竖式计算。其函数原型如下:

/*计算一个一位整数与存储在链表L-H中的数的乘积,结果存储在result的引用返回*/
void multiple_one(list<int>& num, int one, list<int>& result); //模拟竖式计算多位乘以一位/*计算存储在链表L-H中的两个数的和,结果存储在num1*/
void add(list<int>& num1, list<int>& num2); //模拟竖式计算两个数相加

完成快速幂计算后,程序的主体已基本完成,最后的步骤就是根据题目要求格式化输出,即去掉前导与末尾的0。此处实现格式化时得出最后结果,对结果进行操作。也可以直接对原数据进行去零操作。对于程序运行时间来说,后者效果更好。

源代码

~~~Code~~~

C++ 11新标准实现POJ No.1001-Exponentiation的更多相关文章

  1. C++ 11新标准实现POJ No.1002-487-3279

    487-3279(重复的电话号码查询)(标签:优先队列,哈希表) 题目描述 企业喜欢用容易被记住的电话号码.让电话号码容易被记住的一个办法是将它写成一个容易记住的单词或者短语.例如,你需要给滑铁卢大学 ...

  2. C++ 11新标准实现POJ No.2195-GoingHome

    Going Home(回家)(标签:二部图,匈牙利算法,KM算法) 题目描述 在网格地图上,有n个男人和n个房屋. 在每个单位时间内,每个小人都可以水平或垂直移动一个单位步长到相邻的点. 对于每个小人 ...

  3. C++11新标准学习

    <深入理解C++11:C++11新特性解析与应用> <华章科技:深入理解C++11:C++11新特性解析与应用>一共8章:第1章从设计思维和应用范畴两个维度对C++11新标准中 ...

  4. C++11新标准:nullptr关键字

    一.nullptr的意义 1.NULL在C中的定义 #define NULL (void*)0 2.NULL在C++中的定义 #ifndef NULL #ifdef __cplusplus #defi ...

  5. C++11新标准:decltype关键字

    一.decltype意义 有时我们希望从表达式的类型推断出要定义的变量类型,但是不想用该表达式的值初始化变量(如果要初始化就用auto了).为了满足这一需求,C++11新标准引入了decltype类型 ...

  6. C++11新标准:auto关键字

    一.auto意义 编程时常常需要把表达式的值赋给变量,这就要求在声明变量的时候清楚地知道表达式的类型,然后要做到这一点并非那么容易.为了解决这个问题,C++11新标准引入了auto类型说明符,用它就能 ...

  7. c++11新标准for循环和lambda表达式

    :first-child { margin-top: 0px; } .markdown-preview:not([data-use-github-style]) h1, .markdown-previ ...

  8. 关注C++细节——C++11新标准之decltype的使用注意

    c++11新特性--decltype decltype是C++11加入的一个新的keyword,目的是选择并返回操作数的数据类型,重要的是,在此过程中编译器分析表达式并得到它的类型,却不实际计算表达式 ...

  9. 基于c++11新标准开发一个支持多线程高并发的网络库

    背景 新的c++11标准出后,c++语法得到了非常多的扩展,比起以往不论什么时候都要灵活和高效,提高了程序编码的效率,为软件开发者节省了不少的时间. 之前我也写过基于ACE的网络server框架,但A ...

随机推荐

  1. 设计模式----行为型模式之命令模式(Command Pattern)

    下面来自head first设计模式的命令模式一章节. 定义 将"请求"封装成对象,以便使用不同的请求.队列或者日志来参数化其他对象.命令模式也支持可撤销的操作. 类图 注: 1. ...

  2. BZOJ [Scoi2015]情报传递

    Description 奈特公司是一个巨大的情报公司,它有着庞大的情报网络.情报网络中共有n名情报员.每名情报员口J-能有 若T名(可能没有)下线,除1名大头日外其余n-1名情报员有且仅有1名上线.奈 ...

  3. Int类的129为什么转成byte就变成-127了?

    作为一个java开发人员,接触的基本都是上层的,都是以应用为主,根据业务实现功能,但今天无意间发现了一个小问题,int类型的129转成byte类型变成了-127,我知道是因为位数截取的原因,但是还没有 ...

  4. Python实现电子邮件的发送

    利用Python smtplib.SMTP类方法来实现电子邮件的发送. 列举SMTP对象常见的方法: sendmail(from, to ,msg[,mopts,ropts]) :将msg从from发 ...

  5. C# .Net Framework

    .NET Framework概述 .NET Framework是由微软开发,一个致力于敏捷软件开发.快速应用开发.平台无关性和网络透明化的软件开发平台. .NET Framework组成 .NET F ...

  6. 清理 Sketch缓存 Storyist 3.5.1中文破解版 for Mac

    Sketch这款软件可以用来做原型和UI设计,很多设计师和产品经理都在用,时间长了以后,Sketch会占据不少缓存空间,这时候可以试试 Sketch Cache Cleaner 这款软件清理历史文件和 ...

  7. 讲真,MySQL索引优化看这篇文章就够了

    本文主要讨论MySQL索引的部分知识.将会从MySQL索引基础.索引优化实战和数据库索引背后的数据结构三部分相关内容,下面一一展开. 一.MySQL——索引基础 首先,我们将从索引基础开始介绍一下什么 ...

  8. LWIP移植文件介绍

    在介绍文件之前首先介绍一下DMA描述符 stm32以太网模块接收/发送FIFO和内存之间的以太网传输是通过以太网DMA使用DMA描述符完成的,一共有两个描述符列表:一个用于接收,一个用于发送, 两个列 ...

  9. 微信小程序之页面引用utils中的js文件

    /* 只可使用相对路径 */const utils = require('../../../utils/util.js') Page({})

  10. 利用hash远程登陆系统

    有的时候当我们拿到系统管理员hash由于密码复杂度过高无法破解时候可以利用hash直接进行远程登录 我们用到Metasploit里面的模块 肉鸡为windwos server 2003 x32系统 1 ...