double mylog10(bigint a) { double l = 0; while (!a.fits_sint_p()) { a /= 10; ++l; } l += log10(a.get_si()); return l; }