[seqfan] Re: Sums of two primes with prime subscripts
Richard Mathar
mathar at strw.leidenuniv.nl
Fri Mar 26 12:19:28 CET 2010
http://list.seqfan.eu/pipermail/seqfan/2010-March/004185.html says
jvp> Half-sums (averages) of two primes with prime subscripts.
jvp>
jvp> {(A006450(i) + A006450(j))/2} = {(A000040(A000040(i)) + A000040(A000040(j)))/2}
jvp>
jvp> n sum(n) halfsum(n) Note
jvp> 1 6 3 3 + 3
jvp> 2 8 4 3 + 5
jvp> 3 10 5 5 + 5
jvp> 4 14 7 3 + 11
jvp> 5 16 8 5 + 11
jvp> 6 20 10 3 + 17
jvp> 7 22 11 11 + 11 = 17 + 5
jvp> 8 28 14 11 + 17
jvp> 9 34 17 17 + 17 = 3 + 31
jvp> 10 36 18 5 + 31
jvp> 11 42 21 11 + 31
jvp> 12 44 22 3 + 41
jvp> 13 46 23 5 + 41
Sorted, the halfsums(n) would continue 24, 26, 29,31 ,32, 35,
36, 38, 39, 41, 42, 43, 44, 45, 47, 49, 50, 54, 56 ... here, not 48 , 26, 29,..:
3, 4, 5, 7, 8, 10, 11, 14, 17, 18, 21, 22, 23, 24, 26, 29, 31, 32, 35, 36, 38,
39, 41, 42, 43, 44, 45, 47, 49, 50, 54, 56, 57, 59, 60, 62, 63, 65, 66, 67,
69, 70, 71, 72, 75, 79, 80, 81, 83, 84, 87, 88, 91, 92, 93, 94, 95, 96, 97,
98, 99, 101, 104, 105, 107, 108, 109, 110, 111, 112
jvp> 14 14 48 17 + 31 = 7 + 41
jvp> 15 52 26 11 + 41
jvp> 16 58 29 17 + 41
<snip>
In Maple:
hfs := {} ;
for i from 1 to 100 do
for j from 1 to i do
ithprime(ithprime(i))+ithprime(ithprime(j)) ;
hfs := hfs union {%/2}
end do:
end do:
sort(hfs) ;
More information about the SeqFan
mailing list