матрица

#include <locale>
#include <vector>
#include <iostream>
#include <iterator>
#include <algorithm>
#include <functional>
#include <numeric>
 
/*
    В прямоугольной матрице найти максимальное среднее арифметическое
    среди средних арифметических каждой строки матрицы
    и вычесть его из всех элементов матриц, кроме элементов одной заданной строки 
*/
 
template <typename T, size_t Rows>
struct avg : public std::unary_function<const T (&)[Rows], double>
{
    T operator()(const T (&arr)[Rows]) const
    {
        return std::accumulate(arr, arr + Rows, T()) / Rows;
    }
};
 
template <typename T, size_t Cols, size_t Rows>
void go(T (&m)[Cols][Rows])
{
    std::vector<T> v(Cols);
    std::transform(
        &m[0][0], &m[0][0] + Cols*Rows, &m[0][0],
        std::bind2nd(std::minus<T>(), *std::max_element(v.begin(),
        std::transform(m, m + Cols, v.begin(), avg<T, Rows>())))
        );
}
 
template <typename T, size_t Rows>
struct X : public std::unary_function<const T (&)[Rows], void>
{
    X(std::ostream& os) : os_(os) {}
    
    void operator()(const T (&arr)[Rows]) const
    {
        typedef std::ostream_iterator<T> O;
        std::copy(arr, arr+Rows, O(os_, "\t"));
        os_ << std::endl;
    }
 
    std::ostream& os_;
};
 
template <typename T, size_t Cols, size_t Rows>
std::ostream& operator<< (std::ostream& os, const T (&m)[Cols][Rows])
{
    std::for_each(m, m+Cols, X<T, Rows>(os));
    return os;
}
 
int main()
{
    setlocale(LC_ALL, "");
 
    double m[3][4] = {
        {1,2,3,4},
        {1,-2,-3,4},
        {-1,-2,3,4}
        };
 
    std::cout << "Before:" << std::endl << m << std::endl;
    
    go(m);
 
    std::cout << "After:" << std::endl << m;
 
    return 0;
}