login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A075770
Let u(1)=u(2)=u(3)=1, u(n+3)=(n*u(n)+(n+1)*u(n+1)+(n+2)*u(n+2))/(n+3); sequence gives values of n such that u(n) is an integer.
1
1, 2, 3, 10, 14, 350, 702, 817, 1053, 1405, 31426, 44770, 86766, 120118, 216918, 284310, 343278, 347646, 426465, 511758, 767637, 1005030, 1282398, 1483326, 1610510, 1930666, 2434318, 2643990, 2883465, 3345210, 3367430, 3420670, 3647278
OFFSET
1,2
MATHEMATICA
a[1] = a[2] = a[3] = 1; a[n_] := a[n] = ((n - 3)*a[n - 3] + (n - 2)*a[n - 2] + (n - 1)*a[n - 1])/n; Flatten[ Position[ Table[a[n], {n, 1, 10^5}], _Integer]]
CROSSREFS
Sequence in context: A268878 A194544 A348475 * A135101 A108065 A187767
KEYWORD
nonn
AUTHOR
Benoit Cloitre, Oct 09 2002
EXTENSIONS
Edited, corrected and extended by Robert G. Wilson v, Oct 10 2002
a(14) omwards from Lambert Klasen (lambert.klasen(AT)gmx.net), Jul 27 2005
STATUS
approved