Ну, как минимум O(n^4) - потому что пятая переменная вычисляется из предыдущих четырех. А 131^4 - это меньше 300 миллионов, уже практически реально :)
Далее, запишем как K..+J..+H.. = S-R..-U.. - ну, вы поняли :)
Левая часть может иметь 131^3 = 2 с небольшим миллиона значений. Складываем в тот же хеш. Правая часть - и того меньше, 17 тысяч. После чего находим пересечения этих таблиц - вот они и есть решениями...
Update
Вот решение на С++
#include <iostream>
#include <iomanip>
#include <map>
using namespace std;
const long long K = 200, J = -300, H = 749, R = 344, U = -265, S = 1000;
long long intersection(const map<long long,int>& l,
const map<long long,int>& r)
{
long long sum = 0;
auto first1 = l.begin(), last1 = l.end();
auto first2 = r.begin(), last2 = r.end();
while (first1 != last1 && first2 != last2)
{
if (first1->first == first2->first)
{
sum += first1->second*first2->second;
++first1; ++first2;
} else if (first1->first < first2->first) ++first1;
else ++first2;
}
return sum;
}
int main()
{
map<long long,int> l,r;
for(long long y1 = -65; y1 <= 65; ++y1)
for(long long y2 = -65; y2 <= 65; ++y2)
for(long long y3 = -65; y3 <= 65; ++y3)
l[K*y1*y1*y1 + J*y2*y2*y2 + H*y3*y3*y3]++;
for(long long y4 = -65; y4 <= 65; ++y4)
for(long long y5 = -65; y5 <= 65; ++y5)
r[S-R*y4*y4*y4*y4-U*y5*y5*y5*y5*y5]++;
cout << intersection(l,r) << endl;
}
У меня работает около 2 секунд. Накинем на то, что это Python, порядок (понятно, что я утрирую) - все равно до часа далеко...