网上找到的高精度求PI的超简短黑科技

还没弄懂什么意思，不过结果应该是对的。有没有巨巨能帮忙讲一下，谢谢了。

#include <stdlib.h>
#include <stdio.h>

#define BITS 2000

int a = 10000, b, c = BITS * 7 / 2, d, e, f[BITS * 7 / 2 + 1], g;

int main()
{
    for ( ; b - c; )
        f[b++] = a / 5;
    for ( ; d = 0, g = c * 2; c -= 14 , printf ("%.4d",e + d / a), e = d % a)
        for (b = c; d += f[b] * a, f[b] = d % -- g, d /= g--, --b; d *= b);
    return 0;
}