void UniDivFC(poly_fc & q,poly_fc & f,poly_fc & g)
{
poly_fc r;
int k,n;
r.resize(f.size());
for(int i=0;i<f.size();i++)mpfc_set(r[i],f[i]);
n=g.size()-1;
k=r.size()-g.size();
if(k<0)
{
q.resize(0);return;
}
q.resize(k+1);
mpfc_t t;
mpfc_init(t);
do
{
mpfc_div(q[k],r[n+k],g[n]);
if(!mpfc_iszero(q[k]))
{
for(int i=0;i<n;i++)
{
uint j=n+k-1-i;
mpfc_mul(t,q[k],g[j-k]);
mpfc_sub(r[j],r[j],t);
}
}
} while (k--);
mpfc_clear(t);
r.resize(0);
}
inline complexAP complexAP::operator /(complexAP c)
{
complexAP tmp;
mpfc_div(&tmp.value, &value, &c.value);
return tmp;
}
void divide(elem &result, elem a, elem b) const { mpfc_div(&result,&a,&b); }
