Re: divisors() for non-integer arguments

[Bill Allombert <allomber@math.u-bordeaux.fr>]

On Fri, Jun 25, 2004 at 02:51:43PM +0200, Jeroen Demeyer wrote:
> Hello all,
> 
> I noticed that the function divisors() only works for integers.
> I'm not a PARI developer, but I'm sure it cannot be hard at all to
> implement it also for polynomials (main difference would probably be
> that sorting doesn't really make any sense).
> I'd do it myself, but the source code is so cryptic.

If you don't mind cryptic code, you can use the GP function below:

div(P)=local(F,L);F=factor(P);L=[];forvec(V=vector(#F[,2],i,[0,F[i,2]]),L=concat(L,[factorback(concat(Mat(F[,1]),Mat(V~)))]));L

Cheers,
Bill.