C++程序设计:原理与实践(进阶篇)16.5 数值算法

简介:

16.5 数值算法


大多数的标准库算法都涉及处理数据管理问题:它们需要对数据进行拷贝、排序、查找等。但是,只有少数算法涉及数值计算。当我们需要进行计算时,这些数值算法就变得十分重要了,并且这些算法为我们在STL框架中编写数值算法提供了范例。

在STL标准库中只有四种数值算法:

数值算法

x = accumulate(b,e,i) 累加序列中的值;例如,对{a, b, c, d}计算i+a+b+c+d。结果x的类型与初始值i的类型一致

x = inner_product(b,e,b2,i) 将两个序列的对应元素相乘并将结果累加。例如,对{a, b, c, d}和{e, f, g, h}计算i+a·e+b·f+c·g+d·h。结果x的类型与初始值i的类型一致

r = partial_sum(b,e,r) 对一个序列的前n个元素进行累加,并将累加结果生成一个序列。例如,对{a, b, c, d}将生成{a, a+b, a+b+c, a+b+c+d}

r = adjacent_difference(b,e,b2,r) 对一个序列的相邻元素进行减操作,并将得到的差生成一个序列。例如,对{a, b, c, d}将生成{a, b-a, c-b, d-c}

 

这些算法可以在<numeric>中找到。我们将介绍前两个,如果你觉得有需要的话,可以自己去查阅其他两个的详细情况。

16.5.1 累积

accumulate()是最简单但最有用的数值算法。在其最简单的形式中,该算法将一个序列中的值进行累加:

 

给定初始值init,该算法将序列[f?irst:last)中的每个值加到init上,并将和返回。init通常被称为累加器。例如:

 

这段代码将打印15,即0+1+2+3+4+5(0是初始值)。显然,accumulate()能够被用于所有类型的序列:

 

结果(和)的类型与accumulate()用来保存累加器的变量的类型一致。这带来了一定的灵活性,这可能是十分重要的,例如:

 

在一些计算机上long的有效位数要比int更多。与int型相比,double能够表示更大范围的数,但可能精度更差。我们将在第24章中再讨论范围和精度在数值计算中所起的作用。

将保存结果的变量用作初始值是一种常见的指明累加器类型的方法:

 

记住:要初始化累加器并将accumulate()的结果保存在这个变量中。在本例中,s2在初始化之前就用作了初始化器,因此算法的结果是未定义的。我们将s3传递给accumulate()(传值方式,参见8.5.3节),但算法结果并未被保存,这只是浪费时间。

16.5.2 泛化accumulate()

基本的三参数accumulate()版本执行累加运算。但是,我们可能还想在序列上执行很多其他有用的运算,例如乘法和减法。为此,STL提供了另一个四参数的accumulate()版本,允许我们指定要执行的运算:

 

任何接受两个累加器类型实参的操作均能用于这一版本的accumulate()。例如:

 

这段代码将打印35.1384,即1.0×1.1×2.2×3.3×4.4(1.0为初始值)。这里提供的二元运算符multiplies<double>()是一个实现乘法运算的标准库函数对象;multiplies<double>实现double的乘法;multiplies<int>实现int的乘法,等等。还有一些其他的二元函数对象:plus(加法),minus(减法),divides,modulus(取余)。这些对象均在<functional>中定义(见附录C.6.2)。

注意,为了计算浮点数的积,初始值显然应设为1.0。

如sort()例子(见16.4.2节)中所示,我们常常对类对象中包含的数据更感兴趣,而不仅仅是普通内置类型。例如,给定物品的单位价格和单位数,我们可能想要计算所有物品的价值总和:

 

我们可以让accumulate的运算符从一个Record元素中抽取units,将其与单位价格相乘并加到累加器中:

 

我们在这里很“懒惰”,使用了函数而不是函数对象——这仅仅是为了展示可以这么做。我们倾向于优先选用函数对象:

如果需要在调用之间保存值;

或者,如果代码很短以致内联化会带来很大不同(至多是几个原语操作)。

基于第二个原因,我们应该在本例中选用函数对象。

试一试

定义一个vector<Record>,用你所选择物品的四个记录将其初始化,并用上面的函数计算物品的总价值。

16.5.3 内积

给定两个向量,将它们对应位置的元素相乘并将结果累加,这一运算称为向量的内积(inner product),内积在很多领域都十分有用(例如物理和线性代数,参见24.6节)。如果你更喜欢代码而不是文字,下面就是STL版本:

 

这段代码将内积的概念推广到任意元素类型的任何序列。以股票市场指数为例。在股票市场中,每一上市公司都会被分配一个“权重”。例如,在道琼斯工业指数中,我们看到Alcoa公司的最新权重为2.4808。为了获得当前的指数值,我们将每个公司的股票价格与其权重相乘,并将所得所有加权价格相加。显然,这就是价值和权重的内积。例如:

 

 

注意,inner_product()处理两个序列。但它只接受三个参数,对于第二个序列只描述了开始位置。算法假设第二个序列包含的元素个数要等于或多于第一个序列。如果这一假设不成立,将产生运行时错误。对于inner_product()而言,第二个序列包含更多元素是没问题的,那些“多余的元素”将简单地不予处理。

两个序列不需要具有相同的类型,元素类型也不必相同。为了展示这一点,我们使用了vector存储价格、用list存储权重。

16.5.4 泛化inner_product()

inner_product()可以像accumulate()那样泛化,但它需要两个额外参数:一个用于将累加器与新值组合起来(与accumulate()完全一样),另一个用于组合元素值对:

 

在16.6.3节中,我们将回到道琼斯的例子,给出一个更优雅的解决方案,其中就使用了这个泛化的inner_product()。

相关文章
|
14天前
|
小程序 编译器 Linux
C++ 异常原理:以一个小程序为例
作者在调查某个 bug 时涉及到 C++ 异常,借此机会以本文把 C++ 异常机制梳理清楚供大家参考。
|
2天前
|
负载均衡 算法 调度
负载均衡原理及算法
负载均衡原理及算法
8 1
|
6天前
|
Arthas 监控 算法
JVM工作原理与实战(二十五):堆的垃圾回收-垃圾回收算法
JVM作为Java程序的运行环境,其负责解释和执行字节码,管理内存,确保安全,支持多线程和提供性能监控工具,以及确保程序的跨平台运行。本文主要介绍了垃圾回收算法评价标准、标记清除算法、复制算法、标记整理算法、分代垃圾回收算法等内容。
18 0
JVM工作原理与实战(二十五):堆的垃圾回收-垃圾回收算法
|
11天前
|
机器学习/深度学习 自然语言处理 算法
机器学习算法原理与应用:深入探索与实战
【5月更文挑战第2天】本文深入探讨机器学习算法原理,包括监督学习(如线性回归、SVM、神经网络)、非监督学习(聚类、PCA)和强化学习。通过案例展示了机器学习在图像识别(CNN)、自然语言处理(RNN/LSTM)和推荐系统(协同过滤)的应用。随着技术发展,机器学习正广泛影响各领域,但也带来隐私和算法偏见问题,需关注解决。
|
12天前
|
机器学习/深度学习 算法 数据挖掘
【Python机器学习专栏】层次聚类算法的原理与应用
【4月更文挑战第30天】层次聚类是数据挖掘中的聚类技术,无需预设簇数量,能生成数据的层次结构。分为凝聚(自下而上)和分裂(自上而下)两类,常用凝聚层次聚类有最短/最长距离、群集平均和Ward方法。优点是自动确定簇数、提供层次结构,适合小到中型数据集;缺点是计算成本高、过程不可逆且对异常值敏感。在Python中可使用`scipy.cluster.hierarchy`进行实现。尽管有局限,层次聚类仍是各领域强大的分析工具。
|
12天前
|
机器学习/深度学习 算法 前端开发
【Python机器学习专栏】集成学习算法的原理与应用
【4月更文挑战第30天】集成学习通过组合多个基学习器提升预测准确性,广泛应用于分类、回归等问题。主要步骤包括生成基学习器、训练和结合预测结果。算法类型有Bagging(如随机森林)、Boosting(如AdaBoost)和Stacking。Python中可使用scikit-learn实现,如示例代码展示的随机森林分类。集成学习能降低模型方差,缓解过拟合,提高预测性能。
|
14天前
|
Linux 程序员 图形学
C++语言在现代软件开发中的应用与实践
C++语言在现代软件开发中的应用与实践
20 2
|
14天前
|
机器学习/深度学习 算法 数据挖掘
【视频】支持向量机算法原理和Python用户流失数据挖掘SVM实例(下)
【视频】支持向量机算法原理和Python用户流失数据挖掘SVM实例(下)
|
14天前
|
机器学习/深度学习 算法 搜索推荐
【视频】支持向量机算法原理和Python用户流失数据挖掘SVM实例(上)
【视频】支持向量机算法原理和Python用户流失数据挖掘SVM实例
|
14天前
|
编解码 JavaScript 前端开发
【专栏】介绍了字符串Base64编解码的基本原理和在Java、Python、C++、JavaScript及Go等编程语言中的实现示例
【4月更文挑战第29天】本文介绍了字符串Base64编解码的基本原理和在Java、Python、C++、JavaScript及Go等编程语言中的实现示例。Base64编码将24位二进制数据转换为32位可打印字符,用“=”作填充。文中展示了各语言的编码解码代码,帮助开发者理解并应用于实际项目。