OFFSET
1,1
COMMENTS
Initial terms are 2, 2*7, 2*13, 2*19, 2*43, 2*47, 2*53, 2*107, 2*109.
MATHEMATICA
psig[n_]:=Sort[FactorInteger[n][[;; , 2]]]; Drop[With[{nn=220}, SequencePosition[psig/@ LinearRecurrence[ {2, 1}, {0, 1}, nn], {x_, x_}]][[;; , 1]]-1, 2] (* The program generates the first 9 terms of the sequence. *) (* Harvey P. Dale, Jul 07 2024 *)
PROG
(PARI) Pell(n) = ([2, 1; 1, 0]^n)[2, 1];
is(n) = vecsort(factor(Pell(n))[, 2]) == vecsort(factor(Pell(n+1))[, 2]);
CROSSREFS
KEYWORD
nonn,more
AUTHOR
Altug Alkan, Dec 31 2016
EXTENSIONS
a(10)-a(14) from Amiram Eldar, May 19 2024
STATUS
approved