c语言里面，怎么求x的5次方？

pupilzeng

#include<math.h>

double pow (double x, double y) /* x^y */

DoDo

最初由 zug 发表
A=X;
X=X*X;
X=X*X;
X=X*A;

对于X^Y，Y为整数的情况，还是这种方法最好，速度最快

把它推广到一般

1. 
   
2. #include <math.h>
   
3. 
   
4. long int mypow(double x, int y)
   
5. {
   
6.   if(y == 0)
   
7.   {
   
8.     return 1;
   
9.   }
   
10.   int t = (int)(log(y) / log(2));
    
11.   int i;
    
12.   int r = x;
    
13.   int d = 1;
    
14.   for(i=0;i<t;i++)
    
15.   {
    
16.     d *= 2;
    
17.     r *= r;
    
18.   }
    
19.   int d = y - d;
    
20.   if(d)
    
21.   {
    
22.     return r * mypow(x, d);
    
23.   }
    
24.   else
    
25.   {
    
26.     return r;
    
27.   }
    
28. }
    

复制代码

程序没有经过调试，不知是否有问题

DoDo

不知道系统的函数是如何实现的，有人能提供源码吗？

pupilzeng

1. 
   
2. double mypow3 (double x, long unsigned y)
   
3. {
   
4.         if (y == 0)
   
5.                 return 1;
   
6.         double z = mypow3 (x, y / 2);
   
7.         if (y % 2)
   
8.                 return z * z * x;
   
9.         else
   
10.                 return z * z;
    
11. }
    
12. 
    
13. double mypow4(double x, long unsigned y)
    
14. {
    
15.         if (y == 0)
    
16.                 return 1;
    
17.         long i = 1;
    
18.         double z = x;
    
19.         while ( i * 2 <= y )
    
20.          {
    
21.                 x *= x;
    
22.                 i *= 2;
    
23.          }
    
24.         return x * mypow4(z, y - i);
    
25. }
    
26. 
    
27. 
    
28. int main()
    
29. {
    
30.         long unsigned i;
    
31.         double z;
    
32.         for(i = 0; i < 9999999; i++)
    
33.                 z = mypow3(1.001,10000);
    
34.         printf("%lf\n", z);
    
35. }
    

复制代码

测试:
>gcc -o mypow3 mypow.c
>time ./mypow3
21916.681339

real    0m3.912s
user    0m3.900s
sys     0m0.010s


把main中的mypow3改成mypow4
>gcc -o mypow4 mypow.c
>time ./mypow4
21916.681339

real    0m6.428s
user    0m6.400s
sys     0m0.030s

我写mypow4的目的是想减少mypow3的递归次数。
可是没想到性能反而下降了，看来递归的性能还是蛮高的。

wide288

$ time ./pow-1
21916.681339

real    0m7.504s
user    0m7.490s
sys     0m0.010s

我机器是不是太差了。后面这三项是什么意思啊？
我用的是  mypow3() 这个。

kj501

raal是实际占用的时间，user是处于用户模式的时间，sys是处于系统模式的时间。real大致上是它们两者之和。

wide288

哦，谢谢了
那后面的单位是
real 0m7.504s
7秒吗

iDay

去看math.h

pupilzeng

一大堆的宏定义。
看不懂。