A242195 -id:A242195 - cp's OEIS frontend

A242193 Least prime p such that B_{2*n} == 0 (mod p) but there is no k < n with B_{2k} == 0 (mod p), or 1 if such a prime p does not exist, where B_m denotes the m-th Bernoulli number.

1, 1, 1, 1, 5, 691, 7, 3617, 43867, 283, 11, 103, 13, 9349, 1721, 37, 17, 26315271553053477373, 19, 137616929, 1520097643918070802691, 59, 23, 653, 417202699, 577, 39409, 113161, 29, 2003, 31, 1226592271, 839, 101, 688531

Offset: 1

Zhi-Wei Sun, May 07 2014

nonn

Conjecture: a(n) is prime for any n > 4.
It is known that (-1)^(n-1)*B_{2*n} > 0 for all n > 0.
See also A242194 for a similar conjecture involving Euler numbers.

			a(14) = 9349 since the numerator of |B_{28}| is 7*9349*362903 with B_2*B_4*B_6*...*B_{26} not congruent to 0 modulo 9349, but B_{14} == 0 (mod 7).

Peter Luschny, Table of n, a(n) for n = 1..103, (a(1)..a(60) from Zhi-Wei Sun)
Z.-W. Sun, New observations on primitive roots modulo primes, arXiv preprint arXiv:1405.0290 [math.NT], 2014.

Cf. A000040, A027641, A242169, A242170, A242171, A242173, A242194, A242195.

b[n_]:=Numerator[Abs[BernoulliB[2n]]]
f[n_]:=FactorInteger[b[n]]
p[n_]:=Table[Part[Part[f[n],k],1],{k,1,Length[f[n]]}]
Do[If[b[n]<2,Goto[cc]];Do[Do[If[Mod[b[i],Part[p[n],k]]==0,Goto[aa]],{i,1,n-1}];Print[n," ",Part[p[n],k]];Goto[bb];Label[aa];Continue,{k,1,Length[p[n]]}];Label[cc];Print[n," ",1];Label[bb];Continue,{n,1,35}]
(* Second program: *)
LPDtransform[n_, fun_] := Module[{d}, d[p_] := AllTrue[Range[n - 1], !Divisible[fun[#], p]&]; SelectFirst[FactorInteger[fun[n]][[All, 1]], d] /. Missing[_] -> 1];
A242193list[sup_] := Table[LPDtransform[n, Function[k, Abs[BernoulliB[2k]] // Numerator]], {n, 1, sup}]
A242193list[35] (* Jean-François Alcover, Jul 27 2019, after Peter Luschny *)

def LPDtransform(n, fun):
    d = lambda p: all(not p.divides(fun(k)) for k in (1..n-1))
    return next((p for p in prime_divisors(fun(n)) if d(p)), 1)
A242193list = lambda sup: [LPDtransform(n, lambda k: abs(bernoulli(2*k)).numerator()) for n in (1..sup)]
print(A242193list(35)) # Peter Luschny, Jul 26 2019

A242194 Least prime divisor of E_{2n} which does not divide any E_{2k} with k < n, or 1 if such a primitive prime divisor of E_{2*n} does not exist, where E_m denotes the m-th Euler number given by A122045.

Original entry on oeis.org

1, 5, 61, 277, 19, 13, 47, 17, 79, 41737, 31, 2137, 67, 29, 15669721, 930157, 4153, 37, 23489580527043108252017828576198947741, 41, 137, 587, 285528427091, 5516994249383296071214195242422482492286460673697, 5639, 53, 2749, 5303, 1459879476771247347961031445001033, 6821509

Offset: 1

Zhi-Wei Sun, May 07 2014

Conjecture: a(n) is prime for any n > 1.
It is known that (-1)^n*E_{2*n} > 0 for all n = 0, 1, ....
See also A242193 for a similar conjecture involving Bernoulli numbers.

			a(4) = 277 since E_8 = 5*277 with 277 not dividing E_2*E_4*E_6, but 5 divides E_4 = 5.

Peter Luschny, Table of n, a(n) for n = 1..82, (a(1)..a(34) from Zhi-Wei Sun, a(35)..a(38) from Jean-François Alcover).
Z.-W. Sun, New observations on primitive roots modulo primes, arXiv preprint arXiv:1405.0290 [math.NT], 2014.

Cf. A000040, A000364, A122045, A242169, A242170, A242171, A242173, A242193, A242195.

e[n_]:=Abs[EulerE[2n]]
f[n_]:=FactorInteger[e[n]]
p[n_]:=p[n]=Table[Part[Part[f[n],k],1],{k,1,Length[f[n]]}]
Do[If[e[n]<2,Goto[cc]];Do[Do[If[Mod[e[i],Part[p[n],k]]==0,Goto[aa]],{i,1,n-1}];Print[n," ",Part[p[n],k]];Goto[bb];Label[aa];Continue,{k,1,Length[p[n]]}];Label[cc];Print[n," ",1];Label[bb];Continue,{n,1,30}]
(* Second program: *)
LPDtransform[n_, fun_] := Module[{}, d[p_, m_] := d[p, m] = AllTrue[ Range[m-1], ! Divisible[fun[#], p]&]; f[m_] := f[m] = FactorInteger[ fun[m]][[All, 1]]; SelectFirst[f[n], d[#, n]&] /. Missing[_] -> 1];
a[n_] := a[n] = LPDtransform[n, Function[k, Abs[EulerE[2k]]]];
Table[Print[n, " ", a[n]]; a[n], {n, 1, 38}]  (* Jean-François Alcover, Jul 28 2019, non-optimized adaptation of Peter Luschny's Sage code *)

# uses[LPDtransform from A242193]
A242194list = lambda sup: [LPDtransform(n, lambda k: euler_number(2*k)) for n in (1..sup)]
print(A242194list(16)) # Peter Luschny, Jul 26 2019

A242174 Least prime divisor of A005260(n) which does not divide any previous term A005260(k) with k < n, or 1 if such a primitive prime divisor of A005260(n) does not exist.

Original entry on oeis.org

2, 3, 41, 5, 7, 349, 61, 75617, 31, 13, 499, 643897693, 17, 19, 1729774061, 101, 2859112064587, 138407, 83, 167, 59, 29, 653, 257, 997540809461453561581, 347, 13679, 37, 160449179727717672892660463, 211, 151, 43, 97, 73, 47

Offset: 1

Zhi-Wei Sun, May 07 2014

nonn

Conjecture: a(n) is prime for any n > 0. In general, for any r > 2, if n is large enough then f_r(n) = sum_{k=0..n}C(n,k)^r has a prime divisor which does not divide any previous terms f_r(k) with k < n.

			a(3) = 41 since A005260(3) = 2^2*41 with 41 dividing none of A005260(1) = 2 and A005260(2) = 2*3^2.

Zhi-Wei Sun, Table of n, a(n) for n = 1..82

Cf. A000040, A005260, A242169, A242170, A242171, A242173, A242193, A242194, A242195, A242207.

u[n_]:=Sum[Binomial[n,k]^4,{k,0,n}]
f[n_]:=FactorInteger[u[n]]
p[n_]:=Table[Part[Part[f[n],k],1],{k,1,Length[f[n]]}]
Do[If[u[n]<2,Goto[cc]];Do[Do[If[Mod[u[i],Part[p[n],k]]==0,Goto[aa]],{i,1,n-1}];Print[n," ",Part[p[n],k]];Goto[bb];Label[aa];Continue,{k,1,Length[p[n]]}];Label[cc];Print[n," ",1];Label[bb];Continue,{n,1,35}]

A242207 Least prime divisor of the n-th Domb number D(n) which does not divide any D(k) with k < n, or 1 if such a primitive prime divisor of D(n) does not exist.

Original entry on oeis.org

2, 7, 1, 97, 11, 23, 19, 643, 659, 1753, 4922329, 613, 341447, 1193, 2213, 2040452101603, 491, 82461839, 733, 113, 1108394340978316050481, 1034497328556150923437, 593, 73, 17117, 804943, 422291, 1559, 858631, 337655751557

Offset: 1

Zhi-Wei Sun, May 07 2014

nonn

Conjecture: a(n) is prime except for n = 3.

			 a(4) = 97 since D(4) = 2^2*7*97 with 97 dividing none of D(1) = 2^2, D(2) = 2^2*7 and D(3) = 2^8.

Zhi-Wei Sun, Table of n, a(n) for n = 1..72

Cf. A000040, A002895, A242169, A242170, A242171, A242173, A242174, A242193, A242194, A242195.

d[n_]:=Sum[Binomial[n,k]^2*Binomial[2k,k]*Binomial[2(n-k),n-k],{k,0,n}]
f[n_]:=FactorInteger[d[n]]
p[n_]:=Table[Part[Part[f[n],k],1],{k,1,Length[f[n]]}]
Do[If[d[n]<2,Goto[cc]];Do[Do[If[Mod[d[i],Part[p[n],k]]==0,Goto[aa]],{i,1,n-1}];Print[n," ",Part[p[n],k]];Goto[bb];Label[aa];Continue,{k,1,Length[p[n]]}];Label[cc];Print[n," ",1];Label[bb];Continue,{n,1,30}]

A242223 Least prime p such that H(n) == 0 (mod p) but H(k) == 0 (mod p) for no 0 < k < n, or 1 if such a prime p does not exist, where H(n) denotes the n-th harmonic number sum_{k=1..n}1/k.

Original entry on oeis.org

1, 3, 11, 5, 137, 7, 1, 761, 7129, 61, 97, 13, 29, 1049, 41233, 17, 37, 19, 7440427, 11167027, 18858053, 23, 583859, 577, 109, 34395742267, 521, 375035183, 4990290163, 31, 2667653736673, 2917, 269, 3583, 397, 1297, 10839223, 199, 737281, 41

Offset: 1

Zhi-Wei Sun, May 08 2014

nonn

Conjecture: a(n) is prime except for n = 1, 7.

			a(4) = 5 since H(4) = 25/12 == 0 (mod 5), but none of H(1) = 1, H(2) = 3/2 and H(3) = 11/6 is congruent to 0 modulo 5.

Zhi-Wei Sun, Table of n, a(n) for n = 1..184
Z.-W. Sun, New observations on primitive roots modulo primes, arXiv preprint arXiv:1405.0290 [math.NT], 2014.

Cf. A000040, A001008, A002805, A242169, A242170, A242171, A242173, A242174, A242193, A242194, A242195, A242207, A242222.

h[n_]:=Numerator[HarmonicNumber[n]]
f[n_]:=FactorInteger[h[n]]
p[n_]:=Table[Part[Part[f[n],k],1],{k,1,Length[f[n]]}]
Do[If[h[n]<2,Goto[cc]];Do[Do[If[Mod[h[i],Part[p[n],k]]==0,Goto[aa]],{i,1,n-1}];Print[n," ",Part[p[n],k]];Goto[bb];Label[aa];Continue,{k,1,Length[p[n]]}];Label[cc];Print[n," ",1];Label[bb];Continue,{n,1,40}]

A242292 Least prime divisor of 2^n - n which does not divide any 2^k - k with 0 < k < n, or 1 if such a primitive prime divisor of 2^n - n does not exist.

Original entry on oeis.org

1, 2, 5, 3, 1, 29, 11, 31, 503, 13, 7, 1021, 8179, 1637, 4679, 1, 8737, 131063, 524269, 262139, 2097131, 349, 131, 773, 271, 197, 457, 1493, 317, 17, 6733, 73, 41, 157109, 83, 53, 1741, 3329, 49977801259, 997, 149, 2199023255531, 61, 4398046511093, 3769453

Offset: 1

Zhi-Wei Sun, May 10 2014

nonn

Conjecture: a(n) = 1 only for n = 1, 5, 16.

In constrast, a classical theorem of Bang asserts that if n > 1 is different from 6 then 2^n - 1 has a prime divisor which does not divide any 2^k - 1 with 0 < k < n.

			a(4) = 3 since 2^4 - 4 = 2^2*3 with 3 dividing none of 2^1 - 1 = 1, 2^2 - 2 = 2 and 2^3 - 3 = 5.

A. S. Bang, Taltheoretiske Undersgelser, Tidsskrift fur Mat. 4(1886), no. 5, 70--80, 130--137.

Zhi-Wei Sun, Table of n, a(n) for n = 1..215
Z.-W. Sun, New observations on primitive roots modulo primes, arXiv preprint arXiv:1405.0290 [math.NT], 2014.

Cf. A000040, A000325, A112927, A242169, A242170, A242171, A242173, A242174, A242193, A242194, A242195, A242207.

u[n_]:=2^n-n
f[n_]:=FactorInteger[u[n]]
p[n_]:=Table[Part[Part[f[n], k], 1], {k, 1, Length[f[n]]}]
Do[If[u[n]<2, Goto[cc]]; Do[Do[If[Mod[u[i], Part[p[n], k]]==0, Goto[aa]], {i, 1, n-1}]; Print[n, " ", Part[p[n], k]]; Goto[bb]; Label[aa]; Continue, {k, 1, Length[p[n]]}]; Label[cc]; Print[n, " ", 1]; Label[bb]; Continue, {n, 1, 45}]

cp's OEIS Frontend

A242193 Least prime p such that B_{2*n} == 0 (mod p) but there is no k < n with B_{2k} == 0 (mod p), or 1 if such a prime p does not exist, where B_m denotes the m-th Bernoulli number.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

Sage

A242194 Least prime divisor of E_{2*n} which does not divide any E_{2*k} with k < n, or 1 if such a primitive prime divisor of E_{2*n} does not exist, where E_m denotes the m-th Euler number given by A122045.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

Sage

A242174 Least prime divisor of A005260(n) which does not divide any previous term A005260(k) with k < n, or 1 if such a primitive prime divisor of A005260(n) does not exist.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

A242207 Least prime divisor of the n-th Domb number D(n) which does not divide any D(k) with k < n, or 1 if such a primitive prime divisor of D(n) does not exist.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

A242223 Least prime p such that H(n) == 0 (mod p) but H(k) == 0 (mod p) for no 0 < k < n, or 1 if such a prime p does not exist, where H(n) denotes the n-th harmonic number sum_{k=1..n}1/k.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

A242292 Least prime divisor of 2^n - n which does not divide any 2^k - k with 0 < k < n, or 1 if such a primitive prime divisor of 2^n - n does not exist.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Mathematica

A242194 Least prime divisor of E_{2n} which does not divide any E_{2k} with k < n, or 1 if such a primitive prime divisor of E_{2*n} does not exist, where E_m denotes the m-th Euler number given by A122045.