#include <iostream>
#include <cmath>
using namespace std;
int main() {
int n = 100;
cout << "n: " << n << endl;
int h=n>>1;
for(int p=0,pp=0;pp<=h;++p,pp=p*p)
{
int mx=(int)sqrt(n-pp+1);
for(int k=mx-1;k<=mx;++k) if(pp+k*k==n) cout<<n<<" = "<<p<<"^2 + "<<k<<"^2"<<endl;
}
return 0;
}
I2luY2x1ZGUgPGlvc3RyZWFtPgojaW5jbHVkZSA8Y21hdGg+Cgp1c2luZyBuYW1lc3BhY2Ugc3RkOwoKaW50IG1haW4oKSB7CglpbnQgbiA9IDEwMDsKCiAgICBjb3V0IDw8ICJuOiAiIDw8IG4gPDwgZW5kbDsKICAKICAgaW50IGg9bj4+MTsKICAgZm9yKGludCBwPTAscHA9MDtwcDw9aDsrK3AscHA9cCpwKQogICAgIHsKICAgICAgaW50IG14PShpbnQpc3FydChuLXBwKzEpOwogICAgICBmb3IoaW50IGs9bXgtMTtrPD1teDsrK2spIGlmKHBwK2sqaz09bikgY291dDw8bjw8IiA9ICI8PHA8PCJeMiArICI8PGs8PCJeMiI8PGVuZGw7CiAgICAgfQogICAgIAoJcmV0dXJuIDA7Cn0=