cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Showing 1-2 of 2 results.

A103200 a(1)=1, a(2)=2, a(3)=11, a(4)=19; a(n) = a(n-4) + sqrt(60*a(n-2)^2 + 60*a(n-2) + 1) for n >= 5.

Original entry on oeis.org

1, 2, 11, 19, 90, 153, 712, 1208, 5609, 9514, 44163, 74907, 347698, 589745, 2737424, 4643056, 21551697, 36554706, 169676155, 287794595, 1335857546, 2265802057, 10517184216, 17838621864, 82801616185, 140443172858, 651895745267, 1105706761003, 5132364345954
Offset: 1

Views

Author

K. S. Bhanu and M. N. Deshpande, Mar 24 2005

Keywords

Comments

The original version of this question was as follows: Let a(1) = 1, a(2) = 2, a(3) = 11, a(4) = 19; for n = 1..4 let b(n) = sqrt(60 a(n)^2 + 60 a(n) + 1); for n >= 5 let a(n) = a(n-4) + b(n-2), b(n) = sqrt(60 a(n)^2 + 60 a(n) +1). Bhanu and Deshpande ask for a proof that a(n) and b(n) are always integers. The b(n) sequence is A103201.
This sequence is also the interleaving of two sequences c and d that can be extended backwards: c(0) = c(1) = 0, c(n) = sqrt(60 c(n-1)^2 + 60 c(n-1) +1) + c(n-2) giving 0,0,1,11,90,712,5609,... d(0) = 1, d(1) = 0, d(n) = sqrt(60 d(n-1)^2 + 60 d(n-1) +1) + d(n-2) giving 1,0,2,19,153,1208,9514,... and interleaved: 0,1,0,0,1,2,11,19,90,153,712,1208,5609,9514,... lim_{n->infinity} a(n)/a(n-2) = 1/(4 - sqrt(15)), (1/(4-sqrt(15)))^n approaches an integer as n -> infinity. - Gerald McGarvey, Mar 29 2005

References

  • K. S. Bhanu (bhanu_105(AT)yahoo.com) and M. N. Deshpande, An interesting sequence of quadruples and related open problems, Institute of Sciences, Nagpur, India, Preprint, 2005.

Crossrefs

Cf. A103201, A177187 (first differences).

Programs

  • Magma
    I:=[1,2,11,19,90]; [n le 5 select I[n] else Self(n-1)+8*Self(n-2)-8*Self(n-3)-Self(n-4)+Self(n-5): n in [1..30]]; // Vincenzo Librandi, Sep 28 2011
  • Maple
    a[1]:=1: a[2]:=2:a[3]:=11: a[4]:=19: for n from 5 to 31 do a[n]:=a[n-4]+sqrt(60*a[n-2]^2+60*a[n-2]+1) od:seq(a[n],n=1..31); # Emeric Deutsch, Apr 13 2005
  • Mathematica
    RecurrenceTable[{a[1]==1,a[2]==2,a[3]==11,a[4]==19,a[n]==a[n-4]+ Sqrt[60a[n-2]^2+60a[n-2]+1]},a[n],{n,40}] (* or *) LinearRecurrence[ {1,8,-8,-1,1},{1,2,11,19,90},40] (* Harvey P. Dale, Sep 27 2011 *)
    CoefficientList[Series[-x*(1 + x + x^2)/((x - 1)*(x^4 - 8*x^2 + 1)), {x, 0, 40}], x] (* T. D. Noe, Jun 04 2012 *)

Formula

Comments from Pierre CAMI and Gerald McGarvey, Apr 20 2005: (Start)
Sequence satisfies a(0)=0, a(1)=1, a(2)=2, a(3)=11; for n > 3, a(n) = 8*a(n-2) - a(n-4) + 3.
G.f.: -x*(1 + x + x^2) / ( (x - 1)*(x^4 - 8*x^2 + 1) ). Note that the 3 = the sum of the coefficients in the numerator of the g.f., 8 appears in the denominator of the g.f. and 8 = 2*3 + 2. Similar relationships hold for other series defined as nonnegative n such that m*n^2 + m*n + 1 is a square, here m=60. Cf. A001652, A001570, A049629, A105038, A105040, A104240, A077288, A105036, A105037. (End)
a(2n) = (A105426(n)-1)/2, a(2n+1) = (A001090(n+2) - 5*A001090(n+1) - 1)/2. - Ralf Stephan, May 18 2007
a(1)=1, a(2)=2, a(3)=11, a(4)=19, a(5)=90, a(n) = a(n-1) + 8*a(n-2) - 8*a(n-3) - a(n-4) + a(n-5). - Harvey P. Dale, Sep 27 2011

Extensions

More terms from Pierre CAMI and Emeric Deutsch, Apr 13 2005

A131709 Number of partitions into "bus routes" of an n X 1 grid.

Original entry on oeis.org

1, 14, 104, 904, 8004, 71004, 630004, 5590004, 49600004, 440100004, 3905000004, 34649000004, 307440000004, 2727910000004, 24204700000004, 214767900000004, 1905632000000004, 16908641000000004, 150030090000000004, 1331214490000000004, 11811844000000000004, 104806295100000000004, 929944511000000000004, 8251382159000000000004, 73214376480000000000004, 649629943210000000000004
Offset: 0

Views

Author

Yasutoshi Kohmoto, Oct 03 2007, revised Nov 20 2007

Keywords

Comments

If we make bus routes on a graph G, the routes should satisfy the following conditions.
1. One and only one route exists on all edges of G
2. Terminals of two different routes don't meet on the same point
This definition is equivalent to a "partition of graph G into undirected strokes". It is defined as follows.
Given an undirected graph G=(V,E), its partition into strokes is a collection of directed edge-disjoint paths (viewed as sets of directed edges) on V such that (i) union of any two paths is not a path; (ii)union of corresponding undirected paths is E.
So the case of undirected paths is the following.
Definition. Given an undirected graph G=(V,E), its partition into strokes is a collection of edge-disjoint paths (viewed as sets of edges) on V such that (i) union of any two paths is not a path; (ii) union of paths is E.
The first differences 90, 800, 7100, 63000, 559000,... are A177187 multiplied by powers of 10. - R. J. Mathar, Nov 02 2016

Crossrefs

Cf. A131518.

Programs

  • PARI
    Vec(-(30*x^3-30*x^2+3*x+1)/((x-1)*(10*x^2-10*x+1)) + O(x^100)) \\ Colin Barker, Feb 11 2015

Formula

a(n) = Product_{v_i} m_i + Sum_{c_j} (se_j - 1)*(Product_{v_k E (G_n-c_j)} m_k - {number of partitions of (G_n-c_i) which has cycles}) where:
v_i E V_n, G_n={V_n,E_n}, "E" means element
m_i means number of matching of incident edges of v_i
c_j means cycles in G_n
se_j means number of start-end points in c_j
v_k E G_n and not(v_k E c_j)
m_k means number of matching of incident edges of v_k
If (G_n-c_j) is empty then Product_{v_k E (G_n-c_j)} m_k = 1.
For n>=3, a(n)=10*(a(n-1)-a(n-2))+4. - Max Alekseyev, Apr 25 2013
G.f.: -(30*x^3-30*x^2+3*x+1) / ((x-1)*(10*x^2-10*x+1)). - Colin Barker, Feb 11 2015

Extensions

Terms a(4) onward from Max Alekseyev, Apr 25 2013
Showing 1-2 of 2 results.