Деление многочленов. Нахождение остатка

[@Manya@]

Добрый вечер!
Вот собственно код.

#include <stdio.h>
int *remainder(int *a, int &len1, int *b, int len2);
int *remainder(int *p1, int &s1, int *p2, int s2){
	    int *ptr1=p1;
        int del;
        int i;
        int *result = (int *)malloc(sizeof(int)*(s1));
        int sresult =0;
        while (s1>=s2) {
                del = (int)(*ptr1 / *p2);
                result[sresult]=del;
                sresult++;
                for (i=0; i<s2; i++) {
                        *(ptr1+i) -= del * p2[i];
                }
                ptr1++;
                s1--;
        }
		return p1;
}
int main() {
	int s1=9;
	int s2=8;
	int *pol1 = new int [s1];
	int *pol2=new int [s2];
	for (int s=0; s<s1; s++)
	{
		printf("%d ", s);
		printf("coefficient ");
		scanf("%d", &pol1[s]);
	}
	for (int s=0; s<s2; s++)
	{
		printf("%d ", s);
		printf("coefficient ");
		scanf("%d", &pol2[s]);
	}
	printf("\n");
	printf("-> ");
	 int *rem=new int[s1];
	 rem=remainder(pol1, s1, pol2, s2);
	 for (int s=0; s<s1; s++)
		 printf("%d", rem[s]);
        printf("\n");
}

Многочлены вводятся так: например, 1+x^3 будет 1001. Когда коэффициент делителя нулевой, то естественно возникает ошибка деления на 0. Как бы это исправить?
Спасибо

[Sweta]

Проверяйте перед делением знаменатель на "0", и выводите соответствующее сообщение с предложением ввести новые коэффициенты или выходите на конец программы.

[@Manya@]

Если бы это было возможно, то мне не пришлось бы спрашивать
Проблема в том, что это лишь одна функция большой программы, многочлены получаются в результате некоторых матричных операций, а там так не пройдет) Все равно надо поделить. Да и к тому же, там ошибка при делении ненулевых

Вот так на 0 не ругается, но остаток не верен . Обыкновенное деление столбиком ТЫК
ЗЫ многочлен записывается начиная с младшей степени
Периодически даже правильно работает. Но не всегда

int *remaind(int *p1, int &len1, int *p2, int len2)
{
	int k= len1-1;
	for (int i=0; i<len1/2; i++)
	{
		int temp = p1[i];
		p1[i] = p1[k-i];
		p1[k-i] = temp;
	}
	k= len2-1;
	for (int i=0; i<len2/2; i++)
	{
		int temp = p2[i];
		p2[i] = p2[k-i];
		p2[k-i] = temp;
	}
	int *ptr1;
	ptr1= p1;
	while (len1>=len2)
	{
		int del = (int)(*ptr1/ *p2);
		for (int i=0; i<len2; i++)
			*(ptr1+i) -= del * p2[i];
		len1--;
		ptr1++;
	}
	p1 = ptr1;
	return p1;
}