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));
        );
    }