How to write this formula in a C program to output the closest approximation of pi. But you cannot use any functions from the "math.h" library.

p(n)=k=0n16k(8k+148k+428k+518k+61)