A264137
Largest prime factor of the n-th Pell number, A000129(n).
Original entry on oeis.org
2, 5, 3, 29, 7, 13, 17, 197, 41, 5741, 11, 33461, 239, 269, 577, 8297, 199, 179057, 59, 45697, 5741, 982789, 1153, 29201, 33461, 146449, 337, 44560482149, 269, 3272609, 665857, 52734529, 15607, 1800193921, 199, 1101341, 9369319, 4605197, 5521, 1746860020068409
Offset: 2
-
Table[FactorInteger[Fibonacci[n, 2]][[-1, 1]], {n, 25}] (* Alonso del Arte, Dec 10 2016 *)
FactorInteger[#][[-1,1]]&/@LinearRecurrence[{2,1},{2,5},60] (* Harvey P. Dale, Jun 08 2019 *)
-
a(n) = vecmax(factor(([2, 1; 1, 0]^n)[2, 1])[,1]); \\ Daniel Suteu, May 26 2022
A364820
a(n) is the smallest prime factor of the n-th Pell number A000129(n).
Original entry on oeis.org
2, 5, 2, 29, 2, 13, 2, 5, 2, 5741, 2, 33461, 2, 5, 2, 137, 2, 37, 2, 5, 2, 229, 2, 29, 2, 5, 2, 44560482149, 2, 61, 2, 5, 2, 13, 2, 593, 2, 5, 2, 1746860020068409, 2, 11437, 2, 5, 2, 3761, 2, 13, 2, 5, 2, 68480406462161287469, 2, 29, 2, 5, 2
Offset: 2
-
FactorInteger[#][[1, 1]] & /@ LinearRecurrence[{2, 1}, {2, 5}, 57] (* Amiram Eldar, Oct 21 2023 *)
A285314
Numbers k such that the k-th term of some (generalized) Lucas sequence has no primitive prime factor.
Original entry on oeis.org
1, 2, 3, 4, 5, 6, 7, 8, 10, 12, 13, 18, 30
Offset: 1
If (P, Q, D) = (1, -1, 5) (giving the Fibonacci sequence), U(12) = 144 = 2^4 * 3^2, while U(4) = 3 and U(6) = 8 = 2^3. Hence U(12) with (P, Q, D) = (1, -1, 5) has no primitive prime factor and 12 belongs to this sequence.
- Y. Bilu, G. Hanrot, P. M. Voutier, Existence of primitive divisors of Lucas and Lehmer numbers J. reine Angew. Math. 539 (2001), 75--122, a preprint version available from here.
- R. D. Carmichael, On the numerical factors of the arithmetic forms a^n +- b^n, Ann. of Math., 15 (1913), 30--70.
- P. M. Voutier, Primitive divisors of Lucas and Lehmer sequences, Math. Comp. 64 (1995), 869--888.
- M. Yabuta, A simple proof of Carmichael's theorem on primitive divisors, Fibonacci Quart., 39 (2001), 439--443.
Cf.
A086597,
A086600 (number of primitive prime factors in Fibonacci(n) and Lucas(n)).
A247250
Indices of Pell numbers having exactly one primitive prime factor.
Original entry on oeis.org
2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 18, 21, 24, 29, 30, 32, 33, 35, 38, 41, 42, 50, 53, 54, 56, 58, 59, 66, 69, 89, 90, 94, 95, 97, 99, 101, 104, 117, 118, 120, 135, 138, 160, 167, 181, 191, 210, 221, 237, 242, 247
Offset: 1
Pell(1) = 1, which has no prime factors, so 1 is not in this sequence.
Pell(4) = 12 = 2^2 * 3, but 2 is not a primitive prime factor, and 3 is the only primitive prime factor of Pell(4), so 4 is in this sequence.
Pell(5) = 29, which is a prime and the only primitive prime factor of itself, so 5 is in this sequence.
Pell(12) = 13860 = 2^2 * 3^2 * 5 * 7 * 11, but none of 2, 3, 5, 7 is a primitive prime factor, and 11 is the only primitive prime factor of Pell(12), so 12 is in this sequence.
Pell(14) = 80782 = 2 * 13^2 * 239, but neither 2 nor 13 is a primitive prime factor, and 239 is the only primitive prime factor of Pell(14), so 14 is in this sequence.
Pell(19) = 6625109 = 37 * 179057, both of which are primitive prime factors of Pell(19), so 19 is not in this sequence.
Cf.
A152012 (for Fibonacci numbers).
-
Select[Range[1000], PrimePowerQ[(1-Sqrt[2])^EulerPhi[#]*Cyclotomic[#, (1+Sqrt[2])/(1-Sqrt[2])]/GCD[Cyclotomic[#, (1+Sqrt[2])/(1-Sqrt[2])], # ]]&] - Eric Chen, Dec 12 2014
pell[n_] := pell[n] = ((1+Sqrt[2])^n-(1-Sqrt[2])^n )/(2*Sqrt[2]) // Round; primitivePrimeFactors[n_] := Cases[FactorInteger[pell[n]][[All, 1]], p_ /; And @@ (GCD[p, #] == 1 & /@ Array[pell, n-1])]; Reap[For[n=2, n <= 200, n++, If[Length[primitivePrimeFactors[n]] == 1, Print[n]; Sow[n]]]][[2, 1]] (* Jean-François Alcover, Dec 12 2014 *)
-
pell(n) = imag((1 + quadgen(8))^n);
isok(pf, vp) = sum(i=1, #pf, vecsearch(vp, pf[i]) == 0) == 1;
lista(nn) = {vp = []; for (n=2, nn, pf = factor(pell(n))[,1]; if (isok(pf, vp), print1(n, ", ")); vp = vecsort(concat(vp, pf),, 8););} \\ Michel Marcus, Nov 29 2014
Two incorrect terms (72 and 110) deleted by
Colin Barker, Nov 29 2014
Showing 1-4 of 4 results.
Comments