A080088 Let r+i*s be the sum of the first-quadrant Gaussian primes dividing prime(n); sequence gives r values.

2, 3, 3, 7, 11, 5, 5, 19, 23, 7, 31, 7, 9, 43, 47, 9, 59, 11, 67, 71, 11, 79, 83, 13, 13, 11, 103, 107, 13, 15, 127, 131, 15, 139, 17, 151, 17, 163, 167, 15, 179, 19, 191, 19, 15, 199, 211, 223, 227, 17, 21, 239, 19, 251, 17, 263, 23, 271, 23, 21, 283, 19, 307, 311, 25, 25

Offset: 1

Author

Vladeta Jovovic, Jan 27 2003

Keywords

• nonn

Links

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

Crossrefs

Cf. A078908, A078909, A080089.

Programs

Mathematica

s[n_] := Module[{f = FactorInteger[n, GaussianIntegers->True]}, p = f[[;;,1]]; e = f[[;;,2]]; Re[Plus @@ ((If[Abs[#]==1, 0, #]& /@ p) * e)]]; s /@ Select[ Range[100], PrimeQ] (* Amiram Eldar, Feb 28 2020 *)