分数化循环小数(纯/混)

网上看到这个问题,查了一下除了字符串查找没看到别的解法,自己写了一个,要考虑一些特殊情况,比如除0,负数等等,时间复杂度n^2,小数点长度N大于一万就慢了,感觉哪里还能再优化一下,最多跑过11111111/59595961,一共591万位,release要跑20秒5G内存,再大内存就不够了。

写的乱七八糟的,不写点注释以后估计自己都看不懂。

#include <iostream>
#include <vector>
using namespace std;
int main() {
    int upperLimit = 10;//小数点位数上限,当达到上限仍未找到循环节时乘以三继续寻找。
    int loopLen = 1;//循环节长度
    __int64 num;
    __int64 quotient;//商的整数部分
    __int64 a, b;
    vector<int> rem;//remainder小数部分
    cin >> a >> b;
    if (b == 0) {
        cout << "Integer division by zero。"; 
        ::system("pause"); 
        return 0;
    }
    if (a == 0) {
        cout << "0"; 
        ::system("pause"); 
        return 0;
    }
    if (float(a) / float(b) < 0) cout << "-";//负数判断
    a = abs(a);
    b = abs(b);
    quotient = a / b;
    a -= quotient * b;
    while (upperLimit *= 3,1) {
        auto HPC = [&]() {//高精度除法,除至第upperLimit位
            while (rem.size() < upperLimit) {
            a *= 10;
            num = a / b;
            rem.push_back(num);
            a -= b * num;
        }};
        auto Validate = [&]() {//结果有可能为0.(11...112)类似情况,至少验证多一倍的长度加上100位
        //以确保不会得出0.(1)这种答案,当然存在有一百多位的0.(11...112)循环,但估计分子分母是写不下了
            upperLimit *= 2 + 100;
            HPC();
            int front = rem.size() - 2 * loopLen;
            int back = rem.size() - loopLen;
            for (int i = upperLimit / 2; i < upperLimit; i++) {
                if (rem[i] != rem[front + (i - front) % loopLen]) return false;
            }
            return true;
        };
        auto SearchRepetation = [&]() {
            HPC();
            for (; loopLen <= rem.size() / 3; loopLen++) {
                auto HasRepeat = [&]() {
                    for (int index = 1; index <= loopLen; index++)
                        if (rem[rem.size() - index] != rem[rem.size() - loopLen - index] 
                        || rem[rem.size()-index-loopLen]!= rem[rem.size() - index - 2*loopLen]) 
                        return false;
                        //查找时确保循环节出现至少三次,其实两次就够了,更多次数也不能保证特殊情况,
                        //主要靠上面的验证步骤。上面的Validate多出来的一倍和100位也可以合并到这里面。
                    return true;
                };
                if (HasRepeat() && Validate()) {//确定结果,输出
                    int front = rem.size() - 2 * loopLen;
                    int back = rem.size() - loopLen;
                    while (front > 0 && rem[front-1] == rem[back-1]) { front--; back--; }
                    if (back - front == 1 && rem[front] == 0) {//无循环
                        if (front != 0)
                            cout << quotient << ".";
                        else//整数情况
                            cout << quotient;
                        for (int i = 0; i < front; i++) {
                            cout << rem[i];
                        }
                    }
                    else {//有循环
                        cout << quotient << ".";
                        for (int i = 0; i < back; i++) {
                            if (i == front)cout << "(";
                            cout << rem[i];
                        }
                        cout << ")";
                    }
                    return true;
                }
            }
            return false;
        };
        if (SearchRepetation()) break;
    }
    ::system("pause");
}
  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值