秒大刀 博客

//1_20
/*
试编写算法求一元多项式Pn(x)=Sum(ai.x^i),(i=0,...n)的值Pn(x0),并确定算法中
每一语句的执行次数和整个算法的时间复杂度。注意选择你认为较好的输入的输出
方法。本题的输入为ai(i=0,1,...,n),x0和n,输出为Pn(x0)。
*/
//完成日期：2005/3/3
#include <iostream>
using namespace std;
////////////////////////////////////////////////////////////////////////////////////
double polynomial(const double* pa,const unsigned a_len,const double x)
//pa为指向系数数组的指针，a_len为数组的长度，返回将x带入多项式的结果(多项式为升幂排列)
//核心为秦九韶(shao2)算法
{
#ifdef _DEBUG
    cout<<pa[0];
    for(int j(1);j<a_len;++j)
    {
           cout<<'+'<<'('<<pa[j]<<'*'<<x<<'^'<<j<<')';
    }
    cout<<'=';
#endif
    double b=*(pa+a_len-1);//
    for(int i(a_len-2);i>=0;--i)
    {
           b=*(pa+i)+x*b;
    }
    return b;
}
//以下为测试代码
int main()
{
    double a[]={1.5,2,3};
    int a_len=3;
    double x=2.0;
    cout<<polynomial(a,a_len,x)<<endl;
    return 0;
}