Магия числ

[spamsink]

Вот простенькая программка без единой операции деления или вычисления остатка, без условных операторов или выражений, и без единого явно упомянутого простого числа, которая печатает простые числа вплоть до 67. А в числах с плавающей точкой типа double значащих бит всего 53, включая скрытый.

#include <stdio.h>
#include <math.h>

double magic(double f) {
    double ff = floor(f);
    double r = log(ff+6)*(log(ff+6)-1);
    return (f-ff)*floor(r)+f;
}

double C = 2.6358597414547913;

int main() {
    int i;
    for (i = 1; i < 20; ++i) {
        printf("%d\n", (int)C);
        C = magic(C);
    }
    return 0;
}

Сможете улучшить результат?

Comments

[juan_gandhi]

Фигасе.

[spamsink]

Фарида Фирузбэхт была голова.

[vak]

И правда магия.

[archaicos]

А нам нужно long double?

[spamsink]

Если с long double, logl() и константой 2.6358597414547912655L, то более или менее ожидаемо дотягивает до 73, но long double, если я правильно понимаю, есть не у всех. Поэтому лучше пользоваться double, чтобы можно было честно сравнивать.