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-5 of 5 results.

A189226 Curvatures in the nickel-dime-quarter Apollonian circle packing, ordered first by generation and then by size.

Original entry on oeis.org

-11, 21, 24, 28, 40, 52, 61, 157, 76, 85, 96, 117, 120, 132, 181, 213, 237, 376, 388, 397, 132, 156, 160, 189, 204, 205, 216, 237, 253, 285, 288, 309, 316, 336, 349, 405, 412, 421, 453, 460, 469, 472, 517, 544, 565, 616, 628, 685, 717, 741, 1084, 1093, 1104, 1125, 1128, 1140
Offset: 1

Views

Author

Jonathan Sondow, Apr 18 2011

Keywords

Comments

For a circle, curvature = 1/radius. The curvatures of a quarter, nickel, and dime are approximately proportional to 21, 24, and 28, respectively. Three mutually tangent circles with curvatures 21, 24, 28 can be inscribed in a circle of curvature 11.
Apollonius's and Descartes's Theorems say that, given three mutually tangent circles of curvatures a, b, c, there are exactly two circles tangent to all three, and their curvatures are a + b + c +- 2*sqrt(ab + ac + bc). (Here negative curvature of one of the two circles means that the three circles are inscribed in it.)
Fuchs (2009) says "An Apollonian circle packing ... is made by repeatedly inscribing circles into the triangular interstices in a Descartes configuration of four mutually tangent circles. Remarkably, if the original four circles have integer curvature, all of the circles in the packing will have integer curvature as well." That is because if a + b + c - 2s*qrt(ab + ac + bc) is an integer, then so is a + b + c + 2*sqrt(ab + ac + bc).
For n > 1, the n-th generation of the packing has 4*3^(n-2) circles.
Infinitely many of the curvatures are prime numbers A189227. In fact, in any integral Apollonian circle packing that is primitive (i.e., the curvatures have no common factor), the prime curvatures constitute a positive fraction of all primes (Bourgain 2012) and there are infinitely many pairs of tangent circles both of whose curvatures are prime (Sarnak 2007, 2011).
Fuchs and Sanden (2012) report on experiments with the nickel-dime-quarter Apollonian circle packing, which they call the coins packing P_C.

Examples

			The 1st-generation curvatures are -11, 21, 24, 28, the 2nd are 40, 52, 61, 157, and the 3rd are 76, 85, 96, 117, 120, 132, 181, 213, 237, 376, 388, 397. The 4th generation begins 132, 156, 160, 189, 204, 205, 216, ....
As 21 + 24 + 28 +- 2*sqrt(21*24 + 21*28 + 24*28) = 157 or -11, the sequence begins -11, 21, 24, 28, ... and 157 is in it.
The primes 157 and 397 are the curvatures of two circles that are tangent.
		

Crossrefs

Programs

  • Mathematica
    root = {-11, 21, 24, 28};
    triples = Subsets[root, {3}];
    a = {root};
    Do[
      ng = Table[Total@t + 2 Sqrt@Total[Times @@@ Subsets[t, {2}]], {t, triples}];
      AppendTo[a, Sort@ng];
      triples = Join @@ Table[{t, r} = tr; Table[Append[p, r], {p, Subsets[t, {2}]}], {tr, Transpose@{triples, ng}}]
      , {k, 3}];
    Flatten@a (* Andrey Zabolotskiy, May 29 2022 *)

Formula

a(n) == 0, 4, 12, 13, 16, or 21 (mod 24).

Extensions

Terms a(28) and beyond from Andrey Zabolotskiy, May 29 2022

A154637 a(n) is the ratio of the sum of squares of the bends of the circles that are added in the n-th generation of Apollonian packing, to the sum of squares of the bends of the initial three circles.

Original entry on oeis.org

1, 2, 66, 1314, 26082, 517698, 10275714, 203961186, 4048396578, 80356048002, 1594975770306, 31658447262114, 628384017931362, 12472705016840898, 247568948283023874, 4913960850609954786, 97536510167350024098, 1935988320795170617602, 38427156885401362279746, 762735172745641733742114
Offset: 0

Views

Author

Colin Mallows, Jan 13 2009

Keywords

Comments

For more references and links, see A189226.

Examples

			Starting with three circles with bends -1,2,2, the ssq is 9. The first derived generation has two circles, each with bend 3. So a(1) = (9+9)/9 = 2.
		

Crossrefs

For starting with four circles, see A137246. For sums of bends, see A135849 and A154636. For three dimensions, see A154638 - A154645.
Cf. also A189226, A189227.

Programs

  • Mathematica
    CoefficientList[Series[(29 z^2 - 18 z + 1)/(3 z^2 - 20 z + 1), {z, 0, 100}], z] (* and *) LinearRecurrence[{20, -3}, {1, 2, 66}, 100] (* Vladimir Joseph Stephan Orlovsky, Jul 03 2011 *)
  • PARI
    Vec((1-18*x+29*x^2)/(1-20*x+3*x^2) + O(x^30)) \\ Colin Barker, Nov 16 2016

Formula

G.f.: (1-18*x+29*x^2) / (1-20*x+3*x^2).
From Colin Barker, Nov 16 2016: (Start)
a(n) = ((133-13*sqrt(97))*(10+sqrt(97))^n - (10-sqrt(97))^n*(133+13*sqrt(97))) / (3*sqrt(97)) for n>0.
a(n) = 20*a(n-1) - 3*a(n-2) for n>2.
(End)

Extensions

More terms from N. J. A. Sloane, Nov 22 2009

A265189 Soddy circles: the two circles tangent to each of three mutually tangent circles.

Original entry on oeis.org

69, 46, 23, 6, 138, 70, 30, 21, 5, 105, 132, 33, 11, 4, -132, 138, 92, 46, 12, 276, 140, 60, 42, 10, 210, 153, 136, 72, 17, 306, 207, 138, 69, 18, 414, 210, 90, 63, 15, 315, 216, 135, 24, 10, -135, 238, 119, 102, 21, 357, 252, 63, 28, 9, 0
Offset: 1

Views

Author

Colin Barker, Dec 04 2015

Keywords

Comments

For any three mutually tangent circles (with radii a, b, and c), one can construct a fourth circle (the inner Soddy circle, with radius d) that is mutually tangent internally to the three circles, and a fifth circle (the outer Soddy circle, with radius e) that is mutually tangent externally to the three circles. For this sequence all five radii have integral lengths.
The sequence is an array of 5-tuples (a,b,c,d,e) ordered by increasing values of a, with a > b > c.
A positive value for the outer Soddy circle indicates that it contains the three circles; a negative value indicates that it is exterior to the three circles; a value of 0 indicates that it has an infinite radius, that is, it is a straight line.

Crossrefs

Cf. A256694.
See also the many sequences arising from Apollonian circle packing: A135849, A137246, A154636, etc.
Also the sequences related to Soddy's circle packings: A046159, A046160, A062536, etc.

Programs

  • PARI
    soddy(amax) = {
      my(L=List(), abc, t, u);
      for(a=1, amax,
        for(b=1, a-1,
          for(c=1, b-1,
            abc=a*b*c;
            if(issquare(abc*(a+b+c), &t),
              u=a*b+a*c+b*c;
              if(abc%(u+2*t) == 0,
                if(u-2*t != 0,
                  if(abc%(u-2*t) == 0,
                    listput(L, [a,b,c,abc\(u+2*t),-abc\(u-2*t)])
                  )
                ,
                  listput(L, [a,b,c,abc\(u+2*t),0])
                )
              )
            )
          )
        )
      );
      Vec(L)
    }
    soddy(253)

A154635 Ratio of the sum of the bends of the 5-dimensional spheres added in the n-th generation of Apollonian packing to the sum of the bends of the initial configuration of seven mutually tangent spheres.

Original entry on oeis.org

1, 2, 15, 108, 774, 5544, 39708, 284400, 2036952, 14589216, 104492016, 748400832, 5360254560, 38391631488, 274971524544, 1969422407424, 14105550112128, 101027866452480, 723589630947072, 5182549848861696, 37118861005211136, 265855588948518912
Offset: 0

Views

Author

Colin Mallows, Jan 13 2009

Keywords

Examples

			Starting with seven 5-dimensional spheres with bends 0,0,1,1,1,1,1 summing to 5, the first derived generation has seven spheres, with bends 1,1,1,1,1,5/2,5/2 summing to 10. So a(1) = 10/5 = 2.
		

Crossrefs

Cf. A135849 for dim=2. A137146 for the sum of squares of bends when dim=2. A154636 and A154637 for starting with three spheres in 2 dimensions. A154638-A154645 for results in the three-dimensional case.

Programs

Formula

G.f. (1-x)*(1-5*x) / (1-8*x+6*x^2).
From Colin Barker, Nov 16 2016: (Start)
a(n) = (((4-sqrt(10))^n*(-8+sqrt(10))+(4+sqrt(10))^n*(8+sqrt(10))))/(12*sqrt(10)) for n>0.
a(n) = 8*a(n-1) - 6*a(n-2) for n>2.
(End)

A256694 The radius of the largest of four circles with different integer radii arranged so that each circle is tangent externally to the other three circles.

Original entry on oeis.org

69, 70, 132, 138, 140, 153, 198, 207, 210, 216, 238, 252, 264, 264, 264, 270, 276, 280, 285, 290, 306, 345, 350, 390, 396, 396, 414, 420, 429, 432, 459, 476, 483, 490, 504, 504, 520, 528, 528, 528, 539, 540, 552, 560, 567, 570, 580, 594, 595, 612, 621, 630
Offset: 1

Views

Author

Colin Barker, Apr 08 2015

Keywords

Crossrefs

See also the many sequences arising from Apollonian circle packing: A135849, A137246, A154636, etc.
Also the sequences related to Soddy's circle packings: A046159, A046160, A062536, etc.

Programs

  • PARI
    soddy(k) = {
      s=[];
      for(a=1, k,
        for(b=1, a-1,
          for(c=1, b-1,
            if(issquare(a*b*c*(a+b+c), &t),
              if(a*b*c % (a*b+a*c+b*c+2*t) == 0,
                s=concat(s, a)
              )
            )
          )
        )
      );
      s
    }
    soddy(500)
Showing 1-5 of 5 results.