cp's OEIS Frontend

A117929 Number of partitions of n into 2 distinct primes.

%I A117929 #47 Feb 16 2024 11:27:19
%S A117929 0,0,0,0,1,0,1,1,1,1,0,1,1,1,1,2,0,2,1,2,1,2,0,3,1,2,0,2,0,3,1,2,1,3,
%T A117929 0,4,0,1,1,3,0,4,1,3,1,3,0,5,1,4,0,3,0,5,1,3,0,3,0,6,1,2,1,5,0,6,0,2,
%U A117929 1,5,0,6,1,4,1,5,0,7,0,4,1,4,0,8,1,4,0,4,0,9,1,4,0,4,0,7,0,3,1,6,0,8,1,5,1
%N A117929 Number of partitions of n into 2 distinct primes.
%C A117929 Number of distinct rectangles with prime length and width such that L + W = n, W < L. For example, a(16) = 2; the two rectangles are 3 X 13 and 5 X 11. - _Wesley Ivan Hurt_, Oct 29 2017
%H A117929 T. D. Noe, <a href="/A117929/b117929.txt">Table of n, a(n) for n = 1..10000</a>
%F A117929 G.f.: Sum_{j>0} Sum_{i=1..j-1} x^(p(i)+p(j)), where p(k) is the k-th prime.
%F A117929 G.f.: A(x)^2/2 - A(x^2)/2 where A(x) = Sum_{p in primes} x^p. - _Geoffrey Critzer_, Nov 21 2012
%F A117929 a(n) = [x^n*y^2] Product_{i>=1} (1+x^prime(i)*y). - _Alois P. Heinz_, Nov 22 2012
%F A117929 a(n) = Sum_{i=2..floor((n-1)/2)} A010051(i) * A010051(n-i). - _Wesley Ivan Hurt_, Oct 29 2017
%e A117929 a(24) = 3 because we have [19,5], [17,7] and [13,11].
%p A117929 g:=sum(sum(x^(ithprime(i)+ithprime(j)),i=1..j-1),j=1..35): gser:=series(g,x=0,130): seq(coeff(gser,x,n),n=1..125);
%p A117929 # alternative
%p A117929 A117929 := proc(n)
%p A117929     local a,i,p ;
%p A117929     a := 0 ;
%p A117929     p := 2 ;
%p A117929     for i from 1 do
%p A117929         if 2*p >= n then
%p A117929             return a;
%p A117929         end if;
%p A117929         if isprime(n-p) then
%p A117929             a := a+1 ;
%p A117929         end if;
%p A117929         p := nextprime(p) ;
%p A117929     end do:
%p A117929 end proc:
%p A117929 seq(A117929(n),n=1..80) ; # _R. J. Mathar_, Oct 01 2021
%t A117929 l = {}; For[n = 1, n <= 1000, n++, c = 0; For[k = 1, Prime[k] < n/2, k++, If[PrimeQ[n - Prime[k]], c = c + 1] ]; AppendTo[l, c] ] l (* _Jake Foster_, Oct 27 2008 *)
%t A117929 Table[Count[IntegerPartitions[n,{2}],_?(AllTrue[#,PrimeQ]&&#[[1]]!= #[[2]] &)],{n,120}] (* Requires Mathematica version 10 or later *) (* _Harvey P. Dale_, Jul 26 2020 *)
%o A117929 (PARI) a(n)=my(s);forprime(p=2,(n-1)\2,s+=isprime(n-p));s \\ _Charles R Greathouse IV_, Feb 26 2014
%o A117929 (Python)
%o A117929 from sympy import sieve
%o A117929 from collections import Counter
%o A117929 from itertools import combinations
%o A117929 def aupton(max):
%o A117929     sieve.extend(max)
%o A117929     a = Counter(c[0]+c[1] for c in combinations(sieve._list, 2))
%o A117929     return [a[n] for n in range(1, max+1)]
%o A117929 print(aupton(105)) # _Michael S. Branicky_, Feb 16 2024
%Y A117929 Cf. A010051, A045917, A061358, A073610, A166081 (positions of 0), A077914 (positions of 2), A080862 (positions of 6).
%Y A117929 Column k=2 of A219180. - _Alois P. Heinz_, Nov 13 2012
%K A117929 nonn,easy
%O A117929 1,16
%A A117929 _Emeric Deutsch_, Apr 03 2006