OFFSET
1,15
COMMENTS
LINKS
Clark Kimberling, Table of n, a(n) for n = 1..1000
EXAMPLE
a(312) counts these 7 products: 3*104, 4*78, 6*52, 8*39, 12*26, 24*13, 156*2
MAPLE
N:= 1000: # to get a(1) to a(N)
A201:= [seq(floor(n*phi), n=1..N)]:
A1950:= [seq(floor(n*phi^2), n=1..N)]:
A:= Vector(N):
for i from 1 to N do
for j from 1 do
m:= A201[i]*A1950[j];
if m > N then break fi;
A[m]:= A[m]+1;
od
od:
convert(A, list); # Robert Israel, Dec 23 2014
MATHEMATICA
r = (1 + Sqrt[5])/2; s = r/(r - 1); t = Flatten[Table[Floor[r*j]*Floor[s*k], {j, 1, 300}, {k, 1, 300}]]; a[n_] := Count[t, n]; u = Table[a[n], {n, 1, 300}]
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Clark Kimberling, Dec 23 2014
STATUS
approved