#include <iostream>;
using namespace std;

double *poly_div(double *a, double  *b, int deg_a, int deg_b)
{
	double *ab = new double[deg_a - deg_b + 1];
	for (int i = 0; i != deg_a - deg_b + 1; i++)
	{
		ab[deg_a - deg_b - i] = a[deg_a - 1 - i] / b[deg_b - 1 ];
		cout << a[deg_a - 1] << endl;
		for (int j = 0; j != deg_b; j++)
		{
			a[deg_a - 1 - j] -= b[deg_b - 1 - j] * ab[deg_a - deg_b - i];
		}
	}
	return (ab);
}

int main ()
{
	const int deg_a = 5;
	const int deg_b = 2;
	double a[deg_a] = { 1, 1, 2, 1, 1};
	double b[deg_b] = { 1, 1 };

	double *ab = poly_div(a, b, deg_a, deg_b);
	for (int i = 0; i != deg_a - deg_b + 1; i++)
	{
		cout << ab[i] << "\t";
	}
	cout << endl;
	system("pause");
	return (0);
}