Added lucas.py, which implements Lucas' theorem for faster combinations