我使用c/c++中的时钟点clock()函数(获取当前时间点)和clock_t变量(容纳clock返回值)和CLK_TCK常量(1秒钟有多少时钟点)进行效率计算。
问题:
使用秦九韶 算法思路
一般地,一元n次多项式的求值需要经过(n+1)*n/2次乘法和n次加法,而秦九韶算法只需要n次乘法和n次加法。在人工计算时,一次大大简化了运算过程。
把一个n次多项式
改写成如下形式:
求多项式的值时,首先计算最内层括号内一次多项式的值,即
然后由内向外逐层计算一次多项式的值,即
这样,求n次多项式f(x)的值就转化为求n个一次多项式的值。
结论:对于一个n次多项式,至多做n次乘法和n次加法
一、代码实现
#include <iostream>
#include <ctime>
#include <cmath>
#include <iomanip>
using namespace std;
/*记录开始和结束时的时钟点*/
clock_t start,stop;
/*方法一:使用循环累加*/
double f1(int N,double x)
{
double p = 0;
for(int i = 1; i <= N;i++)
{
p+= i * pow(x,i);
}
return p;
};
/*方法二:使用秦九韶算法*/
double f2(int N,double x)
{
double p = N;
for(int i = N; i > 0; i--)
{
p = i-1 + x * p;
}
return p;
};
/*根据传参选择运行函数*/
void the_clock(int N, double x,double temp(int ,double))
{
cout << setprecision(0) << fixed;
/*记录开始时钟点*/
start = clock();
/*因为一次太快了,如果小于一,int无法捕捉*/
for(int i = 0; i < 10000000; i++)
{
temp(N,x);
}
/*记录结束时钟点*/
stop = clock();
/*打印计算结果*/
cout << temp(N,x) << endl;
/*计算所用时间*/
double time = (double)(stop - start)/CLK_TCK/10000000;
/*打印时间*/
cout << setprecision(20) << fixed;
cout << time << endl;
}
int main(int argc, const char * argv[]) {
the_clock(20, 1, f1);
the_clock(20, 1, f2);
return 0;
}
二、输出
210
0.00121658700000000009
210
0.00050331699999999998
三、效率差距
从输出结果可以看出秦九韶算法至少比普通算法高一倍的速度。
要点:编译器在运行+-的速度比*/的速度快很多,秦九韶算法每次循环只有一次乘法运算,而普通循环算法使用pow算法,使得使用乘法运算的次数以几何增长,所以秦九韶算法快很多。