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

以下是C语言函数计算N!的代码：

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;
    }
}

以下是调用该函数计算S的代码：

#include <stdio.h>
#include <stdlib.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 sum(int n) //求和函数
{
    int i = 0; //循环变量
    double sum = 0; //和
    int temp = 1; //用来存储1+4!+7!+...+19!
    for (i = 1; i <= n; i++) {
        temp += fac(3 * i - 2); //计算1+4!+7!+...+19!
        sum += 1.0 / fac(temp); //计算1/（1+4!+7!+...+19!）
    }
    return sum;
}

int main()
{
    int n = 0;
    printf("请输入n:>");
    scanf("%d", &n);
    printf("sum = %lf\n", sum(n));
    system("pause");
    return 0;
}