3

%
�a��@s�dZd
dl
Z
Gdd�de�Zeed�dd�Zeed�d	d
�Zeeed�dd
�Zeee
jeeefd�dd�Z	eeed�dd�Z
e
jee
jeed�dd�Ze
dkr�d
dlZej�
dS)z/Common functionality shared by several modules.�Ncs,eZ
dZdeeeedd��f
dd�
Z�ZS)�NotRelativePrimeError�N)�
a�
b�
d�msg�returncs0t�j
|pd
|
||f�

|
|_||_||_dS)Nz.%d and %d are not relatively prime, divider=%i)�super�__init__rrr)�selfrrrr)
�	__class__��/usr/lib/python3.6/common.pyr
s



zNotRelativePrimeError.__init__)
r)�__name__�
__module__�__qualname__�int�strr
�
__classcell__r
r
)
rrrs
r)�numrc
CsBy|j�St
k
r<
}

ztd
t|�
�
|
�WYdd}
~
XnXdS)a�

    Number of bits needed to represent a integer excluding any prefix
    0 bits.

    Usage::

        >>> bit_size(1023)
        10
        >>> bit_size(1024)
        11
        >>> bit_size(1025)
        11

    :param num:
        Integer value. If num is 0, returns 0. Only the absolute value of the
        number is considered. Therefore, signed integers will be abs(num)
        before the number's bit length is determined.
    :returns:
        Returns the number of bits in the integer.
    z,bit_size(num) only supports integers, not %rN)�
bit_length�AttributeError�	TypeError�type)rZexr
r
r�bit_sizes


r)�numberrc

Cs|d
krdStt
|�
d�S)a�

    Returns the number of bytes required to hold a specific long number.

    The number of bytes is rounded up.

    Usage::

        >>> byte_size(1 << 1023)
        128
        >>> byte_size((1 << 1024) - 1)
        128
        >>> byte_size(1 << 1024)
        129

    :param number:
        An unsigned integer
    :returns:
        The number of bytes required to hold a specific long number.
    r
�
�)�ceil_divr)
rr
r
r�	byte_size8s

r)r�divrcCst||
�\}}|r|d
7}|S)av

    Returns the ceiling function of a division between `num` and `div`.

    Usage::

        >>> ceil_div(100, 7)
        15
        >>> ceil_div(100, 10)
        10
        >>> ceil_div(1, 4)
        1

    :param num: Division's numerator, a number
    :param div: Division's divisor, a number

    :return: Rounded up result of the division between the parameters.
    r)
�divmod)rr Zquanta�modr
r
rrQs


r)rrrc	Cs�d
}d}d}d
}|}|
}xF|
d
kr^||
}|
||
}}
||||}}||||}}qW|d
krp||7}|d
kr�||7}|||fS)z;Returns a tuple (r, i, j) such that r = gcd(a, b) = ia + jbr
rr
)	rr�
x�
yZlxZlyZoa�ob�
qr
r
r�extended_gcdis 















r')r#�
nrcCs(t||
�\}}}|d
kr$t
||
|��
|S)z�Returns the inverse of x % n under multiplication, a.k.a x^-1 (mod n)

    >>> inverse(7, 4)
    3
    >>> (inverse(143, 4) * 143) % 4
    1
    r)r'r)r#r(Zdivider�inv�
_r
r
r�inverse�s	
r+)�a_values�
modulo_valuesrc	Cs`d
}d}x|
D]}||9}qWx<t|
|�D].\}}||}t
||�}|||||}q*W|S)a�
Chinese Remainder Theorem.

    Calculates x such that x = a[i] (mod m[i]) for each i.

    :param a_values: the a-values of the above equation
    :param modulo_values: the m-values of the above equation
    :returns: x such that x = a[i] (mod m[i]) for each i


    >>> crt([2, 3], [3, 5])
    8

    >>> crt([2, 3, 2], [3, 5, 7])
    23

    >>> crt([2, 3, 0], [7, 11, 15])
    135
    rr
)�zipr+)	r,r-�
mr#ZmoduloZm_iZa_iZM_ir)r
r
r�crt�s





r0�__main__)�__doc__Ztyping�
ValueErrorrrrrrZTupler'r+�Iterabler0rZdoctestZtestmodr
r
r
r�<module>s#