A259548 -id:A259548 - cp's OEIS frontend

A340576 Decimal expansion of Product_{primes p == 5 (mod 6)} 1/(1-1/p^2).

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

Offset: 1

Jean-François Alcover, Jan 12 2021

The four similar sequences for products of primes mod 6 are these:
  A175646 for Product_{primes p == 1 (mod 6)} 1/(1-1/p^2),
  A340576 for Product_{primes p == 5 (mod 6)} 1/(1-1/p^2),
  A340577 for Product_{primes p == 1 (mod 6)} 1/(1+1/p^2),
  A340578 for Product_{primes p == 5 (mod 6)} 1/(1+1/p^2).

			1.06054829316911072817412636430987203493077130204487163127994372...

R. J. Mathar, Table of Dirichlet L-series and prime zeta modulo functions for small moduli, arXiv:1008.2547 [math.NT], 2010-2015, Zeta_{6,5}(2) in section 3.2.

Similar constants: A175646, A175647, A248930, A248938, A301429, A333240, A334826, A335963, A340577, A340578.

Cf. A259548.

a := n -> 3^(2^(-n-2))*((1-3^(-2^(n+1)))/2)^(2^(-n-1)):
b := n -> Zeta(n)/Im(polylog(n, (-1)^(2/3))):
c := n -> a(n)*b(2^(n+1))^(1/2^(n+1)):
Digits := 107: evalf((3/4)*mul(c(n), n=0..9)); # Peter Luschny, Jan 14 2021

digits = 105;
precision = digits + 10;
prodeuler[p_, a_, b_, expr_] := Product[If[a <= p <= b, expr, 1], {p, Prime[Range[PrimePi[a], PrimePi[b]]]}];
Lv3[s_] := prodeuler[p, 1, 2^(precision/s), 1/(1 - KroneckerSymbol[-3, p]*p^-s)] // N[#, precision] &;
Lv4[s_] := 2*Im[PolyLog[s, Exp[2*I*Pi/3]]]/Sqrt[3];
Lv[s_] := If[s >= 10000, Lv3[s], Lv4[s]];
gv[s_] := (1 - 3^(-s))*Zeta[s]/Lv[s];
pB = (3/4)*Product[gv[2^n*2]^(2^-(n+1)), {n, 0, 11}] // N[#, precision]&;
RealDigits[pB, 10, digits][[1]] (* Most of this code is due to Artur Jasinski *)
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}];
Z[m_, n_, s_] := (w = 1; sumz = 0; difz = 1; While[Abs[difz] > 10^(-digits - 5), difz = P[m, n, s*w]/w; sumz = sumz + difz; w++]; Exp[sumz]);
$MaxExtraPrecision = 1000; digits = 121; RealDigits[Chop[N[Z[6,5,2], digits]], 10, digits-1][[1]] (* Vaclav Kotesovec, Jan 15 2021 *)

g = A143298 = (9 - PolyGamma(1, 2/3) + PolyGamma(1, 4/3))/(4 sqrt(3));
h = A301429;
Equals (3*sqrt(3)*h^2)/2.
Equals (3/4)*A333240.
A340577 = Pi^4/(243*g*h^2);
A340578 = (45*g*h^2)/(2*Pi^2).
Equals Pi^2/(9*A175646). - Artur Jasinski, Jan 11 2021
Equals Sum_{k>=1} 1/A259548(k)^2. - Amiram Eldar, Jan 24 2021

A343431 Part of n composed of prime factors of the form 6k-1.

Original entry on oeis.org

1, 1, 1, 1, 5, 1, 1, 1, 1, 5, 11, 1, 1, 1, 5, 1, 17, 1, 1, 5, 1, 11, 23, 1, 25, 1, 1, 1, 29, 5, 1, 1, 11, 17, 5, 1, 1, 1, 1, 5, 41, 1, 1, 11, 5, 23, 47, 1, 1, 25, 17, 1, 53, 1, 55, 1, 1, 29, 59, 5, 1, 1, 1, 1, 5, 11, 1, 17, 23, 5, 71, 1, 1, 1, 25, 1, 11, 1, 1, 5, 1, 41, 83, 1, 85, 1, 29, 11, 89, 5

Offset: 1

Peter Munn, Apr 15 2021

Completely multiplicative with a(p) = p if p is of the form 6k-1 and a(p) = 1 otherwise.

Largest term of A259548 that divides n.

Equivalent sequence for distinct prime factors: A170825.
Equivalent sequences for prime factors of other forms: A000265 (2k+1), A343430 (3k-1), A170818 (4k+1), A097706 (4k-1), A248909 (6k+1), A065330 (6k+/-1), A065331 (<= 3), A355582 (<= 5).
Range of terms: A259548.

f[p_, e_] := If[Mod[p, 6] == 5, p^e, 1]; a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 100] (* after Amiram Eldar at A248909 *)

a(n) = {my(f = factor(n)); for (i=1, #f~, if ((f[i, 1] + 1) % 6, f[i, 1] = 1); ); factorback(f); } \\ after Michel Marcus at A248909

from math import prod
from sympy import factorint
def A343431(n): return prod(p**e for p, e in factorint(n).items() if not (p+1)%6) # Chai Wah Wu, Dec 26 2022

a(n) = n / A065331(n) / A248909(n) = A065330(n) / A248909(n).

cp's OEIS Frontend

A340576 Decimal expansion of Product_{primes p == 5 (mod 6)} 1/(1-1/p^2).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

Formula