50分，求助大佬！

#include <iostream>
using namespace std;
int j(int a){
    unsigned long long b=1;
    for(int i=1;i<=a;i++){
        b*=i;
    }
    return b;
}
int main(){
    int n;
    unsigned long long S=0,a;
    cin>>n;
    for(int i=1;i<=n;i++){
        a=j(i);
        S+=a;
    }
    cout<<S;
}