HDU 5226 Tom and matrix

此题的关键在于一个神奇的结论 即一列中组合数的和 可以转化为两个数的差 (具体公式去看官方题解吧 公式太难打)  

而这个神奇的结论还是很好证明的 按照C（n,m）的定义等于n!/m!/(n-m)! 2边展开 就可以得出相等 而比赛时我也没有想出来 (＞﹏＜)  （感觉很难想到）

有了这个结论 再套一个卢卡斯的版子即可 （n,m是10的5次方 感觉不用卢卡斯也一样的）

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cctype>
#include<cmath>
#include<vector>
#include<queue>
#include<map>
#include<algorithm>
#include<set>
#define scnaf scanf
#define cahr char
#define bug puts("bugbugbug");
using namespace std;
typedef long long ll;
const int maxn=1e5+15;
const int inf=1e9;
int mod;
ll powMod(ll a, ll b){
    ll ans =1;
    for( a%=mod; b; b>>=1, a = a * a % mod)
        if(b&1)  ans = ans * a % mod;
    return  ans;
}
ll fac[maxn],inv[maxn];
void init(){
    int len=min(100010,mod);
    fac[0] = 1;
    for(int i = 1; i < len; i++)
        fac[i] = fac[i-1] * i % mod;
    inv[len - 1]=powMod( fac[len - 1], mod - 2);
    for(int i = len-2; i >= 0 ; i--)
        inv[i] = inv[i+1] * (i+1) % mod;
}
ll C(ll n, ll m){
    return  m > n ? 0 : fac[n] * inv[m] % mod * inv[n-m] % mod;
}
ll lucas(ll n,ll m){
   return m? (C(n%mod, m%mod) * lucas(n/mod, m/mod)) % mod : 1;
}
int main()
{
    int x1,x2,y1,y2;
    while(~scanf("%d%d%d%d%d",&x1,&y1,&x2,&y2,&mod))
    {
        init();
        ll ans=0;
        for(int i=y1;i<=y2;i++){
            ans=((ans+lucas(x2+1,i+1)-lucas(x1,i+1))%mod+mod)%mod;
        }
        printf("%I64d\n",ans);
    }
}