Mbarek Haynoux

Problème:

m=[79,431,751,1039,1423,1759,1871,2143,2399,2447,3023,3343,3631,3727,4127];  
{
  for(i=1,#m,
    bnf=bnfinit(x^4-m[i],1);
    print(m[i],":", bnfisintnorm(bnf,2),":",bnfisintnorm(bnf,-2));
    );
}