OFFSET
0,2
COMMENTS
The Berndt-type sequence number 7 for the argument 2Pi/7 defined by the relation: sqrt(7)*a(n) = s(1)*c(4)^(2*n) + s(2)*c(1)^(2*n) + s(4)*c(2)^(2*n), where c(j):=2*cos(2*Pi*j/7) and s(j):=2*sin(2*Pi*j/7). If we additionally defined the following sequences:
sqrt(7)*b(n) = s(2)*c(4)^(2*n) + s(4)*c(1)^(2*n) + s(1)*c(2)^(2*n),
sqrt(7)*c(n) = s(4)*c(4)^(2*n) + s(1)*c(1)^(2*n) + s(2)*c(2)^(2*n), and
sqrt(7)*a1(n) = s(1)*c(4)^(2*n+1) + s(2)*c(1)^(2*n+1) + s(4)*c(2)^(2*n+1),
sqrt(7)*b1(n) = s(2)*c(4)^(2*n+1) + s(4)*c(1)^(2*n+1) + s(1)*c(2)^(2*n+1),
sqrt(7)*c1(n) = s(4)*c(4)^(2*n+1) + s(1)*c(1)^(2*n+1) + s(2)*c(2)^(2*n+1), then the following simple relationships between elements of these sequences hold true: a(n)=c1(n), c(n+1)=a1(n), -a(n)-b(n)=b1(n), which means that the sequences a1(n), b1(n), and c1(n) are completely and in very simple way determined by the sequences a(n), b(n) and c(n). However the last one's satisfy the following system of recurrence equations: a(n+1) = 2*a(n) + b(n), b(n+1) = a(n) + 2*b(n) - c(n), c(n+1) = c(n) - b(n). We have b(n)=A215694(n) and c(n)=A215695(n).
From general recurrence relation: a(n) = 5*a(n-1) - 6*a(n-2) + a(n-3), i.e. a(n) = 5*(a(n-1)-a(n-2)) + (a(n-3)-a(n-2)) the following summation formula can be easily obtained: sum{k=3,..,n} a(k) = 5*a(n-1)-a(n-2)+a(0)-5*a(1). Hence in discussed sequence it follows that: sum{k=3,..,n} a(k) = 5*a(n-1) - a(n-2) - 14.
LINKS
G. C. Greubel, Table of n, a(n) for n = 0..1000
Roman Witula and Damian Slota, New Ramanujan-Type Formulas and Quasi-Fibonacci Numbers of Order 7, Journal of Integer Sequences, Vol. 10 (2007), Article 07.5.6
Index entries for linear recurrences with constant coefficients, signature (5, -6, 1).
FORMULA
G.f.: (1-2*x-x^2)/(1-5*x+6*x^2-x^3).
EXAMPLE
We have a(6) = 10*a(4)+a(1), a(5) = 11*(a(3)-a(1)), a(10)-a(4)+a(3)+a(1)+a(0) = 77*10^3, and a(11)-a(4)+a(3)-a(2)+a(0) = 25*10^4 = (5^6)*(2^4).
MATHEMATICA
LinearRecurrence[{5, -6, 1}, {1, 3, 8}, 50]
PROG
(PARI) x='x+O('x^30); Vec((1-2*x-x^2)/(1-5*x+6*x^2-x^3)) \\ G. C. Greubel, Apr 23 2018
(Magma) I:=[1, 3, 8]; [n le 3 select I[n] else 5*Self(n-1) - 6*Self(n-2) + Self(n-3): n in [1..30]]; // G. C. Greubel, Apr 23 2018
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Roman Witula, Aug 14 2012
STATUS
approved