第三章 pascal算法讲义-第二讲.pdf

百度Pascal 吧公开培训教材-Pascal 培训课程算法讲义-第二讲 第二讲 预备算法 一、高精度算法 高精度，在OI 中是极为常见的一种运算极大数的算法。 高精度运算，是指参与运算的数(加数，减数，因子……）范围 大大超出了标准数据类型（longint 2^31，int64 2^63 等）能表示的范 围的运算。例如，求两个200 位的数的和。 这时，就要用到高精度算法了。 首先让我们来回忆一下小学时学习 的加减乘竖式，如图 高精度运算的原理就是以数组一位一位的储存两个操作数，模拟 各种竖式，最后得出用数组储存的结果数。 Part1.基础高精度 首先是高精度加法 其中a,b 是字符串类型的读入数，后接两个for 把他们转成数组 第1 页，共25 页 百度Pascal 吧公开培训教材-Pascal 培训课程算法讲义-第二讲 abn 和bbn。在数组中储存可以使用顺序或倒序，顺序的话进位的时 候会很难办，所以推荐使用倒序,就像这样 因为是倒序，所以这个数表示的就是3211 。 当然为了模块化我们可以把这个大数数组类型定义成一个新类 型 type bigint=array [0..max] of longint;// max 是可能的最大位数 知道了以上这些，那么高精度加法的基本规则也就很简单了 设c 为结果，a ，b 为操作数 c[i]=a[i]+b[i] 然而并没有这么简单，进位的问题也是需要考虑的 因为是十进制，c 中每一位只有可能是0~9 ，出现超过9 的很明 显需要进位。 Inc(c[i+1],c[i] div 10); c[i]:=c[i] mod 10; 这样就完成了这一位的进位操作：将整 10 送给高位，其余留下 来。 这样我们就得到了c 这一个结果，倒序输出即可。 最后我们还得考虑最高位是否有进位，如最高位进位了那么结果 的长度应该+1 高精度减法和高精度加法的规则是非常相似的 (c[i]=a[i]-b[i],设a 为被减数b 为减数)，这里就不再阐述了。 第2 页，共25 页 百度Pascal 吧公开培训教材-Pascal 培训课程算法讲义-第二讲 Ww 与ff 为被减数与减数, L1 为最终结果长度 主要讲一下借位的问题： 当a[i]b[i]时，很明显是不够减的。那么我们就要在a 的上一位 中借出一位 dec(a[i+1]);inc(a[i],10); 这样就够减了，请读者自己思考 为什么。 至于高精度乘法呢，观察一下竖式规律即可得到，a,b 中每两位 都要互相乘,a[i]*b[j]最终是要加到c[i+j-1]这个位。 第3 页，共25 页 百度Pascal 吧公开培训教材-Pascal 培训课程算法讲义-第二讲 这样的话最后就处理一下c 的进位问题，然后输出即可。 c 的长度请读者自己思考。 高精度除法难度不大，却较为复杂，鉴于其编写调试均需很长时 间，在竞赛范围不可能出现，故这里不做详细讨论。有兴趣的同学可 以自己试试编写。 Part2.时间复杂度分析 加法减法自然就是O(len); 乘法O(lena*lenb) Div 和mod 也是特别慢的 假如有两个1 万位的大数相乘，时间就要接近1s，这么大的 时间开销对于一般题目来说是无法承受的,so 时间超限。 所以我们就要对他们做上一些优化 Part3.压位 用数组储存高精度数，一个能存9 位数的longint 只存一位 数是不是有点浪费呢？ 现在我们就让longint 存下4 位数。 (即把longint 当成一个10000进制数) 如果存5 位的话就有可能会215 爆掉，请读者自己