#!/usr/bin/env python
""" factorial.py -- n! and log2( n! ) """

from math import log
from gamma import Lgamma, Gamma

LOG2 = log(2.0)


def factorial( n ):
   return Gamma( n + 1 )


def log2_factorial( n ):
   return Lgamma( n + 1 ) / LOG2