cp's OEIS Frontend

A354709 Decimal expansion of Sum_{p prime} 3*(2p-1)*log(p)/(p^3 + p^2 - 3p + 1).

%I A354709 #16 Jun 06 2022 02:49:14
%S A354709 2,5,2,9,0,6,6,1,7,3,5,8,0,9,2,9,9,2,9,2,5,9,5,8,7,1,2,9,3,0,1,8,9,4,
%T A354709 5,9,2,3,0,0,0,9,2,2,3,9,9,4,4,3,9,9,7,6,1,1,8,8,9,9,2,5,6,2,7,0,1,3,
%U A354709 5,7,8,0,0,6,6,2,8,6,4,7,7,4,9,6,1,5,1,7,2,2,4,6,7,7,6,3,3,2,0,4,4,3,2,6,5
%N A354709 Decimal expansion of Sum_{p prime} 3*(2p-1)*log(p)/(p^3 + p^2 - 3p + 1).
%C A354709 Also logarithmic derivative of A(s,w) at (0,0), where A(s,w) = Product_{p prime} (1 - (1 - (p*(1 - p^(-1-s))^3)/(-1+p))*(1 - (p*(1 - p^(-1-w))^3)/(-1+p))), with A(0,0) = A256392.
%e A354709 2.52906617358092992925958712930189459230009223994439976118899256270135780066...
%t A354709 Block[{$MaxExtraPrecision = 1000},
%t A354709 Do[CC = Join[{0},
%t A354709     Series[(3 (-1 + 2 p))/(1 - 3 p + p^2 + p^3) //. p -> 1/x, {x, 0,
%t A354709        t}][[3]]];
%t A354709   Print[N[-Sum[
%t A354709         CC[[k]]*(PrimeZetaP'[k] + Log[2]/2^k), {k, 1, Length[CC]}] + (
%t A354709       3 (-1 + 2 p) Log[p])/(1 - 3 p + p^2 + p^3) //. p -> 2, 75]], {t,
%t A354709     1000, 1500, 100}]]
%t A354709 ratfun = 3*(2*p - 1)/(p^3 + p^2 - 3*p + 1); zetas = 0; ratab = Table[konfun = Simplify[ratfun + c/(p^power - 1)] // Together; coefs = CoefficientList[Numerator[konfun], p]; sol = Solve[Last[coefs] == 0, c][[1]]; zetas = zetas + c*Zeta'[power]/Zeta[power] /. sol; ratfun = konfun /. sol, {power, 2, 20}]; Do[Print[N[Sum[Log[p]*ratfun /. p -> Prime[k], {k, 1, m}] + zetas, 120]], {m, 2000, 20000, 2000}] (* _Vaclav Kotesovec_, Jun 04 2022 *)
%Y A354709 Cf. A256392.
%K A354709 nonn,cons
%O A354709 1,1
%A A354709 _David Nguyen_, Jun 03 2022
%E A354709 More digits from _Vaclav Kotesovec_, Jun 04 2022