rnfpolredbest(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))
K=nfinit(quadpoly(29,y)); T=quadray(29,17)
rnfpolredabs(K,T)
rnfpolredbest(K,T)
rnfpolredabs(K,[T,100],16)
rnfpolredbest(K,[T,100])
K=nfinit(y^3-y-1); T=x^3-x-1;
rnfpolredabs(K,T)
rnfpolredabs(K,T,2)
[P,a]=rnfpolredbest(K,T,1);
subst(K.pol,y,a)
[P,a,b]=rnfpolredbest(K,T,3);
subst(K.pol,y,a)
substvec(T,[x,y],[a,b])
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)
\\#1228
polredabs(x^5-13*x^3+3*x^2+5*x-1)
\\#1229
polredabs(x^6+21471450*x^2+71643071500)
polredbest(x^6+21471450*x^2+71643071500)
polredbest(4*x^4+146077*x^2+2629386)

polredabs(2*x+1)
polredabs(2*x+1,1)
polred([x,[1]])
polred(2*x+1)
polred(2*x+1,1)
polred(4*x^2+1,2)
polredbest(2*x+1)
polredbest(2*x+1,1)
polredbest(4*x^2+1,1)
polredord(2*x+1)
polredord(4*x^2+1)
\\ #1519, test T->unscale != 1
polredbest(x^8+24*x^6+80*x^5+1040*x^4-4288*x^3-12736*x^2-61952*x+311296,1)
\\ #1511
polredabs(x^16-4*x^15-334*x^14+264*x^13+32231*x^12+57392*x^11-1031422*x^10-3628868*x^9+7185297*x^8+42417784*x^7+11283472*x^6-137773504*x^5-127243504*x^4+69059728*x^3+56307944*x^2-6264432*x+6436)
