OFFSET
1,2
COMMENTS
This sequence is the union of {1} and of three infinite and disjoint subsequences.
-> Numbers k divisible by 4 but not of the form 8q^2 or 8q(q+1) = {4, 12, 20, 24, 28, ...} (see A182834). For these numbers, the corresponding unique m = k/2 (see example for k = 4).
-> Even numbers k not divisible by 4 and of the form k = 2*A055792 = 2*q^2, q>1 in A001541 = {18, 578, ...}. For these numbers, the corresponding unique m = k/2 - 2 = q^2-2 (see example for k = 18)
-> Even numbers k not divisible by 4, that are in A060626 but not of the form k=2q^2-4 with q>1 in A001541 = {2, 34, 62, 98, 142, 194, ...} (A349496). For these numbers, the corresponding unique m = k/2 + 1 (see example for k = 2).
See A348692 for further information.
LINKS
Rick Mabry and Laura McCormick, Square products of punctured sequences of factorials, Gaz. Aust. Math. Soc., 2009, pages 346-352.
EXAMPLE
For k = 2, 2$ / 2! = 1^2, hence 2 is a term.
For k = 4, 4$ /1! = 288, 4$ / 3! = 48, 4$ / 4! = 12 but for m = 2, 4$ / 2! = 12^2, hence 4 is a term.
For k = 18 and m = 7, we have 18$ / 7! = 29230177671473293820176594405114531928195727360000000000000^2 and there is no other solution m, hence 18 is a term.
MATHEMATICA
q[n_] := Count[BarnesG[n + 2]/Range[n]!, _?(IntegerQ@Sqrt[#] &)] == 1; Select[Range[100], q] (* Amiram Eldar, Nov 20 2021 *)
PROG
(PARI) sf(n) = prod(k=2, n, k!); \\ A000178
isok(m) = my(s=sf(m)); #select(issquare, vector(m, k, s/k!), 1) == 1; \\ Michel Marcus, Nov 20 2021
CROSSREFS
KEYWORD
nonn
AUTHOR
Bernard Schott, Nov 20 2021
STATUS
approved