usqrt
stop = usqrt(val) + 1;
fact_lim = usqrt(start + TABSIZE + TABSIZE + 1);
fact_lim = usqrt(stop + 1);
static u_int64_t usqrt(u_int64_t);