login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A283929 Number of ways of writing n as a sum of a twin prime (A001097) and a squarefree semiprime (A006881). 1

%I #12 Mar 29 2017 19:44:32

%S 0,0,0,0,0,0,0,0,0,1,0,1,0,2,0,1,0,3,1,2,1,2,1,2,1,3,2,3,2,3,0,2,2,3,

%T 2,2,1,3,3,4,3,4,2,3,3,4,4,2,1,3,3,5,4,4,2,3,3,4,4,1,2,1,5,4,5,6,2,4,

%U 5,5,4,2,3,2,5,5,6,5,2,4,5,5,6,2,3,4,4,6,5,4,3,3,5,6,8,3,7,4,9,6,6,3,3,3,5,6,7,4,5,3,5

%N Number of ways of writing n as a sum of a twin prime (A001097) and a squarefree semiprime (A006881).

%C Conjecture: a(n) > 0 for all n > 30.

%H Robert Israel, <a href="/A283929/b283929.txt">Table of n, a(n) for n = 0..10000</a>

%H Ilya Gutkovskiy, <a href="/A283929/a283929.pdf">Extended graphical example</a>

%H Eric Weisstein's World of Mathematics, <a href="http://mathworld.wolfram.com/Semiprime.html">Semiprime</a>

%H Eric Weisstein's World of Mathematics, <a href="http://mathworld.wolfram.com/Squarefree.html">Squarefree</a>

%H Eric Weisstein's World of Mathematics, <a href="http://mathworld.wolfram.com/TwinPrimes.html">Twin Primes</a>

%F G.f.: (Sum_{k>=1} x^A001097(k))*(Sum_{k>=1} x^A006881(k)).

%e a(17) = 3 because we have [14, 3], [11, 6] and [10, 7].

%p N:= 200: # to get a(0) to a(N)

%p V:= Vector(N):

%p Primes:= select(isprime,[2,seq(i,i=3..N+2)]):

%p PS:= convert(Primes,set);

%p Twins:= PS intersect map(`-`,PS,2):

%p Twins:= Twins union map(`+`,Twins,2):

%p Twins:= sort(convert(Twins,list)):

%p for i from 1 to nops(Twins) do

%p for j from 1 to nops(Primes) while Twins[i]+2*Primes[j] <= N do

%p for k from 1 to j-1 do

%p v:= Twins[i]+Primes[k]*Primes[j];

%p if v > N then break fi;

%p V[v]:= V[v]+1;

%p od od od:

%p 0, seq(V[i],i=1..N); # _Robert Israel_, Mar 29 2017

%t nmax = 110; CoefficientList[Series[Sum[Boole[PrimeQ[k] && (PrimeQ[k - 2] || PrimeQ[k + 2])] x^k, {k, 1, nmax}] Sum[MoebiusMu[k]^2 Floor[2/PrimeOmega[k]] Floor[PrimeOmega[k]/2] x^k, {k, 2, nmax}], {x, 0, nmax}], x]

%o (PARI) concat([0, 0, 0, 0, 0, 0, 0, 0, 0], Vec(sum(k=1, 110, (isprime(k) && (isprime(k - 2) || isprime(k + 2)))* x^k) * sum(k=2, 110, moebius(k)^2 * floor(2/bigomega(k)) * floor(bigomega(k)/2) * x^k) + O(x^111))) \\ _Indranil Ghosh_, Mar 18 2017

%Y Cf. A001097, A006881, A098983, A100949, A282192, A282318.

%K nonn

%O 0,14

%A _Ilya Gutkovskiy_, Mar 18 2017

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified March 28 20:05 EDT 2024. Contains 371254 sequences. (Running on oeis4.)