霍纳规则计算多项式

霍纳规则

假设有n＋2个实数a0，a1，…，an,和x的序列，要对多项式Pn(x)= anxn ＋an－１xn－１＋…＋a１x＋a０求值，直接方法是对每一项分别求值，并把每一项求的值累加起来，这种方法十分低效，它需要进行n＋(n－1)＋…＋1＝n(n＋１)/2次乘法运算和n次加法运算。有没有更高效的算法呢?答案是肯定的。通过如下变换我们可以得到一种快得多的算法，即Pn(x)= anxn ＋an－１xn－１＋…＋a１x＋a０＝((…(((anx ＋an－１)x＋an－２)x+ an－3)…)x＋a1)x＋a０，这种求值的安排我们称为霍纳法则。

例如，当x=3时，计算p(x)=2x4－x3＋3x2＋x－5的值。对于多项式p(x)=2x4－x3＋3x2＋x－5，我们按霍纳法则进行变换，有：

p(x)=2x4－x3＋3x2＋x－5

=x(2x3－x2＋3x＋1)－5

=x(x(2x2－x＋3)＋1)－5

=x(x(x(2x－1)＋3)＋1)－5

#include <iostream>
#include <ctime>
using namespace std;

int horner(int *a,int n,int x)
{
	int y=0;
	for(int i=n-1;i>=0;--i){
		y=a[i]+x*y;
	}
	return y;
}

int naive_poly(int *a,int n,int x)
{
	int y=0,k=1;
	for(int i=0;i<n;++i){
		/********************/
		/* 解法一           */
		/*********************/
		/*int k=1;                    
		for(int j=0;j<i;++j){
		k=x*k;
		}
		y=a[i]*k+y
		*/
		/*****************/
		/* 解法二      */
		/**************/
		k*=x;
		y=a[i]*k/x+y;
	}
	return y;
}

int main()
{
	srand(time(NULL));
	int num;
	while ((num=rand()%8)<3);
	//num=4;
	int *a=new int[num];
	for(int i=0;i<num;++i){
		a[i]=rand()%10;
	}
	int x=2;

	for(int i=0;i<num;++i){
		cout<<a[i]<<" ";
	}
	cout<<endl;

	int result=horner(a,num,x);
	int result2=naive_poly(a,num,x);
	
	cout<<result<<endl;
	cout<<result2<<endl;
}