Stirling estimation for binomial coefficient
ptime add year, conversion to and from struct tm