Lightoj 1067 - Combinations

http://lightoj.com/volume_showproblem.php?problem=1067

1. 
   
2. #include <bits/stdc++.h>
3. using namespace std;
4. #define ll long long int
5. #define M 1000003
6. ll fact[M];
7. ll bigmod(ll p, ll q){
8.     if(q == 1) return p%M;
9.     if(q%2 == 0){
10.         ll a = bigmod(p, q/2);
11.         return ((a%M)*(a%M))%M;
12.     }
13.     else{
14.         ll a = bigmod(p, q-1);
15.         return (a*p)%M;
16.     }
17. }
18. void pre(){
19.     fact[0] = 1;
20.     for(int i = 1; i <M; i++){
21.         fact[i] = (i*(fact[i-1]%M))%M;
22.         //if(fact[i]<=0)printf(" fact[%d] =  %lld \n", i, fact[i]);
23.     }
24. }
25. int main(){
26.     pre();
27.     int ts;
28.     cin>>ts;
29.     for(int k = 1; k <= ts; k++){
30.         int n, r;
31.         scanf("%d%d", &n, &r);
32.         if(n<r){
33.             printf("Case %d: 0\n", k);
34.             continue;
35.         }
36.         else if(n == r || r == 0){
37.             printf("Case %d: 1\n", k);
38.             continue;
39.         }
40.         ll s = bigmod((fact[n-r]*fact[r])%M,M-2)%M;
41.         ll p = (fact[n]*s)%M;
42.         printf("Case %d: %lld\n", k,p);
43.     }
44.     return 0;
45. }