多项式求值

最新推荐文章于 2022-05-10 17:45:53 发布

原创最新推荐文章于 2022-05-10 17:45:53 发布 · 1k 阅读

CC 4.0 BY-SA版权

文章标签：

该博客展示了一个Polynomial类模板的代码实现，它继承自vector容器。类中包含构造函数用于初始化，还实现了秦九韶算法的求值函数，可根据不同情况计算多项式的值，代码使用了迭代器和模板技术。

// 版权所有(C) 梁意, 2004 // 最后修改: 2004.7.杭州

#ifndef Polynomial_HEADER

#define Polynomial_HEADER

#include<vector>

using namespace std;

template<typename T,typename Container=vector<T> >

class Polynomial: public Container

public:

   Polynomial() { }

   Polynomial(const T *pT,size_t N)

      copy(pT,pT+N,insert_iterator<Polynomial<T,Container> >(*this,begin()));

   Polynomial(const Polynomial<T,Container> &other)
   {

      copy(other.begin(),other.end(),insert_iterator<Polynomial<T,Container> >(*this,begin()));

   template<typename _In>

   inline _In QinJiuShao_evaluate(_In x,_In identity)//identity 为单位元,如: 整数1、浮点数1.0、单位矩阵

      size_t N=size();

      if(N==0) return identity;

      iterator iter=begin();

      if(N==1) return *iter;

      erator iter2=end();

      _In Res=(*iter)*x;

      iter++; iter2--;

      for(;iter!=iter2;iter++)

         Res=(Res+(*iter)*identity)*x;

      return Res+(*iter2)*identity;

   template<typename _In>

   inline  _In QinJiuShao_evaluate(_In x) //这个更快

   size_t N=size();

   if(N==0) return 1.0;

   iterator iter=begin();

   if(N==1) return *iter;

   iterator iter2=end();

   _In Res=(*iter)*x;

   iter2--;

   iter++;

   for(;iter!=iter2;iter++)

      Res=(Res+(*iter))*x;

   return Res+(*iter2);

};

#endif