cj.h

#ifndef  LVV_CODEJAM_H
#define  LVV_CODEJAM_H

#include <limits>
#include <vector>
#include <map>
#include <math.h>


//#ifndef LVV_MATRIX_H
//#include "scc/matrix.h"
//#endif

///////////////////////////////////////////////////////////////////// UTILS



/*
	template<typename T>
T gcd(T a, T b) {
    if (a < 0) a = -a;
    if (b < 0) b = -b;
    if (a == 0) return b;
    if (a > b) return gcd(b, a);
    return gcd(b % a, a);
}
*/
	// to check http://stackoverflow.com/questions/110344/algorithm-to-calculate-the-number-of-divisors-of-a-given-number
	//


////// HISTOGRAM

	template <typename Rg, typename T = typename Rg::value_type>
	std::vector<unsigned>
histogram_vector(Rg const& rg) {
	std::vector<unsigned>    result(std::numeric_limits<T>::max() - std::numeric_limits<T>::min());
	for(const auto& el : rg)   ++result[el - std::numeric_limits<T>::min()];
	return result;
}

	template <typename Rg, typename T = typename Rg::value_type>
	std::map<T, unsigned>
histogram(Rg const& rg) {
	std::map<T, unsigned> result;
	for(const auto& el : rg)  ++result[el];
	return result;
}


#define countof(array) (sizeof (array) / sizeof(array[0]))


#define scc_RO	1
#include <ro/ro.h>		// http://github.com/lvv/ro
#include <scc/simple.h>		// http://github.com/lvv/scc

#endif  // LVV_CODEJAM_H