OFFSET
1,2
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
EXAMPLE
a(1) counts this sum: 6+9.
a(2) counts these sums: 6+15, 9+15, 10+15, 12+15.
a(3) counts these: 4+21, 6+21, 9+21, 12+21, 14+21, 15+21, 18+21.
MAPLE
Comps:= remove(isprime, [$4..10^5]):
OComps:= select(type, Comps, odd):
f:= proc(n) local d, m;
d:= OComps[n];
m:= ListTools:-BinarySearch(Comps, d);
nops(remove(c -> isprime(c+d), Comps[1..m-1]))
end proc:
map(f, [$1..100]); # Robert Israel, Jun 07 2026
MATHEMATICA
z = 400; p = Prime[Range[z]];
c = Select[Range[2, z], ! PrimeQ@# &]; (* A002808 *)
d = Select[Range[2, z], ! PrimeQ@# && OddQ@# &]; (* A014076 *)
f[n_] := Select[c, # < d[[n]] &];
g[n_] := d[[n]] + Select[c, # < d[[n]] &];
q[n_] := Length[Intersection[p, g[n]]];
tq = Table[q[n], {n, 1, 120}] (* A336406 *)
tc = Table[Length[f[n]], {n, 1, 120}] (* A336407 *)
m = Min[Length[tq], Length[tc]]; Take[tc, m] - Take[tq, m] (* A336408 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Clark Kimberling, Jul 20 2020
EXTENSIONS
Definition corrected by Robert Israel, Jun 07 2026
STATUS
approved
