Variational reduced-density-matrix method using three-particle N-representability conditions with application to many-electron molecules