// JiPP.lab3.cpp : Defines the entry point for the console application.
//

#include "stdafx.h"
#define _CRT_SECURE_NO_WARNINGS
#include <iostream>
#include <conio.h>
#include <math.h>
using namespace std;
/*public void podaj_dane(int *a, int *b, int *c, int *d)
{
printf("Podaj a: ");
scanf("%d", a);
printf("Podaj b: ");
scanf("%d", b);
printf("Podaj c ");
scanf("%d", c);
printf("Podaj d: ");
scanf("%d", d);

}*/


struct l_zespolone {
	float cz_r;
	float cz_u;
};

class Rownanie
{

protected:
	float suma, roznica, iloraz;
	float p_delta, delta;
public:
//	Rownanie();
//	~Rownanie();
	int *dane_wejsciowe;
	virtual void podaj_dane()=0;
	virtual void formatuj_dane()=0;
	virtual void dodaj()=0;
	virtual void odejmij()=0;
	virtual void pomnoz()=0;
};

class Rzeczywista :public Rownanie
{
protected:
	Rownanie *r;
	l_zespolone *l1;
	l_zespolone *l2;
	l_zespolone *l3;
	l_zespolone *l4;
	void podaj_dane();
	void formatuj_dane();
	void dodaj();
	void odejmij();
	void pomnoz();
	Rzeczywista();
	~Rzeczywista();
public:

	float oblicz_d();
	float oblicz_p_d(float delta);
};

Rzeczywista::Rzeczywista()
{
	r = new Rzeczywista;
	l1 = new l_zespolone;
	l2 = new l_zespolone;
	l3 = new l_zespolone;
	l4 = new l_zespolone;
	r->dane_wejsciowe = new int[4];
	l1->cz_r = 0;	l1->cz_u = 0;
	l2->cz_r = 0;	l2->cz_u = 0;
	l3->cz_r = 0;	l3->cz_u = 0;
	l4->cz_r = 0;	l4->cz_u = 0;
	suma = 0;	roznica = 0;	iloraz = 0;
	p_delta = 0;	delta = 0;


}
class Urojona:Rzeczywista
{
	Urojona();
	~Urojona();
public:
	float oblicz_d();
	float oblicz_p_d(float delta);
	void oblicz_rne(Rownanie *r, l_zespolone *a, l_zespolone *b, l_zespolone *c, l_zespolone *d);
};
Urojona::Urojona()
{
	r = new Urojona;
	l1 = new l_zespolone;
	l2 = new l_zespolone;
	l3 = new l_zespolone;
	l4 = new l_zespolone;
	r->dane_wejsciowe = new int[4];
	l1->cz_r = 0;	l1->cz_u = 0;
	l2->cz_r = 0;	l2->cz_u = 0;
	l3->cz_r = 0;	l3->cz_u = 0;
	l4->cz_r = 0;	l4->cz_u = 0;
	suma = 0;	roznica = 0;	iloraz = 0;
	p_delta = 0;	delta = 0;


}
void Rzeczywista::podaj_dane()
{
	for (int i = 0; i < 4; i++)
	{
		cin >> dane_wejsciowe[i];
	}
}
void Rownanie::formatuj_dane()
{
	if (dane_wejsciowe[0] != 0)
	{
		if (dane_wejsciowe[0] == 1)
			printf("xx");
		else if (dane_wejsciowe[0] == -1)
			printf("-xx");
		else
			printf("%dxx", dane_wejsciowe[0]);
	}
	if (dane_wejsciowe[1] != 0)
	{
		if (dane_wejsciowe[1] == 1 && dane_wejsciowe[0] == 0)
			printf("x");
		else if (dane_wejsciowe[1] == 1 && dane_wejsciowe[0] != 0)
			printf("+x");
		else if (dane_wejsciowe[1] == -1)
			printf("-x");
		else if (dane_wejsciowe[1] > 1 && dane_wejsciowe[0] == 0)
			printf("%dx", dane_wejsciowe[1]);
		else if (dane_wejsciowe[1] > 1 && dane_wejsciowe[0] != 0)
			printf("+%dx", dane_wejsciowe[1]);
		else
			printf("%dx", dane_wejsciowe[1]);
	}
	if (dane_wejsciowe[2] != 0)
	{
		if (dane_wejsciowe[2] == 1 && dane_wejsciowe[0] == 0 && dane_wejsciowe[1] == 0)
			printf("1");
		else if (dane_wejsciowe[2] == 1)
			printf("+1");
		else if (dane_wejsciowe[2] == -1)
			printf("-1");
		else if (dane_wejsciowe[2] > 1 && dane_wejsciowe[0] == 0 && dane_wejsciowe[1] == 0)
			printf("%d", dane_wejsciowe[2]);
		else if (dane_wejsciowe[2] > 1)
			printf("+%d", dane_wejsciowe[2]);
		else
			printf("%d", dane_wejsciowe[2]);
	}
	if (dane_wejsciowe[3] != 0)
	{
		if (dane_wejsciowe[3] == 1 && dane_wejsciowe[0] == 0 && dane_wejsciowe[1] == 0 && dane_wejsciowe[2] == 0)
			printf("i");
		else if (dane_wejsciowe[3] == 1)
			printf("+i");
		else if (dane_wejsciowe[3] == -1)
			printf("-i");
		else if (dane_wejsciowe[3] > 1 && dane_wejsciowe[0] == 0 && dane_wejsciowe[1] == 0 && dane_wejsciowe[2] == 0)
			printf("%di", dane_wejsciowe[3]);
		else if (dane_wejsciowe[3] > 1)
			printf("+%di", dane_wejsciowe[3]);
		else
			printf("%di", dane_wejsciowe[3]);
	}
	printf("\n");
}

float Rzeczywista::oblicz_d()
{
	this->delta = (r->dane_wejsciowe[1])*(r->dane_wejsciowe[1]) - (4.0*(r->dane_wejsciowe[0])*(r->dane_wejsciowe[2]));
	return  delta;
}
float Urojona::oblicz_d()
{
	this->delta= (-4.0 * (r->dane_wejsciowe[0] * r->dane_wejsciowe[3]));
	return delta;
}





float Rzeczywista::oblicz_p_d(float delta)
{
	return (sqrt((sqrt(oblicz_d()*oblicz_d() + (-4.0 * (r->dane_wejsciowe[0] * r->dane_wejsciowe[3]))*(-4.0 * (r->dane_wejsciowe[0] * r->dane_wejsciowe[3]))) + oblicz_d()) / 2));
}

float Urojona::oblicz_p_d(float delta)
{
	return (sqrt((sqrt(Rzeczywista::oblicz_d()*Rzeczywista::oblicz_d() + delta*delta) - Rzeczywista::oblicz_d()) / 2));
}

void Urojona::oblicz_rne(Rownanie *r, l_zespolone *x1, l_zespolone *x2, l_zespolone *x3, l_zespolone *x4)
{
	if (r->dane_wejsciowe[0] != 0 && r->dane_wejsciowe[3] == 0)
	{
		//CASE 1.1
		if (Rzeczywista::oblicz_d() > 0)
		{
			x1->cz_r = (-1.0*r->dane_wejsciowe[1] - sqrt(Rzeczywista::oblicz_d())) / 2;
			x2->cz_r = (-1.0*r->dane_wejsciowe[1] + sqrt(Rzeczywista::oblicz_d())) / 2;
		}

		//CASE 1.2
		else if (Rzeczywista::oblicz_d() == 0)
		{
			x1->cz_r = (-1.0*r->dane_wejsciowe[1]) / (2.0 * a);
		}
		//CASE 1.3
		else if (Rzeczywista::oblicz_d() < 0)
		{
			x1->cz_r = (-1.0*r->dane_wejsciowe[1]) / (2.0 * r->dane_wejsciowe[0]);
			x1->cz_u = -(sqrt((Rzeczywista::oblicz_d()*Rzeczywista::oblicz_d()) / 2) / (2.0 * r->dane_wejsciowe[0]));
			x2->cz_r = x1->cz_r;
			x2->cz_u = -x1->cz_u;
		}
	}
	//CASE 2
	else if (r->dane_wejsciowe[0] == 0 && r->dane_wejsciowe[1] != 0 && r->dane_wejsciowe[3] == 0)
	{
		x1->cz_r = (-1.0*r->dane_wejsciowe[2]) / (1.0*r->dane_wejsciowe[1]);
	}
	//CASE 3
	else if (r->dane_wejsciowe[0] == 0 && b == 0 && (c != 0 || d != 0))
	{
		printf("rownanie sprzeczne \n");
	}
	//CASE 4
	else if (a == 0, b == 0, c == 0, d == 0)
	{
		printf("rownanie tozsamosciowe \n");
	}
	//CASE 5
	else if (a == 0 && b != 0 && d != 0)
	{
		*x1r = (-1.0*c) / (1.0*b);
		*x1u = (-1.0*d) / (1.0*b);
	}
	else if (a != 0 && d != 0)
	{

		*x1r = ((-1.0*b - oblicz_p_dr(oblicz_d_r(a, b, c), oblicz_d_u(a, d))) / (2.0 * a));
		*x2r = ((-1.0*b - oblicz_p_dr(oblicz_d_r(a, b, c), oblicz_d_u(a, d))) / (2.0 * a));
		*x3r = ((-1.0*b + oblicz_p_dr(oblicz_d_r(a, b, c), oblicz_d_u(a, d))) / (2.0 * a));
		*x4r = ((-1.0*b + oblicz_p_dr(oblicz_d_r(a, b, c), oblicz_d_u(a, d))) / (2.0 * a));
		*x1u = (oblicz_p_du(oblicz_d_r(a, b, c), oblicz_d_u(a, d)) / (2 * a));
		*x2u = (-oblicz_p_du(oblicz_d_r(a, b, c), oblicz_d_u(a, d)) / (2 * a));
		*x3u = (oblicz_p_du(oblicz_d_r(a, b, c), oblicz_d_u(a, d)) / (2 * a));
		*x4u = (-oblicz_p_du(oblicz_d_r(a, b, c), oblicz_d_u(a, d)) / (2 * a));
	}
}
void dodaj(float deltar, float deltau, float *x1r, float *x2r, float *x3r, float *x4r, float *x1u, float *x2u, float *x3u, float *x4u, float *sr, float *su)
{
	if (deltar > 0 && deltau == 0)
	{
		*sr = *x1r + *x2r;
	}
	else if (deltar < 0 && deltau == 0)
	{
		*sr = *x1r + *x2r;
		*su = *x1u + *x2u;
	}
	else if (deltau != 0)
	{
		*sr = *x1r + *x2r + *x3r + *x4r;
		*su = *x1u + *x2u + *x3u + *x4u;
	}
}
void odejmij(float deltar, float deltau, float *x1r, float *x2r, float *x3r, float *x4r, float *x1u, float *x2u, float *x3u, float *x4u, float *rr, float *ru)
{
	if (deltar > 0 && deltau == 0)
	{
		*rr = *x1r - *x2r;
	}
	else if (deltar < 0 && deltau == 0)
	{
		*rr = *x1r - *x2r;
		*ru = *x1u - *x2u;
	}
	else if (deltau != 0)
	{
		*rr = *x1r - *x2r - *x3r - *x4r;
		*ru = *x1u - *x2u - *x3u - *x4u;
	}
}
void pomnoz(float deltar, float deltau, float *x1r, float *x2r, float *x3r, float *x4r, float *x1u, float *x2u, float *x3u, float *x4u, float *ilr, float *ilu)
{
	if (deltar > 0 && deltau == 0)
	{
		*ilr = *x1r * *x2r;
	}
	else if (deltar < 0 && deltau == 0)
	{
		*ilr = (*x1r * *x2r) - (*x1u * *x2u);
		*ilu = (*x1r * *x2u) - (*x1u * *x2r);
	}
	else if (deltau != 0)
	{
		*ilr = ((*x1r * *x2r) - (*x1u * *x2u))*((*x3r * *x4r) - (*x3u * *x4u)) -
			((*x1u * *x2u) + (*x2r * *x1r))* (((*x3r * *x4u) + (*x4r * *x3u)));
		*ilu = ((*x1r * *x2r) - (*x1u * *x2u))*((*x3r * *x4u) - (*x3u * *x4r)) +
			((*x1r * *x2u) + (*x2r * *x1u))* (((*x3r * *x4r) + (*x4u * *x3u)));
	}
}
void wyswietl_wynik(int a, int b, int c, int d, float x1r, float x2r, float x3r, float x4r, float x1u, float x2u, float x3u, float x4u, float sr, float su, float rr, float ru, float ilr, float ilu)
{
	printf("a=%d b=%d c=%d d=%d \n", a, b, c, d);
	if (a != 0 && d == 0)
	{
		//CASE 1.1
		if (oblicz_d_r(a, b, c) > 0)
		{
			printf("x1r=%f, x2r=%f, \n", x1r, x2r);
		}
		else if (oblicz_d_r(a, b, c) == 0)
		{
			printf("x1r=%f, \n", x1r);
		}
		else if (oblicz_d_r(a, b, c) < 0)
		{
			printf("x1r=%f, x2r=%f, x1u=%fi, x2u=%fi \n", x1r, x2r, x1u, x2u);
		}
	}
	else if (a == 0 && b != 0 && d == 0)
	{
		printf("x1r=%f, \n", x1r);
	}
	else if (a == 0 && b == 0 && (c != 0 || d != 0))
	{
		printf("rownanie sprzeczne \n");
	}
	else if (a == 0, b == 0, c == 0, d == 0)
	{
		printf("rownanie tozsamosciowe \n");
	}
	else if (a == 0 && b != 0 && d != 0)
	{
		printf("x1r=%f, x1u=%fi \n", x1r, x1u);
	}
	else
		printf("x1r=%f, x2r=%f, x3r=%f, x4r=%f, x1u=%fi, x2u=%fi, x3u=%fi, x4u=%fi \n", x1r, x2r, x3r, x4r, x1u, x2u, x3u, x4u);
	if (oblicz_d_r(a, b, c) > 0 && oblicz_d_u(a, d) == 0)
	{
		printf("sr=%f, rr=%f, ilr=%f, \n", sr, rr, ilr);
	}
	else
		printf("sr=%f, su=%f, rr=%f, ru=%f, ilr=%f, ilu=%f \n", sr, su, rr, ru, ilr, ilu);

}


int main()
{
	int a = 0, b = 0, c = 0, d = 0;
	float x1r = 0, x2r = 0, x3r = 0, x4r = 0, x1u = 0, x2u = 0, x3u = 0, x4u = 0;
	float sr = 0, su = 0, rr = 0, ru = 0, ilr = 0, ilu = 0;
	float deltar = oblicz_d_r(a, b, c);
	float deltau = oblicz_d_u(a, d);
	oblicz_rne(a, b, c, d, &x1r, &x2r, &x3r, &x4r, &x1u, &x2u, &x3u, &x4u);
	dodaj(deltar, deltau, &x1r, &x2r, &x3r, &x4r, &x1u, &x2u, &x3u, &x4u, &sr, &su);
	odejmij(deltar, deltau, &x1r, &x2r, &x3r, &x4r, &x1u, &x2u, &x3u, &x4u, &rr, &ru);
	pomnoz(deltar, deltau, &x1r, &x2r, &x3r, &x4r, &x1u, &x2u, &x3u, &x4u, &ilr, &ilu);
	wyswietl_wynik(a, b, c, d, x1r, x2r, x3r, x4r, x1u, x2u, x3u, x4u, sr, su, rr, ru, ilr, ilu);

	getchar();
	return 0;
}