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

A309132 a(n) is the denominator of F(n) = A027641(n-1)/n + A027642(n-1)/n^2.

Original entry on oeis.org

1, 1, 1, 16, 1, 36, 1, 64, 27, 100, 1, 144, 1, 196, 75, 256, 1, 324, 1, 400, 49, 484, 1, 576, 125, 676, 243, 784, 1, 900, 1, 1024, 363, 1156, 1225, 1296, 1, 1444, 169, 1600, 1, 1764, 1, 1936, 135, 2116, 1, 2304, 343, 2500, 867, 2704, 1, 2916, 3025, 3136, 361, 3364, 1, 3600, 1, 3844, 1323, 4096, 845, 4356, 1
Offset: 1

Views

Author

Thomas Ordowski, Jul 14 2019

Keywords

Comments

It seems that the numerator of F(n) is the numerator of (B(n-1) + 1/n), where B(k) is the k-th Bernoulli number; if so, for n > 2, the numerator of F(n) is A174341(n-1). How to prove it?
Conjecture: for n > 1, a(n) = 1 if and only if n is prime.
Is this conjecture equivalent to the Agoh-Giuga conjecture?
Theorem 1. If p is prime, then a(p) = 1. Proof. a(2) = 1, so let p be an odd prime. By the von Staudt-Clausen theorem, if k is even, then B(k) = A(k) - Sum_{prime q, q-1 | k} 1/q, where A(k) is an integer and the sum is over all primes q such that q-1 divides k. Thus B(k) = N(k)/D(k) with D(k) = Product_{prime q, q-1 | k} q. Now let k = p-1. Then N(p-1)/D(p-1) = B(p-1) = A(p-1) - 1/p - Sum_{prime q < p, q-1 | p-1} 1/q (*). Add 1/p to both sides of (*) and multiply by p*D(p-1) to get p*N(p-1) + D(p-1) = p*D(p-1)*(A(p-1) - Sum_{prime q < p, q-1 | p-1} 1/q) (**). Now p | D(p-1), so p^2 | p*D(p-1) in (**). The denominators on the right side of (**) are all of the form q < p. Therefore, p^2 divides both sides of (**). Hence F(p) = N(p-1)/p + D(p-1)/p^2 is an integer, so a(p) = 1. - Jonathan Sondow, Jul 14 2019
Conjecture: composite numbers n such that a(n) is squarefree are only the Carmichael numbers A002997. Cf. A309235. - Thomas Ordowski, Jul 15 2019
Conjecture checked up to n = 101101. - Amiram Eldar, Jul 16 2019
Theorem 2. If n is a prime or a Carmichael number, then a(n) = A326690(n) = denominator of (Sum_{prime p | n} 1/p - 1/n). The proof is a generalization of that of Theorem 1. (Note that Theorem 2 implies Theorem 1, since if n is prime, then (Sum_{prime p | n} 1/p - 1/n) = 1/n - 1/n = 0/1, so a(p) = A326690(n) = 1.) For n a prime or a Carmichael number, an application of Theorem 2 is computing a(n) without calculating Bernoulli(n-1) which may be huge; see A309268 and A326690. - Jonathan Sondow, Jul 19 2019
The values of F(n) when n is prime are A327033. - Jonathan Sondow, Aug 16 2019

Examples

			F(n) = 2/1, 0/1, 1/1, 1/16, 1/1, 1/36, 1/1, 1/64, 7/27, 1/100, 1/1, 1/144, -37/1, 1/196, 37/75, 1/256, -211/1, 1/324, 2311/1, 1/400, -407389/49, ...
		

Crossrefs

Programs

  • Magma
    [Denominator(Numerator(Bernoulli(n-1))/n + Denominator(Bernoulli(n-1))/n^2): n in [1..70]]; // Vincenzo Librandi, Jul 14 2019
  • Mathematica
    Table[Denominator[Numerator[BernoulliB[n - 1]] / n + Denominator[ BernoulliB[ n - 1]] / n^2], {n, 70}] (* Vincenzo Librandi, Jul 14 2019 *)
  • PARI
    a(n) = denominator(numerator(bernfrac(n-1))/n + denominator(bernfrac(n-1))/n^2); \\ Michel Marcus, Jul 14 2019
    

Formula

a(p) = 1 for prime p.
a(2k) = (2k)^2 for k > 1.
Conjecture: for k > 0, a(2k+1) = (2k+1)^2 iff 2k+1 is in A121707.
Denominator(F(p)/p) = 1 for the primes p = 2 and p = 1277 but for no other prime p < 1.5 * 10^4. Does denominator(F(p)/p) = 1 for any prime p > 1.5 * 10^4? - Jonathan Sondow, Jul 14 2019
Similarly, Sum_{k=1..p-1} k^(p-1) == -1 (mod p^2) for the prime p = 1277. - Thomas Ordowski, Jul 15 2019
a(n) = denominator(Sum_{prime p | n} 1/p - 1/n) if n is a prime or a Carmichael number. - Jonathan Sondow, Jul 19 2019

A174341 a(n) = Numerator of Bernoulli(n, 1) + 1/(n+1).

Original entry on oeis.org

2, 1, 1, 1, 1, 1, 1, 1, 7, 1, 1, 1, -37, 1, 37, 1, -211, 1, 2311, 1, -407389, 1, 37153, 1, -1181819909, 1, 76977929, 1, -818946931, 1, 277930363757, 1, -84802531453217, 1, 90219075042851, 1, -711223555487930419, 1, 12696640293313423, 1, -6367871182840222481, 1, 35351107998094669831, 1, -83499808737903072705023, 1, 12690449182849194963361, 1
Offset: 0

Views

Author

Paul Curtz, Mar 16 2010

Keywords

Comments

a(n) is numerator of (A164555(n)/A027642(n) + 1/(n+1)).
1/(n+1) and Bernoulli(n,1) are autosequences in the sense that they remain the same (up to sign) under inverse binomial transform. This feature is kept for their sum, a(n)/A174342(n) = 2, 1, 1/2, 1/4, 1/6, 1/6, 1/6, 1/8, 7/90, 1/10, ...
Similar autosequences are also A000045, A001045, A113405, A000975 preceded by two zeros, and A140096.
Conjecture: the numerator of (A164555(n)/(n+1) + A027642(n)/(n+1)^2) is a(n) and the denominator of this fraction is equal to 1 if and only if n+1 is prime or 1. Cf. A309132. - Thomas Ordowski, Jul 09 2019
The "if" part of the conjecture is true: see the theorems in A309132 and A326690. The values of the numerator when n+1 is prime are A327033. - Jonathan Sondow, Aug 15 2019

Crossrefs

Programs

  • Magma
    [2,1] cat [Numerator(Bernoulli(n)+1/(n+1)): n in [2..40]]; // Vincenzo Librandi, Jul 18 2019
  • Maple
    A174341 := proc(n) bernoulli(n,1)+1/(n+1); numer(%) end proc: # R. J. Mathar, Nov 19 2010
  • Mathematica
    a[n_] := Numerator[BernoulliB[n, 1] + 1/(n + 1)];
    Table[a[n], {n, 0, 47}] (* Peter Luschny, Jul 13 2019 *)
  • PARI
    B(n)=if(n!=1, bernfrac(n), -bernfrac(n));
    a(n)=numerator(B(n) + 1/(n + 1));
    for(n=0, 50, print1(a(n),", ")) \\ Indranil Ghosh, Jun 19 2017
    
  • PARI
    a(n)=numerator(bernpol(n, 1) + 1/(n + 1)); \\ Michel Marcus, Jun 26 2025
    
  • Python
    from sympy import bernoulli, Integer
    def a(n): return (bernoulli(n) + 1/Integer(n + 1)).numerator # Indranil Ghosh, Jun 19 2017
    

Extensions

Reformulation of the name by Peter Luschny, Jul 13 2019
Showing 1-2 of 2 results.