C语言高精度阶乘的问题

网上看到的这个算法，我在中间注释了不懂的地方，求大神梳理思路

#include <stdio.h>
#include <string.h>

int main()
{
    int a[3000];//用来存储阶乘后的数
    int n,i,j;
    int c = 0;//起初将进位设置为0
    int s;//每个位上的数
    scanf("%d",&n);
    //先将数组a所有的空间都初始化为0
    memset(a,0,sizeof(a));
    a[0] = 1;//1的阶乘
    //计算阶乘
        //*********这一段计算阶乘不懂啊***********
    for(i = 2; i <= n; i++)
    {
        for(j = 0; j < 3000; j++)
        {
            s = a[j] * i + c;
            a[j] = s % 10;
            c = s / 10;
        }
    }
    //去掉前面为0的
    for(i = 3000 - 1; i >= 0; i--)
    {
        if(a[i] != 0)
            break;
    }
    for(j = i; j >= 0; j--)
    {
        printf("%d",a[j]);
    }
}