OpenJ_Bailian - 2810 完美立方 【枚举】

题目连接：https://cn.vjudge.net/problem/OpenJ_Bailian-2810

题目描述：

形如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。

Input

一个正整数N (N≤100)。

Output

每行输出一个完美立方。输出格式为： 
Cube = a, Triple = (b,c,d) 
其中a,b,c,d所在位置分别用实际求出四元组值代入。 

请按照a的值，从小到大依次输出。当两个完美立方等式中a的值相同，则b值小的优先输出、仍相同则c值小的优先输出、再相同则d值小的先输出。

Sample Input

24

Sample Output

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)

AC代码： 

#include <iostream>
#include <cstdio>
using namespace std;
int FastPow(int a,int n) {      //快速幂
    int ans=1;
    while(n != 0) {
        if(n%2)
            ans *= a;
        a *= a;
        n /=2;
    }
    return ans;
}
int main()
{
    int N,a,b,c,d;
    cin >> N;
    for(int a=2;a<=N;a++) {
        for(int b=2;b<a;b++) {
            for(int c=b;c<a;c++) {
                for(int d=c;d<a;d++) {
                    if(FastPow(a,3) == FastPow(b,3)+FastPow(c,3)+FastPow(d,3)) {
                        printf("Cube = %d, Triple = (%d,%d,%d)\n",a,b,c,d);
                    }
                }
            }
        }
    }
}