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

A189227 Primes among the curvatures in the nickel-dime-quarter Apollonian circle packing A189226.

Original entry on oeis.org

-11, 61, 157, 181, 349, 373, 397, 421, 541, 661, 709, 733, 829, 853, 877, 997, 1021, 1069, 1093, 1213, 1237, 1381, 1429, 1597, 1621, 1669, 1693, 1741, 1861, 2029, 2221, 2293, 2341, 2389, 2557, 2677, 2749, 2917, 3037, 3061, 3109, 3181, 3229, 3253, 3301, 3373
Offset: 1

Views

Author

Jonathan Sondow, Apr 22 2011

Keywords

Comments

See A189226 for comments, references, links, examples, and crossrefs.

Crossrefs

Programs

  • Mathematica
    (* terms < 10^4 *) t = Range[9999]*0; w = {-11, 21, 24, 28}; s[1] = {{-1,2,2,2}, {0,1,0,0}, {0,0,1,0}, {0,0,0,1}}; s[2] = {{1,0,0,0}, {2,-1,2,2}, {0,0,1,0}, {0,0,0,1}}; s[3] = {{1,0,0,0}, {0,1,0,0}, {2,2,-1,2}, {0,0,0,1}}; s[4] = {{1,0,0,0}, {0,1,0,0}, {0,0,1,0}, {2,2,2,-1}}; r[m_, j_, p_] := Block[{v = (m.w)[[p]]}, If[v < 9999, t[[v]] = 1; Do[ If[i != j, r[m.s[i], i, p]], {i, 4}]]]; Do[ r[s[i], i, i], {i, 4}]; Prepend[ Select[ Flatten@ Position[t,1], PrimeQ], -11] (* Giovanni Resta, Jan 02 2014 *)

Formula

a(n) == 13 (mod 24) (because a(n) is prime, a(n) = A189226(k) for some k, and all terms of A189226 are == 0, 4, 12, 13, 16, or 21 (mod 24)).

Extensions

Corrected and extended by Steven Finch, Jan 02 2014
a(16)-a(46) from Giovanni Resta, Jan 02 2014

A248930 Decimal expansion of c = 2*Product_{prime p == 3 (mod 4)} (1 - 2/(p*(p-1)^2)), a constant related to the problem of integral Apollonian circle packings.

Original entry on oeis.org

1, 6, 4, 9, 3, 3, 7, 6, 8, 9, 0, 9, 8, 0, 3, 0, 7, 0, 1, 0, 2, 5, 9, 4, 2, 9, 3, 3, 3, 6, 0, 1, 7, 8, 9, 6, 3, 6, 6, 9, 2, 3, 5, 7, 6, 6, 2, 5, 6, 6, 1, 1, 4, 4, 9, 0, 5, 7, 7, 2, 4, 8, 8, 3, 8, 4, 2, 5, 6, 4, 5, 1, 8, 9, 4, 8, 0, 7, 7, 2, 5, 2, 0, 6, 9, 0, 2, 0, 4, 2, 4, 8, 5, 2, 5, 3, 6, 0, 1, 0, 2, 7, 0, 1, 7
Offset: 1

Views

Author

Jean-François Alcover, Oct 17 2014

Keywords

Examples

			1.64933768909803...
		

Crossrefs

Programs

  • Mathematica
    kmax = 25; Do[ P[k] = Product[p = Prime[n]; If[Mod[p, 4] == 3, 1 - 2/(p*(p - 1)^2) // N[#, 40]&, 1], {n, 1, 2^k}]; Print["P(", k, ") = ", P[k]], {k, 10, kmax}]; c = 2*P[kmax]; RealDigits[c, 10, 15] // First
    (* -------------------------------------------------------------------------- *)
    $MaxExtraPrecision = 1000; digits = 121;
    f[p_] := (1 - 2/(p*(p - 1)^2));
    coefs = Rest[CoefficientList[Series[Log[f[1/x]], {x, 0, 1000}], x]];
    S[m_, n_, s_] := (t = 1; sums = 0; difs = 1; While[Abs[difs] > 10^(-digits - 5) || difs == 0, difs = (MoebiusMu[t]/t) * Log[If[s*t == 1, DirichletL[m, n, s*t], Sum[Zeta[s*t, j/m]*DirichletCharacter[m, n, j]^t, {j, 1, m}]/m^(s*t)]]; sums = sums + difs; t++]; sums);
    P[m_, n_, s_] := 1/EulerPhi[m] * Sum[Conjugate[DirichletCharacter[m, r, n]] * S[m, r, s], {r, 1, EulerPhi[m]}] + Sum[If[GCD[p, m] > 1 && Mod[p, m] == n, 1/p^s, 0], {p, 1, m}];
    m = 2; sump = 0; difp = 1; While[Abs[difp] > 10^(-digits - 5) || difp == 0, difp = coefs[[m]]*P[4, 3, m]; sump = sump + difp; m++];
    RealDigits[Chop[N[2*Exp[sump], digits]], 10, digits - 1][[1]] (* Vaclav Kotesovec, Jan 16 2021 *)

Extensions

More digits from Vaclav Kotesovec, Jun 27 2020

A248938 Decimal expansion of beta = G^2*(2/3)*Product_{prime p == 3 (mod 4)} (1 - 2/(p*(p-1)^2)) (where G is Catalan's constant), a constant related to the problem of integral Apollonian circle packings.

Original entry on oeis.org

4, 6, 1, 2, 6, 0, 9, 0, 8, 6, 1, 3, 8, 6, 1, 3, 0, 3, 3, 2, 8, 5, 2, 9, 8, 4, 6, 4, 2, 4, 6, 0, 7, 5, 1, 5, 8, 0, 1, 3, 8, 3, 4, 4, 3, 7, 6, 5, 8, 8, 2, 0, 6, 3, 0, 0, 7, 0, 3, 9, 7, 7, 5, 1, 9, 0, 7, 1, 2, 8, 1, 6, 0, 7, 2, 2, 0, 7, 4, 9, 8, 3, 7, 9, 1, 0, 4, 2, 6, 0, 7, 2, 6, 2, 1, 4, 8, 0, 7, 2, 3, 1, 6, 3, 1, 6
Offset: 0

Views

Author

Jean-François Alcover, Oct 17 2014

Keywords

Examples

			0.4612609086138613...
		

Crossrefs

Programs

  • Mathematica
    kmax = 25; Clear[P]; Do[P[k] = Product[p = Prime[n]; If[Mod[p, 4] == 3 , 1 - 2/(p*(p - 1)^2) // N[#, 40]&, 1], {n, 1, 2^k}]; Print["P(", k, ") = ", P[k]], {k, 10, kmax}]; beta = Catalan^2*(2/3)*P[kmax]; RealDigits[beta, 10, 16] // First
    (* -------------------------------------------------------------------------- *)
    $MaxExtraPrecision = 1000; digits = 121;
    f[p_] := (1 - 2/(p*(p - 1)^2));
    coefs = Rest[CoefficientList[Series[Log[f[1/x]], {x, 0, 1000}], x]];
    S[m_, n_, s_] := (t = 1; sums = 0; difs = 1; While[Abs[difs] > 10^(-digits - 5) || difs == 0, difs = (MoebiusMu[t]/t) * Log[If[s*t == 1, DirichletL[m, n, s*t], Sum[Zeta[s*t, j/m]*DirichletCharacter[m, n, j]^t, {j, 1, m}]/m^(s*t)]]; sums = sums + difs; t++]; sums);
    P[m_, n_, s_] := 1/EulerPhi[m] * Sum[Conjugate[DirichletCharacter[m, r, n]] * S[m, r, s], {r, 1, EulerPhi[m]}] + Sum[If[GCD[p, m] > 1 && Mod[p, m] == n, 1/p^s, 0], {p, 1, m}];
    m = 2; sump = 0; difp = 1; While[Abs[difp] > 10^(-digits - 5) || difp == 0, difp = coefs[[m]]*P[4, 3, m]; sump = sump + difp; m++];
    RealDigits[Chop[N[2*Catalan^2/3 * Exp[sump], digits]], 10, digits - 1][[1]] (* Vaclav Kotesovec, Jan 16 2021 *)

Formula

beta = (G^2/3)*A248930, where G is Catalan's constant A006752.

Extensions

More digits from Vaclav Kotesovec, Jun 27 2020

A135849 a(n) is the ratio of the sum of the bends (curvatures) of the circles in the n-th generation of an Apollonian packing to the sum of the bends in the initial four-circle configuration.

Original entry on oeis.org

1, 5, 39, 297, 2259, 17181, 130671, 993825, 7558587, 57487221, 437222007, 3325314393, 25290849123, 192350849805, 1462934251071, 11126421459153, 84622568920011, 643601286982629, 4894942589100999, 37228736851860105, 283145067047577843, 2153474325825042429
Offset: 1

Views

Author

Colin Mallows, Mar 06 2008

Keywords

Comments

These ratios are independent of the starting configuration.
For more comments, references and links, see A189226.

Examples

			Starting with the configuration with bends (-1,2,2,3) with sum(bends) = 6, the next generation contains four circles with bends 3,6,6,15. The sum is 30 = 6*a(2). The third generation has 12 circles with sum(bends) = 234 = 6*a(3).
		

Crossrefs

Programs

  • Magma
    I:=[1, 5, 39]; [n le 3 select I[n] else  8*Self(n-1) - 3*Self(n-2): n in [1..30]]; // Vincenzo Librandi, Dec 25 2012
  • Mathematica
    CoefficientList[Series[(2 z^2 - 3 z + 1)/(3 z^2 - 8 z + 1), {z, 0, 100}], z] (* and *) LinearRecurrence[{8, -3}, {1, 5, 39}, 100] (* Vladimir Joseph Stephan Orlovsky, Jul 03 2011 *)
  • PARI
    Vec((2*x^3 - 3*x^2 + x)/(3*x^2 - 8*x + 1)+O(x^99)) \\ Charles R Greathouse IV, Jul 03 2011
    

Formula

For n >= 4, a(n) = 8*a(n-1) - 3*a(n-2).
For n>2, [a(n+2), a(n+3)] = the 2 X 2 matrix [0,1; -3,8]^n * [5,39]. Example: [0,1; -3,8]^3 * [5,39] = [a(5), a(6)] = [2259, 17181]. - Gary W. Adamson, Mar 09 2008 (typo corrected by Jonathan Sondow, Dec 24 2012)
a(n) = floor(C * A138264(n)), where C = 1.057097576... = (1/2)*((1/9) + sqrt((1/81) + 4)). Example: a(7) = 130671 = floor(C * A138264(7)) = floor(C * 123613). A135849(n)/A138264(n) tends to C. - Gary W. Adamson, Mar 09 2008
O.g.f.: 2*x/3 +7/9 +(59*x-7)/(9*(1-8*x+3*x^2)). - R. J. Mathar, Apr 24 2008
a(n) = 31*sqrt(13)*(A^n - B^n)/234 - 7*(A^n + B^n)/18 for n>1 where A=3/(4-sqrt(13)) and B=3/(4+sqrt(13)). - R. J. Mathar, Apr 24 2008

A137246 a(n) is the ratio of the sum of the squares of the bends (curvatures) of the n-th generation of an Apollonian packing to the sum of the squares of the bends of the initial four-circle configuration.

Original entry on oeis.org

1, 17, 339, 6729, 133563, 2651073, 52620771, 1044462201, 20731381707, 411494247537, 8167690805619, 162119333369769, 3217883594978523, 63871313899461153, 1267772627204287491, 25163838602387366361, 499473454166134464747, 9913977567515527195857
Offset: 1

Views

Author

Colin Mallows, Mar 09 2008

Keywords

Comments

These ratios are independent of the starting configuration. Similar ratios of third and higher moments are not so independent.
See A189226 for additional comments, references and links.

Examples

			Starting with the configuration with bends (-1,2,2,3) with sum(bends^2) = 18, the next generation contains four circles with bends 3,6,6,15. The sum of their squares is 306 = 18*a(2). The third generation has 12 circles with sum(bends^2) = 6102 = 18*a(3).
		

Crossrefs

Programs

  • GAP
    a:=[1,17,339];; for n in [4..30] do a[n]:=20*a[n-1]-3*a[n-2]; od; a; # G. C. Greubel, May 24 2019
  • Magma
    R:=PowerSeriesRing(Integers(), 30); Coefficients(R!(x*(1-x)*(1-2*x)/(1-20*x+3*x^2)));  // Bruno Berselli, Jul 04 2011
    
  • Mathematica
    CoefficientList[Series[(2z^2-3z+1)/(3z^2-20z+1), {z, 0, 30}], z] (* and *) LinearRecurrence[{20, -3}, {1, 17, 339}, 30] (* Vladimir Joseph Stephan Orlovsky, Jul 03 2011 *)
  • PARI
    Vec(x*(1-2*x)*(1-x)/(1-20*x+3*x^2)+O(x^30)) \\ Charles R Greathouse IV, Jul 03 2011
    
  • Sage
    a=(x*(1-x)*(1-2*x)/(1-20*x+3*x^2)).series(x, 30).coefficients(x, sparse=False); a[1:] # G. C. Greubel, May 24 2019
    

Formula

For n >= 4, a(n) = 20*a(n-1) - 3*a(n-2).
O.g.f.: x*(1-x)*(1-2*x)/(1-20*x+3*x^2). - R. J. Mathar, Mar 31 2008
a(n) = ((41+sqrt(97))*(10+sqrt(97))^(n-1) - (41-sqrt(97))*(10-sqrt(97))^(n-1))/(6*sqrt(97)) for n>1. - Bruno Berselli, Jul 04 2011

A154636 a(n) is the ratio of the sum of the bends of the circles that are drawn in the n-th generation of Apollonian packing to the sum of the bends of the circles in the initial configuration of 3 circles.

Original entry on oeis.org

1, 2, 18, 138, 1050, 7986, 60738, 461946, 3513354, 26720994, 203227890, 1545660138, 11755597434, 89407799058, 679995600162, 5171741404122, 39333944432490, 299156331247554, 2275248816682962, 17304521539721034, 131610425867719386, 1000969842322591986
Offset: 0

Views

Author

Colin Mallows, Jan 13 2009

Keywords

Comments

For comments and more references and links, see A189226.

Examples

			Starting from three circles with bends -1,2,2 summing to 3, the first derived generation consists of two circles, each with bend 3. So a(1) is (3+3)/3 = 2.
		

Crossrefs

Other sequences relating to the two-dimensional case are A135849, A137246, A154637. For the three-dim. case see A154638 - A154645. Five dimensions: A154635.
Cf. also A189226, A189227.

Programs

  • Mathematica
    CoefficientList[Series[(5 z^2 - 6 z + 1)/(3 z^2 - 8 z + 1), {z, 0, 100}], z] (* and *) LinearRecurrence[{8, -3}, {1, 2, 18}, 100] (* Vladimir Joseph Stephan Orlovsky, Jul 03 2011 *)
  • PARI
    Vec((1 - x)*(1 - 5*x) / (1 - 8*x + 3*x^2) + O(x^30)) \\ Colin Barker, Jul 15 2017

Formula

G.f.: (1 - x)*(1 - 5*x) / (1 - 8*x + 3*x^2).
From Colin Barker, Jul 15 2017: (Start)
a(n) = ((-(-7+sqrt(13))*(4+sqrt(13))^n - (4-sqrt(13))^n*(7+sqrt(13)))) / (3*sqrt(13)) for n>0.
a(n) = 8*a(n-1) - 3*a(n-2) for n>2.
(End)

Extensions

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

A060790 Inscribe two circles of curvature 2 inside a circle of curvature -1. Sequence gives curvatures of the smallest circles that can be sequentially inscribed in such a diagram.

Original entry on oeis.org

-1, 2, 2, 3, 15, 38, 110, 323, 927, 2682, 7754, 22403, 64751, 187134, 540822, 1563011, 4517183, 13054898, 37729362, 109039875, 315131087, 910745750, 2632104062, 7606921923, 21984412383, 63536130986, 183622826522, 530679817859, 1533693138351, 4432455434478
Offset: 0

Views

Author

Brian Galebach, Apr 26 2001

Keywords

Comments

The ratio of successive terms approaches the constant phi+sqrt(phi) ~= 2.89005363826396..., where phi is the golden ratio (sqrt(5)+1)/2. The ratio between the curvatures of two successively smaller circles approaches this constant in any apollonian packing as the curvatures increase.
For more comments, references and links, see A189226.

Examples

			After circles of 2, 2, 3, 15 have been inscribed in the diagram, the next smallest circle that can be inscribed has a curvature of 38.
		

Crossrefs

Programs

  • GAP
    a:=[-1,2,2,3];; for n in [5..30] do a[n]:=2*a[n-1]+2*a[n-2]+2*a[n-3]-a[n-4]; od; a; # Muniru A Asiru, Sep 12 2018
  • Maple
    seq(coeff(series((x-1)*(1-3*x-3*x^2)/(1-2*x-2*x^2-2*x^3+x^4),x,n+1), x, n), n = 0 .. 30); # Muniru A Asiru, Sep 12 2018
  • Mathematica
    CoefficientList[Series[(-3 z^3 + 4 z - 1)/(z^4 - 2 z^3 - 2 z^2 - 2 z + 1), {z, 0, 100}], z] (* and *) LinearRecurrence[{2, 2, 2, -1}, {-1, 2, 2, 3}, 100] (* Vladimir Joseph Stephan Orlovsky, Jul 03 2011 *)
  • PARI
    { for (n=0, 200, if (n>3, a=2*a1 + 2*a2 + 2*a3 - a4; a4=a3; a3=a2; a2=a1; a1=a, if (n==0, a=a4=-1, if (n==1, a=a3=2, if (n==2, a=a2=2, a=a1=3)))); write("b060790.txt", n, " ", a); ) } \\ Harry J. Smith, Jul 12 2009
    

Formula

a(n) = 2a(n-1) + 2a(n-2) + 2a(n-3) - a(n-4).
G.f.: -(1-x)*(1 - 3*x - 3*x^2)/(1 - 2*x - 2*x^2 - 2*x^3 + x^4). - Colin Barker, Apr 22 2012
Lim_{n -> inf} a(n)/a(n-1) = A318605. - A.H.M. Smeets, Sep 12 2018

Extensions

Corrected by T. D. Noe, Nov 08 2006

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

A045673 Curvatures in diagram constructed by inscribing 2 circles of curvature 0 and 1 inside circle of curvature 0, continuing indefinitely to inscribe circles wherever possible.

Original entry on oeis.org

0, 1, 4, 9, 12, 16, 24, 25, 28, 33, 36, 40, 49, 52, 57, 60, 64, 72, 73, 76, 81, 84, 88, 96, 97, 100, 105, 108, 112, 121, 124, 129, 136, 144, 145, 148, 153, 156, 160, 168, 169, 172, 177, 180, 184, 192, 193, 196, 201, 204, 216, 217, 220, 225, 228, 232
Offset: 0

Views

Author

Keywords

Comments

Appears to be a superset of {A008784 - 1}. - Ralf Stephan, Jan 26 2005
See A189226 for additional comments, references, links, examples, and crossrefs. - Jonathan Sondow, Aug 24 2012

Crossrefs

A294510 Residues modulo 24 of curvatures in the Apollonian circle packing A042944.

Original entry on oeis.org

2, 3, 6, 11, 14, 15, 18, 23
Offset: 1

Views

Author

Jonathan Sondow, Nov 16 2017

Keywords

Comments

Fuchs and Sanden proved that all curvatures in the Apollonian circle packing -1, 2, 3, 6 are congruent mod 24 to either 2, 3, 6, 11, 14, 15, 18, or 23.
Additional comments and links are in A189226.

Crossrefs

Showing 1-10 of 10 results.