xtoymn.c:

/* Modular exponentiation */
long xtoymn(x,y,n) long x,y,n; {
   long result;

   if (n<=0)  return 0;
   result = 1;

   while (y>0) {
     if (y & 1) {
       result = xbyymn(result,x,n);
     } 
     x = xbyymn (x,x,n) ;
     y = y>>1 ;
   }

   return result;
}

Go to ...

• Factoring Page
• John Halleck's Home Page

This page is http://www.cc.utah.edu/~nahaj/factoring/xtoymn.c.html
© Copyright 2003 by John Halleck, All Rights Reserved.
This snapshot was last modified on January 12th, 2009
And the underlying file was last modified on December 12th, 2005