function f(c,x)
b=(1:c).^x;d=bsxfun(@plus,b,b');d(d>c^x)=0;[~,r]=max(d(:));sprintf('%i^%i + %i^%i < %i^%i',[mod(r-1,c)+1 ceil(r/c) c;x x x])
end % needed in online interpreter only
f(8,3)
ZnVuY3Rpb24gZihjLHgpCmI9KDE6YykuXng7ZD1ic3hmdW4oQHBsdXMsYixiJyk7ZChkPmNeeCk9MDtbfixyXT1tYXgoZCg6KSk7c3ByaW50ZignJWleJWkgKyAlaV4laSA8ICVpXiVpJyxbbW9kKHItMSxjKSsxIGNlaWwoci9jKSBjO3ggeCB4XSkKCmVuZCAlIG5lZWRlZCBpbiBvbmxpbmUgaW50ZXJwcmV0ZXIgb25seQoKZig4LDMp