src/integer/algebraic/cl_I_rootp_I.cc

   1 // rootp().
   2
   3 // General includes.
   4 #include "cl_sysdep.h"
   5
   6 // Specification.
   7 #include "cln/integer.h"
   8
   9
  10 // Implementation.
  11
  12 #include "cl_I.h"
  13
  14 namespace cln {
  15
  16 // Methode:
  17 // Falls x=0 oder x=1: x = x^n -> JA, x als Ergebnis.
  18 // Hier also x>1. Suche ein Integer y > 1 mit x=y^n.
  19 // Falls n >= integer_length(x): NEIN. (Da y>=2, müßte x>=2^n gelten.)
  20 // Hier also n>0 klein...
  21
  22 cl_boolean rootp (const cl_I& x, const cl_I& n, cl_I* w)
  23 {
  24         if (eq(x,0) || eq(x,1)) // x=0 oder x=1 ?
  25           { *w = x; return cl_true; } // ja -> x als Ergebnis
  26         if (n >= (cl_I)(unsigned long)integer_length(x))
  27           { return cl_false; }
  28         // Nun ist n < (integer-length x). Also paßt n in ein uintC.
  29         return cl_rootp_aux(x,cl_I_to_ulong(n),w);
  30 }
  31
  32 }  // namespace cln