Если бы он полноценно заменял, тогда результат бы сходился. А он не сходится.
Ну асимптотически-то сходится. Полноценно -- это перебор только.
Для того, чтобы понять, где примерно наступает максимум плотности - достаточно наверное и этой формулы. А чтобы получить количество кортежей с погрешностью меньше процента - надо больше констант.
-- добавлено через 3 минуты --А вы понимаете почему интересующая область примерно от 0 до 1e21 — 1e23 ?
Я не хочу играть с вами в загадки. :)
Считайте, что на ваши вопросы "а вы понимаете..." ответ по умолчанию "нет, не понимаю".
-- добавлено через 3 минуты --Так о чём и речь. Претендуйте.
Зачем?
Я вроде бы всё рассказал и показал, ни одного секрета за пазухой не утаил. Или всё-таки какой-то момент остался непонятен?
Ну, хотелось бы увидеть изложение вычисления что вы там вычисляете не в виде "читайте тему, там всё есть" или скриптов pari/gp, а в виде ясных формул, от начала до конца. Но это так, если вдруг вам делать нечего будет :) Теорминимум, короче.