OFFSET
0,2
MATHEMATICA
a7913[n_]:=a7913[n]=Times@@(#[[1]]^Mod[#[[2]], 2])&[Transpose[FactorInteger[n]]];
Map[a7913[Total[Map[a7913, Binomial[#, Range[0, #]]]]]&, Range[0, 50]] (* Peter J. C. Moses, Oct 28 2014 *)
PROG
(PARI) a(n) = core(sum(i=0, n, core(binomial(n, i)))); \\ Michel Marcus, Nov 13 2014
CROSSREFS
KEYWORD
nonn
AUTHOR
Vladimir Shevelev, Oct 28 2014
EXTENSIONS
More terms from Peter J. C. Moses, Oct 28 2014
STATUS
approved