励志用少的代码做高效的表达


题目(提交)链接→UVA-11809


算是个数学题吧,虽然在AOAPC上面给放到象征水题的第三章里面了。

这个题基本就是帮着你复习了一遍浮点数的存储方式了。浮点数在计算机里是分三部分表示的,最前面一位表示符号,后面一部分是尾数,最后一部分是阶码,表示方法类似于科学记数法,不过是二进制的,尾数是M阶码是E的话那么表示起来就是M × 2^E了。然后对于M还有一个要求,就是1/2 ≤ M < 1,所以用二进制表示M的话就应该是0.1XX……,用计算机表示的时候就把最前面的“0.1”这个永远不变的部分给省略掉,只表示可能变化的部分。阶码部分则是只用二进制表示E。
  
    上面的图就给出了一个例子,前面的0表示是正数。后面8位表示尾数m,这里是0.111111111(注意后面是9个1,因为头一个省略了)。之后那个0表示分割,最后面6位表示e的二进制为111111。所以这个数就是,用十进制表示就是。

在计算机中用二进制表示M和E的时候如果位数不同,那么它们所能表示的最大值也不同。现在给你所能够表示的最大的浮点数的值,让你倒回去求M和E分别有多少位。输入格式为AeB,表示最大浮点数为,并且0 < A < 10,并且保证输出的结果中0 ≤ M ≤ 9且1 ≤ E ≤ 30。输入以0e0表示结束,0e0本身不计算。

这个如果直接去算的话相当麻烦,当E很大的时候数会直接超出上限。这个时候可以反过来想,最大的时候M和E的每一位肯定都是1,并且又有0 ≤ M ≤ 9且1 ≤ E ≤ 30的限定,所以一共只有300种情况,自然就想到了打表,先用二重循环枚举M和E可能出现位数的所有情况打一张表,然后输入的时候倒回去找即可。

假设当前一层M和E的值为m和e,它们的位数分别为i和j。

首先计算m的值,用二进制表示的话,m的值为0.11…,也就是m = 2^(-1) + 2^(-2) + … + 2^(-1 - i)(i比实际1的个数少1个),也就是m = 1 - 2^(-1 - i)。

接下来就是计算e的值,不难得出,e = 2^j - 1。

那么也就有m * 2^e = A * 10^B,似乎可以直接计算了。然而,直接这样算的话是不行的,因为当e太大的话(e最大可以是1073741823,注意这还只是2的指数),等号左边的数就会超出上限,所以要想继续算下去,就得自己去想办法再写出满足要求的类来,这显然太麻烦了。所以,这个时候我们对等式两边同时取对数,这个时候就有 log10(m) + e × log10(2) = log10(A) + B。因为此时m和e的值都是确定的,所以不妨令等式左边为t,也就有t = log10(A) + B。

这个时候就有问题了,A和B怎么算。

写题解的时候突然意识到了这个问题,读题的时候很多人,包括我,都把AeB默认为了科学记数法,在ACM协会群里面讨论的时候很多人也都说这是科学计数法。先来看如果是科学记数法的时候应该怎么办。

如果是科学记数法的话,那么对于A,就有1 ≤ A < 10。那么0 < log10(A) < 1。所以t的小数部分就是log10(A),整数部分就是B,即B = ⌊t⌋,A = 10^(t - B)。那么接下来,我们只需要开出两个二维数组来,分别记录对应i和j下A和B的大小,之后从输入里提取出A和B的大小,去二维数组里面查找对应的i和j即可。

这种办法在UVA上面是可以直接AC的,但是我却感觉这题这样A了有点数据太水的感觉,秉着处女座+强迫症死磕到底的精神,我们看下哪里有问题。

其实回头读下题,我们发现科学记数法1 ≤ A < 10的条件是我们脑补出来的,题目里面根本没有提及,只是简单交待0 < A < 10。也就是说,对于确定的M和E的位数,十进制的表示可以有多种,例如样例中的5.699141892149156e76,下面的数据应当也是完全可能的,而且结果应当与样例的结果是相同的(当然是在保证精度可以计算出结果的前提下):

0.569914189214915e77
0.056991418921491e78
0.005699141892149e79
0.000569914189214e80

带着这个想法我分别拿着上面的数据去UVA toolkit和uDebug上试了试,UVA toolkit依旧能够输出“5 8”的结果来,但是uDebug告诉我我的输入不合法……果真是我想多了么……
  不过这个问题也好办,还是看上面的数据,忽略掉后面几位精度丢失的问题的话,上面的几个数完全可以通过“A *= 10, B -= 1”或者“A /= 10, B += 1”的操作来相互转化。那么对于0 < A < 1的A的值,我们就可以通过“A *= 10, B -= 1”的操作来使其满足科学记数法的条件。

另外,在查表的时候还应该注意精度的问题,15位有效数字对于double来说精度似乎也不够,而且计算出所需要的整数值其实需要的精度也没有那么高,所以这里的精度就只用到了1e-4的程度。


代码:

#include <bits/stdc++.h>
using namespace std;
int main() {double M[20][40];long long E[20][40];// 打表for(int i = 0; i <= 9; ++i) for(int j = 1; j <= 30; ++j) {double m = 1 - pow(2, -1 - i), e = pow(2, j) - 1;double t = log10(m) + e * log10(2);E[i][j] = t, M[i][j] = pow(10, t - E[i][j]);}// 输入并输出结果string in;while(cin >> in && in != "0e0") {// 处理输入for(string::iterator i = in.begin(); i != in.end(); ++i) if(*i == 'e') *i = ' ';istringstream ss(in);double A; int B;ss >> A >> B;// 在打好的表中寻找答案for(int i = 0; i <= 9; ++i) for(int j = 1; j <= 30; ++j) if(B == E[i][j] && (fabs(A - M[i][j]) < 1e-4) {cout << i << ' ' << j << endl;break;}}
}

拨云见日,未来可期

28行代码AC——习题3-12 浮点数(UVA 11809 - Floating-Point Numbers)——解题报告相关推荐

  1. 11行代码AC——习题2-4 子序列的和(subsequence)——解题报告

    励志用尽量少的代码做高效的表达. 题目描述: 输入两个正整数n<m<106,输出1/(n²)+1/((n+1)²)+--+1/(m²),保留5位小数.输入包含多组数据,结束标记为n=m=0 ...

  2. 31行代码AC——PTA 求二叉树的叶子结点个数 (20分)——解题报告

    励志用尽量少的代码做高效的表达. 以二叉链表作为二叉树的存储结构,求二叉树的叶子结点个数. 输入格式: 输入二叉树的先序序列. 提示:一棵二叉树的先序序列是一个字符串,若字符是'#',表示该二叉树是空 ...

  3. 21行代码AC——习题5-1 代码对齐(Alignment of Code, UVa1593)——解题报告

    题意: 输入若干行代码,要求各列单词的左边界对齐且尽量靠左,单词之间至少要空一格,每个单词不超过80个字符,每行不超过180个字符,一共最多1000行. 思路: 1.输入内容存入二维数组 2.找出每列 ...

  4. 26行代码AC——习题3-2 分子量 (UVa1586,Molar Mass)——解题报告

    大意: 给出分子式,式中只包含以下四种元素.求分子量. C:12.01 H: 1.008 O: 16.00 N: 14.01 Sample Input 4 C C6H5OH NH2CH2COOH C1 ...

  5. 28行代码AC——Minimum Sum LCM UVA - 10791(最大质因子)

    励志用尽量少的代码做高效表达 题目(提交)链接-->UVA - 10791 题意 输入正整数n,找至少两个数,使得他们的最小公倍数(Least Common Multiple)为n且要输出最小的 ...

  6. 21行代码AC——习题3-7 DNA序列(UVa-1368)_解题报告

    励志用尽量少的代码做高效表达. 题目(提交)链接→UVa-1368 思路: DNA序列:按列遍历,记录每一列出现次数最多(若同样多,则字典序最小)的字母,录入s串累加. 距离:重新遍历,录入出现次数比 ...

  7. 22行代码AC——习题5-6 对称轴(Symmetry,UVa1595)——解题报告

    励志用尽量少的代码做高效的表达. 题目(提交)链接→UVa-1595 思路: 此题本质是一道笛卡尔坐标系上的对称性问题. 判定性问题:由于只要能判别图像是否左右对称即可,无需确认关于哪条垂直线对称,那 ...

  8. 15行代码AC——习题3-3 数数字 (UVa1225,Digit Counting)

    大意: 把n(n<=10000)个整数顺序写在一起,求0~9分别出现多少次 Sample Input 2 3 13 Sample Output 0 1 1 1 0 0 0 0 0 0 1 6 2 ...

  9. 15行代码AC——习题3-1 得分 (UVa1585,Score)

    题意概述 给出一个由O和X组成的串,统计得分,O得分为目前连续出现的O的个数,X的得分为0. 如:OOXXOXXOOO得分为:1+2+0+0+1+0+0+1+2+3=10 Sample Input 5 ...

最新文章

  1. C#/Net代码精简优化技巧(3)
  2. Android 博客园客户端 (八) 下拉刷新、分页,AsyncTask
  3. 【NetApp】ROOT VOLUME NOT WORKING PROPERLY: RECOVERY REQUIRED
  4. [BZOJ]3436: 小K的农场
  5. 用学生编程记录预测学习成果,第二届计算机教育数据挖掘大赛, 赢取现金奖励+顶刊发表机会!...
  6. as3.0横向渐变发光字
  7. Word2013 设置默认缩进
  8. php 使用redis队列简单实用
  9. cad计算机试题及答案,CAD考试试题库及参考答案整理版.doc
  10. 燕东微通过注册:预计年营收超20亿 亦庄国投与京东方是股东
  11. device-side assert triggered原因和解决方法
  12. fastjson 是反射吗_6种超声检测灵敏度,你都知道吗?
  13. 电脑右下角的扬声器出现了红叉
  14. python考研成绩查询_2020-09-08考研成绩预测模型
  15. @Value为null问题排查及解决方法
  16. LDAP认证服务客户端配置
  17. (转)QQ在线客服代码
  18. 理解pandas的groupby().apply()
  19. unity3D---3Dmaxs给带骨骼蒙皮的模型绑定bip
  20. 环模制粒机设计(说明书+CAD)

热门文章

  1. Python 中的url,Base64和MD5编码解码的使用
  2. Django框架实现支付宝第三方支付
  3. linux常用安装命令集锦
  4. 经典|深入理解 Go高性能网络编程模型
  5. 过滤请求绝技 — 布隆过滤器与布谷鸟过滤器
  6. LiveVideoStack Meet | 苏州:视频会议研发中心一日游
  7. WebRTC/Chromium在2020年的更新
  8. 音视频技术开发周刊 83期
  9. c/c++比较灵活的方法:回调函数和函数指针
  10. CentOS 7 的 Python2 升级成 Python 3.7