fpdiv(n):={ //renvoie la liste des diviseurs de n en utilisant factprem local l1,l2,l3,d,ex,s; l3:=factprem(n); s:=size(l3); l1:=[1]; for (k:=0;k<s-1;k:=k+2) { l2:=[]; d:=l3[k]; ex:=l3[k+1]; for (j:=1;j<=ex;j++) { l2:=concat(l2,l1*(d^j)); } l1:=concat(l1,l2); } return(l1); }