MyException - 我的异常网
当前位置:我的异常网» Java相关 » 求立方根算法-个人反抗方根算法的穷举和优化

求立方根算法-个人反抗方根算法的穷举和优化

www.MyException.Cn  网友分享于:2013-08-11  浏览:0次
求立方根算法--个人对立方根算法的穷举和优化

在hpe实训中心学习,遇到了求立方根的题目,在此做一下算法笔记,

分析过程:

数n的立方根就是n=i*i**i;所以我们会优先想到一下方法.

    static double g32(double n){ //简易版
        double i = 0, k = 0.0005f;
        if (n < 0) {    //输入负数判断
            k /= -1;
        }
        do{
            i+=k;
        }while(abs(i*i*i)<abs(n)); //abs为自己写的求绝对值方法
        return i;
    }

可以看出此方法的求解精度为0.001;且当输入数据过大时效率堪忧,所以就有了以下优化

    static double g33(double n){    //优化2
        double i = 0, k = 5f;
        if (n < 0) {
            k /= -1;
        }for (int t = 0; t < 15; t++) {//精度到小数点后15位
            do {
                i += k;  //开始时每次加5快速逼近正确值
            } while (abs(i * i * i) < abs(n));  //当i^3>=n时退出
            i-=k;k /= 9.2;  //i还原到退出前的值,k缩小,进入下一次逼近
        } 
        return i;
    }

此方法可以快速求得立方根,输入数值n不太大时使用,当n太大在逼近过程中i^3与(i+k)^3差距太大,循环次数剧增,进入死循环状态.在我电脑上当n=9999999999 (10个9)时就会进入死循环

所以想到一种解决方法,设置一循环次数计数器w,使k值随循环次数增加而增加,在一定程度上解决死循环问题.下面是代码

    static double g34(double n){    //最终优化
        double i = 0, k = 5f;
        if (n < 0) {
            k /= -1;
        }
        int w=0;
        for (int t = 0; t < 15; t++) {
            do {
                i += k;
                k=k+w*k/50000;
                w++;  //循环计数器
            } while (abs(i * i * i) < abs(n));
            i-=k;k /= 9.5;
        }
        System.out.println(w);
        return i;
    }

更改后n等于20个9也不会死循环.

最后附上所有程序代码.

package com.gfuzan.test;

import java.util.Scanner;

public class Test {

    /**
     * 开发: GFuZan
     * 时间: 2017.08.08
     * 功能: 立方根
     */
    public static void main(String[] args) {
        System.out.print("请输入一个数: ");
        Scanner sc = new Scanner(System.in);
        double n = sc.nextDouble();
        sc.close();
        System.out.println("Math:    \t"+ Math.pow(n, 1.0/3));
        System.out.println("My简易版:  \t"+g32(n)); 
        System.out.println("My优化2:  \t"+g33(n)); 
        System.out.println("My最终优化:\t"+g34(n)); 
        
    }
    
    static double g32(double n){ //简易版
        double i = 0, k = 0.0005f;
        if (n < 0) {
            k /= -1;
        }
        do{
            i+=k;
        }while(abs(i*i*i)<abs(n));
        return i;
    }
    static double g34(double n){    //最终优化
        double i = 0, k = 5f;
        if (n < 0) {
            k /= -1;
        }
        int w=0;
        for (int t = 0; t < 15; t++) {
            do {
                i += k;
                k=k+w*k/50000;
                w++;
            } while (abs(i * i * i) < abs(n));
            i-=k;k /= 9.5;
        }
        return i;
    }
    static double g33(double n){    //优化2
        double i = 0, k = 5f;
        if (n < 0) {
            k /= -1;
        }
        for (int t = 0; t < 15; t++) {
            do {
                i += k;
            } while (abs(i * i * i) < abs(n));
            i-=k;k /= 9.2;
        }
        return i;
    }
    static double abs(double f) {
        if (f < 0) {
            return 0 - f;
        }
        return f;
    }
}

运行结果

 

文章评论

程序员眼里IE浏览器是什么样的
程序员眼里IE浏览器是什么样的
做程序猿的老婆应该注意的一些事情
做程序猿的老婆应该注意的一些事情
当下全球最炙手可热的八位少年创业者
当下全球最炙手可热的八位少年创业者
漫画:程序员的工作
漫画:程序员的工作
十大编程算法助程序员走上高手之路
十大编程算法助程序员走上高手之路
为什么程序员都是夜猫子
为什么程序员都是夜猫子
聊聊HTTPS和SSL/TLS协议
聊聊HTTPS和SSL/TLS协议
我的丈夫是个程序员
我的丈夫是个程序员
程序猿的崛起——Growth Hacker
程序猿的崛起——Growth Hacker
Web开发者需具备的8个好习惯
Web开发者需具备的8个好习惯
为啥Android手机总会越用越慢?
为啥Android手机总会越用越慢?
科技史上最臭名昭著的13大罪犯
科技史上最臭名昭著的13大罪犯
“懒”出效率是程序员的美德
“懒”出效率是程序员的美德
编程语言是女人
编程语言是女人
如何区分一个程序员是“老手“还是“新手“?
如何区分一个程序员是“老手“还是“新手“?
2013年中国软件开发者薪资调查报告
2013年中国软件开发者薪资调查报告
程序员都该阅读的书
程序员都该阅读的书
总结2014中国互联网十大段子
总结2014中国互联网十大段子
不懂技术不要对懂技术的人说这很容易实现
不懂技术不要对懂技术的人说这很容易实现
程序员必看的十大电影
程序员必看的十大电影
我是如何打败拖延症的
我是如何打败拖延症的
2013年美国开发者薪资调查报告
2013年美国开发者薪资调查报告
 程序员的样子
程序员的样子
每天工作4小时的程序员
每天工作4小时的程序员
程序员应该关注的一些事儿
程序员应该关注的一些事儿
Java程序员必看电影
Java程序员必看电影
10个调试和排错的小建议
10个调试和排错的小建议
程序员的一天:一寸光阴一寸金
程序员的一天:一寸光阴一寸金
5款最佳正则表达式编辑调试器
5款最佳正则表达式编辑调试器
程序员周末都喜欢做什么?
程序员周末都喜欢做什么?
程序员和编码员之间的区别
程序员和编码员之间的区别
60个开发者不容错过的免费资源库
60个开发者不容错过的免费资源库
那些性感的让人尖叫的程序员
那些性感的让人尖叫的程序员
一个程序员的时间管理
一个程序员的时间管理
写给自己也写给你 自己到底该何去何从
写给自己也写给你 自己到底该何去何从
团队中“技术大拿”并非越多越好
团队中“技术大拿”并非越多越好
代码女神横空出世
代码女神横空出世
旅行,写作,编程
旅行,写作,编程
老程序员的下场
老程序员的下场
程序员最害怕的5件事 你中招了吗?
程序员最害怕的5件事 你中招了吗?
10个帮程序员减压放松的网站
10个帮程序员减压放松的网站
如何成为一名黑客
如何成为一名黑客
Google伦敦新总部 犹如星级庄园
Google伦敦新总部 犹如星级庄园
中美印日四国程序员比较
中美印日四国程序员比较
Web开发人员为什么越来越懒了?
Web开发人员为什么越来越懒了?
那些争议最大的编程观点
那些争议最大的编程观点
Java 与 .NET 的平台发展之争
Java 与 .NET 的平台发展之争
什么才是优秀的用户界面设计
什么才是优秀的用户界面设计
程序员的鄙视链
程序员的鄙视链
软件开发程序错误异常ExceptionCopyright © 2009-2015 MyException 版权所有