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”).

A144437
Period 3: repeat [3, 3, 1].
15
3, 3, 1, 3, 3, 1, 3, 3, 1, 3, 3, 1, 3, 3, 1, 3, 3, 1, 3, 3, 1, 3, 3, 1, 3, 3, 1, 3, 3, 1, 3, 3, 1, 3, 3, 1, 3, 3, 1, 3, 3, 1, 3, 3, 1, 3, 3, 1, 3, 3, 1, 3, 3, 1, 3, 3, 1, 3, 3, 1, 3, 3, 1, 3, 3, 1, 3, 3, 1, 3, 3, 1, 3, 3, 1, 3, 3, 1, 3, 3, 1, 3, 3, 1, 3, 3
OFFSET
1,1
COMMENTS
The sequence is generated from numerators in the energy differences of the hydrogen spectrum: A005563(1), A061037(4), A061039(6), A061041(8), A061043(10), A061045(12), A061047(14), A061049(16), ...
Conjecture: a(n) is the separatix. See A045944.
Also the decimal expansion of the constant 3310/999. - R. J. Mathar, May 21 2009
Continued fraction expansion of A171417.
Greatest common divisor of (n+1)^2-1 and (n+1)^2+2. - Bruno Berselli, Mar 08 2017
FORMULA
a(n) = (7-4*cos(2*Pi*n/3))/3. - Jaume Oliver Lafont, Nov 23 2008
G.f.: x*(3 + 3*x + x^2)/((1 - x)*(1 + x + x^2)). - R. J. Mathar, May 21 2009
a(n) = 3/gcd(n,3). - Reinhard Zumkeller, Oct 30 2009
a(n) = denominator(n^k/3), where k>0 is an integer. - Enrique Pérez Herrero, Oct 05 2011
a(n) = gcd(T(n+1), T(2)) = A256095(n+1, 2), with the triangular numbers T = A000217, for n >= 1. - Wolfdieter Lang, Mar 17 2015
a(n) = a(n-3) for n>3; a(n) = A169609(n) for n>0. - Wesley Ivan Hurt, Jul 02 2016
E.g.f.: (1/3)*(7*exp(x) - 4*exp(-x/2)*cos(sqrt(3)*x/2) - 3). - G. C. Greubel, Aug 24 2017
From Nicolas Bělohoubek, Nov 11 2021: (Start)
a(n) = 9/(a(n-2)*a(n-1)).
a(n) = 7 - a(n-2) - a(n-1). See also A052901 or A069705. (End)
MAPLE
seq(op([3, 3, 1]), n=1..50); # Wesley Ivan Hurt, Jul 02 2016
MATHEMATICA
A144437[n_]:=Denominator[n/3]; Array[A144437, 100] (* Enrique Pérez Herrero, Oct 05 2011 *)
CoefficientList[Series[(3 + 3 x + x^2)/(1 - x^3), {x, 0, 120}], x] (* Michael De Vlieger, Jul 02 2016 *)
Table[Mod[2*n^2 + 1, 3, 1], {n, 1, 50}] (* G. C. Greubel, Aug 24 2017 *)
PROG
(PARI) a(n)=if(n%3, 3, 1) \\ Charles R Greathouse IV, Sep 28 2015
(Magma) &cat [[3, 3, 1]^^30]; // Wesley Ivan Hurt, Jul 02 2016
CROSSREFS
Numerators in the energy differences of the hydrogen spectrum: A005563(1), A061037(4), A061039(6), A061041(8), A061043(10), A061045(12), A061047(14), A061049(16), ...
Sequence in context: A088420 A103585 A154595 * A169609 A353527 A220670
KEYWORD
nonn,easy
AUTHOR
Paul Curtz, Oct 05 2008
EXTENSIONS
Edited by R. J. Mathar, May 21 2009
STATUS
approved