login
This site is supported by donations to The OEIS Foundation.

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A083557 a(n) is the greatest prime factor of 3*a(n-1)+2. 1
3, 11, 7, 23, 71, 43, 131, 79, 239, 719, 127, 383, 1151, 691, 83, 251, 151, 13, 41, 5, 17, 53, 23, 71, 43, 131, 79, 239, 719, 127, 383, 1151, 691, 83, 251, 151, 13, 41, 5, 17, 53, 23, 71, 43, 131, 79, 239, 719, 127, 383, 1151, 691, 83, 251, 151, 13, 41, 5, 17, 53, 23 (list; graph; refs; listen; history; text; internal format)
OFFSET

1,1

COMMENTS

Conjecture: if a(1)=m then the sequence becomes cyclic, for any m.

Conjecture verified up to 25000000 by Jud McCranie, Jun 11 2003

LINKS

Table of n, a(n) for n=1..61.

FORMULA

G.f.: x*(3 + 11*x + 7*x^2 + 23*x^3 + 71*x^4 + 43*x^5 + 131*x^6 + 79*x^7 + 239*x^8 + 719*x^9 + 127*x^10 + 383*x^11 + 1151*x^12 + 691*x^13 + 83*x^14 + 251*x^15 + 151*x^16 + 13*x^17 + 41*x^18 + 2*x^19 + 6*x^20 + 46*x^21) / (1 - x^19) (conjectured). - Colin Barker, Jul 15 2017

MATHEMATICA

f[n_] := Flatten[Table[ #[[1]], {1}] & /@ FactorInteger[ 3n + 2 ]][[ -1]]; NestWhileList[f, 3, UnsameQ, All]

NestList[FactorInteger[3#+2][[-1, 1]]&, 3, 70] (* Harvey P. Dale, Feb 21 2013 *)

PROG

(PARI) lista(nn) = {print1(a = 3, ", "); for (n=1, nn, a = vecmax(factor(3*a+2)[, 1]); print1(a, ", "); ); } \\ Michel Marcus, Jul 15 2017

CROSSREFS

Cf. A031439, A031440, A031442, A082021, A082132, A084923.

Sequence in context: A288832 A164808 A153285 * A119324 A322364 A250034

Adjacent sequences:  A083554 A083555 A083556 * A083558 A083559 A083560

KEYWORD

nonn

AUTHOR

Yasutoshi Kohmoto, Jun 05 2003

STATUS

approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recent
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified April 22 22:38 EDT 2019. Contains 322380 sequences. (Running on oeis4.)