a(n)= N(n)/D(n) with N(n) := sfac(n-1)*sfac(2*n+1)/sfac(n+1)= A089501(n) and D(n) := gcd([seq((2*m+1)*((m*(m+1))^n)*N(n)/((n+m+1)!*(n-m)!), m=1..n)]), where sfac(n)=A000178(n) (superfactorials) and gcd(L) is the greatest common divisor >1 of a list of numbers L.