OFFSET
1,1
LINKS
Charles R Greathouse IV, Table of n, a(n) for n = 1..10000
EXAMPLE
MAPLE
N:= 10^6:
P:= select(isprime, [$1..N]):
S:= select(t -> numtheory:-bigomega(t)=2, [$1..N]):
P10:= {seq(add(P[i], i=m..m+9), m=1..nops(P)-9)}:
S10:= {seq(add(S[i], i=m..m+9), m=1..nops(S)-9)}:
sort(convert(P10 intersect S10, list)); # Robert Israel, Mar 20 2017
MATHEMATICA
With[{nn = 12600}, Intersection[Total /@ Partition[Prime@ Range@ PrimePi@ nn, 10, 1], Total /@ Partition[Select[Range@ nn, PrimeOmega@ # == 2 &], 10, 1]]] (* Michael De Vlieger, Mar 20 2017 *)
PROG
(PARI) list(lim)=if(lim<6504, return([])); my(v=List(), u=v, P=primes(9), x=(lim+10*log(lim))\1, t); forprime(p=2, x\2, forprime(q=2, min(x\p, p), listput(u, p*q))); u=Set(u); while(u[#u]+1+(t=sum(i=0, 8, u[#u-i]))<=lim, for(n=x+1, lim-t, if(issemi(n), u=concat(u, n); next(2))); break); for(i=1, #u-9, u[i]+=sum(j=1, 9, u[i+j])); t=vecsum(P); forprime(p=P[#P]+1, , t+=p; if(t>lim, break); if(setsearch(u, t), listput(v, t)); t-=P[1]; P=concat(P[2..9], p)); Vec(v) \\ Charles R Greathouse IV, Mar 20 2017
CROSSREFS
KEYWORD
nonn
AUTHOR
Zak Seidov, Mar 20 2017
STATUS
approved