#include // for lround, round long lround(double x) { return round(x); }