完美立方(POJ 2810)
Note:==
1.为避免重复计算。用一个数组cube[101]存储事先计算好的立方值。
2.题目要求2<=b<=c<=d<a所以循环由外至内为b,c,d它们的取值上限都是a-1,
下限分别是2,b,c==
问题描述:
形如a3= b3 + c3 + d3的等式被称为完美立方等式。例如123= 63 + 83 + 103 。编写一个程序,对任给的正整数N (N≤100),寻找所有的四元组(a, b, c, d),使得a3 = b3 + c3 + d3,其中a,b,c,d 大于 1, 小于等于N,且b<=c<=d。
输入
一个正整数N (N≤100)。
输出
每行输出一个完美立方。输出格式为:
Cube = a, Triple = (b,c,d)
其中a,b,c,d所在位置分别用实际求出四元组值代入。
请按照a的值,从小到大依次输出。当两个完美立方等式中a的值相同,则b值小的优先输出、仍相同则c值小的优先输出、再相同则d值小的先输出。
样例输入
24
样例输出
Cube = 6, Triple = (3,4,5)
Cube = 12, Triple = (6,8,10)
Cube = 18, Triple = (2,12,16)
Cube = 18, Triple = (9,12,15)
Cube = 19, Triple = (3,10,18)
Cube = 20, Triple = (7,14,17)
Cube = 24, Triple = (12,16,20)
#include<cstdio>
#include<iostream>
using namespace std;
int main(){
int a, b, c, d, N, i;
int cube[101];
cin >> N;
for(i = 2; i <= N; i++)//预处理2~N的立方
cube[i] = i * i * i;
for(a = 6; a <= N; a++)//枚举a
for(b = 2; b < a; b++)//枚举b
for(c = b; c < a; c++)//枚举c
for(d = c; d < a; d++)//枚举d
if(cube[a] == cube[b] + cube[c] + cube[d])
printf("Cube = %d, Triple = (%d,%d,%d)\n",a,b,c,d);
return 0;
}