OFFSET
1,2
COMMENTS
MAPLE
c:= proc(n) option remember; local k; if n=1 then 4
else for k from 1+c(n-1) while isprime(k) do od; k fi
end:
a:= n-> mul(mul(c(i)+c(j), i=1..j-1), j=2..n):
seq(a(n), n=1..10); # Alois P. Heinz, Jul 23 2017
MATHEMATICA
t = Table[If[PrimeQ[k], 0, k], {k, 1, 100}];
composite = Rest[Rest[Union[t]]] (* A002808 *)
f[j_] := composite[[j]]; z = 20;
v[n_] := Product[Product[f[k] + f[j], {j, 1, k - 1}], {k, 2, n}]
d[n_] := Product[(i - 1)!, {i, 1, n}] (* A000178 *)
Table[v[n], {n, 1, z}] (* A203530 *)
Table[v[n + 1]/v[n], {n, 1, z - 1}] (* A203532 *)
Table[v[n]/d[n], {n, 1, 20}] (* A203533 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Clark Kimberling, Jan 03 2012
EXTENSIONS
Name edited by Alois P. Heinz, Jul 23 2017
STATUS
approved