Или заучите полдесятка (кажется их не больше практически полезных) формул преобразования логарифмов: сложение/вычитание, умножение/деление на число/константу, замена основания, определение. Вроде и всё.
Для оценки погрешностей полезно знать производные логарифма и экспоненты.
primepi считает очень медленно, для PARI/GP в разы быстрее будет тот финт что я применил выше: посчитать количество простых в диапазоне командой
forprime(p=a,b,n++) (для миллионного интервала после
![$10^{100}$ $10^{100}$](https://dxdy-01.korotkov.co.uk/f/4/0/f/40f92fd1553ca05a6465c8d79275a86a82.png)
выполнялся 6 секунд) и сравнить с разностью
![$Sf(a+b)-Sf(a)$ $Sf(a+b)-Sf(a)$](https://dxdy-04.korotkov.co.uk/f/7/f/d/7fd4852cb5f0496dbc3421dd40af5fb882.png)
(которую тоже неплохо бы математически упростить). Правда тут надо аккуратнее с количеством значащих цифр, если их не хватит можно получить бред (или отрицательное число, или ноль). Но PARI/GP позволяет указывать желаемую точность вычислений (командой
\p КоличествоЦифр), для чисел
![$\approx 10^{100}$ $\approx 10^{100}$](https://dxdy-03.korotkov.co.uk/f/2/a/e/2ae6fe32f39a215d947cf01cfe177e8882.png)
пришлось ставить не меньше 100 цифр (т.к.
![$Sf(10^{100})\approx 4\cdot 10^{97}$ $Sf(10^{100})\approx 4\cdot 10^{97}$](https://dxdy-04.korotkov.co.uk/f/3/3/7/3379b20b000a1d79d032e927f1accdae82.png)
, а надо чтобы была точной младшая цифра целой части). Ну это чисто от лени преобразовать выражение к нормальному виду чтобы не вычитать два близких стозначных числа.