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

A209885 Numbers such that A172470(n) = 9.

Original entry on oeis.org

1, 3, 6, 8, 12, 14, 17, 19, 23, 25, 28, 30, 32, 35, 37, 41, 43, 46, 48, 52, 54, 57, 59, 63, 65, 68, 70, 72, 75, 77, 81, 83, 86, 88, 92, 94, 97, 99, 103, 105, 109, 111, 114, 116, 120, 122, 125, 127, 131, 133, 136, 138, 142, 144, 148, 150, 153, 155, 159, 161
Offset: 1

Views

Author

Stephen Crowley, Mar 14 2012

Keywords

Crossrefs

Cf. A172470.

A209886 Numbers such that A172470(n) = 52.

Original entry on oeis.org

2, 7, 13, 18, 24, 29, 31, 36, 42, 47, 53, 58, 64, 69, 71, 76, 82, 87, 93, 98, 104, 110, 115, 121, 126, 132, 137, 143, 149, 154, 160, 165, 171, 176, 178, 183, 189, 194, 200, 205, 211, 217, 222, 228, 233, 239, 244, 246, 251, 257, 262, 268, 273, 279, 285, 290
Offset: 1

Views

Author

Stephen Crowley, Mar 14 2012

Keywords

Crossrefs

Cf. A172470.

A209887 Numbers such that A172470(n) = 61.

Original entry on oeis.org

4, 5, 9, 10, 11, 15, 16, 20, 21, 22, 26, 27, 33, 34, 38, 39, 40, 44, 45, 49, 50, 51, 55, 56, 60, 61, 62, 66, 67, 73, 74, 78, 79, 80, 84, 85, 89, 90, 91, 95, 96, 100, 101, 102, 106, 107, 108, 112, 113, 117, 118, 119, 123, 124, 128, 129, 130, 134, 135, 139, 140
Offset: 1

Views

Author

Stephen Crowley, Mar 14 2012

Keywords

Crossrefs

Cf. A172470.

A172468 Numbers such that A166986(n)/2 - A167389(n) = 1.

Original entry on oeis.org

50, 59, 111, 120, 181, 242, 251, 303, 312, 373, 434, 495, 504, 556, 565, 626, 687, 696, 748, 757, 818, 879, 940, 949, 1001, 1010, 1071, 1132, 1141, 1193, 1202, 1254, 1263, 1324, 1385, 1394, 1446, 1455, 1516, 1577, 1638, 1647, 1699, 1708, 1769, 1830, 1839
Offset: 1

Views

Author

Stephen Crowley, Feb 03 2010

Keywords

Comments

It is conjectured that the successive differences of this sequence, A172470, are limited to three numbers: 9, 52 and 61 where it is noted that 61 - 52 = 9, (52 + 2)/9 = 6, (61 + 2)/9 = 7 and we have lcm(9, 52, 61) = 28548 = 13^4 - 13 and 6 + 7 = 13.
From Travis Scott, Oct 16 2022: (Start)
Given sequences S(n), T(n) such that S'(n), T'(n) both ~ r for some real number r, if S(n) - T(n) converges to c then I(n) = floor(S(n)) - floor(T(n)) - floor(c) converges to the indicator [(r*n) mod 1 < c mod 1]. Take a new sequence 1(k) from the k-th n indicated by I(n). If S(n) - T(n) - floor(c) is [I] nonincreasing and [II] < 1 for all n > m, it is easy to see that the first differences of 1(k) for all values > m are capped at max(a, b) by the earliest unordered pair of positive and negative residues {r*a, r*b} == {-x, y} (mod 1)_[-1/2, 1/2) satisfying x + y < c mod 1, since they are jointly sufficient to map any interval [0, c mod 1 <= c_n < 1] back to itself modulo 1.
Rearranging terms, an equivalent statement of A172468 is 1(k) indicated by [((n + 2)/log(2)) mod 1 < c_n] with [III] c_n = (n + 2 - Im(W(n,-log(2)/2))/(2*Pi))/log(2) - 5/2. By the asymptotics of W_n [see for example Corless et al., s. 4] we have Im(W(n,z)) ~ 2*Pi*n + arg(z) - Pi/2 as n->oo, taking the entire negative real axis to [IV] Im(W(n,-r)) ~ 2*Pi*n + Pi/2 and [ibid., from Eq. 4.20] Im(W(n,-r)) - 2*Pi*n = C strictly increasing from n >= 0 for the slice [V] (-1/e, 0). [Empirically this holds to -0.93568951... due to the fact that W_n is a discrete assembly of branches.]
The conjectured closure of {9,52,61} then follows from plugging [IV] into [III] to get c = lim_{n->oo} c_n = 7/log(16) - 5/2 and by inspection of (n/log(2)) mod 1 for n <= 61, noting that -log(2)/2 satisfies [V] so that c_n satisfies [I] and [II]. Some basic pigeonholing further restricts the differences to three admissible runs [{52,9}, {61,61,9}, {61,61,61,9}] and issuing the sequence simplifies to evaluating ((n + 2)/log(2)) mod 1 once per run, comparing it to two constants, and measuring a low watermark for c_n as infrequently as possible.
Asymptotically, with {a, b} = {9/log(2), 52/log(2)} mod 1_[-1/2, 1/2), the three differences shuffle [0, c] by sending [c+a, c] to [0, -a] at n+9, [0, b] to [c-b, c] at n+52, and (b, c+a) to (-a, c-b) at n+61. 1254, 9892, 111768, 137237, 3194660, 11530771, 47096480, 208252803, 2084612060, 2581695828, 8931808997, 29473399808, 36320596745... is the subsequence of terms that transgress the asymptote, in that c < ((n + 2)/log(2)) mod 1 < c_n and a run that would be {61,61,61,9} splits to {52,9}, {61,61,9}.
Note finally that A172468 runs naturally from 52, 61, ... and that the offset to 50, 59, ... is an artifact of defining A166986 in terms of (n+2)/log(2) instead of n/log(2). (End)

Crossrefs

Programs

  • Maple
    [ListTools[SearchAll](1, [seq(round(evalf(floor((n+2)/ln(2))-2-(argument(exp(-(ln(2)+LambertW(n, -ln(sqrt(2))))/ln(2)))*ln(2)+Im(LambertW(n,-ln(sqrt(2)))))/(2*Pi*ln(2)))), n = 1 .. 10000)])]
  • Mathematica
    Select[Range@1888,Floor[(#+2)/Log@2]==Floor[Im@LambertW[#,-Log@2/2]/Log@4/Pi+7/2]&] (* Or, accelerated per comment: *) Module[{x=0,m,z=1},Flatten@Table[m=Mod[x/Log@2,1];If[m<145/2-201/Log@16||(mTravis Scott, Oct 16 2022 *)

Formula

a(n) ~ c*n, where c = 4*log(2)/(7 - 10*log(2)) = 40.4590949.... - Travis Scott, Oct 16 2022
Showing 1-4 of 4 results.