A070191 -id:A070191 - cp's OEIS frontend

A069040 Numbers k that divide the numerator of B(2k) (the Bernoulli numbers).

1, 5, 7, 11, 13, 17, 19, 23, 25, 29, 31, 35, 37, 41, 43, 47, 49, 53, 59, 61, 65, 67, 71, 73, 77, 79, 83, 85, 89, 91, 95, 97, 101, 103, 107, 109, 113, 115, 119, 121, 125, 127, 131, 133, 137, 139, 143, 145, 149, 151, 155, 157, 161, 163, 167, 169, 173, 175, 179, 181

Offset: 1

Benoit Cloitre, Apr 03 2002

nonn

Equivalently, k is relatively prime to the denominator of B(2k). Equivalently, there are no primes p such that p divides k and p-1 divides 2k. These equivalences follow from the von Staudt-Clausen and Sylvester-Lipschitz theorems.

The listed terms are the same as those in A070191, but the sequences are not identical. (The similarity is mostly explained by the absence of multiples of 2, 3 and 55 from both sequences.) See A070192 and A070193 for the differences.

G. H. Hardy and E. M. Wright, An Introduction to the Theory of Numbers, 3rd ed., Oxford Univ. Press, 1954.

Amiram Eldar, Table of n, a(n) for n = 1..10000 (terms 1..300 from Seiichi Manyama)
I. Sh. Slavutskii, A note on Bernoulli numbers, Journal of Number Theory, Vol. 53, No. 2 (1995), pp. 309-310.

Cf. A027641, A070191, A070192, A070193.

A069040 := proc(n)
    option remember;
    if n=1 then
        1;
    else
        for k from procname(n-1)+1 do
            if numer(bernoulli(2*k)) mod k = 0 then
                return k;
            end if;
        end do:
    end if;
end proc: # R. J. Mathar, Jan 06 2013

testb[n_] := Select[First/@FactorInteger[n], Mod[2n, #-1]==0&]=={}; Select[Range[200], testb]

isok(k) = {my(p = factor(k)[,1]); for(i = 1, #p, if(!((2*k) % (p[i]-1)), return(0))); 1;} \\ Amiram Eldar, Apr 24 2025

More information from Dean Hickerson, Apr 26 2002

A070192 Numbers k such that k divides the numerator of B(2k) (the Bernoulli numbers), but gcd(3k, 8^k+1) > 3.

Original entry on oeis.org

301, 737, 1505, 1655, 2107, 3197, 3311, 3913, 5117, 5159, 5219, 5719, 6275, 6923, 7385, 7513, 7525, 8107, 8275, 8729, 9331, 9581, 9835, 10535, 10849, 11137, 11585, 12341, 12529, 12943, 13301, 14003, 14147, 14749, 15953, 15985, 17759, 18361

Offset: 1

Benoit Cloitre and Dean Hickerson, Apr 26 2002

nonn

Equivalently, numbers that are in A069040 but not in A070191.

Amiram Eldar, Table of n, a(n) for n = 1..10000

Cf. A027641, A069040, A070191, A070193.

testb[n_] := Select[First/@FactorInteger[n], Mod[2n, #-1]==0&]=={}; test8[n_] := GCD[3n, PowerMod[8, n, 3n]+1]==3; Select[Range[19000], testb[ # ]&&!test8[ # ]&]

isA070191(k) = gcd(3*k, Mod(8, 3*k)^k + 1) == 3;
isok(k) = {my(p = factor(k)[,1]); for(i = 1, #p, if(!((2*k) % (p[i]-1)), return(0))); !isA070191(k);} \\ Amiram Eldar, Apr 24 2025

A070193 Numbers k such that gcd(3k,8^k+1) = 3 but k does not divide the numerator of B(2k) (the Bernoulli numbers).

Original entry on oeis.org

253, 1081, 1771, 2485, 2783, 3289, 4301, 4807, 5405, 5819, 7337, 7567, 7843, 9361, 10373, 10879, 11891, 12397, 12425, 13409, 13861, 14053, 14927, 15433, 17395, 17963, 18145, 18377, 18469, 19481, 19987, 20539, 20999, 22517, 23023, 24541

Offset: 1

Benoit Cloitre and Dean Hickerson, Apr 26 2002

nonn

Equivalently, numbers is in A070191 but not in A069040.

Amiram Eldar, Table of n, a(n) for n = 1..10000

Cf. A027641, A069040, A070191, A070192.

testb[n_] := Select[First/@FactorInteger[n], Mod[2n, #-1]==0&]=={}; test8[n_] := GCD[3n, PowerMod[8, n, 3n]+1]==3; Select[Range[25000], test8[ # ]&&!testb[ # ]&]

isA070191(k) = gcd(3*k, Mod(8, 3*k)^k + 1) == 3;
isok(k) = if(!isA070191(k), 0, my(p = factor(k)[,1]); for(i = 1, #p, if(!((2*k) % (p[i]-1)), return(1))); 0); \\ Amiram Eldar, Apr 24 2025

cp's OEIS Frontend

A069040 Numbers k that divide the numerator of B(2k) (the Bernoulli numbers).

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Extensions

A070192 Numbers k such that k divides the numerator of B(2k) (the Bernoulli numbers), but gcd(3k, 8^k+1) > 3.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

PARI

A070193 Numbers k such that gcd(3k,8^k+1) = 3 but k does not divide the numerator of B(2k) (the Bernoulli numbers).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

PARI