nov15-euler-constant

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

int main() {
    long long N = 1;
    N <<= 34;

    long long int border[35];
    double part[34];

    for (int i = 0; i <= 33; ++i)
        part[i] = 0.0;

    border[0] = 1;
    for (int i = 1; i <= 34; ++i)
        border[i] = border[i - 1] << 1;

    long long k = N;
    for (int i = 33; i >= 0; --i) {
        for (k; k > border[i]; --k)
            part[i] += 1.0 / (double) k;
    }

    double S = 0.0;

    for (int i = 33; i >= 0; --i)
        S += part[i];

    S += 1.0;

    double logn = log(N);

    double C = S - logn;

    printf("%1.10f\n", C);

    return 0;
}