Richard Heylen on Sat, 16 Nov 2013 22:43:37 +0100
|
[Date Prev] [Date Next] [Thread Prev] [Thread Next] [Date Index] [Thread Index]
Shorter way to convert polynomial in x^3 to x?
|
- To: pari-users@pari.math.u-bordeaux.fr
- Subject: Shorter way to convert polynomial in x^3 to x?
- From: Richard Heylen <richardheylen@gmail.com>
- Date: Sat, 16 Nov 2013 21:43:26 +0000
- Delivery-date: Sat, 16 Nov 2013 22:43:37 +0100
- Dkim-signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=gmail.com; s=20120113; h=mime-version:date:message-id:subject:from:to:content-type; bh=Hu4L1r6pvV6d+Tu90UuOEBvfhRF9Z1VF7lyPlcUHJcU=; b=yfIc8b9ouXyszbuxIUNKRnyWo2tYkgJyawXrcWDcDS1kVNJgMGnMUmHAfI3Hr8eiOq PY6U7p8PjgSECQ16ab4O7WI6b0wF1PDc3tf3Mq7At9lHzt514vWKMY7ApiMURpp9ilQ0 BcAKs1CAZKsYCxYr6ZqA+hIMrqRQtZaRP3bEKqrSuyWTRACw+LZfbJsMYw8Z3N7shYem 6F6Ky5IGZxJegkH+IfuuFSgU6RzyT8sWuNRhu5/8wAnpJwwhaasT7cpbLbwQ9tFM0NQ7 tb0s+4Wt/swcWqf8WHLeaMGl1VrNsOtJ3NkMPCS/5StqS6jUmMMUl5gkiDr1VOL+BbCW 6hVA==
? f=factor(elldivpol(ellinit([0,1]),13))[1,1]
%32 = 13*x^12 + 52*x^9 + 1536*x^6 + 1024*x^3 + 256
? Pol(vector(poldegree(f)/3+1,i,polcoeff(f,poldegree(f)-3*i+3)))
%33 = 13*x^4 + 52*x^3 + 1536*x^2 + 1024*x + 256
I tried using subst(f,x,y^(1/3)) but to no avail. Is there a short way
of doing this?
Richard