- 注册时间
- 2007-12-27
- 最后登录
- 1970-1-1
- 威望
- 星
- 金币
- 枚
- 贡献
- 分
- 经验
- 点
- 鲜花
- 朵
- 魅力
- 点
- 上传
- 次
- 下载
- 次
- 积分
- 40155
- 在线时间
- 小时
|
发表于 2019-11-7 08:34:18
|
显示全部楼层
(08:28) gp > y=x^2+x+1
%1 = x^2 + x + 1
(08:28) gp > Mod(y,x^3+t*x+1)^2
%2 = Mod((-t + 3)*x^2 + (-2*t + 1)*x - 1, x^3 + t*x + 1)
(08:28) gp > Mod(y,x^3+t*x+1)^3
%3 = Mod((t^2 - 6*t + 3)*x^2 + (3*t^2 - 5*t - 3)*x + (3*t - 5), x^3 + t*x + 1)
(08:28) gp > [0,0,1;1,1,1;(-t+3),(-2*t+1),-1;(t^2 - 6*t + 3),(3*t^2 - 5*t - 3),(3*t - 5)]
%4 =
[ 0 0 1]
[ 1 1 1]
[ -t + 3 -2*t + 1 -1]
[t^2 - 6*t + 3 3*t^2 - 5*t - 3 3*t - 5]
(08:31) gp > %4~
%5 =
[0 1 -t + 3 t^2 - 6*t + 3]
[0 1 -2*t + 1 3*t^2 - 5*t - 3]
[1 1 -1 3*t - 5]
(08:31) gp > matker(%5)
%6 =
[-t^2 + 2*t - 4]
[ t^2 - 3*t + 6]
[ 2*t - 3]
[ 1]
(08:31) gp > y^3+(2*t-3)*y^2+(t^2 - 3*t + 6)*y+(-t^2 + 2*t - 4)
%7 = x^6 + 3*x^5 + (2*t + 3)*x^4 + (4*t + 1)*x^3 + (t^2 + 3*t + 3)*x^2 + (t^2 + t + 3)*x + t
(08:32) gp > polroots(x^3+2*x+1)
%8 = [-0.45339765151640376764474653900019218887 + 0.E-38*I, 0.22669882575820188382237326950009609443 - 1.4677115087102242702017782875332674014*I, 0.22669882575820188382237326950009609443 + 1.4677115087102242702017782875332674014*I]~
(08:32) gp > kill(y);
(08:32) gp > y^3+(2*t-3)*y^2+(t^2 - 3*t + 6)*y+(-t^2 + 2*t - 4)
%10 = (y - 1)*t^2 + (2*y^2 - 3*y + 2)*t + (y^3 - 3*y^2 + 6*y - 4)
(08:32) gp > subst(%,t,2)
%11 = y^3 + y^2 + 4*y - 4
(08:32) gp > polroots(%11)
%12 = [0.75217177888418654405728207878363123751 + 0.E-38*I, -0.87608588944209327202864103939181561876 - 2.1331684598630377469079688991623451133*I, -0.87608588944209327202864103939181561876 + 2.1331684598630377469079688991623451133*I]~
(08:33) gp > %8[1]^2+%8[1]+1
%13 = 0.75217177888418654405728207878363123751 + 0.E-38*I
(08:33) gp > %8[2]^2+%8[2]+1
%14 = -0.87608588944209327202864103939181561877 - 2.1331684598630377469079688991623451133*I
(08:33) gp > %8[3]^2+%8[3]+1
%15 = -0.87608588944209327202864103939181561877 + 2.1331684598630377469079688991623451133*I |
|