A283461 Second-largest prime factor of 2^n - 1, if composite, or 1 otherwise.

1, 1, 3, 1, 3, 1, 5, 7, 11, 23, 7, 1, 43, 31, 17, 1, 19, 1, 31, 127, 89, 47, 17, 601, 2731, 73, 113, 1103, 151, 1, 257, 89, 43691, 127, 73, 223, 174763, 8191, 41, 13367, 337, 9719, 683, 631, 178481, 4513, 257, 127, 1801, 11119, 2731

Offset

2,3

Comments

For clarification: if the largest prime factor occurs more than once, then that prime factor is selected.

Links

Charles R Greathouse IV, Table of n, a(n) for n = 2..1206

Formula

a(n) = A006530(A000225(n)/A005420(n)).

Mathematica

a[n_] := If[PrimeQ[2^n-1], 1, Block[{f = FactorInteger[2^n-1]}, If[f[[-1, 2]] == 1, f[[-2, 1]], f[[-1, 1]]]]]; a /@ Range[2, 52] (* Giovanni Resta, Mar 08 2017 *)

Prog

(PARI) a(n)=my(f=factor(2^n-1), t=#f~); if(f[t, 2]>1, f[t, 1], if(t>1, f[t-1, 1], 1))

Crossrefs

Cf. A005420, A193615, A006530, A000225, A108974.

Sequence in context: A346241 A328573 A369066 * A146907 A322865 A238694

Adjacent sequences: A283458 A283459 A283460 * A283462 A283463 A283464

Keyword

nonn

Author

Charles R Greathouse IV, Mar 08 2017

Status

approved