A096431 Denominator of (9*(n^4 - 2*n^3 + 2*n^2 - n) + 2)/(2*(2*n - 1)).
1, 3, 1, 7, 9, 11, 13, 3, 17, 19, 21, 23, 5, 27, 29, 31, 33, 7, 37, 39, 41, 43, 9, 47, 49, 51, 53, 11, 57, 59, 61, 63, 13, 67, 69, 71, 73, 15, 77, 79, 81, 83, 17, 87, 89, 91, 93, 19, 97, 99, 101, 103, 21, 107, 109, 111, 113, 23, 117, 119, 121, 123, 25, 127, 129, 131, 133
Offset: 1
Examples
1, 28/3, 38, 703/7, 1891/9, 4186/11, ...
Links
- G. C. Greubel, Table of n, a(n) for n = 1..2000
- Eric Weisstein's World of Mathematics, Magic Hexagon
- Index entries for linear recurrences with constant coefficients, signature (0,0,0,0,2,0,0,0,0,-1).
Programs
-
Magma
A096431:= func< n | Numerator((2*n-1)/5) >; [A096431(n): n in [1..60]]; // G. C. Greubel, Oct 14 2024
-
Mathematica
Table[ (9(n^4-2n^3+2n^2-n)+2)/(2(2n-1)),{n,80}]//Denominator (* or *) LinearRecurrence[{0,0,0,0,2,0,0,0,0,-1},{1,3,1,7,9,11,13,3,17,19},80] (* Harvey P. Dale, Aug 25 2021 *) -
PARI
a(n) = numerator((2*n-1)/5); \\ Altug Alkan, Apr 13 2018
-
PARI
first(n) = my(res = vector(n, i, 2*i - 1)); forstep(i = 3, n, 5, res[i]/=5); res \\ David A. Corneth, Apr 15 2018
-
SageMath
def A096431(n): return numerator((2*n-1)/5) [A096431(n) for n in range(1,61)] # G. C. Greubel, Oct 14 2024
Formula
Satisfies a linear recurrence with characteristic polynomial (1-x^5)^2.
G.f.: x*(1+x)*(1+2*x-x^2+8*x^3+x^4+8*x^5-x^6+2*x^7+x^8)/((1-x)^2*(1+x+x^2+x^3+x^4)^2). - R. J. Mathar, Mar 11 2011
Comments