Sony GO FOR IT 2)実数の階乗

# coding: utf-8

# ライセンス１：
# このプログラムのライセンス２以前の部分は下記のURLに記載されているプログラムです。
# 下記のサイトのライセンス（CC-BY-SA）に従ってください。
# http://en.wikipedia.org/wiki/Lanczos_approximation

from cmath import *
 
# Coefficients used by the GNU Scientific Library
g = 7
p = [0.99999999999980993, 676.5203681218851, -1259.1392167224028,
     771.32342877765313, -176.61502916214059, 12.507343278686905,
     -0.13857109526572012, 9.9843695780195716e-6, 1.5056327351493116e-7]
 
def gamma(z):
    z = complex(z)
    # Reflection formula
    if z.real < 0.5:
        return pi / (sin(pi*z)*gamma(1-z))
    else:
        z -= 1
        x = p[0]
        for i in range(1, g+2):
            x += p[i]/(z+i)
        t = z + g + 0.5
        return sqrt(2*pi) * t**(z+0.5) * exp(-t) * x

# ライセンス２：このプログラムの以降の部分は、好きに使ってください。
print "Input a."
a = input()
print "%s! = %s"%(a,gamma(a+1).real)