@@ -0,0 +1,436 @@
+'''%0A%0Afermits little theorem%0AnCr %25 p = (fac%5Bn%5D* modinverse(fac%5Br%5D) %25 p *  modinverse(fac%5Bn-r%5D) %25 p) %25 p%0A%0Amodinverse(x) %25 p = x%5Ep-2%0A%0A'''%0A%0A%0Adef ncr(n, r, p): %0A    numerator = denominator = 1 %0A    for i in range(r): %0A        numerator = (numerator * (n - i)) %25 p %0A        denominator = (denominator * (i + 1)) %25 p     %0A    return (numerator * pow(denominator , p - 2 , p)) %25 p %0A%0A%0An, r, p = 10**5, 20000, 1000000007%0A%0Aprint(ncr(n, r, p)) %0A%0A%0A%0A

Create nCr.py