#include <algorithm>
#include <cstdlib>
#include <boost/accumulators/accumulators.hpp>
#include <boost/accumulators/statistics.hpp>
#include <boost/algorithm/minmax_element.hpp>
#include <boost/array.hpp>
#include <boost/timer.hpp>

using namespace boost;
using namespace std;

typedef array<int,100000000> dataarray;

inline void do_tests(dataarray& data)
{

  {
    timer t;
    dataarray::iterator resmin = std::min_element(data.begin(),data.end());
    dataarray::iterator resmax = std::max_element(data.begin(),data.end());
    cout<<"min_element,max_element:"<<t.elapsed()<<endl;
    cout<<"minVal="<<*resmin<<",maxVal="<<*resmax<<endl;
  }

  {
    boost::accumulators::accumulator_set<int, boost::accumulators::features<boost::accumulators::tag::min,boost::accumulators::tag::max> > acc;

    timer t;
    acc = std::for_each(data.begin(),data.end(),acc);
    cout<<"min_max: accumulator:"<<t.elapsed()<<endl;
    int minVal = boost::accumulators::extract_result< boost::accumulators::tag::min >(acc);
    int maxVal = boost::accumulators::extract_result< boost::accumulators::tag::max >(acc);
    cout<<"minVal="<<minVal<<",maxVal="<<maxVal<<endl;
  }

  {
    timer t;
    pair< vector<int>::iterator, vector<int>::iterator > result = boost::minmax_element(data.begin(), data.end());
    cout<<"min_max: minmax_element:"<<t.elapsed()<<endl;   
    cout<<"minVal="<<*result.first<<",maxVal="<<*result.second<<endl;
  }

  {
    int minVal = (std::numeric_limits<int>::max)();
    int maxVal = (std::numeric_limits<int>::min)();
    timer t;
    for(dataarray::iterator it = data.begin(); it!=data.end();it++)
    {
        minVal = (std::min)(minVal,*it);
        maxVal = (std::max)(maxVal,*it);
    }
    cout<<"min/max handcoded:"<<t.elapsed()<<endl;   
    cout<<"minVal="<<minVal<<",maxVal="<<maxVal<<endl;
  }
}

int main(int argc, char const* const* argv)
{
  dataarray data;

  cout << "sorted array:" << endl;

  for(unsigned int i = 0; i< data.size();i++)
  {
    data[i] = i;
  }

  do_tests(data);

  cout << "random array:" << endl;

  std::srand(0);

  for(unsigned int i = 0; i< data.size();i++)
  {
    data[i] = std::rand();
  }

  do_tests(data);

  return 0;
}