Posté par
vincprof vincprofc'est bon j'ai trouvé une solution
> derivk:=proc(a)
> #k:=1;
> #while (deriv(k)<>a and k<=(a^2)/4) do
> # k:=k+1;
> # if (deriv(k)=a)=true then
> # RETURN(false);
> # fi;
> # true;
> # od;
> for j from 1 to floor((a^2)/4) do
> if deriv(j)=a then RETURN(false); fi;
> true;
> od;
> end:
>
>
Warning, `j` is implicitly declared local to procedure `derivk`
> derivk(10);
false
> lissolut2:=proc(n)
> L := [];
> k:=3;
> for a from 2 to n do
> if evalb((a mod 2)=1)=true then
> if evalb(isprime(a-2))=false then
> if (derivk(a))=true then L:=append1([L,[a]]); print(a);
> fi
> fi;
> fi;
> od;
> RETURN(L)
> end:
Warning, `L` is implicitly declared local to procedure `lissolut2`
Warning, `k` is implicitly declared local to procedure `lissolut2`
Warning, `a` is implicitly declared local to procedure `lissolut2`
> lissolut2(1000);
>
3
11
17
23
29
Warning, computation interrupted
mais il a mis plus de 5 minutes pour trouver ca...