(PARI) bc=[[5, 0, 19, 24], [5, -16, 75, 24], [5, -18, 64, 24], [5, -18, 88, 24], [5, -34, 24, 24], [5, -34, 192, 24], [5, -34, 216, 24], [5, -34, 0, 24]];
m=[[1, 6, 1], [5, 6, 1], [3, 6, 2], [2, 12, 3], [10, 12, 3], [4, 12, 4], [8, 12, 4], [0, 60, 5], [6, 60, 6], [18, 60, 6], [42, 60, 6], [54, 60, 6], [12, 60, 7], [24, 60, 7], [36, 60, 7], [48, 60, 7], [30, 60, 8]];
ix(n)=for(i=1, length(m), x=m[i]; if(n%x[2]==x[1], return(x[3]))); -1
a(n)=x=bc[ix(n)]; (x[1]*n^2+x[2]*n+x[3])/x[4]
vector(200, x, a(x+2))