STL algorithm (1) : numeric part

STL algorithm (1) : numeric part
includes functions: accumulate, partial_sum, adjacent_difference, inner_product

1. accumulate

用途: 对于序列 X[n], 计算其中某个区段[first, last) 的统计特性.

具体来说对于计算 ∑fun(X_i)   i∈[first, lase) 很有用。

例如：X[n] = {1， 2， 3， 4， 5}；

我们需要计算  ∑（3X_i²-5X_i+10）  i∈[1, 5]

我们这样写代码：

//代码经过vc6.0 调试， 使用SGI-STL3.0 (此版本代码可阅读性比较好，侯先生的<<STL源码剖析>>就是用的此版 本)

#include <numeric>
#include<iostream.h>

template<class T>
struct myFun
{

 T operator()(const T& x, const T& y)const
 { return x + 3*y*y -5*y+10; }
};

using namespace  std;
void main()
{
 int data[5] = {1, 2, 3, 4, 5};
 cout<< accumulate(data, data+5, 0, myFun<int>())<<endl; //140
}

函数声名如下：

template <class InputIterator, class T>
T accumulate(InputIterator first, InputIterator last, T init);

template <class InputIterator, class T, class BinaryFunction>
T accumulate(InputIterator first, InputIterator last, T init,
             BinaryFunction binary_op);

 2.partial_sum 部分和

用途: 对于序列 X[n], 计算其中区段[first, last) 中逐个位置的累积特性.

该函数和 accumulate 有相似之处在于都是计算累积特性， 但是 accumulate 函数只有一个返回值， 而partial_sum 可以计算逐个位置的累积特性也。

例如：X[n] = {1， 2， 3， 4， 5}；

我们需要计算5个值：

∑（3X_i²-5X_i+10）  i∈[1, k] , k ∈[1, 5]

我们这样写代码：

#include <numeric>
#include<iostream>
#include <algorithm>

template<class T>
struct myFun
{

 T operator()(const T& x, const T& y)const
 { return x + 3*y*y -5*y+10; }
};

using namespace  std;
void main()
{
 int data[5] = {1, 2, 3, 4, 5};
 partial_sum (data, data+5, data, myFun<int>());
 copy (data, data+5, ostream_iterator<int>(cout, "* "));  //1* 13* 35* 73* 133*
 cout<<endl;
}

函数声名如下：

template <class InputIterator, class OutputIterator>
OutputIterator partial_sum(InputIterator first, InputIterator last,
                           OutputIterator result);

template <class InputIterator, class OutputIterator, class BinaryOperation>
OutputIterator partial_sum(InputIterator first, InputIterator last,
                           OutputIterator result, BinaryOperation binary_op);

 3. adjacent_difference 相邻差

用途: 对于序列 X[n], 计算区段[first, last) 中相邻位置相关的特性.

该函数和 partial_sum 为互逆运算（两个函数的第一个声明）

例如：X[n] = {1， 2， 3， 4， 5}；

我们需要计算：

       y_i =x_i²-5x_i-1+10 ,    ( i∈[2, 5] , y₁=x₁)

我们这样写代码：

#include <numeric>
#include<iostream>
#include <algorithm>

template<class T>
struct myFun
{
 T operator()(const T& x, const T& y)const
 { return x*x - 5*y + 10; }
};

using namespace  std;
void main()
{
 int data[5] = {1, 2, 3, 4, 5};
 adjacent_difference (data, data+5, data, myFun<int>());
 copy (data, data+5, ostream_iterator<int>(cout, "* "));  //1* 9* 9* 11* 15*
 cout<<endl;
}
函数声名如下：

template <class InputIterator, class OutputIterator>
OutputIterator adjacent_difference(InputIterator first, InputIterator last, 
                                   OutputIterator result);

template <class InputIterator, class OutputIterator, class BinaryFunction>
OutputIterator adjacent_difference(InputIterator first, InputIterator last,
                                   OutputIterator result,
                                   BinaryFunction binary_op);

4.inner_product 内积

 用途: 计算两个序列的统计特性。

例如：

X[n] = {1, 2, 3, 4, 5}；

Y[n] = {5, 4, 3, 2 ,1};

我们需要计算：

       ∑(x_i²-5Y_i+10 ),    ( i∈[1, 5] )

我们这样写代码：

#include <numeric>
#include<iostream>
#include <algorithm>

template<class T>
struct myFun
{
 T operator()(const T& x, const T& y)const
 { return x*x - 5*y + 10; }
};

using namespace  std;
void main()
{
 int x[5] = {1, 2, 3, 4, 5};
 int y[5] = {5, 4, 3, 2 ,1};
 cout<<inner_product (x, x+5, y, 0, plus<int>(),myFun<int>()) //30
  <<endl;
}

简单小结如下：

1.  accumulate , inner_product 计算序列的一个统计特性数据;

     partial_sum, adjacent_difference 计算已知序列的若干特性。

2。partial_sum, adjacent_difference 中的 result 可以是 beging

3。pattial_sum , adjacent_difference 的第一个模板函数互逆