factorion[n, base] := sum[map["factorial", integerDigits[n, base]]]

for base = 9 to 12
{
   for n = 1 to 1_499_999
      if n == factorion[n, base]
         println["$base\t$n"]
}