- 注册时间
- 2022-3-21
- 最后登录
- 1970-1-1
- 威望
- 星
- 金币
- 枚
- 贡献
- 分
- 经验
- 点
- 鲜花
- 朵
- 魅力
- 点
- 上传
- 次
- 下载
- 次
- 积分
- 8890
- 在线时间
- 小时
|
发表于 2024-5-20 20:49:05
|
显示全部楼层
感谢提醒。
以下是gp/pari代码
test09() = {
a = 16557825695384603350011256367;
b = 7402106235083715173833501130;
c = 22633927789404195625906842766836083950159361724001877;
d = 10122181481681239360668340460448686303713447708473285;
t5 = quadunit(20) - 2;
d2 = (a-b*t5)^2 - 2*11^4*(c-d*t5);
print(d2);
d2 = (a+b*t5)^2 - 2*11^4*(c+d*t5);
print(d2);
default(realprecision,100);
c = 8669126517449449573024696780490669721303755676965485;
d = 3876911040311425488549116190343499942853937801760782;
T1 = (-a+b*t5+115*sqrt(c-d*t5)*I)/2;
T2 = (-a+b*t5-115*sqrt(c-d*t5)*I)/2;
T3 = (-a-b*t5+115*sqrt(c+d*t5)*I)/2;
T4 = (-a-b*t5-115*sqrt(c+d*t5)*I)/2;
T = [T1,T2,T3,T4];
s = prod(n=1,4,x-T[n]);
print(round(s));
w = exp(2*Pi*I/5);
W = vector(5,n,w^(n-1));
Ls = [
[0,0,0,0],[1,4,2,3],[2,3,4,1],[3,2,1,4],[4,1,3,2]
];
X = List();
for(n = 1,5,
x0 = -209341+11^(3/5)*sum(k=1,4,T[k]^0.2*W[Ls[n][k]+1]);
listput(X,x0/10);
);
eqn = prod(n=1,5,x-X[n])*32;
print(eqn);
print(round(eqn));
};
test09(); |
|