CPP开方

简弘

已于 2022-04-08 00:21:00 修改

阅读量802

点赞数

文章标签： c++ 算法

于 2022-04-07 19:32:59 首次发布

本文链接：https://blog.youkuaiyun.com/AA209935/article/details/124024461

版权

C++ 算法平方根迭代法高精度计算

关键词由优快云通过智能技术生成

该博客主要展示了两个用于计算平方根的C++函数实现，`mySqrt` 和 `mySqrt1`。这两个函数采用迭代法逐步逼近目标值，精确到指定的小数位数。通过输入一个数字，程序将输出其平方根的高精度表示。

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


#include "mySqrt.hpp"
#include <iostream>
#include <string>
#include <assert.h>

using namespace std;

int *mySqrt(double, int max = 150);
double mySqrt1(double num, int max);

int main()
{

    int num;
    cout << "输入一个数：";
    cin >> num;
    int i = 1;
    int *a = mySqrt(num);
    cout << "√(" << num << ") = " << a[0];
    if (a[i] != -1)
    {
        cout << ".";
    }
    while (a[i] != -1)
    {
        cout << a[i];
        i++;
    }
    cout << endl;

    // system("pause");

    return 0;
    // cout << mySqrt(25) << endl;
}

/*
 num：需要开方的数字
 max：精确到多少位
 */
int *mySqrt(double num, int max)
{
    int *a = new int[max];
    /*for (int w = 0; w < max; w++) {
        a[w] = 0.0;
    }*/
    assert(num > 0);
    int i = 1;
    int d = 1;
    double divisor;  // 除数
                     // int divisor_Next; // 传给下一个后
    double dividend; // 被除数
    double down;     // 下边的数
    double trading;  // 商
                     // double finalTrad;

    while (i * i < num)
    {
        trading = i;
        i++;
    }

    if (i * i == num)
    {
        a[0] = i;
        a[1] = -1;
        return a;
    }
    else
    {
        a[0] = (i - 1);
    }

    // cout << trading << endl;
    dividend = num - (i - 1) * (i - 1);
    // cout << dividend << endl;
    divisor = trading * 20;
    // cout << divisor << endl; //40
    for (int j = 0; j < max; j++)
    {
        int k = 0;
        for (k = 0; k <= 9; k++)
        {
            if ((divisor + k) * k >= dividend * 100)
            {
                //cout << (k - 1) << " ";

                break;
            }
        }
        a[d] = k - 1;
        //cout << a[d] << " ";
        // cout << "k" << k << endl;
        dividend *= 100; // 除的数字 300 上边的数字
        down = (divisor + (k - 1)) * (k - 1);
        dividend = dividend - down;
        // cout << "下边" << dividend << endl;
        // cout << "下" << xia << endl;
        // cout << dividend << endl;
        trading = trading * 10 + (k - 1); // 最上边的最后的得数
        // cout << trading << endl;
        // cout << j << endl;
        divisor = trading * 20; // 左边的数
        d++;
    }
    // cout << d << endl;
    // while ((d - 1) > 0)
    // {
    //     trading /= 10.0;
    //     d--;
    // }
    // return trading;
    a[d - 1] = -1;

    return a;

    // cout << "   " << divisor;
}

double mySqrt1(double num, int max)
{
    // assert(num > 0);
    int i = 1;
    int d = 1;
    double divisor; // 除数
    // int divisor_Next; // 传给下一个后
    double dividend; // 被除数
    double down;     // 下边的数
    double trading;  // 商
    // double finalTrad;

    while (i * i < num)
    {
        trading = i;
        i++;
        if (i * i == num)
        {
            return i;
        }
    }
    // cout << trading << endl;
    dividend = num - (i - 1) * (i - 1);
    // cout << dividend << endl;
    divisor = trading * 20;
    // cout << divisor << endl; //40
    for (int j = 0; j < max; j++)
    {
        int k = 0;
        for (k = 0; k <= 9; k++)
        {
            if ((divisor + k) * k >= dividend * 100)
            {
                break;
            }
        }
        // cout << "k" << k << endl;
        dividend *= 100; // 除的数字 300 上边的数字
        down = (divisor + (k - 1)) * (k - 1);
        dividend = dividend - down;
        // cout << "下边" << dividend << endl;
        // cout << "下" << xia << endl;
        // cout << dividend << endl;
        trading = trading * 10 + (k - 1); // 最上边的最后的得数
        // cout << trading << endl;
        // cout << j << endl;
        divisor = trading * 20; // 左边的数
        d++;
    }

    while ((d - 1) > 0)
    {
        trading /= 10.0;
        d--;
    }
    return trading;

    // cout << "   " << divisor;
}