OFFSET
1,3
COMMENTS
a(n) = (A117991(n) - 1)/2.
LINKS
Robert Israel, Table of n, a(n) for n = 1..6000
EXAMPLE
n=4, the first 4 composite numbers are 4,6, 8 and 9:
a(3)=#{4/6,4/8,4/9,6/8,6/9,8/9}=#{4/9,1/2,2/3,3/4,8/9}=5.
MAPLE
N:= 1000: # to use composites up to N
comps:= remove(isprime, [$4..100]):
nc:= nops(comps):
S:= {}: A[1]:= 0:
for n from 2 to nc do
S:= S union {seq(comps[j]/comps[n], j=1..n-1)};
A[n]:= nops(S);
od:
seq(A[i], i=1..nc); # Robert Israel, Jan 30 2018
MATHEMATICA
M = 100; (* to use composites up to M *)
comps = Select[Range[4, M], CompositeQ];
nc = Length[comps]; S = {}; A[1] = 0;
For[n = 2, n <= nc, n++, S = S ~Union~ Table[comps[[j]]/comps[[n]], {j, 1, n - 1}]; A[n] = Length[S]];
Array[A, nc] (* Jean-François Alcover, Mar 10 2019, after Robert Israel *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Reinhard Zumkeller, Apr 08 2006
STATUS
approved