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.

User: Stephen Crowley

Stephen Crowley's wiki page.

Stephen Crowley has authored 14 sequences. Here are the ten most recent ones:

A182063 Leading coefficient (divided by n+1) of the denominator of the factorization of a certain integral of a truncated series converging to the Riemann zeta function.

Original entry on oeis.org

1, 1, 1, 2, 1, 2, 1, 2, 3, 6, 1, 6, 3, 2, 3, 12, 2, 12, 3, 4, 6, 12, 1, 12, 6, 4, 3, 12, 2, 12, 3, 20, 30, 12, 5, 60, 30, 20, 3, 60, 10, 60, 15, 4, 30, 60, 5, 60, 6, 20, 15, 60, 10, 12, 15, 20, 30, 60, 1, 60, 30, 20, 15, 12, 10, 60, 15, 20, 6, 60, 5, 60, 30, 4, 15
Offset: 1

Author

Stephen Crowley, Jun 19 2012

Keywords

Programs

  • Maple
    G:=N->1+N*(Li(N+1)-Li((N+1)^2))/(N+1)+sum(n/ln(n+1), n = 1 .. N-1); GFact:=proc (N) options remember, operator, arrow; return op(1, sort(convert(denom(factor(G(N))), list)))/(N+1) end proc; [seq(GFact(m), m = 1 .. 200)]

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

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

Author

Stephen Crowley, Mar 14 2012

Keywords

Crossrefs

Cf. A172470.

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

Author

Stephen Crowley, Mar 14 2012

Keywords

Crossrefs

Cf. A172470.

A172470 First differences of A172468.

Original entry on oeis.org

9, 52, 9, 61, 61, 9, 52, 9, 61, 61, 61, 9, 52, 9, 61, 61, 9, 52, 9, 61, 61, 61, 9, 52, 9, 61, 61, 9, 52, 9, 52, 9, 61, 61, 9, 52, 9, 61, 61, 61, 9, 52, 9, 61, 61, 9, 52, 9, 61, 61, 61, 9, 52, 9, 61, 61, 9, 52, 9, 61, 61, 61, 9, 52, 9, 61, 61, 9, 52, 9, 52, 9, 61, 61, 9, 52, 9, 61, 61, 61, 9
Offset: 1

Author

Stephen Crowley, Feb 04 2010

Keywords

Comments

The numbers in the sequence are conjectured to be restricted to 9,52, or 61. Note that 61-52 = 9, (52+2)/(61-52) = (52+2)/9 = 6, (61+2)/(61-52) = (61+2)/9 = 7 and we have lcm(9,52,61) =28548 = 13^4 - 13 and (52+2)/9 + (61+2)/9 = 6+7 = 13.

Crossrefs

Formula

a(n) = A172468(n+1) - A172468(n).

A172513 Complement of A167389.

Original entry on oeis.org

1, 4, 7, 11, 14, 17, 20, 24, 27, 30, 33, 37, 40, 43, 46, 50, 53, 56, 59, 63, 66, 69, 73, 76, 79, 82, 86, 89, 92, 95, 99, 102, 105, 108, 112, 115, 118, 121, 125, 128, 131, 134, 138, 141, 144, 147, 151, 154, 157, 161, 164, 167, 170, 174, 177, 180, 183, 187, 190, 193
Offset: 1

Author

Stephen Crowley, Feb 05 2010

Keywords

Comments

It appears that the sequence of first differences (A172515) consists of only 3's and 4's. - M. F. Hasler, Apr 11 2019

Crossrefs

Programs

  • Maple
    sort(convert((convert([seq(n, n = 1 .. 1000)], set) minus convert([seq(round(evalf((argument(exp(-(ln(2)+LambertW(n, -(1/2)*ln(2)))/ln(2)))*ln(2)+Im(LambertW(n, -(1/2)*ln(2))))/(2*Pi*ln(2)))), n = 1 .. 1000)], set)), list))
  • PARI
    A172513_upto(LIM=200)=setminus([1..LIM],A167389_upto(200)) \\ M. F. Hasler, Apr 11 2019

Extensions

Keyword "hard" removed by M. F. Hasler, Apr 11 2019

A172515 First differences of A172513.

Original entry on oeis.org

3, 3, 4, 3, 3, 3, 4, 3, 3, 3, 4, 3, 3, 3, 4, 3, 3, 3, 4, 3, 3, 4, 3, 3, 3, 4, 3, 3, 3, 4, 3, 3, 3, 4, 3, 3, 3, 4, 3, 3, 3, 4, 3, 3, 3, 4, 3, 3, 4, 3, 3, 3, 4, 3, 3, 3, 4, 3, 3, 3, 4, 3, 3, 3, 4, 3, 3, 3, 4, 3, 3, 3, 4, 3, 3, 3, 4, 3, 3, 4, 3, 3, 3, 4, 3, 3, 3, 4, 3, 3, 3, 4, 3, 3, 3, 4, 3, 3, 3, 4, 3, 3, 3, 4, 3
Offset: 1

Author

Stephen Crowley, Feb 05 2010

Keywords

Comments

This sequence appears to have only '3's and '4's. Indices of '4's seem to be given by A322408, but this is true only for n < 95: then the next term in A322408 is 95, but only a(96) = 4. - M. F. Hasler, Apr 19 2019

Crossrefs

Programs

Extensions

Keyword "hard" removed by M. F. Hasler, Apr 11 2019

A160273 Successive differences (divided by 3) of the average of twin prime pairs divided by 2 (A040040).

Original entry on oeis.org

1, 1, 2, 2, 3, 2, 5, 1, 5, 2, 5, 2, 1, 5, 2, 5, 2, 5, 6, 12, 2, 5, 10, 8, 5, 3, 4, 3, 25, 2, 1, 5, 4, 23, 2, 3, 2, 5, 10, 13, 8, 2, 2, 3, 18, 4, 5, 1, 20, 2, 8, 5, 4, 11, 14, 1, 9, 3, 8, 5, 9, 1, 4, 3
Offset: 1

Author

Stephen Crowley, May 07 2009

Keywords

Comments

From Eric Snyder, Jul 02 2021: (Start)
Also the successive differences of A002822.
Instances of a(n) = 1 correspond to prime quadruples. (End)

Crossrefs

Programs

  • Maple
    ZL := []; for p to 1000000 do if `and`(isprime(p), isprime(p+2)) then ZL := [op(ZL), ((p+2)^2-p^2)*(1/8)] end if end do; a := [seq((ZL[i+1]-ZL[i])*(1/3), i = 2 .. nops(ZL)-1)]
  • Mathematica
    Rest[(Differences[(Mean/@Select[Partition[Prime[Range[500]],2,1], #[[2]]- #[[1]]==2&])/2])/3] (* Harvey P. Dale, Jun 10 2014 *)

Extensions

Keyword "hard" deleted by Harvey P. Dale, Jun 10 2014

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

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

A164360 Period 3: repeat [5, 4, 3].

Original entry on oeis.org

5, 4, 3, 5, 4, 3, 5, 4, 3, 5, 4, 3, 5, 4, 3, 5, 4, 3, 5, 4, 3, 5, 4, 3, 5, 4, 3, 5, 4, 3, 5, 4, 3, 5, 4, 3, 5, 4, 3, 5, 4, 3, 5, 4, 3, 5, 4, 3, 5, 4, 3, 5, 4, 3, 5, 4, 3, 5, 4, 3, 5, 4, 3, 5, 4, 3, 5, 4, 3, 5, 4, 3, 5, 4, 3, 5, 4, 3, 5, 4, 3, 5, 4, 3, 5, 4, 3, 5, 4, 3, 5, 4, 3, 5, 4, 3, 5, 4, 3, 5, 4, 3, 5, 4, 3
Offset: 0

Author

Stephen Crowley, Aug 14 2009

Keywords

Comments

From Klaus Brockhaus, May 29 2010: (Start)
Continued fraction expansion of (32+sqrt(1297))/13.
Decimal expansion of 181/333. (End)

Crossrefs

Cf. A007877 (repeat 0,1,2,1), A068073 (repeat 1,2,3,2), A028356 (repeat 1,2,3,4,3,2), A130784 (repeat 1,3,2), A158289 (repeat 0,1,2,3,4,5,6,7,8,9,8,7,6,5,4,3,2,1).
Cf. A178566 (decimal expansion of (32+sqrt(1297))/13). [Klaus Brockhaus, May 29 2010]

Programs

Formula

a(n) = 4+(-1)^n*((1/2+I*sqrt(3)/6)*((1+I*sqrt(3))/2)^n+(1/2-I*sqrt(3)/6)*((1-I*sqrt(3))/2)^n). [Corrected by Klaus Brockhaus, Sep 17 2009]
a(n) = 4+(1/3)*sqrt(3)*sin(2*n*Pi/3)+cos(2*n*Pi/3). [Corrected by Klaus Brockhaus, Sep 17 2009]
a(n) = a(n-3) for n > 2, with a(0) = 5, a(1) = 4, a(2) = 3.
G.f.: (5+4*x+3*x^2)/((1-x)*(1+x+x^2)). [Klaus Brockhaus, Sep 17 2009]
E.g.f.: 4*exp(x)+(1/3)*sqrt(3)*exp(-(1/2)*x)*sin((1/2)*x*sqrt(3))+exp(-(1/2)*x)*cos((1/2)*x*sqrt(3)).
a(n) = 4 + A057078(n). - Wesley Ivan Hurt, Jul 01 2016

Extensions

Edited by Klaus Brockhaus, Sep 17 2009
Offset changed to 0 and formulas adjusted by Klaus Brockhaus, May 18 2010