This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.
%I A115644 #2 Mar 30 2012 17:40:35 %S A115644 6,9,25,121,841,871,5041,5767,363721,368761,409111,3633841,3992431, %T A115644 3992551,4032121,4037791,39962281,39962311,39963031,40279711,40279801, %U A115644 43585921,43591687,43909207,479047801,479365321,479370271,482631271 %N A115644 Brilliant numbers (A078972) that are sums of distinct factorials. %e A115644 39962281 = 11! + 8! + 7! + 5!+ 1! = 4861*8221. %t A115644 brillQ[n_] := Block[{d = FactorInteger[n]}, Plus@@Last/@d==2 && (Last/@d=={2} || Length@IntegerDigits@((First/@d)[[1]])==Length@IntegerDigits@((First/@d)[[2]]))]; fac=Range[20]!;lst={}; Do[ n = Plus@@(fac*IntegerDigits[k, 2, 20]); If[brillQ[n], AppendTo[lst, n]], {k, 2^20-1}]; lst %Y A115644 Cf. A078972, A025494, A115645, A115646, A089359. %K A115644 nonn %O A115644 1,1 %A A115644 _Giovanni Resta_, Jan 27 2006