mpi.h

Go to the documentation of this file.

/* Part of the C++ Multi-Precision arithmetic library

   Written by Brian Sipos <bjs5075@rit.edu> or <brian@sipos.yi.org>
   Copyright 2001,2002,2003 Brian Sipos

   The MP++ library is free software; you can redistribute it and/or modify
   it under the terms of the GNU General Public License as published by
   the Free Software Foundation; either version 2 of the License, or
   (at your option) any later version.

   The MP++ library is distributed in the hope that it will be useful,
   but WITHOUT ANY WARRANTY; without even the implied warranty of
   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
   GNU General Public License for more details.

   You should have received a copy of the GNU General Public License
   along with the MP++ library; if not, write to the Free Software
   Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA
*/

#ifndef MPI_H
#define MPI_H

#include "mpn.h"

namespace mp {

class mpi_ : public mpn_ {
 public:
  typedef mp::int_type int_type; 

  static const byte capabilities();

  static void setDefaultBase(radix_type b);

  static radix_type getDefaultBase();

  mpi_() {}
  mpi_(int_type n, size_type size = 0) {
    init(n, size);
  }
  mpi_(const mpi_ &n, size_type size = 0) {
    init(n, size);
  }

  void init(int_type n, size_type size = 0);
  void init(const mpi_ &n, size_type size = 0);

  void clear(); 

  void resize(size_type size);

  void swap(mpi_ &n);

  int_type toInt() const;

  bool isNeg() const;

  size_type size() const;

  void setString(const std::string &s, radix_type b = 0);

  std::string getString(radix_type b = string_base) const;

  void shift_l(size_type n);
  void shift_r(size_type n);

  void print() const;

  void add(const mpi_ &n);

  void sub(const mpi_ &n);

  void mul(const mpi_ &n);

  void div(const mpi_ &n);

  void mod(const mpi_ &n);

  void negate(); 
  void abs(); 

  void sqr(); 

  mpi_ pow(int_type e) const;

  mpi_ pow_mod(const mpi_ &e, const mpi_ &m) const;

  int compare(const mpi_ &n) const;

  void checkNeg(); 

 protected:
  static radix_type string_base; 

  bool neg_; 

 private:
  void add_sub(const mpi_ &y, bool subtract); 
};

class mpi : public mpi_ {
 public:
  mpi(const mpi_ &n) : mpi_(n) {} // Casting up

  mpi() {}
  mpi(int_type n, size_type size = 0) : mpi_(n, size) {}
  mpi(const mpi &n, size_type size = 0) : mpi_(n, size) {}
  mpi(const std::string &s, radix_type b = 0) {
    setString(s, b);
  }

  operator int_type() const {
    return toInt();
  }
  operator std::string() const {
    return getString();
  }

  const mpi & operator = (int_type n);
  const mpi & operator = (const mpi &n);
  const mpi & operator = (const std::string &s);

  const mpi & operator <<= (size_type y);
  const mpi & operator >>= (size_type y);

  const mpi & operator &= (const mpi &n);
  const mpi & operator |= (const mpi &n);
  const mpi & operator ^= (const mpi &n);
  mpi operator ~ ();

  mpi & operator ++ ();  // ++var
  const mpi operator ++ (int); //   var++

  mpi & operator -- ();  // --var
  const mpi operator -- (int); //   var--

  const mpi & operator += (const mpi &n);
  const mpi & operator -= (const mpi &n);
  const mpi & operator *= (const mpi &n);
  const mpi & operator /= (const mpi &n);
  const mpi & operator %= (const mpi &n);

  void prime(size_type n);

  bool isPrime(int n = _MP_RM_TESTS);

};

std::istream & operator >> (std::istream &is, mpi &n);

std::ostream & operator << (std::ostream &os, const mpi &n);

bool operator == (const mpi &x, int_type y);
bool operator != (const mpi &x, int_type y);
bool operator <  (const mpi &x, int_type y);
bool operator <= (const mpi &x, int_type y);
bool operator >  (const mpi &x, int_type y);
bool operator >= (const mpi &x, int_type y);

bool operator == (int_type x, const mpi &y);
bool operator != (int_type x, const mpi &y);
bool operator <  (int_type x, const mpi &y);
bool operator <= (int_type x, const mpi &y);
bool operator >  (int_type x, const mpi &y);
bool operator >= (int_type x, const mpi &y);

bool operator == (const mpi &x, const mpi &y);
bool operator != (const mpi &x, const mpi &y);
bool operator <  (const mpi &x, const mpi &y);
bool operator <= (const mpi &x, const mpi &y);
bool operator >  (const mpi &x, const mpi &y);
bool operator >= (const mpi &x, const mpi &y);

mpi operator << (const mpi &x, mpi::size_type y);
mpi operator >> (const mpi &x, mpi::size_type y);

mpi operator & (const mpi &x, const mpi &y);
mpi operator | (const mpi &x, const mpi &y);
mpi operator ^ (const mpi &x, const mpi &y);

mpi operator + (const mpi &x);
mpi operator - (const mpi &x);

mpi operator + (const mpi &x, int_type y);
mpi operator + (int_type x, const mpi &y);
mpi operator + (const mpi &x, const mpi &y);

mpi operator - (const mpi &x, int_type y);
mpi operator - (int_type x, const mpi &y);
mpi operator - (const mpi &x, const mpi &y);

mpi operator * (const mpi &x, int_type y);
mpi operator * (int_type x, const mpi &y);
mpi operator * (const mpi &x, const mpi &y);

mpi operator / (const mpi &x, int_type y);
mpi operator / (int_type x, const mpi &y);
mpi operator / (const mpi &x, const mpi &y);

mpi operator % (const mpi &x, int_type y);
mpi operator % (int_type x, const mpi &y);
mpi operator % (const mpi &x, const mpi &y);

mpi abs(const mpi &n);
mpi sqr(const mpi &n);
mpi pow(const mpi &base, const mpi &p);
mpi pow_mod(const mpi &base, const mpi &p, const mpi &mod);

mpi euclid(const mpi &x, const mpi &y);
mpi factorial(const mpi &n);
mpi gcd(const mpi &x, const mpi &y);

} // End namespace mp
#endif

---