matrix_of.hpp

Go to the documentation of this file.

#ifndef _MPOLY_MATRIX_OF_HPP_
#define _MPOLY_MATRIX_OF_HPP_

template<class Polynomial>
struct matrix_of {
  static Polynomial choose(unsigned n, unsigned d);
};

template<class Polynomial>
Polynomial matrix_of<Polynomial>::choose(unsigned n, unsigned d) {
  typedef typename Polynomial::monom_t Monomial;

  Polynomial L(1);


   for(int i=0;i<n;i++)
    L += Monomial(i,1);
    Polynomial res(L);
   for (unsigned i=1;i<d; ++i)
    res*= L;
   for(typename Polynomial::iterator it=res.begin();it!=res.end();it++)
    it->SetCoeff(1);


  return res;
}

#endif //_MPOLY_MATRIX_OF_HPP_