puimod(a,p,n):={
//calcule de facon recursive la puissance rapide a^p modulo n
 if (p==0){
    return(1);
 } 	
 if (irem(p,2)==0){
    return(puimod(irem(a*a,n),iquo(p,2),n));
 } 
 return(irem(a*puimod(a,p-1,n),n));
};