OFFSET
1,2
LINKS
Alois P. Heinz, Table of n, a(n) for n = 1..20000
EXAMPLE
From Jon E. Schoenfield, Aug 13 2017: (Start)
a(1) = 1 (which has no prime factors);
a(2) = 2 (a prime);
a(3) cannot be 3, because 2 and 3 have the same prime signature (each is a prime); however, the prime signature of 4 (the square of a prime) differs from that of 2, so a(3) = 4. (End)
MAPLE
s:= n-> sort(map(i-> i[2], ifactors(n)[2])):
a:= proc(n) option remember; local k; for k from
1+a(n-1) while s(k)=s(a(n-1)) do od; k
end: a(1):=1:
seq(a(n), n=1..80); # Alois P. Heinz, Mar 09 2018
MATHEMATICA
{1, 2} ~Join~ Select[Range[2, 80], Last /@ FactorInteger[#] != Last /@ FactorInteger[# - 1] &] (* Giovanni Resta, Aug 14 2017 *)
CROSSREFS
KEYWORD
easy,nonn
AUTHOR
Amarnath Murthy, May 07 2003
EXTENSIONS
More terms from James A. Sellers, May 19 2003
Incorrect term 76 removed by Alois P. Heinz, Mar 09 2018
STATUS
approved