|
|
A096362
|
|
Order in which prime factors first occur in the Lucas sequence.
|
|
6
|
|
|
2, 3, 7, 11, 29, 47, 19, 41, 199, 23, 521, 281, 31, 2207, 3571, 107, 9349, 2161, 211, 43, 307, 139, 461, 1103, 101, 151, 90481, 5779, 14503, 59, 19489, 2521, 3010349, 1087, 4481, 9901, 67, 63443, 71, 911, 103681, 54018521, 29134601, 79, 859, 1601, 3041
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,1
|
|
LINKS
|
|
|
MAPLE
|
L:= gfun:-rectoproc({a(n)=a(n-1)+a(n-2), a(0)=2, a(1)=1}, a(n), remember):
S:= {}: Res:= NULL:
for n from 0 to 100 do
P:= numtheory:-factorset(L(n)) minus S;
Res:= Res, op(sort(convert(P, list)));
S:= S union P;
od:
|
|
MATHEMATICA
|
PrimeFactors[n_Integer] := Flatten[ Table[ # [[1]], {1}] & /@ FactorInteger[n]]; L[n_] := Fibonacci[n - 1] + Fibonacci[n + 1]; pf = {}; f[n_] := Block[{p = PrimeFactors[L[n]]}, l = Length[p]; k = 1; While[k <= l, If[ Position[ pf, p[[k]]] == {}, AppendTo[ pf, p[[k]] ]]; k++ ]]; Do[ f[n], {n, 40}]; pf (* Robert G. Wilson v, Jul 01 2004 *)
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
|