#include <iostream>
#include <cmath>
 
using namespace std;
 
int main()
{
    double a = 1, b = 10000, c = 1;
 
    double delta = b*b - 4*a*c;
 
    cout << "Rownanie: x^2 + 10000x + 1 = 0" << endl;
    cout << "Delta = " << delta << endl;
 
    if (delta >= 0)
    {
        double sqrtDelta = sqrt(delta);
 
        // poprawny numerycznie algorytm
        double x1 = (-b - sqrtDelta) / (2*a);
        double x2 = c / (a * x1);
 
        cout << "x1 = " << x1 << endl;
        cout << "x2 = " << x2 << endl;
    }
 
    return 0;
}
 

I2luY2x1ZGUgPGlvc3RyZWFtPgojaW5jbHVkZSA8Y21hdGg+Cgp1c2luZyBuYW1lc3BhY2Ugc3RkOwoKaW50IG1haW4oKQp7CiAgICBkb3VibGUgYSA9IDEsIGIgPSAxMDAwMCwgYyA9IDE7CgogICAgZG91YmxlIGRlbHRhID0gYipiIC0gNCphKmM7CgogICAgY291dCA8PCAiUm93bmFuaWU6IHheMiArIDEwMDAweCArIDEgPSAwIiA8PCBlbmRsOwogICAgY291dCA8PCAiRGVsdGEgPSAiIDw8IGRlbHRhIDw8IGVuZGw7CgogICAgaWYgKGRlbHRhID49IDApCiAgICB7CiAgICAgICAgZG91YmxlIHNxcnREZWx0YSA9IHNxcnQoZGVsdGEpOwoKICAgICAgICAvLyBwb3ByYXdueSBudW1lcnljem5pZSBhbGdvcnl0bQogICAgICAgIGRvdWJsZSB4MSA9ICgtYiAtIHNxcnREZWx0YSkgLyAoMiphKTsKICAgICAgICBkb3VibGUgeDIgPSBjIC8gKGEgKiB4MSk7CgogICAgICAgIGNvdXQgPDwgIngxID0gIiA8PCB4MSA8PCBlbmRsOwogICAgICAgIGNvdXQgPDwgIngyID0gIiA8PCB4MiA8PCBlbmRsOwogICAgfQoKICAgIHJldHVybiAwOwp9Cg==