a(n) = (48*n^3 + 64*n^4 + 44*n^6 + 84*n^7 + 56*n^8 + 12*n^9 + 5*n^12 + 36*n^13 + 18*n^14 + 12*n^15 + 4*n^18 + n^24) / 384.
a(n) = C(n,1) + 49125*C(n, 2) + 740212980*C(n, 3) + 730815102166*C(n, 4) + 151600044933990*C(n, 5) + 11420034970306170*C(n, 6) + 415777158607920585*C(n, 7) + 8643499341510394200*C(n, 8) + 113988734942055623055*C(n, 9) + 1023002477284840979850*C(n, 10) + 6559265715033958749900*C(n, 11) + 31097943476763200314200*C(n, 12) + 111710751446923209781200*C(n, 13) + 309231173588248964052000*C(n, 14) + 666846649590586048584000*C(n, 15) + 1126625898539640346848000*C(n, 16) + 1492173541849975272288000*C(n, 17) + 1541987122059614438208000*C(n, 18) + 1229356526029003532160000*C(n, 19) + 741102367008078915840000*C(n, 20) + 326583680209195368960000*C(n, 21) + 99234043419574103040000*C(n, 22) + 18581137031073576960000*C(n, 23) + 1615751046180311040000*C(n, 24), where the coefficient of C(n,k) is the number of colorings using exactly k colors.