问题描述:
a^3 = b^3 + c^3 + d^3为完美立方等式。
编写程序对任意给的正整数N,寻找所有的四元组(a,b,c,d)使得a^3 = b^3 + c^3 + d^3,其中1
#include <iostream>
using namespace std;
void f(int n)
{
long int temp[101] = {0};
int a, b, c, d, i = 0;
while(i++ <= n)
temp[i] = i * i * i;
for (a = 6; a < n + 1; a++){
for (b = 2; b < a - 2; b++){
for (c = b + 1; c < a - 1; c++){
if (temp[a] < temp[b] + temp[c] + temp[c + 1])
break;
for (d = c + 1; d < a; d++){
if (temp[a] == temp[b] + temp[c] + temp[d])
cout << a <<"\t" << b << " " << c << " " << d << endl;
}}}}}
int main()
{
int n;
cin >> n;
f(n);
return 0;
}
分析:
- 可从数学角度分析,和的最小可能是多少?1,2,3,加下来36,不为某数的立方;2,3,4加下来99,不为立方,所以和的最小可能是5,但是5的立方为125,大于99,所以a应该从6开始;
- b的话,有数学分析可知,为1是不成立,所以最小从2开始;
- 因为是从小到大排列输出,所以d > c, c > b,
所以有 d = c + 1, d < a;
c = b + 1,c < a - 1;
b = 2, b < a - 2; - 为避免穷举过程中不断计算立方值,所以在一开始直接计算好存入数组,这样每个数字只需计算一次,而不是在循环中和循环次数一样多;