快速pow和sqrt的小技巧 hdu4282

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的方法~~


如下:

[cpp] view plaincopy
  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的代码:

[cpp] view plaincopy
  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  

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/382106.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

466. 统计重复个数 golang[转]

转载&#xff08;mark&#xff09; https://blog.by24.cn/archives/leetcode-count-the-repetitions.html func getMaxRepetitions(s1 string, n1 int, s2 string, n2 int) int {len1, len2 : len(s1), len(s2)index1, index2 : 0, 0 // 注意此处直接使用 Ra Rb 的下标&#…

pow(x,y)函数

实现浮点类型的幂运算&#xff0c;函数原型为: double pow(double x, int n) 在求解这个问题的时候是一个很挣扎的过程&#xff0c;因为它不是报错而是一直提示你超出时间&#xff0c;那么必须一次次的考虑怎样降低时间复杂度。 首先最直接的思路是下面这样的&#xff0c;就跟直…

LeetCode 303,560,1248 (前缀求和 )

303. 区域和检索 - 数组不可变 给定一个整数数组 nums&#xff0c;求出数组从索引 i 到 j (i ≤ j) 范围内元素的总和&#xff0c;包含 i, j 两点。 示例&#xff1a; 给定 nums [-2, 0, 3, -5, 2, -1]&#xff0c;求和函数为 sumRange() sumRange(0, 2) -> 1 sumRange…

ZOJ 2060----Fibonacci Again

没什么好说的&#xff0c;就是要求斐波那契&#xff0c;问该数是否能够被三整除&#xff0c;写出了答案的个位便可发现规律。code&#xff1a;#include <iostream>using namespace std;int main() { int n; while (cin>>n) { if ((n - 2) % 4!0) cout<…

分布式是写出来的(一)

分布式对象存储笔记 实现一个单机版本的对象存储 package mainimport ("io""log""net/http""os""strings" )func Handler(w http.ResponseWriter, r *http.Request) {m : r.Methodif m http.MethodPut {put(w, r)return}…

分布式是写出来的(二)

从单机存储进化为接口和存储的分离 概述 接口服务层对外提供REST服务&#xff0c;数据服务层提供数据存储功能。两者之间通过消息队列进行通信&#xff0c;数据服务层的所有数据服务注册dataServer Exchange&#xff0c; 以便client给接口服务层发消息后&#xff0c;接口服务…

ZOJ 1295——Reverse Text

在没学STL之前&#xff0c;可能做这道题需要开数组&#xff0c;在出现空格的地方外加判断&#xff0c;但是学完string&#xff0c;一条reverse语句便将输入逆转&#xff0c;c语音可能需要2-30行&#xff0c;可见c的优化。注意&#xff1a;输入有空格&#xff0c;要用getline(ci…

分布式是写出来的(三)

添加元数据服务 元数据服务就是对元数据提供存取功能的服务。元数据就是系统定义的基本信息&#xff0c;比如一张相片的名字&#xff0c;版本&#xff0c;拍摄时间&#xff0c;散列值等。客户端和接口服务之间根据对象的名字来引用一个对象&#xff0c;一个对象可以有多个版本…

ZOj 2104——Let the Balloon Rise

老师在上课的时候讲过这道题&#xff0c;不过当时做这道题时还是纠结了许久&#xff0c;那时stl不熟&#xff0c;老是想着用数组&#xff0c;去重很麻烦&#xff0c;学了STL后&#xff0c;用map就简单多了。code &#xff1a;#include <iostream>#include <string>…

分布式是写出来的(四)

数据去重 为了防止同一份数据上传多次&#xff0c;我们需要进行hash值校验。 启动服务时&#xff0c;扫描所有文件&#xff0c;计算sha256值&#xff0c;存储在hash表中。上传文件时&#xff0c;在put请求中添加自己计算的sha256值&#xff0c;hash表内的值进行比较&#xff…

Zoj 2947——Abbreviation

呵呵哒&#xff0c;一开始觉的很难&#xff0c;是因为题长&#xff0c;以后不能被长题再坑了&#xff01;后来耐下性子读&#xff0c;才发现就是让比对缩写是否相同&#xff0c;我用字符数组写的&#xff0c;当时还开了二维数组&#xff0c;现在想来&#xff0c;实在是笨&#…

分布式是写出来的(五)

数据冗余策略 RS(Reed Solomon Coding)纠删码 在存储系统中&#xff0c;需要采用数据冗余技术来保证数据的可靠性&#xff0c;相比使用多副本复制机外&#xff0c;使用纠删码能够以更小的数据冗余度获得更高的数据可靠性。 RS纠删码将原文件分成n个数据块&#xff0c;同时为这…

ZOJ 1151——Word Reversal

有是一个字符串翻转问题&#xff0c;唉&#xff0c;自以为处理的很巧妙&#xff0c;因为如果整句输入可能真的不好处理&#xff0c;呵呵&#xff0c;反正没这样简单&#xff01;code&#xff1a;#include <cstring>#include <cstdio>#include <iostream>#inc…

分布式是写出来的(六)

断点下载 GET 当服务端收到GET请求&#xff0c;服务端不会把整个对象返回给客户端&#xff0c;服务端首先做SEEK&#xff0c;查找客户端提供的Range: bytesfirst的字节数&#xff0c;从0-first的内容服务端直接丢弃&#xff0c;那么服务端从first开始传递数据 POST 如果客户…

【转载】输出二进制 C

想知道怎样输出二进制。我们知道C输出十六进制是cout〈〈hex〈〈 a&#xff1b;而八进制是cout〈〈 ocx〈〈 a;二进制则没有默认的输出格式&#xff0c;需要自己写函数进行转换&#xff0c;于是上网搜索了一下。网上思路真是广泛啊。下面列出一些方法。 #include 〈iostream…

6.824 Raft lesson4 2020(一)

raft实现 距离上一篇文章一个月&#xff0c;因为6.824的课程看不懂&#xff0c;基础知识薄弱。现在了解一点Raft算法&#xff08;自己动手实现一遍&#xff09;还需要其他分布式相关的基础知识&#xff08;实现一个分布式对象存储系统&#xff09;&#xff0c;然后再去继续学习…

#define和const的区别

(1) 编译器处理方式不同 define宏是在预处理阶段展开。 const常量是编译运行阶段使用。 (2) 类型和安全检查不同 define宏没有类型&#xff0c;不做任何类型检查&#xff0c;仅仅是展开。 const常量有具体的类型&#xff0c;在编译阶段会执行类型检查。 (3) 存储方式不同 defin…

1190. 反转每对括号间的子串 golang反转字符串

1190. 反转每对括号间的子串 557. 反转字符串中的单词 III 344. 反转字符串 golang 这个题比较新&#xff0c;但是面试竟然被问到了。10分钟之内没写出来 给出一个字符串 s&#xff08;仅含有小写英文字母和括号&#xff09;。 请你按照从括号内到外的顺序&#xff0c;逐层…

【转载】define与typedef区别与联系

define与typedef区别与联系 1) #define是预处理指令&#xff0c;在编译预处理时进行简单的替换&#xff0c;不作正确性检查&#xff0c;不关含义是否正 确照样带入&#xff0c;只有在编译已被展开的源程序时才会发现可能的错误并报错。例如&#xff1a; #define PI 3.1415926 程…

LeetCode(合集)两数之和总结 (1,167,1346)

1.两数之和 167. 两数之和 II - 输入有序数组 golang 1346. 检查整数及其两倍数是否存在 一个数组是否存在两个数是target的整数倍&#xff1f; 如果有人知道这是LeetCode哪个题&#xff0c;给我说一下 给每个数做整除&#xff0c;然后在遍历