(* b stands for
A001501 *) b[n_] := n!^2 Sum[2^(2k-n) 3^(k-n) (3(n-k))! HypergeometricPFQ[{k-n, k-n}, {3(k-n)/2, 1/2 + 3(k-n)/2}, -9/2]/(k! (n-k)!^2), {k, 0, n}]/6^n;
(* c stands for
A246599 *) c[n_] := c[n] = Binomial[2n-1, n] b[n] - Sum[ Binomial[2n-1, 2k] Binomial[2k, k] b[k] c[n-k], {k, 1, n-1}];
a[n_] := a[n] = c[n] + Sum[Binomial[2n-1, 2k-1] c[k] a[n-k], {k, 1, n-1}];