快速pow和sqrt的小技巧 hdu4282

本文探讨了解决复杂数学竞赛题目时所遇到的挑战,并提供了优化算法的策略来避免时间限制错误。通过引入快速幂运算和快速求逆平方根的方法,成功解决了难题。文章还分享了最终通过AC验证的代码实现,包括快速幂算法和牛顿迭代法的使用,为解决类似问题提供了一套实用的解决方案。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

http://acm.hdu.edu.cn/showproblem.php?pid=4282

今年网络赛。。天津赛区。。有道题。。是这样的。。。X^Z + Y^Z + XYZ = K  给出K ,求XYZ,我思路很明确。。。枚举其二,然后二分其一,但是始终TLE。。。。晚上回去之后,看了人家报告,。。。才发现。。。原来是微软的函数pow 惹的祸。。。我本来以为微软的函数写的都很好。。效率很高。。但是我忘了一件事。。就是 pow(int a, int b) 这个函数可以理解为 pow(double a, double b)  之所以TLE  。。是因为微软提高的兼容性。。而导致时间效率的底下。。(但是令我郁闷的就是。。网上有位仁兄。。也用的pow。。居然984ms 压线过了这道题。。。汗!= =||)。。所以我就百度了一下。。找到了快速pow 和sqrt的方法~~


如下:

  1. **=============================================**  
  2. ||快速pow(多次使用时及其有效)               ||  
  3. **=============================================**  
  4. __int64 qpow(int a, int b){  
  5.     __int64 c, d; c = 1; d = a;  
  6.     while (b > 0){  
  7.         if (b & 1)  
  8.           c *= d;  
  9.         b = b >> 1;  
  10.         d = d * d;  
  11.     }  
  12.     return c;  
  13. }  
  14. **=============================================**  
  15. ||快速1/sqrt(x)(牛顿迭代法)                    ||  
  16. **=============================================**  
  17. float InvSqrt (float x) {  
  18.      float xhalf = 0.5f * x;  
  19.      long i = *(long*)&x;  
  20.      i = 0x5f3759df - (i >> 1);  
  21.      x = *(float *)&i;  
  22.      x = x * (1.5f - xhalf * x * x);  
  23.      return x;  
  24. }  

再列出我最后AC的代码:

  1. #include <stdio.h>  
  2. #include <math.h>  
  3.   
  4. __int64 k;  
  5. //其实我也试着写了一个pow。。只不过弱爆了。。不够快。。  
  6. /* 
  7. __int64 myqpow(__int64 x, __int64 y) 
  8. { 
  9.     if (y == 1) 
  10.     { 
  11.         return x; 
  12.     } 
  13.     __int64 tmp = qpow(x, y/2); 
  14.     if (y%2 == 1) 
  15.     { 
  16.         return (tmp * tmp * x); 
  17.     } 
  18.     else 
  19.     { 
  20.         return (tmp * tmp); 
  21.     } 
  22. */  
  23.   
  24. __int64 qpow(int a, int b)  
  25. {  
  26.     __int64 c, d; c = 1; d = a;  
  27.     while (b > 0)  
  28.     {  
  29.         if (b & 1)  
  30.           c *= d;  
  31.         b = b >> 1;  
  32.         d = d * d;  
  33.     }  
  34.     return c;  
  35. }  
  36.   
  37. __int64 solve(__int64 x, __int64 z)  
  38. {  
  39.     __int64 l = x + 1, r = 32768, y = (l + r) >> 1;//r = k - qpow(x, z)  
  40.     while (l <= r)  
  41.     {  
  42.         __int64 tmp = x * y * z + qpow(x, z) + qpow(y, z);  
  43.         if (tmp == k)  
  44.         {  
  45.             return y ;  
  46.         }  
  47.         else if (tmp > k || tmp < 0)  
  48.         {  
  49.             r = y - 1;  
  50.         }  
  51.         else  
  52.         {  
  53.             l = y + 1;  
  54.         }  
  55.         y = (l + r) >> 1;  
  56.     }  
  57. /* 
  58.     if (x * y * z + qpow(x, z) + qpow(y, z) == k) 
  59.     { 
  60.         return y ; 
  61.     } 
  62. */  
  63.     return 0 ;  
  64. }  
  65.   
  66. int main()  
  67. {  
  68.     while (scanf("%I64d", &k), k)      
  69.     {  
  70.         __int64 i, j, ans=0;  
  71.         for (i = 2; i <= 31; i++)  
  72.         {  
  73.             for (j = 1;  ; j++)  
  74.             {  
  75.                 __int64 tmp1 = qpow(j, i);  
  76.                 if (tmp1*2 > k || tmp1 < 0)  
  77.                 {  
  78.                     break ;  
  79.                 }  
  80.   
  81. //                __int64 a = solve(j, i);  
  82.                 if (solve(j, i))  
  83.                 {  
  84.                     ans++;  
  85. //                    printf("%I64d %I64d %I64d\n", j, a, i);  
  86.                 }  
  87.             }  
  88.         }  
  89.         printf("%I64d\n", ans);  
  90.     }  
  91.     return 0;  
  92. }  
  93.   
  94. //  5 6 5  
  95. //  3125 7776 150  
  96. //  11051  
  97. //  40 48 5  
  98. //  102400000  254803968 9600  
  99. //  357213568  
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值