expand {c(4k)(c^2-16k^2), c(3k^2-2k-1)(c^2-(3k^2-2k-1)^2), (3k^2+2k-1)c((3k^2+2k-1)^2-c^2)}, where c=3k^2+1