Blame view

Giac_maj/giac-1.4.9/examples/arit/estpremier 631 Bytes
6663b6c9   adorian   projet complet av...
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
  estprem(n):={
  //teste si n est premier
    local prem,j,k;
    if ((irem(n,2)==0) or (irem(n,3)===0) or (n==1)) {
       return(false);
    }
    if ((n==2) or (n==3)) {
       return(true);
    }
    prem:=true;
    k:=5;
    while ((k*k<=n) and prem) {
       if (irem(n,k)==0 or irem(n,k+2)==0) {
          prem:=false;
       }
       else {
          k:=k+6;
       }
    }
    return(prem);
  };
  estpremc(n):={
  //utilise la fonction crible (fichier erathostene) pour tester si n est premier
  local prem,s;
  prem:=crible(floor(sqrt(n)));
  s:=size(prem);
  if (n==1) return(false);
  for (k:=0;k<s;k++){
   if (irem(n,prem[k])==0){
      return(false);
   } 
  }
  return(true);
  };