A194262 Largest prime that divides the n-th partition number p(n) but does not divide p(1)*p(2)*...*p(n-1), or 1 if none.

1, 2, 3, 5, 7, 11, 1, 1, 1, 1, 1, 1, 101, 1, 1, 1, 1, 1, 1, 19, 1, 167, 251, 1, 89, 29, 43, 13, 83, 467, 311, 23, 1, 1231, 41, 17977, 281, 1, 1, 127, 193, 2417, 71, 97, 1087, 241, 67, 7013, 631, 9283, 661, 53, 5237, 59, 227, 1019, 102359, 3251, 199, 409, 971

Offset

1,2

Comments

It appears that a(n) is prime for all n > 97. See A194259 and A194260 for additional comments and links.

Links

Alois P. Heinz and Giovanni Resta, Table of n, a(n) for n = 1..10000 (first 2000 terms from Alois P. Heinz)

Maple

with(combinat): with(numtheory):
b:= proc(n) option remember;
      `if`(n=1, {}, b(n-1) union factorset(numbpart(n)))
    end:
a:= n-> `if`(n=1, 1, max(1, (b(n) minus b(n-1))[])):
seq(a(n), n=1..120);  # Alois P. Heinz, Aug 21 2011

Mathematica

a[n_] := Complement[FactorInteger[PartitionsP[n]][[All, 1]], FactorInteger[Product[PartitionsP[k], {k, 1, n-1}]][[All, 1]]] /. {} -> {1} // Last; Table[a[n], {n, 1, 100}] (* Jean-François Alcover, Jan 28 2014 *)

Crossrefs

Cf. A000041, A071963, A194259, A194260, A194261.

Sequence in context: A075657 A039711 A194261 * A039710 A087173 A179335

Adjacent sequences: A194259 A194260 A194261 * A194263 A194264 A194265

Keyword

nonn,look

Author

Jonathan Sondow, Aug 21 2011

Status

approved