BZOJ 1213 [HNOI2004]高精度开根 牛顿法
来源:互联网 发布:脸书注册网络错误 编辑:程序博客网 时间:2024/05/17 02:30
Description
晓华所在的工作组正在编写一套高精度科学计算的软件,一些简单的部分如高精度加减法、乘除法早已写完了,现在就剩下晓华所负责的部分:实数的高精度开m次根。因为一个有理数开根之后可能得到一个无理数,所以这项工作是有较大难度的。现在要做的只是这项工作的第一步:只对自然数进行开整数次根,求出它的一个非负根,并且不考虑结果的小数部分,只要求把结果截断取整即可。程序需要根据给定的输入,包括需要开根的次数,以及被开根的整数;计算出它的非负根取整后的结果。
Input
共有两行,每行都有一个整数,并且输入中没有多余的空格:第一行有一个正整数m (1 <= m <= 50),表示要开的根次;第二行有一个整数n (0<=n <= 10^10000),表示被开根的数。
Output
只有一行,包括一个数,即为开根取整后的结果。
Sample Input
3
1000000000
1000000000
Sample Output
1000
HINT
传送门
此题求x^m=a的x的下取整整数解。
另f(x)=x^m-a,则f'(x)=m*x^(m-1)
接下来代进入式子里进行牛顿迭代,,
学了下python。。然后就python水过了。。
python免了高精度花式AC啊日= =
注意精度问题……一开始WA了几发。。最后发现精度差挺大的,暴力逼近就好了。。
好像都写的二分啊= =
m=int(input())n=int(input())if n==0: print 0else: len=0;t=n while t>0:t/=10;len+=1 ws=len/m now=len%m while (now+m<=309)and(ws>0): ws-=1;now+=m MAX=1 for i in range(ws*m):MAX*=10 f=float(n/MAX)**(1.0/m) x=int(f) for i in range(ws):x*=10 x0=-1; while x0!=x: x0=x x=x+x*(n-x**m)/(n*m) while (x+1)**m<=n:x=x+1 print x
阅读全文
0 0
- BZOJ 1213 [HNOI2004]高精度开根 牛顿法
- [牛顿法 || 二分] BZOJ 1213 [HNOI2004]高精度开根 & 51Nod 1166 大数开平方
- BZOJ 1213: [HNOI2004]高精度开根
- 【BZOJ】1213 [HNOI2004]高精度开根
- BZOJ 1213 HNOI2004 高精度开根
- BZOJ 1213 HNOI2004 高精度开根 二分+高(Py)精(thon)度
- 【HNOI2004】【BZOJ1213】高精度开根
- bzoj1213 [HNOI2004]高精度开根
- 二分+高精度——BZOJ1213/Luogu2293 [HNOI2004]高精度开根
- [BZOJ]1210: [HNOI2004]邮递员 插头DP+高精度
- 【BZOJ】1213 高精度开根
- BZOJ 1213, 高精度开根
- Java 高精度开根!(BZOJ 1213)
- 【BZOJ】1211 [HNOI2004]树的计数 prufer编码+组合数学+高精度
- BZOJ 1210: [HNOI2004]邮递员
- BZOJ 1213 Python... 二分+高精度
- [置顶] Project Euler Problem 80-高精度开方-牛顿逼近法
- bzoj 1207 [HNOI2004]打鼹鼠
- AngularJS路由的简单配置
- 楠叔的Oracle DBA进阶之路-起航
- vue2.0子路由配置和跳转 vue-router: 嵌套路由
- nginx+ffmpeg搭建rtmp转播rtsp流的flash服务器
- jdbc笔记
- BZOJ 1213 [HNOI2004]高精度开根 牛顿法
- Spring的面向切面编程在redis中应用
- rasterio
- 第五章 开始Unity Shader学习之旅
- 两个信号量实现生产者消费者模型
- 【Leetcode-Easy-53】Maximum Subarray
- jQuery -- Swiper插件使用ajax动态加载数据出现无法滑动的问题
- final、finally、finalize区别
- 嗯哇和您一起学python