OFFSET
1,2
COMMENTS
Odd terms of A051178.
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
EXAMPLE
The odd number 27 is a term because A000005(27!) = 2^9*3^3*7^2 is divisible by 27.
MAPLE
N:= 2000: # to get all terms <= N
P:= select(isprime, [2, seq(i, i=3..N, 2)]):
d:= Vector(nops(P), 1):
A:= 1:
for n from 2 to N do
f:= ifactors(n)[2];
for t in f do
if member(t[1], P, 'k') then
d[k]:= d[k] + t[2]
fi
od:
if n::odd and convert(d, `*`) mod n = 0 then A:= A, n fi;
od:
A; # Robert Israel, Aug 05 2016
MATHEMATICA
A275250Q = OddQ[#] && Divisible[DivisorSigma[0, #!], #] &; Select[Range[500], A275250Q] (* JungHwan Min, Jul 29 2016 *)
PROG
(PARI) isok(n) = (n % 2) && !(numdiv(n!) % n); \\ Michel Marcus, Jul 26 2016
CROSSREFS
KEYWORD
nonn
AUTHOR
Altug Alkan, Jul 21 2016
STATUS
approved