#include <synaps/arithm/gmp.h>
#include <synaps/mpol.h>       // multivariate polynomials

typedef MPol<QQ>    mpol_t;

mpol_t reduce(const mpol_t & p , const mpol_t & q)
{
  mpol_t r(p);
  ...
  return r;
}

int main(int argc, char** argv) {
  
  mpol_t p("2*x0^3*x1-34*x0^2*x2+23;"),q("3*x0^2-1;"); 
  cout<<p<<endl<<q<<endl;
  cout <<reduce(p,q)<<endl;

}