#include <iostream>
#include <cmath>
#include <stdlib.h>
using namespace std;

int main()
{
    int N, m = 0, p = 0, n = 0, count_pos = 0, count_neg = 0;
    double a, b, c, e, f, step;
    cout << "N="; cin >> N;

    double Y[30], eps = 1e-9;
    cout << "a="; cin >> a;
    cout << "b="; cin >> b;
    cout << "c="; cin >> c;
    cout << "e="; cin >> e;
    cout << "f="; cin >> f;

    step = (f - e) / (N - 1);
    cout << "step=" << step << endl;
    cout << endl;

    cout << "   Y" << "       " << "  Ypozitive" << "    " << "      Ynegative" << endl;

    for (double x = e; x <= f + eps; x += step)
    {
        Y[m] = a * x*x + b * x + c;
        cout << Y[m] << " " << endl;
        if (Y[m] > eps) count_pos++;
        else if (Y[m] < -eps) count_neg++;
        m++;
    }
    cout << endl;

    if (count_pos)
    {
        double*Ypozitive = new double[count_pos];

        for (int i = 0; i < N; i++)
            if (Y[i] > eps) Ypozitive[p++] = Y[i];

        for (int i = 0; i < count_pos; i++)
            cout << Ypozitive[i] << " " << endl;
        cout << endl;
        delete[] Ypozitive;
    }

    if (count_neg)
    {
        double*Ynegative = new double[count_neg];

        for (int i = 0; i < N; i++)
            if (Y[i] < -eps) Ynegative[n++] = Y[i];

        for (int i = 0; i < n; i++)
            cout << Ynegative[i] << " " << endl;
        cout << endl;
        delete[] Ynegative;
    }

    system("pause");
    return 0;
}