Выяснить, можно ли представить данное число в виде суммы трех квадратов натуральных чисел

#include <iostream>
#include <cmath>
using namespace std;
 
const int MAX_ARRAY = 100;
int array [MAX_ARRAY];
 
void init_array ()
{
    for (int i = 1; i < MAX_ARRAY; i ++)
         array [i-1] = (i*i);
}
 
bool func (int n)
{
    for (int i1 = 0; i1 < MAX_ARRAY-1; i1 ++)
         for (int i2 = 0; i2 < MAX_ARRAY-1; i2 ++)
              for (int i3 = 0; i3 < MAX_ARRAY-1; i3 ++)
                   if (array [i1] + array [i2] + array [i3] == n)
                   {
                        cout << sqrt (array [i1]) << "^2 + " << sqrt (array [i2]) << "^2 + " << sqrt (array [i3]) << "^2 = " << n;
                        return 1;
                   }
    return 0;
}
 
int main (int argc, char *argv [])
{
    init_array ();
    // 9 + 16 + 25 = 50
    // 4 + 100 + 600 = 704
    int n = 704;
    if (!func (n))
         cout << "Error!";
    getchar ();
    return 0;
}

const int N=10000;
int a[N];
 
bool f(int n, int k, int j) {  
    if (k==0) return n==0; else if (n<a[j]) return false;
    for (int i=j; i<N; i++) {
        if (f(n-a[i], k-1, i)) {
            cout << (k==1 ? "" : " + ") << i << "^2";
            return true;
        }
    }
    return false;
}
 
int main() {
    for (int i=0; i<N; i++) a[i]=i*i;
    int n, k; cin >> n >> k;    
    if (f(n, k, 1)) cout << " = " << n; else cout << "no solution";
}

 int count = 3; //число слагаемых x^2
    int limit = 100; //предел вычисления x^2+x^2+x^2=30000
    int value = 12; //заданное число для поиска
 
    int *length = new int[count] {};
    for (int k = 0; k < std::ceil(std::pow(limit, count)); k++) {
        for (int i = 0; i < count; i++) {
            length[i]++;
            if (length[i] > limit || length[i] == 1) length[i] = 1;
            else break;
        }
        if (value == std::inner_product(length, length + count, length, 0, std::plus<int>(), std::multiplies<int>())) {
            for (int c = 0; c < count; c++)cout << length[c] << ((c + 1 == count) ? "^2" : "^2+");
            return 0;
        }
    }
    delete[] length;
    cout << "no solution";