勾股数查询的现场计算版

输入一个数，然后现场计算包含输入数的勾股数，最大可以到4294967295（大概率有bug），数字越大，计算越慢。

源码：

#include<iostream>
using namespace std;
int main()
{
    while (true)
    {
        register unsigned long a, b, c;
        cout << ">>> ";
        cin >> a;
        if (a == 0) { break; }
        register unsigned long ta1[3] = { a / 100000000,(a % 100000000) / 10000,a % 10000 };
        register unsigned long ta2[6] = { ta1[0] * ta1[0],ta1[0] * ta1[1],ta1[0] * ta1[2],ta1[1] * ta1[1],ta1[1] * ta1[2],ta1[2] * ta1[2] };
        register unsigned long ta3[5] = { ta2[5],ta2[4] + ta2[4],ta2[2] + ta2[3] + ta2[2],ta2[1] + ta2[1],ta2[0] };
        register unsigned long ta4[6] = { (ta3[0]) % 10000,(ta3[1] + ta3[0] / 10000) % 10000,(ta3[2] + (ta3[1] + ta3[0] / 10000) / 10000) % 10000,(ta3[3] + (ta3[2] + (ta3[1] + ta3[0] / 10000) / 10000) / 10000) % 10000,
            (ta3[4] + (ta3[3] + (ta3[2] + (ta3[1] + ta3[0] / 10000) / 10000) / 10000) / 10000) % 10000,(ta3[4] + (ta3[3] + (ta3[2] + (ta3[1] + ta3[0] / 10000) / 10000) / 10000) / 10000) / 10000 };
        cout << "if b<c<a:" << endl;
        bool none = true;
        for (b = 1; b < a / 1.4; b++)
        {
            register unsigned long tb1[3] = { b / 100000000,(b % 100000000) / 10000,b % 10000 };
            register unsigned long tb2[6] = { tb1[0] * tb1[0],tb1[0] * tb1[1],tb1[0] * tb1[2],tb1[1] * tb1[1],tb1[1] * tb1[2],tb1[2] * tb1[2] };
            register unsigned long tb3[5] = { tb2[5],tb2[4] + tb2[4],tb2[2] + tb2[3] + tb2[2],tb2[1] + tb2[1],tb2[0] };
            register unsigned long tb4[6] = { (tb3[0]) % 10000,(tb3[1] + tb3[0] / 10000) % 10000,(tb3[2] + (tb3[1] + tb3[0] / 10000) / 10000) % 10000,(tb3[3] + (tb3[2] + (tb3[1] + tb3[0] / 10000) / 10000) / 10000) % 10000,
                (tb3[4] + (tb3[3] + (tb3[2] + (tb3[1] + tb3[0] / 10000) / 10000) / 10000) / 10000) % 10000,(tb3[4] + (tb3[3] + (tb3[2] + (tb3[1] + tb3[0] / 10000) / 10000) / 10000) / 10000) / 10000 };
            for (c = b + 1; c < a; c++)
            {
                register unsigned long tc1[3] = { c / 100000000,(c % 100000000) / 10000,c % 10000 };
                register unsigned long tc2[6] = { tc1[0] * tc1[0],tc1[0] * tc1[1],tc1[0] * tc1[2],tc1[1] * tc1[1],tc1[1] * tc1[2],tc1[2] * tc1[2] };
                register unsigned long tc3[5] = { tc2[5],tc2[4] + tc2[4],tc2[2] + tc2[3] + tc2[2],tc2[1] + tc2[1],tc2[0] };
                register unsigned long tc4[6] = { (tc3[0]) % 10000,(tc3[1] + tc3[0] / 10000) % 10000,(tc3[2] + (tc3[1] + tc3[0] / 10000) / 10000) % 10000,(tc3[3] + (tc3[2] + (tc3[1] + tc3[0] / 10000) / 10000) / 10000) % 10000,
                    (tc3[4] + (tc3[3] + (tc3[2] + (tc3[1] + tc3[0] / 10000) / 10000) / 10000) / 10000) % 10000,(tc3[4] + (tc3[3] + (tc3[2] + (tc3[1] + tc3[0] / 10000) / 10000) / 10000) / 10000) / 10000 };
                if (ta4[5] != tb4[5] + tc4[5]) { continue; }
                if (ta4[4] != tb4[4] + tc4[4]) { continue; }
                if (ta4[3] != tb4[3] + tc4[3]) { continue; }
                if (ta4[2] != tb4[2] + tc4[2]) { continue; }
                if (ta4[1] != tb4[1] + tc4[1]) { continue; }
                if (ta4[0] != tb4[0] + tc4[0]) { continue; }
                cout << "    " << b << " " << c << " " << a << " ";
                if (ta4[5] != 0) { cout << ta4[5]; }
                if (ta4[4] != 0) { cout << ta4[4]; }
                if (ta4[3] != 0) { cout << ta4[3]; }
                if (ta4[2] != 0) { cout << ta4[2]; }
                if (ta4[1] != 0) { cout << ta4[1]; }
                cout << ta4[0] << endl;
                none = false;
            }
        }
        if (none)
        {
            cout << "    " << "none" << endl;
        }
        cout << "if b<a<c:" << endl;
        none = true;
        for (b = 1; b < a; b++)
        {
            register unsigned long tb1[3] = { b / 100000000,(b % 100000000) / 10000,b % 10000 };
            register unsigned long tb2[6] = { tb1[0] * tb1[0],tb1[0] * tb1[1],tb1[0] * tb1[2],tb1[1] * tb1[1],tb1[1] * tb1[2],tb1[2] * tb1[2] };
            register unsigned long tb3[5] = { tb2[5],tb2[4] + tb2[4],tb2[2] + tb2[3] + tb2[2],tb2[1] + tb2[1],tb2[0] };
            register unsigned long tb4[6] = { (tb3[0]) % 10000,(tb3[1] + tb3[0] / 10000) % 10000,(tb3[2] + (tb3[1] + tb3[0] / 10000) / 10000) % 10000,(tb3[3] + (tb3[2] + (tb3[1] + tb3[0] / 10000) / 10000) / 10000) % 10000,
                (tb3[4] + (tb3[3] + (tb3[2] + (tb3[1] + tb3[0] / 10000) / 10000) / 10000) / 10000) % 10000,(tb3[4] + (tb3[3] + (tb3[2] + (tb3[1] + tb3[0] / 10000) / 10000) / 10000) / 10000) / 10000 };
            for (c = a + 1; c < a * 1.5; c++)
            {
                register unsigned long tc1[3] = { c / 100000000,(c % 100000000) / 10000,c % 10000 };
                register unsigned long tc2[6] = { tc1[0] * tc1[0],tc1[0] * tc1[1],tc1[0] * tc1[2],tc1[1] * tc1[1],tc1[1] * tc1[2],tc1[2] * tc1[2] };
                register unsigned long tc3[5] = { tc2[5],tc2[4] + tc2[4],tc2[2] + tc2[3] + tc2[2],tc2[1] + tc2[1],tc2[0] };
                register unsigned long tc4[6] = { (tc3[0]) % 10000,(tc3[1] + tc3[0] / 10000) % 10000,(tc3[2] + (tc3[1] + tc3[0] / 10000) / 10000) % 10000,(tc3[3] + (tc3[2] + (tc3[1] + tc3[0] / 10000) / 10000) / 10000) % 10000,
                    (tc3[4] + (tc3[3] + (tc3[2] + (tc3[1] + tc3[0] / 10000) / 10000) / 10000) / 10000) % 10000,(tc3[4] + (tc3[3] + (tc3[2] + (tc3[1] + tc3[0] / 10000) / 10000) / 10000) / 10000) / 10000 };
                if (tc4[5] != tb4[5] + ta4[5]) { continue; }
                if (tc4[4] != tb4[4] + ta4[4]) { continue; }
                if (tc4[3] != tb4[3] + ta4[3]) { continue; }
                if (tc4[2] != tb4[2] + ta4[2]) { continue; }
                if (tc4[1] != tb4[1] + ta4[1]) { continue; }
                if (tc4[0] != tb4[0] + ta4[0]) { continue; }
                cout << "    " << b << " " << a << " " << c << " ";
                if (tc4[5] != 0) { cout << tc4[5]; }
                if (tc4[4] != 0) { cout << tc4[4]; }
                if (tc4[3] != 0) { cout << tc4[3]; }
                if (tc4[2] != 0) { cout << tc4[2]; }
                if (tc4[1] != 0) { cout << tc4[1]; }
                cout << tc4[0] << endl;
                none = false;
            }
        }
        if (none)
        {
            cout << "    " << "none" << endl;
        }
        none = true;
        cout << "if a<b<c:" << endl;
        if (a < 92680)
        {
            register unsigned long d = (a / 2 + 1) * a;
            for (b = a + 1; b < d; b++)
            {
                register unsigned long tb1[3] = { b / 100000000,(b % 100000000) / 10000,b % 10000 };
                register unsigned long tb2[6] = { tb1[0] * tb1[0],tb1[0] * tb1[1],tb1[0] * tb1[2],tb1[1] * tb1[1],tb1[1] * tb1[2],tb1[2] * tb1[2] };
                register unsigned long tb3[5] = { tb2[5],tb2[4] + tb2[4],tb2[2] + tb2[3] + tb2[2],tb2[1] + tb2[1],tb2[0] };
                register unsigned long tb4[6] = { (tb3[0]) % 10000,(tb3[1] + tb3[0] / 10000) % 10000,(tb3[2] + (tb3[1] + tb3[0] / 10000) / 10000) % 10000,(tb3[3] + (tb3[2] + (tb3[1] + tb3[0] / 10000) / 10000) / 10000) % 10000,
                    (tb3[4] + (tb3[3] + (tb3[2] + (tb3[1] + tb3[0] / 10000) / 10000) / 10000) / 10000) % 10000,(tb3[4] + (tb3[3] + (tb3[2] + (tb3[1] + tb3[0] / 10000) / 10000) / 10000) / 10000) / 10000 };
                for (c = b + 1; c < b * 1.5; c++)
                {
                    register unsigned long tc1[3] = { c / 100000000,(c % 100000000) / 10000,c % 10000 };
                    register unsigned long tc2[6] = { tc1[0] * tc1[0],tc1[0] * tc1[1],tc1[0] * tc1[2],tc1[1] * tc1[1],tc1[1] * tc1[2],tc1[2] * tc1[2] };
                    register unsigned long tc3[5] = { tc2[5],tc2[4] + tc2[4],tc2[2] + tc2[3] + tc2[2],tc2[1] + tc2[1],tc2[0] };
                    register unsigned long tc4[6] = { (tc3[0]) % 10000,(tc3[1] + tc3[0] / 10000) % 10000,(tc3[2] + (tc3[1] + tc3[0] / 10000) / 10000) % 10000,(tc3[3] + (tc3[2] + (tc3[1] + tc3[0] / 10000) / 10000) / 10000) % 10000,
                        (tc3[4] + (tc3[3] + (tc3[2] + (tc3[1] + tc3[0] / 10000) / 10000) / 10000) / 10000) % 10000,(tc3[4] + (tc3[3] + (tc3[2] + (tc3[1] + tc3[0] / 10000) / 10000) / 10000) / 10000) / 10000 };
                    if (tc4[5] != tb4[5] + ta4[5]) { continue; }
                    if (tc4[4] != tb4[4] + ta4[4]) { continue; }
                    if (tc4[3] != tb4[3] + ta4[3]) { continue; }
                    if (tc4[2] != tb4[2] + ta4[2]) { continue; }
                    if (tc4[1] != tb4[1] + ta4[1]) { continue; }
                    if (tc4[0] != tb4[0] + ta4[0]) { continue; }
                    cout << "    " << a << " " << b << " " << c << " ";
                    if (tc4[5] != 0) { cout << tc4[5]; }
                    if (tc4[4] != 0) { cout << tc4[4]; }
                    if (tc4[3] != 0) { cout << tc4[3]; }
                    if (tc4[2] != 0) { cout << tc4[2]; }
                    if (tc4[1] != 0) { cout << tc4[1]; }
                    cout << tc4[0] << endl;
                    none = false;
                }
            }
        }
        else
        {
            for (b = a + 1; b < 2863311531; b++)
            {
                register unsigned long tb1[3] = { b / 100000000,(b % 100000000) / 10000,b % 10000 };
                register unsigned long tb2[6] = { tb1[0] * tb1[0],tb1[0] * tb1[1],tb1[0] * tb1[2],tb1[1] * tb1[1],tb1[1] * tb1[2],tb1[2] * tb1[2] };
                register unsigned long tb3[5] = { tb2[5],tb2[4] + tb2[4],tb2[2] + tb2[3] + tb2[2],tb2[1] + tb2[1],tb2[0] };
                register unsigned long tb4[6] = { (tb3[0]) % 10000,(tb3[1] + tb3[0] / 10000) % 10000,(tb3[2] + (tb3[1] + tb3[0] / 10000) / 10000) % 10000,(tb3[3] + (tb3[2] + (tb3[1] + tb3[0] / 10000) / 10000) / 10000) % 10000,
                    (tb3[4] + (tb3[3] + (tb3[2] + (tb3[1] + tb3[0] / 10000) / 10000) / 10000) / 10000) % 10000,(tb3[4] + (tb3[3] + (tb3[2] + (tb3[1] + tb3[0] / 10000) / 10000) / 10000) / 10000) / 10000 };
                for (c = b + 1; c < b * 1.5 && c < 4294967295; c++)
                {
                    register unsigned long tc1[3] = { c / 100000000,(c % 100000000) / 10000,c % 10000 };
                    register unsigned long tc2[6] = { tc1[0] * tc1[0],tc1[0] * tc1[1],tc1[0] * tc1[2],tc1[1] * tc1[1],tc1[1] * tc1[2],tc1[2] * tc1[2] };
                    register unsigned long tc3[5] = { tc2[5],tc2[4] + tc2[4],tc2[2] + tc2[3] + tc2[2],tc2[1] + tc2[1],tc2[0] };
                    register unsigned long tc4[6] = { (tc3[0]) % 10000,(tc3[1] + tc3[0] / 10000) % 10000,(tc3[2] + (tc3[1] + tc3[0] / 10000) / 10000) % 10000,(tc3[3] + (tc3[2] + (tc3[1] + tc3[0] / 10000) / 10000) / 10000) % 10000,
                        (tc3[4] + (tc3[3] + (tc3[2] + (tc3[1] + tc3[0] / 10000) / 10000) / 10000) / 10000) % 10000,(tc3[4] + (tc3[3] + (tc3[2] + (tc3[1] + tc3[0] / 10000) / 10000) / 10000) / 10000) / 10000 };
                    if (tc4[5] != tb4[5] + ta4[5]) { continue; }
                    if (tc4[4] != tb4[4] + ta4[4]) { continue; }
                    if (tc4[3] != tb4[3] + ta4[3]) { continue; }
                    if (tc4[2] != tb4[2] + ta4[2]) { continue; }
                    if (tc4[1] != tb4[1] + ta4[1]) { continue; }
                    if (tc4[0] != tb4[0] + ta4[0]) { continue; }
                    cout << "    " << a << " " << b << " " << c << " ";
                    if (tc4[5] != 0) { cout << tc4[5]; }
                    if (tc4[4] != 0) { cout << tc4[4]; }
                    if (tc4[3] != 0) { cout << tc4[3]; }
                    if (tc4[2] != 0) { cout << tc4[2]; }
                    if (tc4[1] != 0) { cout << tc4[1]; }
                    cout << tc4[0] << endl;
                    none = false;
                }
            }
        }
        if (none)
        {
            cout << "    " << "none" << endl;
        }
        cout << "over" << endl;
    }
}