OFFSET
1,3
COMMENTS
Differs from A047318 (not congruent to 3 (mod 7)) starting at a(n=22) = 24 = 3*7 + 3 which is in this sequence but not in A047318. - M. F. Hasler, Jun 30 2025
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
Eric Weisstein's World of Mathematics, Refactorable Number
MAPLE
notref:= proc(n) option remember; n mod numtheory:-tau(n) <> 0 end proc:
filter:= proc(n)
andmap(t -> notref(t) or notref(n-t) or notref(n-2*t), [$1 .. (n-1)/2])
end proc:
select(filter, [$0..100]); # Robert Israel, Jul 29 2025
MATHEMATICA
Flatten[Table[If[Sum[(1 - Ceiling[(n - 2 i)/DivisorSigma[0, n - 2 i]] + Floor[(n - 2 i)/DivisorSigma[0, n - 2 i]]) (1 - Ceiling[i/DivisorSigma[0, i]] + Floor[i/DivisorSigma[0, i]]) (1 - Ceiling[(n - i)/DivisorSigma[0, n - i]] + Floor[(n - i)/DivisorSigma[0, n - i]]), {i, Floor[(n - 1)/2]}] == 0, n, {}], {n, 0, 100}]]
CROSSREFS
KEYWORD
nonn
AUTHOR
Wesley Ivan Hurt, Jul 27 2019
STATUS
approved
