b[i] = b[i-1] + d ； b[1] = 1
b[i] = d*(i-1) + 1
a[i] = a[i-1] + b[i] ； a[0] = 1
a[i] = 1 + for(k=1,k<=i,k++)b[k]
a[i] = 1 + i + d*i(i-1)/2
S(a[0-i]) = (i+1) + i(i+1)/2 + d/2*[i(i+1)(2i+1)/6 + i(i+1)/2]

這邊看起來有點怪，因為我是從第0巷開始計算，或許可以改好看一些

在本題，i = 49

經過一連串化簡，1275 + d*19600