OFFSET
1,1
COMMENTS
If squarefree numbers were accepted, then Carmichael numbers would be included. - Michel Marcus, Mar 13 2018
MAPLE
isA292815 := proc(n)
local pf, pfs;
pfs := ifactors(n)[2] ;
if nops(pfs) = 1 or issqrfree(n) then
return false;
end if;
for pf in pfs do
if modp(n-1, op(1, pf)^op(2, pf)-1) > 0 then
return false;
end if;
end do:
true ;
end proc:
for n from 1 do
if isA292815(n) then
print(n) ;
end if;
end do: # R. J. Mathar, May 02 2018
MATHEMATICA
fa[n_] := fa[n] = FactorInteger[n];
free[n_] := Union[Table[fa[n]〚i, 2〛, {i, Length[fa[n]]}]] == {1}
tes1[n_] := Union@Table[IntegerQ[(n - 1)/(fa[n][[i, 1]]^fa[n][[i, 2]] - 1)], {i, Length[fa[n]]}] == {True};
Select[1 + Range[3300200], ! free[#] && Length@fa[#] > 1 && tes1[#] &]
PROG
(PARI) isok(n) = {my(f = factor(n)); if ((#f~ > 1) && ! issquarefree(n), for (k=1, #f~, if ((n-1) % (f[k, 1]^f[k, 2] -1), return (0)); ); return (1); ); return (0); } \\ Michel Marcus, Mar 05 2018
CROSSREFS
KEYWORD
nonn
AUTHOR
José María Grau Ribas, Sep 24 2017
STATUS
approved