r=(3a4+16a2b64ac+256d)/256.