#include <stdio.h>
#include <time.h>//引入clock()计时函数
#include <math.h>
clock_t start, stop;//clock_t是变量类型(typedef long clock_t);
double duration;//用于计算秒数
#define MAXN 10//多项式的最大项系数,注意从0开始
#define MAXK 1e7//设置重复调用次数,从而有时间积累,看到函数所花费的时间
double f1(int n, double a[], double x);
double f2(int n, double a[], double x);
int main(void) {
int i;
double a[MAXN];//存储多项式的系数的数组
for (i = 0; i < MAXN; i++)
a[i] = (double)i;
//测试第一个函数运行时间
start = clock();
for (i = 0; i < MAXK; i++)
f1(MAXN - 1, a, 1.1);
stop = clock();
duration = ((double)(stop - start)) / CLK_TCK;//(clock_t)1000
printf("ticks1 = %f\n", (double)(stop - start));
printf("duration1 = %6.2e\n", duration);
//测试第二个函数运行时间
start = clock();
for (i = 0; i < MAXK; i++)
f2(MAXN - 1, a, 1.1);
stop = clock();
duration = ((double)(stop - start)) / CLK_TCK;
printf("ticks2 = %f\n", (double)(stop - start));
printf("duration2 = %6.2e\n", duration);
return 0;
}
//------------------------
//
//计算给定多项式在给定点处的值--探究算法设计效率问题
// n:阶数n
// a[]:存储多项式系数的数组
// x:待计算的这个点
// f(x) = a0+a1*x+a2*x^2+...+an-1*x^(n-1)+an*x^n
//
//-----------------------
double f1(int n, double a[], double x) {
int i;
double p = a[0];
for (i = 1; i <= n; i++)
p += (a[i] * pow(x, i));
return p;
}
//秦九昭算法
double f2(int n, double a[], double x) {
int i;
double p = a[n];
for (i = n; i > 0; i--)
p += a[i - 1] + x * p;
return p;
}
运行结果:
总结:
计算多项式:
用秦九昭算法比普通计算快了将近一个数量级的样子,所以说程序运行的效率还和算法设计有关