【C++ 原创】拉马努金（的士数）_python实现求“的士数”如1729-CSDN博客

本文链接：https://blog.csdn.net/Z02012341234/article/details/139569883

【题目描述】
拉马努金是印度历史上最著名的数学家之一。他没受过正规的高等数学教育，沉迷数论，尤爱牵涉 $eq?%5Cpi$ 、质数等数学常数的求和公式，以及整数分拆。惯以直觉（或者是跳步）导出公式，不喜作证明（事后往往证明他是对的）。他留下的那些没有证明的公式，引发了后来的大量研究。
拉马努金病重，哈代前往探望。哈代说：“我乘出租车来，车牌号码是1729，这数真没趣，希望不是不祥之兆。”拉马努金答道：“不，那是个有趣得很的数。可以用两个立方之和来表达而且有两种表达方式的数之中，1729是最小的。”（即1729 = $eq?1%5E%7B3%7D$ + $eq?12%5E%7B3%7D$ = $eq?9%5E%7B3%7D$ + $eq?10%5E%7B3%7D$ ，后来这类数称为的士数；事实上1729是最小的可以表为两种正立方之和的正整数。

请编写程序，输入一个正整数n，判断其能否用两个立方之和来表达而且有两种表达方式。如果能，请输出分解形式。如果不能，输出“不能”。（注：能以3个不同的方法表示成两个正立方数之和的数最小是87539319）

【输入格式】输入一个正整数n（1000 $eq?%5Cleqslant$ n $eq?%5Cleqslant$ 100000）。

【输出格式】如果能用两个立方之和来表达而且有两种表达方式，输出两种分解方式，每种分解方式中间用空格分开。如果不能，输出“不能”。（输出的第二种分解方式的第一个数大于第一种的第一个数，且每种方式的两个数后一个大于前一个。）

输入样例1:
1729
输出样例1:
1 12
9 10
输入样例2:
1008
输出样例2:
不能

【参考代码】

#include<iostream>
#include<cmath>
using namespace std;
bool isShishu(int i)
{
    int temp = 0,c = 0, a5, a6, a7; 
    for (int j = 1; j <= (pow(i, 1.0 / 3)); j++)//开三次根号节省时间
        for (int t = j + 1; t <= (pow(i, 1.0 / 3)); t++)
        {
            int a = floor(pow(j, 3));
            int b = floor(pow(t, 3));
            if ((a + b) == i)
            {
                if (c == 0)//一次性存储，发现第一个能分解成两个正立方和的数
                {
                    a5 = i, a6 = j, a7 = t;
                    c++;
                    break;//终止
                }
                if (i != a5)//又发现一个能分解的数，比对和上次发现的是否一样
                {
                    a5 = i, a6 = j, a7 = t;//不一样把它更新
                    break;
                }
                else//一样的话说明是另一种分解方式
                {
                    cout <<a6<< " " << a7<< endl;
                    cout <<j << " " << t;
                    return true;
                }
            }
            else if ((a + b) > i)
                break;
            else;
        }
    return false;
}
int main()
{
    int m;
    cin >> m;
    if (isShishu(m))
        ;
    else
        cout << "不能";
}

【运行示例】

【最初版本】

最开始写的是找一个区间中能分解成两个正立方和且有两种分解方式的数

【代码】

#include<iostream>
#include<cmath>
using namespace std;
int main()
{
    int temp = 0, c = 0, a5, a6, a7;
    for (int i = 1000; i <= 10000000; i++)
        for (int j = 1; j <= floor(pow(i, 1.0 / 3)); j++)
            for (int t = j + 1; t <= floor(pow(i, 1.0 / 3)); t++)
            {
                double a = floor(pow(j, 3));
                double b = floor(pow(t, 3));
                if ((a + b) == i)
                {
                    if (c == 0)
                    {
                        a5 = i, a6 = j, a7 = t;
                        c++;
                        break;
                    }
                    if (i != a5)
                    {
                        a5 = i, a6 = j, a7 = t;
                        break;
                    }
                    else
                    {
                        cout << a5 << " " << a6 << " " << a7 << endl;
                        cout << i << " " << j << " " << t << endl<<endl;
                        a5 = i, a6 = j, a7 = t;
                    }
                }
                else if ((a + b) > i)
                    break;
                else;
            }
}

【运行示例】