总时间限制: 1000ms 内存限制: 65536kB

描述

形如a= b + c + d的等式被称为完美立方等式。例如12= 6 + 8 + 10 。编写一个程序,对任给的正整数N (N ≤ 100),寻找所有的四元组(a, b, c, d),使得a = b + c + d,其中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值小的先输出。

样例输入

  1. 24

样例输出

  1. Cube = 6, Triple = (3,4,5)
  2. Cube = 12, Triple = (6,8,10)
  3. Cube = 18, Triple = (2,12,16)
  4. Cube = 18, Triple = (9,12,15)
  5. Cube = 19, Triple = (3,10,18)
  6. Cube = 20, Triple = (7,14,17)
  7. Cube = 24, Triple = (12,16,20)

思路

  1. 题目要求 a 的值从小到大输出, 所以 a 从小到大枚举.
  2. 这是不可避免的四重循环, 但不是每一个元素都在答案空间. 通过思考我们可以得到各个变量的枚举空间:
    • POJ 2810 完美立方 - 图1
    • POJ 2810 完美立方 - 图2
    • POJ 2810 完美立方 - 图3

代码

  1. #include <iostream>
  2. using namespace std;
  3. int main() {
  4. int N;
  5. cin >> N;
  6. for(int a = 2; a <= N; a++) {
  7. for(int b = 2; b <= a - 1; b++) {
  8. for(int c = b; c <= a - 1; c++) {
  9. for(int d = c; d <= a - 1; d++) {
  10. if(a*a*a == b*b*b + c*c*c + d*d*d)
  11. printf("Cube = %d, Triple = (%d,%d,%d)\n", a, b, c, d);
  12. }
  13. }
  14. }
  15. }
  16. return 0;
  17. }