Дано натуральное число N. Определить, если это возможно, пару x,y таких натуральных чисел, что N = x^2+y^2

bool coupleXY(int num, int *x, int *y)
{ int xx, yy;
    for(xx=1 ; xx*xx < num; xx++) {
        for( yy=1; x*x+y*y <=num; yy++) {
          if  (xx*xx + yy*yy == num )    {
            *x = xx:
            *y = yy;
            return true;
          }
       }
    }
    return false;  
}

#include <iostream>
int main()
{
    int N;
    std::cout << "N = "; std::cin >> N;
    for (int x = 1; x * x <= N; x++)
        for (int y = 1; y * y <= N; y++)
            if (x * x + y * y == N)
                std::cout << "x = " << x << "; y = " << y << std::endl;
    return 0;
}

#include <iostream>
#include <iterator>
#include <map>
std::map<size_t, size_t> coupleXY(size_t N)
{
    std::map<size_t, size_t> res;
    for (size_t x = 1; x * x <= N; x++)
        for (size_t y = 1; y * y <= N; y++)
            if (x * x + y * y == N) 
                res.insert(std::pair<size_t, size_t>(x, y));
    return res;
}
int main()
{
    size_t N;
    std::cout << "N = "; std::cin >> N;
    std::map<size_t, size_t> res = coupleXY(N);
    for (std::map<size_t, size_t>::iterator itr= res.begin(); itr != res.end(); itr++)
    std::cout << "x = " << itr->first << "; y = " << itr->second << std::endl;
    return 0;
}

typedef struct (
  imt x, y;
} XY;
XY coupleXY(int num)
{ XY xy;
    for(xy.x=1 ; xy.x*xy.x < num; xy.x++) {
        for( xy.y=xy.x; xy.x*xy.x+xy.y*xy.y <=num; xy.y++) {
          if  (xy.x*xy.x + xy.y*xy.y == num )    {
            return xy;
          }
       }
    }
    xy.x = xy.y= 0; 
    return r;  
}
int main()
{
    int N, XY result;
    
    cout << "Set natural number N: ";
    cin >> N;
    resultY = coupleXY(N);
    cout<<"\nx^2+y^2 = N "<<endl;
    cout<<"x = "<<result.x<<endl;
    cout<<"y = "<<result.y<<endl;
    return 0;
}

XY coupleXY(int num)
{ XY r;
    for(r.x=1 ; r.x*r.x < num; r.x++) {
        for( r.y=r.x; r.x*r.x+xy.y*r.y <=num; r.y++) {
          if  (r.x*r.x + r.y*r.y == num )    {
            return r;
          }
       }
    }
    r.x = r.y= 0; 
    return r;  
}

#include <iostream>
using namespace std;
 
int main(){
    int n, m = 1;
 
    cin >> n;
    while (m * m + (m*m - 1) / 2 * (m*m - 1) / 2 <= n) {
        if (m * m + (m*m - 1) / 2 * (m*m - 1) / 2 == n) {
            cout << m << "  " << (m*m - 1) / 2;
            return 0;
        }
        m++;
    }
    cout << "no";
}

        for(x = 1; x*x < N; x++)
            for(y = 1; y*y < N; y++)

#include <iostream>
#include <cmath>
 
uint8_t decompose(const uint64_t n, uint32_t& a, uint32_t& b)
{
   for (uint32_t i=1;i<sqrt(n);i++)
    {
     double _t(sqrt(n-i*i));
     if (_t==int(_t)) {
        a = i;
        b = int(_t);
        return 1;
     }
    }
    return 0;
}
int main()
{
    uint32_t a,b;
    for (uint64_t n=1; n<100; n++)
    if (decompose(n,a,b))
        std::cout<<n<<" = "<<a<<"^2 + "<<b<<"^2"<<std::endl;
 
    return 0;
}

{b = true; goto met_10;}