请问为什么在电脑上秦九韶算法快不了多少啊？

sunnyrubik

不是都说秦九韶比普通的写法快几倍吗？

1. #include<stdio.h>
   
2. #include<stdlib.h>
   
3. #include<math.h>
   
4. #include<time.h>
   
5. 
   
6. #define MAXK 1e8 //被测函数最大重复调用次数
   
7. 
   
8. double f1(int n,double a[],double x);
   
9. double f2(int n,double a[],double x);//秦九韶算法
   
10. 
    
11. clock_t start,stop;
    
12. double duration;
    
13. 
    
14. int main()
    
15. {
    
16.         int i;
    
17.         double a[9];
    
18.        
    
19.         for(i=0;i<10;i++)
    
20.         {
    
21.                 a[i]=(double)i;
    
22.         }
    
23.          
    
24.         start = clock();
    
25.         for(i = 0;i < MAXK;i++)
    
26.         f1(9,a,1.1);
    
27.         stop = clock();
    
28.         duration = ((double)(stop-start))/CLK_TCK/MAXK;//计算函数单次运行时间
    
29.         printf("ticks1 = %f\n",(double)(stop-start));
    
30.         printf("duration = %6.2e\n",duration);
    
31.        
    
32.         start = clock();
    
33.         for(i = 0;i < MAXK;i++)
    
34.         f2(9,a,1.1);
    
35.         stop = clock();
    
36.         duration = ((double)(stop-start))/CLK_TCK/MAXK;
    
37.         printf("ticks2 = %f\n",(double)(stop-start));
    
38.         printf("duration = %6.2e\n",duration);
    
39.        
    
40.         return 0;
    
41. }
    
42. 
    
43. double f1(int n,double a[],double x)
    
44. {
    
45.         int i;
    
46.         double p = a[0];
    
47.         for(i=1;i<=n;i++)
    
48.         {
    
49.                 p += (a[i] * pow(x,i));
    
50.         }
    
51.         return p;
    
52. }
    
53. 
    
54. double f2(int n,double a[],double x)
    
55. {
    
56.         int i;
    
57.         double p = a[n];
    
58.         for(i=n;i>0;i--)
    
59.         {
    
60.                 p = a[i-1] + x*p;
    
61.         }
    
62.         return p;
    
63. }

复制代码

sunnyrubik

已解决。。。。。。