sum_p=1^a (n^2-(2p-1)n)