Пример 11.22 реализует прямое решение, которое показывает, как следует писать простую обобщенную функцию в стиле STL. Для расчета расстояний между векторами я мог бы использовать функцию inner_product
, однако я не стал использовать функтор, потому что это неоправданно усложнило бы решение. Пример 11.23 показывает, как можно рассчитывать расстояние между векторами, применяя функтор и функцию inner_product
из заголовочного файла
.
#include
#include
#include
#include
using namespace std;
template
struct DiffSquared {
Value_T operator()(Value_T x, Value_T y) const {
return (x - y) * (x - y);
}
};
template
double vectorDistance(Iter_T first, Iter_T last, Iter2_T first2) {
double ret = inner_product(first, last, first2, 0.0L,
plus
return ret > 0.0 ? sqrt(ret) : 0.0;
}
int main() {
int v1[] = { 1, 5 };
int v2[] = { 4, 9 };
cout << "distance between vectors (1,5) and (4,9) is ";
cout << vectorDistance(v1, v1 + 2, v2) << endl;
}
Поскольку реализация функции inner_product()
может быть специально оптимизирована для вашей платформы и вашего компилятора, я предпочитаю ее использовать везде, где это возможно.
11.13. Реализация итератора с шагом
Имеются смежные числовые ряды и требуется обеспечить доступ к каждому n-му элементу.
В примере 11.24 представлен заголовочный файл, реализующий класс итератора с шагом.
#ifndef STRIDE_ITER_HPP
#define STRIDE_ITER_HPP
#include
#include
template
class stride_iter {
public:
// открытые имена, вводимые typedef
typedef typename std::iterator_traits
typedef typename std::iterator_traits
typedef typename std::iterator_traits
difference_type;
typedef typename std::iterator_traits
typedef std::random_access_iterator_tag iterator_category;
typedef stride_iter self;
// конструкторы
stride_iter() : m(NULL), step(0) {};
stride_iter(const self& x) : m(x.m), step(x.step) {}
stride_iter(Iter_T x, difference_type n) : m(x), step(n) {}
// операторы
self& operator++() { m += step; return *this; }
self operator++(int) { self tmp = *this; m += step; return tmp; }
self& operator+=(difference_type x) { m += x * step; return *this; }
self& operator--() { m -= step; return *this; }
self operator--(int) { self tmp = *this; m -= step; return trap; }
self& operator--(difference type x) { m -= x + step; return *this; }
reference operator[](difference_type n) { return m[n * step]; }
reference operator*() { return *m; }
// дружественные операторы
friend bool operator==(const self& x, const self& y) {
assert(x.step == y.step);
return x.m == y.m;
}
friend bool operator!=(const self& x, const self& y) {
assert(x.step == y.step);
return x.m != y.m;
}
friend bool operator<(const self& x, const self& y) {
assert(x.step == y.step);
return x.m < y.m;
}