Получить квадратную матрицу по формулам - C++

#include <iostream>
#include <cmath>
using namespace std;
void main()
{
    const int n=5;
    unsigned long A[n][n];
    int i, j, k, f;
    for (i=0; i<n; i++)
        for (int j=0; j<n; j++)
        {
            f=1;
            k=j+1;
            do
            {
                f*=k;
                k--;
            }
            while (k>1);
                A[i][j]= pow (f, i+1);
        }
    for (i=0; i<n; i++)
    {
        for (int j=0; j<n; j++)
            cout<<A[i][j]<<"\t";
        cout<<"\n";
    }
}