C语言编写一个计算s=1+1/（1+4！）+1/（1+4！+7!)+...+1/(1+4!+7!+...+19!)的程序

以下是C语言编写的计算s的程序：

#include <stdio.h>
#include <stdlib.h>

int main()
{
    int i, j, k;
    double s = 0, fac = 1;
    for (i = 0; i < 5; i++) {
        k = 1;
        for (j = 0; j <= i; j++) {
            k *= 3 * j + 1;
        }
        fac *= k;
        s += 1.0 / fac;    }
    printf("s = %lf\n", s + 1);
    return 0;
}