A217795 -id:A217795 - cp's OEIS frontend

A337608 Decimal expansion of Lal's constant: the Hardy-Littlewood constant for A217795.

7, 9, 2, 2, 0, 8, 2, 3, 8, 1, 6, 7, 5, 4, 1, 6, 6, 8, 7, 7, 5, 4, 5, 5, 5, 6, 6, 5, 7, 9, 0, 2, 4, 1, 0, 1, 1, 2, 8, 9, 3, 2, 2, 5, 0, 9, 8, 6, 2, 2, 1, 1, 1, 7, 2, 2, 7, 9, 7, 3, 4, 5, 2, 5, 6, 9, 5, 1, 4, 1, 5, 4, 9, 4, 4, 1, 2, 4, 9, 0, 6, 6, 0, 2, 9, 5, 3, 8, 8, 3, 9, 8, 0, 2, 7, 5, 2, 9, 2, 7, 8, 7, 3, 9, 7, 3

Offset: 0

Amiram Eldar, Sep 04 2020

Shanks (1967) conjectured that the number of primes of the form (m + 1)^4 + 1 such that (m - 1)^4 + 1 is also a prime (A217795 plus 1), with m <= x, is asymptotic to c * li_2(x), where li_2(x) = Integral_{t=2..n} (1/log(t)^2) dt, and c is this constant. He defined c as in the formula section, evaluated it by 0.79220 and named it after the mathematician Mohan Lal, who conjectured the asymptotic formula without evaluating this constant.

The first 100 digits of this constant were calculated by Ettahri et al. (2019).

			0.792208238167541668775455566579024101128932250986221...

Steven R. Finch, Mathematical Constants, Cambridge University Press, 2003, pp. 90-91.

Keith Conrad, Hardy-Littlewood constants in: Mathematical properties of sequences and other combinatorial structures, Jong-Seon No et al. (eds.), Kluwer, Boston/Dordrecht/London, 2003, pp. 133-154, alternative link.
Salma Ettahri, Olivier Ramaré, Léon Surel, Fast multi-precision computation of some Euler products, arXiv:1908.06808 [math.NT], 2019 (Corollary 1.9).
Mohan Lal, Primes of the form n^4 + 1, Mathematics of Computation, Vol. 21, No. 98 (1967), pp. 245-247.
Daniel Shanks, Lal's constant and generalizations, Mathematics of Computation, Vol. 21, No. 100 (1967), pp. 705-707.
Eric Weisstein's World of Mathematics, Lal's Constant.

Cf. A000068, A037896, A088367, A210630, A217795, A217796.

Similar constants: A005597, A331941, A337606, A337607.

$MaxExtraPrecision = 1000; digits = 121;
f[p_] := (p-8)*(p+1)^4/((p-1)^4*p);
coefs = Rest[CoefficientList[Series[Log[f[1/x]], {x, 0, 1000}], x]];
S[m_, n_, s_] := (t = 1; sums = 0; difs = 1; While[Abs[difs] > 10^(-digits - 5) || difs == 0, difs = (MoebiusMu[t]/t) * Log[If[s*t == 1, DirichletL[m, n, s*t], Sum[Zeta[s*t, j/m]*DirichletCharacter[m, n, j]^t, {j, 1, m}]/m^(s*t)]]; sums = sums + difs; t++]; sums);
P[m_, n_, s_] := 1/EulerPhi[m] * Sum[Conjugate[DirichletCharacter[m, r, n]] * S[m, r, s], {r, 1, EulerPhi[m]}] + Sum[If[GCD[p, m] > 1 && Mod[p, m] == n, 1/p^s, 0], {p, 1, m}];
m = 2; sump = 0; difp = 1; While[Abs[difp] > 10^(-digits - 5) || difp == 0, difp = coefs[[m]]*(P[8, 1, m] - 1/17^m); sump = sump + difp; m++];
RealDigits[Chop[N[f[17] * Pi^4/(2^7 * Log[1+Sqrt[2]]^2) * Exp[sump], digits]], 10, digits - 1][[1]] (* Vaclav Kotesovec, Jan 16 2021 *)

Equals (Pi^4/(2^7 * log(1+sqrt(2))^2)) * Product_{primes p == 1 (mod 8)} (1 - 4/p)^2 * ((p + 1)/(p - 1))^4 * p*(p-8)/(p-4)^2 = (Pi^2/32) * A088367^2 * A334826^2 * A210630 = 2 * A337607^2 * A210630.

More terms from Vaclav Kotesovec, Jan 16 2021

A217796 Primes of the form n^4+1 such that (n+2)^4+1 is also prime.

Original entry on oeis.org

17, 257, 4477457, 8503057, 40960001, 59969537, 384160001, 5802782977, 58594980097, 94197431057, 102627966737, 114733948177, 283982410001, 330123790097, 381671897617, 405519334417, 691798081537, 741637881857, 1700843738897, 1749006250001, 2073600000001

Offset: 1

Michel Lagneau, Oct 12 2012

nonn

The corresponding n are in A217795.

			257 is in the sequence because  4^4+1 = 257 and (4+2)^4+1 = 1297 are both prime.

Vincenzo Librandi, Table of n, a(n) for n = 1..1000

Cf. A000068, A002523, A037896, A217795.

for n from 0 by 2 to 3500 do: if type(n^4+1,prime)=true and type((n+2)^4+1,prime)=true then printf(`%d, `, n^4+1):else fi:od:

lst={}; Do[p=n^4+1; q=(n+2)^4+1;If[PrimeQ[p] && PrimeQ[q], AppendTo[lst, p]], {n, 0, 3500}];lst
Select[Partition[Table[n^4+1,{n,1300}],3,1],AllTrue[{#[[1]],#[[3]]}, PrimeQ]&][[All,1]] (* Requires Mathematica version 10 or later *) (* Harvey P. Dale, Aug 17 2020 *)

cp's OEIS Frontend

A337608 Decimal expansion of Lal's constant: the Hardy-Littlewood constant for A217795.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Mathematica

Formula

Extensions