polynomial operations

**dev7060** · Oct 15 '20, 05:13 PM

I remember working with polynomial arithmetic. I didn't use a particular named algo to achieve the results, just the common sense math. Linked lists serve the purpose and would represent the polynomial better than the array. Prepare a separate list for every expression. Assign three fields to every node; coefficient part, exponent part, and the pointer to the next expression. You need to know how the traversal in the linked list works to implement the operations.

**kristineg** · Oct 16 '20, 02:21 AM

could you please explain how to get into the implementation part?
thanks a lot for the previous explanation but ugh this is bugging me out for long now LOL. i tried reading online and am stuck in one way or the other. there's only code but no clear explanation to me. i understand the list preparation part but everything after that seems just out of the context.

Code:


Node *last = *head_ref;
new_node->c = c;
new_node->e = e; 
new_node->next = NULL;   
if (*head_ref == NULL)   
{   
*head_ref = new_node;
return;
}  
while (last->next != NULL)   
last = last->next;   
last->next = new_node;
return;

i'd really appreciate an example or snippet of your approach.

kristine

**Banfa** · Oct 16 '20, 02:36 PM

What language are you using? If you are using C++ I would say a std::map lends itself well to working with polynomials in a single variable with an unstated maximum degree.

Alternatively in C a variable sized array is probably better than a linked list because order is important but there isn't that much data to carry and being able to index straight into the array for constant for a specific degree seems to be to have some value in efficiency. (You could apply this argument to C++ too actually)

The code snippet you have posted appears to be part of a function to add a new node to a linked list although it appears to be missing the actual allocation of the node.

**dev7060** · Oct 16 '20, 03:55 PM

i tried reading online and am stuck in one way or the other. there's only code but no clear explanation to me. i understand the list preparation part but everything after that seems just out of the context.

You didn't specify what exact part you are having difficulty understanding from what you're reading.

could you please explain how to get into the implementation part?
i'd really appreciate an example or snippet of your approach.

I was building a scientific calculator back in the days as a hobby project. I've pulled up some code from it related to polynomial operations. Keep in mind that this involves a few bad practices since I was new to algorithms and programming at that time. More like I was trying to be the jack of all trades and ended up mixing c & cpp concepts for some reason. But everything does work as one lang is the superset of the other. Hopefully, it would give you the idea of the algo and implementation and prove helpful. Here are few things you may need to consider to change (for c++); use new instead of malloc, use class data type for the list, assign proper names to the variables, free the memory wherever required, plus other patches that you feel like returning and accepting relevant pointers from/to the functions, etc. Also, the code isn't prepared for the boundary and invalid inputs. For the addition part; if any of the two lists is empty, the other will be the result. If both are non-empty; the coefficients of the nodes having the same exponents are added. The next step is to check what list is completely traversed, the rest nodes of the other are then added to the final expression as it is. In the multiplication one, nested loops are used to multiply the fields of one expression with another and to cover every node. Add the coefficients of the terms having the same exponents to produce the final expression. I've used a separate list to keep track of the unique coefficients. Other checks can also be included that were there in the addition part.

Code:

/*
*
*
* Author: deV
* Last modified: Jan 28, 2018
* License: GPL-3.0
** Master branch src **
*/
#include <iostream>
#include <stdlib.h>
#include "poly.h"
using namespace std;
Poly::Poly()
{
	start = NULL;
	begin1 = NULL;
	add = NULL;
	global_ch = 0;	//global_ch=1 -> subtraction
}

Poly::~Poly() {}

void Poly::first()
{
	int cont;
	struct POLYS *node, *temp1;
	cout << "Insertion of first expression:\n";
	do {
		node = (struct POLYS *) malloc(sizeof(struct POLYS));
		cout << "Enter coef: ";
		cin >> node->coef;
		cout << "Enter exponent:";
		cin >> node->exp;
		node->next = NULL;
		temp1 = start;
		if (start == NULL)
		{
			start = node;
		}
		else
		{
			while (temp1->next != NULL)
			{
				temp1 = temp1->next;
			}

			temp1->next = node;
		}

		cout << "First expression is: ";
		struct POLYS *view1 = start;
		int a = 0, i = 0;
		while (view1 != NULL)
		{
			a++;
			view1 = view1->next;
		}

		view1 = start;
		while (view1 != NULL)
		{
			if (view1->coef != 0)
				cout << "[" << view1->coef << "x^" << view1->exp << "]";
			if (i < (a - 1))
				cout << "+";
			i++;
			view1 = view1->next;
		}

		cout << "\nPress 5 to continue the first expression insertion: ";
		cin >> cont;
	} while (cont == 5);
}

void Poly::second()
{
	int cont;
	struct POLYS *node, *temp2;
	cout << "Insertion of second expression:\n";
	do {
		node = (struct POLYS *) malloc(sizeof(struct POLYS));
		cout << "Enter coef: ";
		cin >> node->coef;
		cout << "Enter exponent:";
		cin >> node->exp;
		node->next = NULL;
		temp2 = begin1;
		if (begin1 == NULL)
		{
			begin1 = node;
		}
		else
		{
			while (temp2->next != NULL)
			{
				temp2 = temp2->next;
			}

			temp2->next = node;
		}

		cout << "Second expression is: ";
		struct POLYS *view2 = begin1;
		int a = 0, i = 0;
		while (view2 != NULL)
		{
			a++;
			view2 = view2->next;
		}

		view2 = begin1;
		while (view2 != NULL)
		{
			if (view2->coef != 0)
				cout << "[" << view2->coef << "x^" << view2->exp << "]";
			if (i < (a - 1))
				cout << "+";
			i++;
			view2 = view2->next;
		}

		cout << "\nPress 5 to continue the second expression insertion: ";
		cin >> cont;
	} while (cont == 5);
}

void Poly::addpoly()
{
	struct POLYS *temp1, *temp2, *temp3, *node;
	if (start == NULL && begin1 == NULL)
	{
		cout << "No expression!";
	}
	else if (start != NULL && begin1 == NULL)
	{
		temp1 = start;
		while (temp1 != NULL)
		{
			cout << temp1->coef << "x^" << temp1->exp;
			temp1 = temp1->next;
		}
	}
	else if (start == NULL && begin1 != NULL)
	{
		temp2 = begin1;
		while (temp2 != NULL)
		{
			cout << temp2->coef << "x^" << temp2->exp;
			temp2 = temp2->next;
		}
	}
	else if (start != NULL && begin1 != NULL)
	{
		temp1 = start;
		temp2 = begin1;
		while (temp1 != NULL && temp2 != NULL)
		{
			node = (struct POLYS *) malloc(sizeof(struct POLYS));
			if (temp1->exp > temp2->exp)
			{
				node->coef = temp1->coef;
				node->exp = temp1->exp;
				node->next = NULL;
				temp1 = temp1->next;
			}
			else if (temp2->exp > temp1->exp)
			{
				node->coef = temp2->coef;
				node->exp = temp2->exp;
				node->next = NULL;
				temp2 = temp2->next;
			}
			else if (temp1->exp == temp2->exp)
			{
				node->coef = temp1->coef + temp2->coef;
				node->exp = temp1->exp;
				node->next = NULL;
				temp1 = temp1->next;
				temp2 = temp2->next;
			}

			temp3 = add;
			if (add == NULL)
			{
				add = node;
			}
			else
			{
				while (temp3->next != NULL)
				{
					temp3 = temp3->next;
				}

				temp3->next = node;
			}
		}	//end of while
		if (temp1 != NULL && temp2 == NULL)
		{
			while (temp1 != NULL)
			{
				node = (struct POLYS *) malloc(sizeof(struct POLYS));
				node->coef = temp1->coef;
				node->exp = temp1->exp;
				node->next = NULL;
				temp3 = add;
				while (temp3->next != NULL)
					temp3 = temp3->next;
				temp3->next = node;
				temp1 = temp1->next;
			}
		}
		else if (temp2 != NULL && temp1 == NULL)
		{
			while (temp2 != NULL)
			{
				node = (struct POLYS *) malloc(sizeof(struct POLYS));
				node->coef = temp2->coef;
				node->exp = temp2->exp;
				node->next = NULL;
				temp3 = add;
				while (temp3->next != NULL)
					temp3 = temp3->next;
				temp3->next = node;
				temp2 = temp2->next;
			}
		}

		if (global_ch == 0)
			cout << "Added expression is: ";
		else if (global_ch == 1)
			cout << "Subtracted expression is: ";
		temp3 = add;
		int a = 0, i = 0;
		while (temp3 != NULL)
		{
			a++;
			temp3 = temp3->next;
		}

		temp3 = add;
		while (temp3 != NULL)
		{
			if (temp3->coef != 0)
				cout << "[" << temp3->coef << "x^" << temp3->exp << "]";
			if (i < (a - 1))
			{
				cout << "+";
			}

			i++;
			temp3 = temp3->next;
		}
	}	//end of outer if else
}

void Poly::subpoly()
{
	struct POLYS *temp = begin1;
	while (temp != NULL)
	{
		temp->coef = -temp->coef;
		temp = temp->next;
	}

	global_ch = 1;
	addpoly();
	global_ch = 0;
}

void Poly::mulpoly()
{
	struct POLYS *temp1 = start;
	struct POLYS *temp2 = begin1;
	struct POLYS * mul_beg;
	mul_beg = NULL;
	struct UN_COE_LIST * unique_coef_list;
	unique_coef_list = NULL;
	struct POLYS * temp4;
	while (temp1 != NULL)
	{
		temp2 = begin1;
		while (temp2 != NULL)
		{
			struct POLYS * node;
			node = (struct POLYS *) malloc(sizeof(struct POLYS));
			node->coef = (temp1->coef) *(temp2->coef);
			node->exp = (temp1->exp) + (temp2->exp);
			node->next = NULL;
			if (mul_beg == NULL)
			{
				mul_beg = node;
			}
			else
			{
				temp4 = mul_beg;
				while (temp4->next != NULL)
				{
					temp4 = temp4->next;
				}

				temp4->next = node;
			}

			temp2 = temp2->next;
		}

		temp1 = temp1->next;
	}

	temp4 = mul_beg;
	int indicator;
	//loop to add the unique exp values to the unique_coef_list
	while (temp4 != NULL)
	{
		struct UN_COE_LIST * node;
		node = (struct UN_COE_LIST *) malloc(sizeof(struct POLYS));
		node->exp = temp4->exp;
		node->next = NULL;
		indicator = 10;
		//10=coef not present
		//11=coef already present
		struct UN_COE_LIST * temp5;
		temp5 = unique_coef_list;
		while (temp5 != NULL)
		{
			if (temp5->exp == node->exp)
			{
				indicator = 11;
			}
			else
			{
				indicator = 10;
			}

			temp5 = temp5->next;
		}

		if (indicator == 10)
		{
			if (unique_coef_list == NULL)
			{
				unique_coef_list = node;
			}
			else
			{
				struct UN_COE_LIST * temp6;
				temp6 = unique_coef_list;
				while (temp6->next != NULL)
				{
					temp6 = temp6->next;
				}

				temp6->next = node;
			}
		}

		temp4 = temp4->next;
	}

	struct UN_COE_LIST * temp6;
	temp6 = unique_coef_list;
	int a = 0, i = 0;
	while (temp6 != NULL)
	{
		a++;
		temp6 = temp6->next;
	}

	temp6 = unique_coef_list;
	/*testing unique_coef_list
	while(temp6!=NULL){
	        cout<<"["<<temp6->exp<<"]";
	        temp6=temp6->next;
	    }

	*/
	/*testing mul_beg list
	temp6=mul_beg;
	while(temp6!=NULL){
	        if(temp6->coef!=0)
	            cout<<"["<<temp6->coef<<"x^"<<temp6->exp<<"]";
	        if(i < (a-1)){
	            cout<<"+";
	        }

	        i++;
	        temp6=temp6->next;
	    }

	    */
	struct UN_COE_LIST * ntemp1;
	struct POLYS * ntemp2;
	struct POLYS *final_list = NULL;
	ntemp1 = unique_coef_list;
	ntemp2 = mul_beg;
	int
	var = 0;
	while (ntemp1 != NULL)
	{
		ntemp2 = mul_beg;
		while (ntemp2 != NULL)
		{
			if (ntemp2->exp == ntemp1->exp)
			{
				var =
				var +(ntemp2->coef);
			}

			ntemp2 = ntemp2->next;
		}

		struct POLYS * node;
		node = (struct POLYS *) malloc(sizeof(struct POLYS));
		node->coef =
			var;
		node->exp = ntemp1->exp;
		node->next = NULL;
		struct POLYS * ptr;
		ptr = final_list;
		//to skip the same exp term
		int p = 0;
		//p=0 if not present
		//p=1 if present
		while (ptr != NULL)
		{
			if (ptr->exp == node->exp)
				p = 1;
			ptr = ptr->next;
		}

		if (p == 0)
		{
			if (final_list == NULL)
			{
				final_list = node;
			}
			else
			{
				struct POLYS * temp7;
				temp7 = final_list;
				while (temp7->next != NULL)
				{
					temp7 = temp7->next;
				}

				temp7->next = node;
			}
		}

		p = 0;
		var = 0;
		ntemp1 = ntemp1->next;
	}

	struct POLYS * temp8;
	temp8 = final_list;
	i = 0;
	a = 0;
	while (temp8 != 0)
	{
		a++;
		temp8 = temp8->next;
	}

	temp8 = final_list;
	cout << "Multiplied expression: ";
	while (temp8 != NULL)
	{
		if (temp8->coef != 0)
			cout << "[" << temp8->coef << "x^" << temp8->exp << "]";
		if (i < (a - 1))
		{
			cout << "+";
		}

		i++;
		temp8 = temp8->next;
	}
}

Structure blueprints:

Code:

struct POLYS{
    int coef;
    int exp;
    struct POLYS *next;
};
struct UN_COE_LIST{
    int exp;
    struct UN_COE_LIST *next;
};

**kristineg** · Oct 17 '20, 03:38 AM

I am using C++. thank you both Banfa and dev7060 for the advice. it starts making sense now when i have a little theory. dev7060, i will have a deep study of your code and reach out for the doubts. thank you so much. your post was very helpful.

polynomial operations

polynomial operations

Comment

Comment

Comment

Comment

Comment