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

A090495 Numbers k such that numerator(Bernoulli(2*k)/(2*k)) is different from numerator(Bernoulli(2*k)/(2*k*(2*k-1))).

Original entry on oeis.org

574, 1185, 1240, 1269, 1376, 1906, 1910, 2572, 2689, 2980, 3238, 3384, 3801, 3904, 4121, 4570, 4691, 4789, 5236, 5862, 5902, 6227, 6332, 6402, 6438, 6568, 7234, 7900, 8113, 8434, 8543, 8557, 8566, 9232, 9611, 9670, 9824, 9891, 9898, 10564, 10587, 10754, 11230, 11247, 11535, 11691, 11896, 12562, 12965, 13019, 13228, 13246, 13355, 13484, 13894, 14560, 14714, 14957, 15176, 15226, 15346, 15892, 16558, 16668, 16944, 17035, 17224, 17387, 17890, 18379, 18406, 18534, 18556, 18761, 19222, 19598, 19888, 20090
Offset: 1

Views

Author

N. J. A. Sloane, Feb 03 2004

Keywords

Comments

Michael Somos (Feb 01 2004) discovered the remarkable fact that A001067 is different from A046968, even though they agree for the first 573 terms.
Numbers n such that A001067 is different from A046968, or alternatively, those n such that gcd(A001067(n),2n-1) is > 1.
If gcd(A000367(n), A000367(n+2)) <>1 then n = A090495(n) - (3*A090496(n) + 1)/2. - Mohammed Bouayoun (bouyao(AT)wanadoo.fr), Feb 08 2004
So far, all terms correspond to irregular primes. Notice that these numbers are generated by n=((2k+1)p+1)/2 where p is an irregular prime and k is some integer = 1,2,... . In the Excel spreadsheet provided at the link, you will notice that much larger firstborn irregular primes p tend to produce smaller values of k. E.g., p = 691, 683, 653, k = 5, 15, 23. So by some guessing we could test a given large irregular prime for the first few values of k. I found ip's 257, 293, 311 this way, but not the index. Also the spreadsheet shows the corresponding irregular primes where the Bacher forecast fails for firstborn irregular prime. - Cino Hilliard, Feb 15 2004

Crossrefs

Programs

  • Maple
    a := n->numer(bernoulli(2*n)/(2*n)): b := n->numer(bernoulli(2*n)/(2*n*(2*n-1))): for n from 1 to 2000 do if a(n)<>b(n) then print(n,a(n)/b(n)); fi; od:
  • Mathematica
    a[n_] := Numerator[BernoulliB[2n]/(2n)] (* A001067 *); b[n_] := Numerator[BernoulliB[2n]/(2n(2n-1))] (* A046968 *); For[n=1, n <= 580, n++, If[ a[n] != b[n], Print[n, " ", a[n]/b[n]] ] ]
    k = 1; lst = {}; While[k < 38001, b = BernoulliB[2 k]; If [Numerator[b/(2 k)] != Numerator[b/(2 k (2 k - 1))], AppendTo[lst, k]; Print[{k}]]; k++ ]; lst (* Robert G. Wilson v, Aug 19 2010 *)
  • PARI
    bern2(c,m1,m2) = { for(n=m1,m2, n2=n+n; a = numerator(bernfrac(n2)/(n2)); \ A001067 b = numerator(a/(n2-1)); if(a <> b,print("A("c") = "n","a/b);c++) ) } \\ Cino Hilliard

Extensions

a(1)-a(7) from Michael Somos and W. Edwin Clark, Feb 03 2004
a(8)-a(9) from Robert G. Wilson v, Feb 03 2004
a(10)-a(12) from Eric W. Weisstein, Feb 03 2004
a(13)-a(39) from Cino Hilliard, Feb 03 2004
a(40) from Eric W. Weisstein, Feb 04 2004
Many further terms from Cino Hilliard, Feb 15 2004

A090496 Ratio of numerator(Bernoulli(2*n)/(2*n)) to numerator(Bernoulli(2*n)/(2*n*(2*n-1))) for n's for which they are different.

Original entry on oeis.org

37, 103, 37, 59, 131, 37, 67, 37, 283, 59, 37, 101, 691, 37, 67, 37, 59, 157, 37, 617, 37, 593, 67, 59, 103, 37, 37, 37, 59, 101, 67, 157, 37, 37, 149, 233, 59, 131, 37, 37, 683, 67, 37, 271, 59, 103, 37, 37, 67, 263, 37, 59, 307, 101, 37, 37, 577, 59, 67, 37, 653, 37, 37, 59, 103, 157, 37, 67, 37, 59, 131, 101
Offset: 1

Views

Author

N. J. A. Sloane, Feb 03 2004

Keywords

Comments

A001067(n) / A046968(n) when they are different, or alternatively, gcd(A001067(n),2n-1) when that number is > 1.
These numbers are always products of irregular primes (A000928).
All values yielding 37 are of the form 574+666*k, k=0,1,2,3,4,... and form thus an arithmetic progression with step 666=18*37=((37-1)/2)*37. All values yielding 59 are of the form 1269+1711*k, k=0,1,2,3 and 1711=28*59=((59-1)/2)*59. The two values yielding 67 are at distance 2211=((67-1)/2)*67. Conjecture: all indices yielding a given prime p form an arithmetic progression of step ((p-1)/2)*p. See A092291. - Roland Bacher, Feb 04 2004
The positions where 37 occurs appear to coincide with A026352. - Mohammed Bouayoun, Feb 05 2004
Roland Bacher conjectures that values of n yielding the same quotient p form an arithmetic progression n0+d*k, where d = p(p-1)/2. Actual and conjectured values of n0 are in the sequence A092291.
Composite values do occur. An example is 2n = 272876, which yields a quotient of 37*59. This was found by tdn using the Kummer congruences and CRT: using the irregular pairs (37,32) and (59,44), we know that the following Diophantine equations must be solved for (k,l,m): 32+36*k = 44+58*l = 1+37*59*m. Some quotients are not possible, e.g., 37*67, 37*103. All quotients are the product of irregular primes A000928. Composite quotients imply there are missing terms in the arithmetic progression conjectured by Bacher. - T. D. Noe, Feb 12 2004

Crossrefs

Programs

  • Mathematica
    A090496 = {}; Do[ r = Numerator[ b = BernoulliB[2n]/(2n) ] / Numerator[ b/(2n-1) ]; If[ r > 1, Print[n, " ", r]; AppendTo[ A090496, r] ], {n, 1, 20000}]; A090496 (* Jean-François Alcover, Jan 24 2012 *)

Extensions

a(1)-a(7) from Michael Somos and W. Edwin Clark, Feb 03 2004
a(8), a(9) from Robert G. Wilson v, Feb 03 2004
a(10)-a(12) from Eric W. Weisstein, Feb 03 2004
a(13)-a(39) from Cino Hilliard, Feb 03 2004
a(40)-a(44) from Eric W. Weisstein, Feb 04 2004
Terms from a(45) onwards from David Wasserman, Dec 06 2005

A091721 Babylonian sexagesimal (base 60) expansion of 1/11.

Original entry on oeis.org

5, 27, 16, 21, 49, 5, 27, 16, 21, 49, 5, 27, 16, 21, 49, 5, 27, 16, 21, 49, 5, 27, 16, 21, 49, 5, 27, 16, 21, 49, 5, 27, 16, 21, 49, 5, 27, 16, 21, 49, 5, 27, 16, 21, 49, 5, 27, 16, 21, 49, 5, 27, 16, 21, 49, 5, 27, 16, 21, 49, 5, 27, 16, 21, 49, 5, 27, 16, 21, 49, 5, 27, 16, 21
Offset: 0

Views

Author

Jeppe Stig Nielsen, Feb 01 2004

Keywords

Comments

Period 5: repeat [5, 27, 16, 21, 49]. - Wesley Ivan Hurt, May 25 2024

Crossrefs

Programs

  • Mathematica
    RealDigits[ 1/11, 60, 75] [[1]] (* Robert G. Wilson v, Feb 02 2004 *)
    CoefficientList[Series[(5 + 27 x + 16 x^2 + 21 x^3 + 49 x^4)/(1 - x^5), {x, 0, 40}], x] (* Wesley Ivan Hurt, May 25 2024 *)

Formula

From Wesley Ivan Hurt, May 25 2024: (Start)
a(n+5) = a(n).
G.f.: (5+27*x+16*x^2+21*x^3+49*x^4)/(1-x^5). (End)
Showing 1-3 of 3 results.