rnfpolredabs(nfinit(quadpoly(5,y)),x^7-14*x^5+56*x^3-56*x+22)
\\ rnfpolred(nfinit(quadpoly(1996,y)),quadray(1996,1))
\\ rnfpolred(nfinit(quadpoly(904,y)),quadray(904,1))
rnfpolredabs(nfinit(quadpoly(29,y)),quadray(29,17))
rnfpolredabs(nfinit(y^3-y-1),x^3-x-1)
rnfpolred(nfinit(quadpoly(1129,y)),quadray(1129,1));

nfinit(Pol([256,-2560,5120,6400,-60320,6976,116320,72560,-456615,270630,-129362]));
polredabs(quadpoly(14586217464))
{
p=x^32 - 1680*x^30 + 1026480*x^28 - 289096080*x^26 + 39933334980*x^24 -
2715474610800*x^22 + 88886222283600*x^20 - 1619990314513200*x^18 +
17928141864081750*x^16 - 125620995771054000*x^14 + 565267786831818000*x^12 -
1629524362237758000*x^10 + 2978275448322310500*x^8 - 3393290168363970000*x^6 +
2319247705779270000*x^4 - 866101453967610000*x^2 + 135345425000900625;
#polredabs(p,4)
}
polredabs(x^12+139968*x^10+24814646784*x^8+1854249948463104*x^6+177954917169813848064*x^4-1827912356210202139164672*x^2+68504919608701082757419237376);
polredabs(x^4+10^21*x^2+1)
polredabs(x^4+146077*x^2+2629386)
\\#1146
polredabs(x^9-4*x^7-3*x^6+9*x^5+8*x^4-6*x^3-9*x^2-4*x-1)
