polspecial.hpp

Go to the documentation of this file.

template<typename mon, typename T,typename Base>
mon choice_naif(const pol<mon,T> & p, const Base &b)
{
  list<mon> tmp;
  mon m; 
  int maxdeg=-1;
  //int maxtmpdeg=0;
  for(int i=0;i<p.size;i++)
    {
      if (!Iszero(p.nf[i]))
    {
      mon tmpmon;
      int2mon(i,tmpmon);
      if(tmpmon.GetDegree()>maxdeg)
        {
          tmp.erase(tmp.begin(),tmp.end());
          maxdeg=tmpmon.GetDegree();
        }
      if(tmpmon.GetDegree()==maxdeg)
        {
          tmpmon.SetCoeff(p.nf[i]);
          tmp.push_back(tmpmon);
        }
    }
    }
  m=tmp.front();
  for(typename list<mon>::iterator itertmp=tmp.begin();itertmp!=tmp.end();itertmp++)
    {
      if(m.GetCoeff()<itertmp->GetCoeff()) {m=*itertmp;};
    }
  //  cout<<"monom choisi "<<m<<endl;
  return m;
}

template<typename mon, typename T,typename Base>
mon choice_grevlex(const pol<mon,T> & p,const Base &b)
{
  list<mon> tmp;
  mon m; 
  int maxdeg=-1;
  //int maxtmpdeg=0;
  for(int i=0;i<p.size;i++)
    {
      if (!Iszero(p.nf[i]))
    {
      mon tmpmon;
      int2mon(i,tmpmon);
      if(tmpmon.GetDegree()>maxdeg)
        {
          tmp.erase(tmp.begin(),tmp.end());
          maxdeg=tmpmon.GetDegree();
        }
      if(tmpmon.GetDegree()==maxdeg)
        {
          tmpmon.SetCoeff(p.nf[i]);
          tmp.push_back(tmpmon);
        }
    }
    }
  m=tmp.front();
  for(typename list<mon>::iterator itertmp=tmp.begin();itertmp!=tmp.end();itertmp++)
    {
      if(my_ord(*itertmp,m)<0) {m=*itertmp;};
    }
  // cout<<"monom choisi "<<m<<endl;
  return m;
}

template<typename mon, typename T,typename Base>
mon choice_dlex(const pol<mon,T> & p,const Base &b)
{
  list<mon> tmp;
  mon m; 
  int maxdeg=-1;
  //int maxtmpdeg=0;
  for(int i=0;i<p.size;i++)
    {
      if (!Iszero(p.nf[i]))
    {
      mon tmpmon;
      int2mon(i,tmpmon);
      if(tmpmon.GetDegree()>maxdeg)
        {
          tmp.erase(tmp.begin(),tmp.end());
          maxdeg=tmpmon.GetDegree();
        }
      if(tmpmon.GetDegree()==maxdeg)
        {
          tmpmon.SetCoeff(p.nf[i]);
          tmp.push_back(tmpmon);
        }
    }
    }
  m=tmp.front();
  for(typename list<mon>::iterator itertmp=tmp.begin();itertmp!=tmp.end();itertmp++)
    {
      if(my_ord(*itertmp,m)>0) {m=*itertmp;};
    }
  //  cout<<"monom choisi "<<m<<endl;
  return m;
}


template<typename mon, typename T,typename Base>
mon choice_mac(const pol<mon,T> & p,const Base &b)
{
  list<mon> tmp;
  mon m; 
  int maxdeg=-1;
  int maxtmpdeg=0;
  for(int i=0;i<p.size;i++)
    {
      if (!Iszero(p.nf[i]))
    {
      mon tmpmon;
      int2mon(i,tmpmon);
      if(tmpmon.GetDegree()>maxdeg)
        {
          tmp.erase(tmp.begin(),tmp.end());
          maxdeg=tmpmon.GetDegree();
        }
      if(tmpmon.GetDegree()==maxdeg)
        {
          tmpmon.SetCoeff(p.nf[i]);
          tmp.push_back(tmpmon);
        }
    }
    }
  m=tmp.front();
  for(typename list<mon>::iterator itertmp=tmp.begin();itertmp!=tmp.end();itertmp++)
    {
      //cout<<"monomes regarde "<<*itertmp<<" Lavr  "<<lvar(itertmp->rep)<<endl;
      for(int i=0;i<=lvar(itertmp->rep);i++)
    {
      if(itertmp->GetDegree(i)>=maxtmpdeg)
        {
          maxtmpdeg=itertmp->GetDegree(i);
          m=*itertmp;
        }
    }
    }
  return m;
}


template<typename mon, typename T,typename Base>
mon choice_mpq(const pol<mon,T> & p,const Base &b)
{
  int i,degp=Degree(p);
  int mini=INT_MAX;
  mon tmpmon;
  mon stockmon;
  for(i=0;i<p.size;i++)
    if(!Iszero(p.nf[i]))
      {
    int2mon(i,tmpmon);
    if(tmpmon.GetDegree()==degp)
      {
        if(mpz_size(mpq_numref(&(p.nf[i].rep())))+
           mpz_size(mpq_denref(&(p.nf[i].rep())))<mini) 
          {
        mini=mpz_size(mpq_numref(&(p.nf[i].rep())))+
          mpz_size(mpq_denref(&(p.nf[i].rep())));
        //cout<<"newmax "<<max<<endl;
        //cout<<" les monomes "<<*iter<<endl;
        stockmon=tmpmon*p.nf[i];
          }
      }
      }
  //cout<<"je retourne "<<stockmon<<endl;
  return stockmon;
}

template<typename mon, typename T, typename Base>
mon choice_fall(const pol<mon,T> & p, const Base &b)
{
  int i;
  int maxi=-1;
  mon tmpmon;
  mon stockmon;
  mon res=choice_mac(p,b);
  typename Base::const_iterator iterb;
  list<mon> tmplist;
  int degp=Degree(p);
  
  for(i=0;i<p.size;i++)
    if(!Iszero(p.nf[i]))
      {
    int2mon(i,tmpmon);
    if(tmpmon.GetDegree()==degp)
      tmplist.push_back(tmpmon*p.nf[i]);
    //cout<<" k "<<i<<endl;
      }
  for(typename list<mon>::iterator iterlist=tmplist.begin()
    ;iterlist!=tmplist.end();iterlist++)
    {
      for(iterb=b.begin();iterb!=b.end();iterb++)
    {
      int isgood=0;
      for(int i=0;i<iterb->taille1;i++)
        if(isdivisible(iterb->accept[i],*iterlist,b))
          {isgood=1;break;}
      for(int i=0;i<iterb->taille2;i++)
        if(isdivisible(iterb->refuse[i],*iterlist,b))
          {isgood=0;break;}
      if(isgood)
        {

          for(int i=0;i<iterb->taille2;i++)
        {
          mon m=lcm(*iterlist,iterb->refuse[i]);
          int dist=0;
          for(int j=0;j<iterb->taille2;j++)
            if ((j!=i) && (isdivisible(iterb->refuse[j],m,b)))
              isgood=0;
          if(isgood)
            {
              for(int j=0;j<b.nbvar();j++)
            if(m.GetDegree(j)>iterb->refuse[i].GetDegree(j))
              dist+=m.GetDegree(j)-iterb->refuse[i].GetDegree(j);
              if(dist>maxi)
            {
              res=*iterlist;
              maxi=dist;
            }
            }
        }
        }
      
    }
    }
  //cout<<"le pol "<<invconv<Poly>(p)<<endl<<"le mon choisi "<<res<<endl;
  return res;
}