OFFSET
1,1
COMMENTS
First differs from A060292 at n = 4.
If neither of the two unordered triples contains a 1, then the corresponding term is also in A337080. - Matej Veselovac, Sep 14 2020
LINKS
Alois P. Heinz, Table of n, a(n) for n = 1..10000
L. F. Meyers and R. See, The Census-Taker Problem, Math. Mag. 63 (1990) 86-88
I. J. L. Garces and M. L. Loyola, Revisiting a Number-Theoretic Puzzle: The Census-Taker Problem, Intersection, 11 (2010), pp. 28-38; arXiv:1204.2071 [math.HO], 2012.
EXAMPLE
36 = 6*6*1 = 9*2*2. 6+6+1 = 9+2+2. So 36 is in the sequence.
This example also explains the name of the sequence:
"A census taker knocks on a door. A mother answers.
The census taker says, "I need to know the number of children you have, and their ages." The woman responds in puzzle-ese, "I have three daughters, the product of their ages is 36, and the sum of their ages is equal to the house number next door."
The census taker, who never wastes questions, computes for a while and then asks, "Does your oldest daughter love dogs?" The mother answers affirmatively. The census taker says, "Thank you. I now know the ages."
What are the ages of the children?" (Garces and Loyola 2010, 28)
3264 is not in the sequence: (34, 32, 3) and (48, 17, 4) give sum 69, but (48, 34, 2), (64, 17, 3) and (68, 12, 4) give sum 84.
MAPLE
b:= proc(n, k, t) option remember; expand(`if`(t=1, `if`(k<n, 0, x^n),
add(`if`(d>k, 0, b(n/d, d, t-1)*x^d), d=numtheory[divisors](n))))
end:
q:= n-> (p-> degree(p)=2 and coeff(p, x, 2)=1)(
add(x^i, i=[coeffs(b(n$2, 3))])):
select(q, [$1..2000])[]; # Alois P. Heinz, May 16 2020
MATHEMATICA
b[n_, k_, t_] := b[n, k, t] = If[t==1, If[k<n, 0, x^n], Sum[If[d>k, 0, b[n/d, d, t-1] x^d], {d, Divisors[n]}]];
q[n_] := Exponent[#, x]==2 && Coefficient[#, x, 2]==1& @ Total[ x^CoefficientList[b[n, n, 3], x]];
Select[Range[2000], q] (* Jean-François Alcover, Nov 16 2020, after Alois P. Heinz *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Thomas Gawlick, May 16 2020
STATUS
approved