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

以下是C语言函数计算N!的代码和调用该函数计算S的代码：

#include <stdio.h>

int fac(int n) // 阶乘函数
{
    int i = 0;
    int y = 1; // 用来存返回值
    if (0 == n || 1 == n) {
        return 1;
    } else {
        for (i = 2; i <= n; i++) {
            y *= i;
        }
        return y;
    }
}

double calc_S() // 计算S的函数
{
    int i = 0, j = 0; // 循环变量
    double S = 0; // 和
    for (i = 0; i < 5; i++) { // 一共有5项
        int n = 1 + i * 3; // 每一项的n
        int f = fac(n); // 每一项的n!
        int denominator = 1; // 每一项的分母
        for (j = i; j >= 0; j--) { // 计算每一项的分母
            denominator += 3 * j;
        }
        S += 1.0 / (denominator * f); // 计算每一项的值并加到和上
    }
    return S;
}

int main()
{
    printf("S = %lf\n", calc_S());
    return 0;
}