Julia: 一般項がn、n²、n³、...、n^10、n^50の数列の和

SymPy で n^10 までと，n^50 をやってみました

julia> summation(k^50, (k, 1, n)) |> factor |> print

n*(n + 1)*(2*n + 1)*(429*n^48 + 10296*n^47 + 75504*n^46 - 118404*n^45 - 3433716*n^44 + 5209776*n^43 + 177855964*n^42 - 269388834*n^41 - 8790556346*n^40 + 13320528936*n^39 + 399690120544*n^38 - 606195445284*n^37 - 16528902129916*n^36 + 25096450917516*n^35 + 617619215226134*n^34 - 938977048297959*n^33 - 20729072134582867*n^32 + 31563096726023280*n^31 + 620996836434449376*n^30 - 947276803014685704*n^29 - 16489621590673441320*n^28 + 25208070787517504832*n^27 + 385039595475766105512*n^26 - 590163428607407910684*n^25 - 7834990698171041513900*n^24 + 12047567761560266226192*n^23 + 137484855642671592627616*n^22 - 212251067344787522054520*n^21 - 2054997867949212429383888*n^20 + 3188622335596212405103092*n^19 + 25783333201061356027601406*n^18 - 40269310969390140243953655*n^17 - 266742505880274126845610963*n^16 + 420248414305106260390393272*n^15 + 2225375820438985887388497936*n^14 - 3548187937811031961277943540*n^13 - 14547150581890312152150257764*n^12 + 23594819841740984208864358416*n^11 + 71655262530103266642836795708*n^10 - 119280303716025392068687372770*n^9 - 251166419012752676999004054970*n^8 + 436389780377141711532849768840*n^7 + 569535901810418645522117212000*n^6 - 1072498742904198824049600702420*n^5 - 679511742411895276079901986300*n^4 + 1555516985069942326144653330660*n^3 + 137940192459690838148194087930*n^2 - 984668781224507420294617797225*n + 328222927074835806764872599075)/43758