In the same vein, here is a blog post from 2011 that uses SAS to define and manipulate "BigIntegers" for computing the factorials of numbers when the factorial value exceeds the value that can be stored in a double (1.798E308). I could write a similar program for combinations. Notice that I use the LFACT function, which is the analog to the LCOMB function that Pierre and Data_Null use here. Like Fried Egg's post, this way of solving the problem is more for amusement than for practical scientidic computing, so I think the post uses all of the features that Art likes.
... View more
