f(p,q) = (-pq^2 -p^2 + pq + p + 1)/2