Кубическое уравнение

#include <iostream>
 
using namespace std;
 
int horner(int *p, int n, int x)
{
    int i,r=0;
    for (i=0; i<n; i++) r=r*x+p[i];
    return r;
}
 
int main()
{
    int n,i,j,t,q,u,ptr,z=0;
    cout << "n=";
    cin >> n;
    
    int *p = new int[n+1];
    int *roots = new int[n];
    
    for (i=0; i<=n; i++)
    {
        cout << "p[" << (n-i) << "]=";
        cin >> p[i];
        if (p[i]==0) z++;
    }
 
    if (z == n)
    {
        cout << "-1" << endl;
        delete [] p;
        delete [] roots;
        return 0;
    }
 
    q=p[n];
    ptr=0;
    i=1;
    
    while(i*i<abs(q))
    {
        if (q%i==0)
        {
            if (horner(p,n+1,i)==0)
            {
                roots[ptr++]=i;
                u=q/i;
                if ((u != i) && (horner(p,n+1,u)==0))
                {
                   roots[ptr++]=i;
                   if (horner(p,n+1,-u)==0) roots[ptr++]=-u;
                }
                if (horner(p,n+1,-i)==0) roots[ptr++]=-i;
            }
        }
        
        i++;
    }    
 
    cout << ptr << endl;
    
    if (ptr>0)
    {
        for (i=0; i<ptr-1; i++)
            for (j=i+1; j<ptr; j++)
                if (roots[i]>roots[j])
                {
                    t=roots[i];
                    roots[i]=roots[j];
                    roots[j]=t;
                }
        for (i=0; i<ptr; i++) cout << roots[i] << endl;        
    }
 
    delete [] p;
    delete [] roots;
    return 0;
}