1 | // [The FreeBSD Licence]
|
---|
2 | //
|
---|
3 | // Copyright (c) 2008
|
---|
4 | // Sebastian Mies, Institute of Telematics, UniversitÀt Karlsruhe (TH)
|
---|
5 | // All rights reserved.
|
---|
6 | //
|
---|
7 | // Redistribution and use in source and binary forms, with or without
|
---|
8 | // modification, are permitted provided that the following conditions are
|
---|
9 | // met:
|
---|
10 | //
|
---|
11 | // * Redistributions of source code must retain the above copyright notice,
|
---|
12 | // this list of conditions and the following disclaimer.
|
---|
13 | // * Redistributions in binary form must reproduce the above copyright
|
---|
14 | // notice, this list of conditions and the following disclaimer in the
|
---|
15 | // documentation and/or other materials provided with the distribution.
|
---|
16 | // * Neither the name of the author nor the names of its contributors may be
|
---|
17 | // used to endorse or promote products derived from this software without
|
---|
18 | // specific prior written permission.
|
---|
19 | //
|
---|
20 | // THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS
|
---|
21 | // IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO,
|
---|
22 | // THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR
|
---|
23 | // PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER OR
|
---|
24 | // CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL,
|
---|
25 | // EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO,
|
---|
26 | // PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS;
|
---|
27 | // OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY,
|
---|
28 | // WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR
|
---|
29 | // OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF
|
---|
30 | // ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
|
---|
31 | //
|
---|
32 | #ifndef VINT_HPP_
|
---|
33 | #define VINT_HPP_
|
---|
34 |
|
---|
35 | #include <boost/mpl/if.hpp>
|
---|
36 |
|
---|
37 | #include "detail/helper.hpp"
|
---|
38 | #include "detail/vint_big.hpp"
|
---|
39 | #include "detail/vint_small.hpp"
|
---|
40 |
|
---|
41 | using namespace _vint::detail;
|
---|
42 |
|
---|
43 | template<size_t l, bool s> class vint;
|
---|
44 | template<size_t l, bool s>
|
---|
45 | std::ostream& operator<<(std::ostream&, const vint<l, s>& v);
|
---|
46 |
|
---|
47 | /**
|
---|
48 | * This class implements an adaptive, scalable integer type.<br />
|
---|
49 | *
|
---|
50 | * Key features include:<br />
|
---|
51 | * <ul>
|
---|
52 | * <li>Transparent type switching according to native types</li>
|
---|
53 | * <li>Optimal memory consumption and optimization properties</li>
|
---|
54 | * <li>Automatic big integer arithmetic using gmp or other libraries</li>
|
---|
55 | * <li>Unconventional integer types (e.g. length not a power of 2)</li>
|
---|
56 | * <li>Additional math, like intervals, log2, divrem, gcd, ...</li>
|
---|
57 | * </ul>
|
---|
58 | *
|
---|
59 | * ... short: a real scalable integer type at compile-time cost!<br />
|
---|
60 | *
|
---|
61 | * @author Sebastian Mies <mies@tm.uka.de>
|
---|
62 | */
|
---|
63 | template<size_t __length = 0, bool __sign = false>
|
---|
64 | class vint {
|
---|
65 | // stream operator
|
---|
66 | friend std::ostream& operator<< <__length, __sign> (std::ostream&,
|
---|
67 | const vint<__length, __sign>&);
|
---|
68 |
|
---|
69 | public:
|
---|
70 | // select proper type
|
---|
71 | typedef typename boost::mpl::if_<
|
---|
72 | boost::mpl::bool_<(__length> 64 || __length == 0)>,
|
---|
73 | vint_big<__length, __sign>, vint_small<__length, __sign>
|
---|
74 | >::type type;
|
---|
75 |
|
---|
76 | // selected value type
|
---|
77 | type value;
|
---|
78 |
|
---|
79 | // internal self-template
|
---|
80 | typedef vint<__length, __sign> _self;
|
---|
81 |
|
---|
82 | public:
|
---|
83 | /* construct an zero integer */
|
---|
84 | finline vint() : value() {}
|
---|
85 |
|
---|
86 | /* construct an integer from a constant */
|
---|
87 | template<typename T> finline vint(const T v, if_integral(T) ) : value(v) {}
|
---|
88 |
|
---|
89 | /* construct a variable integer from a given string */
|
---|
90 | finline vint(const char* text, int base = 10) : value(text,base) {}
|
---|
91 |
|
---|
92 | /* copy constructor */
|
---|
93 | finline vint( const type& v ) {value = v;}
|
---|
94 |
|
---|
95 | /* copy constructor */
|
---|
96 | finline vint( const _self& v ) {value = v.value;}
|
---|
97 |
|
---|
98 | //-- conversion -----------------------------------------------------------
|
---|
99 |
|
---|
100 | template<size_t ___length, bool ___sign>
|
---|
101 | finline void convert_to(vint<___length, ___sign>& v ) const {
|
---|
102 | return value.convert_to(v.value);
|
---|
103 | }
|
---|
104 |
|
---|
105 | template<class T>
|
---|
106 | finline void convert_to( T& v, if_integral(T) ) const {
|
---|
107 | return value.convert_to(v);
|
---|
108 | }
|
---|
109 |
|
---|
110 | finline std::string to_string(int base = 10,
|
---|
111 | int str_size = 0, char str_fill = '0') const {
|
---|
112 | return value.to_string(base,str_size,str_fill);
|
---|
113 | }
|
---|
114 |
|
---|
115 | std::string to_debug_string(int base = 10,
|
---|
116 | int str_size = 0, char str_fill = '0') const {
|
---|
117 | return value.to_debug_string(base,str_size,str_fill);
|
---|
118 | }
|
---|
119 |
|
---|
120 | //-- sub integers ---------------------------------------------------------
|
---|
121 |
|
---|
122 | finline bool get_bit(size_t index) const {
|
---|
123 | return value.get_bit(index);
|
---|
124 | }
|
---|
125 |
|
---|
126 | finline void set_bit(size_t index, bool v) {
|
---|
127 | value.set_bit(index, v);
|
---|
128 | }
|
---|
129 |
|
---|
130 | finline void set_subint(_self& v, size_t index) {
|
---|
131 | value.set_subint(v.value,index);
|
---|
132 | }
|
---|
133 |
|
---|
134 | finline void set_subint(uintmax_t v, size_t index) {
|
---|
135 | value.set_subint(v,index);
|
---|
136 | }
|
---|
137 |
|
---|
138 | finline _self get_subint(size_t index, size_t length) const {
|
---|
139 | return value.get_subint(index,length);
|
---|
140 | }
|
---|
141 |
|
---|
142 | finline uintmax_t get_subint(size_t index) const {
|
---|
143 | return value.get_subint(index);
|
---|
144 | }
|
---|
145 |
|
---|
146 | //-- assignment -----------------------------------------------------------
|
---|
147 |
|
---|
148 | public:
|
---|
149 |
|
---|
150 | template<size_t _len, bool _sign>
|
---|
151 | finline void assign( const vint<_len,_sign>& v ) {
|
---|
152 | value.assign(v.value);
|
---|
153 | }
|
---|
154 |
|
---|
155 | template<typename T>
|
---|
156 | finline void assign(const T& v, if_integral(T) ) {
|
---|
157 | value.assign(v);
|
---|
158 | }
|
---|
159 |
|
---|
160 | finline void assign(const std::string& text, int base = 10) {
|
---|
161 | value.assign( text, base );
|
---|
162 | }
|
---|
163 |
|
---|
164 | void assign(const char*& text, int base = 10) {
|
---|
165 | assign(std::string(text), base);
|
---|
166 | }
|
---|
167 |
|
---|
168 | //-- compare operations ---------------------------------------------------
|
---|
169 |
|
---|
170 | finline int compare_to(const _self& v) const {
|
---|
171 | return value.compare_to(v.value);
|
---|
172 | }
|
---|
173 |
|
---|
174 | template<class T>
|
---|
175 | finline int compare_to(const T& v, if_integral(T) ) const {
|
---|
176 | return value.compare_to(v);
|
---|
177 | }
|
---|
178 |
|
---|
179 | //-- arithmetic operations ------------------------------------------------
|
---|
180 |
|
---|
181 | // addition
|
---|
182 | template<class T>
|
---|
183 | finline void add(const T& rhs, if_integral(T) ) {value.add(rhs.value);}
|
---|
184 | finline void add(const _self& rhs) {value.add(rhs.value);}
|
---|
185 |
|
---|
186 | // subtraction
|
---|
187 | template<class T>
|
---|
188 | finline void sub(const T& rhs, if_integral(T) ) {value.sub(rhs.value);}
|
---|
189 | finline void sub(const _self& rhs) {value.sub(rhs.value);}
|
---|
190 |
|
---|
191 | // multiplication
|
---|
192 | template<class T>
|
---|
193 | finline void mul( const T& rhs, if_integral(T) ) {value.mul(rhs);}
|
---|
194 | finline void mul( _self& rhs) {value.mul(rhs.value);}
|
---|
195 | finline void mul( _self& res, const _self& rhs ) const {value.mul(res.value,rhs.value);}
|
---|
196 |
|
---|
197 | // division
|
---|
198 | template<class T>
|
---|
199 | finline void div( const T& rhs, if_integral(T) ) {value.div(rhs);}
|
---|
200 | finline void div(_self& rhs) {value.div(rhs.value);}
|
---|
201 | finline void div(_self& res, const _self& rhs ) const {value.div(res.value, rhs);}
|
---|
202 |
|
---|
203 | // modulo
|
---|
204 | template<class T>
|
---|
205 | finline void mod( const T& rhs, if_integral(T) ) {value.mod(rhs);}
|
---|
206 | finline void mod(_self& rhs) {value.mod(rhs.value);}
|
---|
207 | finline void mod(_self& res, const _self& rhs ) const {value.mod(res.value, rhs);}
|
---|
208 |
|
---|
209 | // exception of modulo: result may be of integral type
|
---|
210 | template<class T>
|
---|
211 | finline void mod( T& res, const T& rhs, if_integral(T) ) const {value.mod(res,rhs);}
|
---|
212 |
|
---|
213 | // division with remainder
|
---|
214 | template<class T>
|
---|
215 | finline void divrem(_self& res_div, T& res_rem, const T& rhs, if_integral(T) ) const {
|
---|
216 | value.divrem(res_div.value, res_rem, rhs );
|
---|
217 | }
|
---|
218 | finline void divrem(_self& res_div, _self& res_rem, const _self& rhs ) const {
|
---|
219 | value.divrem(res_div.value, res_rem.value, rhs.value);
|
---|
220 | }
|
---|
221 | template<class T>
|
---|
222 | finline void divrem( T& res_rem, const T& rhs, if_integral(T) ) {
|
---|
223 | value.div(res_rem, rhs.value);
|
---|
224 | }
|
---|
225 | finline void divrem( _self& res_rem, const _self& rhs ) {
|
---|
226 | value.div(res_rem.value, rhs.value);
|
---|
227 | }
|
---|
228 |
|
---|
229 | // logarithm to the base of 2
|
---|
230 | finline uintmax_t log2() const {return value.log2();}
|
---|
231 |
|
---|
232 | //-- logical bit operations -----------------------------------------------
|
---|
233 |
|
---|
234 | finline void or_ (const _self& v) {value.or_(v.value);}
|
---|
235 | finline void and_(const _self& v) {value.and_(v.value);}
|
---|
236 | finline void xor_(const _self& v) {value.xor_(v.value);}
|
---|
237 | finline void complement() {value.complement();}
|
---|
238 | finline void lshift(size_t steps) {value.lshift(steps);}
|
---|
239 | finline void rshift(size_t steps) {value.rshift(steps);}
|
---|
240 |
|
---|
241 | //-- integer dynamic ------------------------------------------------------
|
---|
242 |
|
---|
243 | inline _self normalized() const {return value.normalized();}
|
---|
244 |
|
---|
245 | //-- general information --------------------------------------------------
|
---|
246 |
|
---|
247 | finline size_t length() const {return value.length();}
|
---|
248 | finline void set_length(size_t length) {value.set_length(length);}
|
---|
249 | finline bool is_unspecified() {return value.is_unspecified();}
|
---|
250 |
|
---|
251 | finline _self max() const {_self v; v.value = value.max(); return v;}
|
---|
252 | finline _self min() const {_self v; v.value = value.min(); return v;}
|
---|
253 | finline _self zero() const {_self v; v.value = value.zero(); return v;}
|
---|
254 |
|
---|
255 | //-- general extensions ---------------------------------------------------
|
---|
256 |
|
---|
257 | finline void convert_to( std::string& str ) const {
|
---|
258 | str = to_string();
|
---|
259 | }
|
---|
260 |
|
---|
261 | /**
|
---|
262 | * Returns true, if the integer is inside the interval [a,b]
|
---|
263 | */
|
---|
264 | finline bool is_between(const _self& a, const _self& b) const {
|
---|
265 | return (a <= b) ?
|
---|
266 | ((*this >= a) && (*this <= b))
|
---|
267 | :
|
---|
268 | ((*this <= a) && (*this >= b))
|
---|
269 | ;
|
---|
270 | }
|
---|
271 |
|
---|
272 | /**
|
---|
273 | * Returns true, if the integer is inside the interval [a,b] on
|
---|
274 | * a ring.
|
---|
275 | */
|
---|
276 | finline bool is_between_ring(const _self& a, const _self& b) const {
|
---|
277 | return (b<a) ?
|
---|
278 | (is_between(a, a.max())||is_between(b.zero(), b))
|
---|
279 | :
|
---|
280 | ((*this >= a) && (*this <= b))
|
---|
281 | ;
|
---|
282 | }
|
---|
283 |
|
---|
284 | //-- convenience operator overloads ---------------------------------------
|
---|
285 |
|
---|
286 | // conversion operators
|
---|
287 | template<typename T> finline T convert() const {T v;convert_to(v);return v;}
|
---|
288 | template<typename T> finline operator T () const {return convert<T>();}
|
---|
289 |
|
---|
290 | // assign operators
|
---|
291 | template<class T>
|
---|
292 | finline _self& operator=(const T& copy) {assign(copy);return *this;}
|
---|
293 |
|
---|
294 | // add operators
|
---|
295 | template<class T>
|
---|
296 | finline _self& operator+=(const T& v) {add(v);return *this;}
|
---|
297 | template<class T>
|
---|
298 | finline _self& operator++() {add(1);return *this;}
|
---|
299 | template<class T>
|
---|
300 | finline _self operator+ (const T& v) const {_self x=*this;x.add(v);return x;}
|
---|
301 |
|
---|
302 | // substract operators
|
---|
303 | template<class T>
|
---|
304 | finline _self& operator-=(const T& v) {sub_model(v);return *this;}
|
---|
305 | finline _self& operator--() {sub(1);return *this;}
|
---|
306 | template<class T>
|
---|
307 | finline _self operator- (const T& v) const {
|
---|
308 | _self x(*this); x.sub_model(v); return x;
|
---|
309 | }
|
---|
310 |
|
---|
311 | // multiply operators
|
---|
312 | template<class T> finline _self& operator*=(const T& v) {
|
---|
313 | value.mul(v);
|
---|
314 | return *this;
|
---|
315 | }
|
---|
316 | template<class T> finline _self operator*(const T& v) const {
|
---|
317 | _self res = *this;res.mul(v); return res;
|
---|
318 | }
|
---|
319 |
|
---|
320 | // bit shift operators convenience
|
---|
321 | finline _self operator<< (size_t steps) {_self x=*this;return(x<<=steps);}
|
---|
322 | finline _self& operator>>=(size_t steps) {rshift(steps);return *this;}
|
---|
323 | finline _self operator>> (size_t steps) {_self x=*this;return(x<<=steps);}
|
---|
324 | finline _self& operator<<=(size_t steps) {lshift(steps);return *this;}
|
---|
325 |
|
---|
326 | // bit operators convenience
|
---|
327 | finline _self& operator&=(const _self& v) {and_(v); return *this;}
|
---|
328 | finline _self operator & (const _self& v) {_self x=*this;return(x&=v);}
|
---|
329 | finline _self& operator|=(const _self& v) {or_(v); return *this;}
|
---|
330 | finline _self operator | (const _self& v) {_self x=*this;return(x|=v);}
|
---|
331 | finline _self& operator^=(const _self& v) {xor_(v);return *this;}
|
---|
332 | finline _self operator^ (const _self& v) {_self x=*this;return(x ^= v);}
|
---|
333 | finline _self operator~() {_self x=*this;return x.complement();}
|
---|
334 |
|
---|
335 | // compare operators convenience
|
---|
336 | template<class T>
|
---|
337 | finline bool operator< (const T& v) const {return compare_to(v)< 0;}
|
---|
338 | template<class T>
|
---|
339 | finline bool operator> (const T& v) const {return compare_to(v)> 0;}
|
---|
340 | template<class T>
|
---|
341 | finline bool operator<=(const T& v) const {return compare_to(v)<=0;}
|
---|
342 | template<class T>
|
---|
343 | finline bool operator>=(const T& v) const {return compare_to(v)>=0;}
|
---|
344 | template<class T>
|
---|
345 | finline bool operator==(const T& v) const {return compare_to(v)==0;}
|
---|
346 | template<class T>
|
---|
347 | finline bool operator!=(const T& v) const {return compare_to(v)!=0;}
|
---|
348 | };
|
---|
349 |
|
---|
350 | template<size_t l, bool s>
|
---|
351 | std::ostream& operator<<(std::ostream& os, const vint<l, s>& v) {
|
---|
352 | return os << v.value;
|
---|
353 | }
|
---|
354 |
|
---|
355 | #endif /* VINT_HPP_ */
|
---|
356 |
|
---|