怎样学好哲学（lucas+费马小定理）

怎样学习哲学

时间限制: 1 Sec  内存限制: 128 MB
提交: 97  解决: 27
[提交][状态][讨论版]

题目描述

OI大师抖儿在夺得银牌之后，顺利保送pku。这一天，抖儿问长者：“虽然我已经保送了，但是我还要参加学考。马上就要考政治了，请问应该怎样学习哲学，通过政治考试？”
 长者回答：“你啊，Too Young Too Simple,Sometimes Naive!哲学这种东西，不是说想懂就能懂的，需要静心撕烤。你去后面的森林里好好想想。”

长者的后院有一片哲♂学森林。由于一些奥妙重重的原因，这片森林构成了一个n*m的矩形，其中每个点就代表了一棵树。此外，由于辣鸡出题人KJDH从中捣鬼，有些树被连根拔起（也就是消失了）。抖儿每天都要到树下撕烤，因此他想要在每一行选择一棵树。但是他非常讨厌走回头路，因此第i行选择的树必须比第i-1行的靠右。现在抖儿想知道，总共有多少种选择的方案。

输入

第一行三个整数n,m,p，分别表示森林的长、宽，以及消失的树的数目。

接下来p行每行两个整数，表示第ai行第bi列的树消失了。

输出

一行一个整数，表示方案数。由于答案可能很大，请对1000003取模。

样例输入

3 5 2
2 3
3 4

样例输出

5

提示

【样例说明】

方案一：选（1,1）（2,2）（3,3）

方案二：选（1,1）（2,2）（3,5）

方案三：选（1,1）（2,4）（3,5）

方案四：选（1,2）（2,4）（3,5）

方案五：选（1,3）（2,4）（3,5）

 

题解，可以将其看成三角形的一个类似的，走法问题，就是半三角形走法，然后就是发现方案数是C（n，m），这个是可以推出来，

然后就是dp，当前节点的方案数总，是它左上部分经过不合法点到达其的方案数之和，相减即为走到该点方案数。

这样可以证明，到该点的方案数是所有，因为任何经过左上的dp[i]方案中，是表示到达dp[i]的合法方案数，因此通过数学归纳法得证，

这个推断是正确的，为了简便，将n+1，m+1这棵树拔掉，然后这个点的方案数，就为结果了。

 #include<cstdio>
 #include<algorithm>
 #include<cstring>
 #include<cmath>
 #include<iostream>
 #define mod 1000003
 #define ll long long
 #define Q 2007
 using namespace std;
 
 int n,m,q;
 ll p[mod+],inv[mod+],dp[Q];
 struct Node
 {
     int x,y;
 }a[Q];
 
 bool cmp(Node x,Node y)
 {
     return x.x<y.x;
 }
 ll ksm(ll a,ll b)
 {
     ll ans=;
     while (b)
     {
         if (b&) ans=a*ans%mod;
         b/=;
         a=a*a%mod;
     }
     return ans;
 }
 ll Lucas_C(int n,int m)
 {
     if (n<m) return ;
     if (m==) return ;
     if (n==m) return ;
     if (n<mod) return p[n]*inv[m]%mod*inv[n-m]%mod;
     else return Lucas_C(n%mod,m%mod)*Lucas_C(n/mod,m/mod)%mod;
 }
 int main()
 {
     p[]=;
     for (int i=;i<=mod;i++)
         p[i]=(p[i-]*i)%mod;
     for (int i=;i<=mod;i++)
         inv[i]=ksm(p[i],mod-);
     scanf("%d%d%d",&n,&m,&q);
 
     for (int i=;i<=q;i++)
         scanf("%d%d",&a[i].x,&a[i].y);
     q++,a[q].x=n+,a[q].y=m+;
     sort(a+,a+q+,cmp);
     for (int i=;i<=q;i++)
     {
         dp[i]=Lucas_C(a[i].y-,a[i].x-);
         for (int j=;j<i;j++)
             if (a[i].x>a[j].x&&a[i].y>a[j].y)
                 dp[i]=(dp[i]-dp[j]*Lucas_C(a[i].y-a[j].y-,a[i].x-a[j].x-)%mod+mod)%mod;
     }
     printf("%lld",dp[q]);
 }