clock()函数调用模版

课程主要讲了对多项式的优化，

给定一个多项式p=a0+a1*x+a2*x^2+...+an*x^n;

正常解法：代码如下

double f1(int n,double a[],double x)
{
int i=0;
double p=a[0];
  for(i=1;i<n;i++)
  {
  p+=(a[i]*pow(x,i));
 }
 return p;
}

从第一项开始一直到第n项；

优化后的解法：代码如下

double f2(int n,double a[],double x)
{
int i;
double p=a[n];
for(i=a[n];i>=1;i--)
{
p=(a[i-1]+x*p);
}
return p;
}

每项提取一个x

a0+x*(a1+x*(a2+...+x*(an-1+x*an))）

经过对算法的优化后时间发生了指数级的减小；

源代码：

#include <cstdio>
#include <math.h>
#include <time.h>
clock_t start,stop;
double duration;
#define MAXN 10
#define MAXK 1e7
double f1(int n,double a[],double x)
{
int i=0;
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=a[n];i>=1;i--)
{
p=(a[i-1]+x*p);
}
return p;
}
int main()
{   
int i;
double a[MAXN];
double k;
for(int i=0;i<MAXN;i++) a[i]=(double)i;
start=clock();
for(i=0;i<MAXK;i++)
{
f1(MAXN-1,a,1.1);
}
// printf("%f ",k);
stop=clock();
duration=((double)((stop-start))/CLK_TCK)/MAXK;
printf("ticks1=%f   ",(double)(stop-start));
printf("duration=%6.2e\n",duration);
//函数2 
start=clock();
for(i=0;i<MAXK;i++)
{
f2(MAXN-1,a,1.1);
}
// printf("%f ",k);
stop=clock();
duration=((double)((stop-start))/CLK_TCK)/MAXK;
printf("ticks2=%f   ",(double)(stop-start));
printf("duration=%6.2e\n",duration);
}

这里由于一个函数跑起来非常快，所以让函数进行循环，得出结果。

得出结论：解决问题方法的效率，跟算法的巧妙程度有关。