#include <vector>
#include "histogram.h"

using namespace std;

static
void
find_minmax(const vector<double>& numbers, double& min, double& max) {
    min = numbers[0];
    max = numbers[0];
    for (size_t i = 0; i < numbers.size(); i++){
        if (numbers[i] < min){
            min = numbers[i];
        }
        if (numbers[i] > max){
            max = numbers[i];
        }
    }
}

vector <double> make_histogram(const vector<double> numbers, size_t bin_count){

    double min, max;
    find_minmax(numbers, min, max);
    double bin_size = (max - min) / bin_count;
    vector <double> bins;
    bins.resize(bin_count);

    for (std::size_t i = 0; i < numbers.size(); i++)
    {
        bool found = false;
        for (std::size_t j = 0; (j < bin_count - 1) && !found; j++)
        {
            auto lo = min + j * bin_size;
            auto hi = min + (j + 1) * bin_size;
            if ((lo <= numbers[i]) && (numbers[i] < hi))
            {
                bins[j]++;
                found = true;
            }
        }
        if (!found)
        {
            bins[bin_count - 1]++;
        }
    }
    double out;
    for (std::size_t i = 0; i < bin_count-1; i++)
    {
        if (bins[i] > bins[i+1])
        {
            out = bins[i];
            bins[i] = bins[i + 1];
            bins[i + 1] = out;
        }
    }
    return bins;
}