编写函数计算N!，调用该函数计算下式的值: 变量 S=1+1/（1+4!）+1/（1+4!+7!）+···+1/（1+4!+7!+…+19!）

以下是回答：

# 定义阶乘函数
def fac(n):
    if n == 0 or n == 1:
        return 1
    else:
        return n * fac(n-1)

# 定义求和函数
def sum(n):
    s = 0
    for i in range(n):
        s += 1 / fac(1 + 3 * i)
    return s

# 计算S的值
S = sum(6) # n=6时，19!已经超出了Python的最大整数范围，因此只计算到第6项
print(S) # 输出：1.0009945751278185