We examine the polynomial form of the scattering equations by means of computational algebraic geometry. The scattering equations are the backbone of the Cachazo-He-Yuan (CHY) representation of the S-matrix. We explain how the Bezoutian matrix facilitates the calculation of amplitudes in the CHY formalism, without explicitly solving the scattering equations or summing over the individual residues. Since for n-particle scattering, the size of the Bezoutian matrix grows only as (n − 3) × (n − 3), our algorithm is very efficient for analytic and numeric amplitude computations.

We develop a systematic procedure for computing maximal unitarity cuts of multiloop Feynman integrals in arbitrary dimension. Our approach is based on the Baikov representation in which the structure of the cuts is particularly simple. We examine several planar and nonplanar integral topologies and demonstrate that the maximal cut inherits IBPs and dimension shift identities satisfied by the uncut integral. Furthermore, for the examples we calculated, we find that the maximal cut functions from different allowed regions, form the Wronskian matrix of the differential equations on the maximal cut.