reduce for p: -p D[b-ap,p]/(b-ap)>1, 0 < p < b/a