static long ldiv_round(long a, long b) {
  ret (a+b/2)/b;
}

static long ldiv_round(double a, double b) {
  ret lround(a/b);
}