|
|
A047550
|
|
Numbers that are congruent to {5, 7} mod 8.
|
|
7
|
|
|
5, 7, 13, 15, 21, 23, 29, 31, 37, 39, 45, 47, 53, 55, 61, 63, 69, 71, 77, 79, 85, 87, 93, 95, 101, 103, 109, 111, 117, 119, 125, 127, 133, 135, 141, 143, 149, 151, 157, 159, 165, 167, 173, 175, 181, 183, 189, 191, 197, 199, 205, 207, 213, 215, 221, 223, 229, 231, 237, 239, 245, 247, 253, 255, 261
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,1
|
|
LINKS
|
|
|
FORMULA
|
a(1)=5, a(2)=7, a(3)=13; for n>3, a(n) = a(n-1)+a(n-2)-a(n-3). - Harvey P. Dale, Jun 04 2012
G.f.: x*(5+2*x+x^2) / ((1-x)^2*(1+x)). - Colin Barker, Aug 26 2016
Sum_{n>=1} (-1)^(n+1)/a(n) = Pi/8 - sqrt(2)*log(sqrt(2)+1)/4. - Amiram Eldar, Dec 19 2021
E.g.f.: 1 + 4*x*exp(x) - exp(-x). - David Lovler, Sep 02 2022
|
|
MAPLE
|
|
|
MATHEMATICA
|
With[{r8=8*Range[0, 40]}, Sort[Join[r8+5, r8+7]]] (* or *) LinearRecurrence[ {1, 1, -1}, {5, 7, 13}, 80] (* Harvey P. Dale, Jun 04 2012 *)
|
|
PROG
|
(PARI) Vec(x*(5+2*x+x^2)/((1-x)^2*(1+x)) + O(x^100)) \\ Colin Barker, Aug 26 2016
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn,easy
|
|
AUTHOR
|
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
|