OFFSET
0,2
COMMENTS
Least splitter is defined for x < y at A227631 as the least positive integer d such that x <= c/d < y for some integer c; the number c/d is called the least splitting rational of x and y. Conjecture: a(n) is the least splitter of s(n) and s(n+1), where s(n) = n*sin(1/n).
LINKS
Clark Kimberling, Table of n, a(n) for n = 0..1000
Leo Tavares, Illustration: Hexagonal Star Rays.
Index entries for linear recurrences with constant coefficients, signature (3,-3,1).
FORMULA
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3).
G.f.: (1 + 4*x + 7*x^2)/(1 - x)^3.
a(n) = A287326(2n, n). - Kolosov Petro, Nov 06 2017
From Amiram Eldar, Jul 15 2020: (Start)
Sum_{n>=0} 1/a(n) = (1 + (Pi/sqrt(6))*coth(Pi/sqrt(6)))/2.
Sum_{n>=0} (-1)^n/a(n) = (1 + (Pi/sqrt(6))*csch(Pi/sqrt(6)))/2. (End)
From Amiram Eldar, Feb 05 2021: (Start)
Product_{n>=0} (1 + 1/a(n)) = sqrt(2)*csch(Pi/sqrt(6))*sinh(Pi/sqrt(3)).
Product_{n>=1} (1 - 1/a(n)) = (Pi/sqrt(6))*csch(Pi/sqrt(6)).(End)
From Leo Tavares, Nov 20 2021: (Start)
a(n) = A033581(n) + 1. (End)
E.g.f.: exp(x)*(1 + 6*x + 6*x^2). - Stefano Spezia, Sep 14 2024
EXAMPLE
The first eight least splitting rationals for {n*sin(1/n), n >=1 } are these fractions: 6/7, 24/25, 54/55, 96/97, 150/151, 216/217, 294/295, 384/385.
MATHEMATICA
z = 40; r[x_, y_] := Module[{c, d}, d = NestWhile[#1 + 1 &, 1, ! (c = Ceiling[#1 x - 1]) < Ceiling[#1 y] - 1 &]; (c + 1)/d]; s[n_] := s[n] = n*Sin[1/n]; t = Table[r[s[n], s[n + 1]], {n, 1, z}] (* least splitting rationals *); fd = Denominator[t] (* Peter J. C. Moses, Jul 15 2013 *)
Array[6 #^2 + 1 &, 45] (* Michael De Vlieger, Nov 08 2017 *)
LinearRecurrence[{3, -3, 1}, {7, 25, 55}, 50] (* Harvey P. Dale, Dec 16 2017 *)
PROG
(PARI) a(n)=6*n^2+1 \\ Charles R Greathouse IV, Jun 17 2017
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Clark Kimberling, Jul 30 2013
EXTENSIONS
a(0) = 1 prepended by Robert P. P. McKone, Oct 09 2023
STATUS
approved