hdu5794A Simple Chess(lucas定理+dp)

题目链接：点这里！！！

题意：

给你一个n*m的方格(1<=n,m<=10^18),方格上面有r(0<=r<=100)个障碍物，问你从1,1走到n,m有多少种方案？走的时候只能走日。

(例如x1,y1->x2,y2满足(x1-x2)*(x1-x2)+(y1-y2)*(y1-y2)=5)

题解：

1、我们可以能够看出他走的其实是杨辉三角，我们可以利用组合数来求解，又因为C(x,y)%p，x,y非常大，我们就用lucas定理来求解组合数。

2、我们从1,1走到n,m。我们先将障碍物排好序，我们假设dp[i]为从1,1走到i节点的步数。我们假设dp[0]=1,dp[r+1]是我们要求解的答案(r+1节点是n,m)。然后利用容斥的原理得到dp[i]=cel(0,i)-dp[j]*cel(j,i)(1<=j<i)。

代码：

#include<cstdio>
#include<cstring>
#include<iostream>
#include<sstream>
#include<algorithm>
#include<vector>
#include<bitset>
#include<set>
#include<queue>
#include<stack>
#include<map>
#include<cstdlib>
#include<cmath>
#define LL long long
#define pb push_back
#define pa pair<int,int>
#define clr(a,b) memset(a,b,sizeof(a))
#define lson lr<<1,l,mid
#define rson lr<<1|1,mid+1,r
#define bug(x) printf("%d++++++++++++++++++++%d\n",x,x)
#define key_value ch[ch[root][1]][0]
#pragma comment(linker, "/STACK:102400000000,102400000000")
const LL  MOD = 110119;
const int N = 2e5+15;
const int maxn = 1e5+15;
const int letter = 130;
const LL INF = 1e18;
const double pi=acos(-1.0);
const double eps=1e-6;
using namespace std;
inline int read()
{
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
    return x*f;
}
LL F[N],Fn[N],inv[N];
LL n,m;
LL f[110];
int r;
struct node{
    LL x,y;
    bool operator < (const node &p)const {
        if(x!=p.x) return x<p.x;
        return y<p.y;
    }
}p[110];
inline void init(LL p){
    F[1]=F[0]=inv[1]=Fn[0]=Fn[1]=1;
    for(LL i =2;i <= p;i++){
        F[i] = F[i-1]*1ll*i % (MOD);
        inv[i]=inv[MOD%i]*(MOD-MOD/i)%MOD;
        Fn[i]=Fn[i-1]*inv[i]%MOD;
    }
}
inline LL Lucas(LL n,LL m,LL p){
    LL ans = 1;
    while(n&&m){
        LL a = n%p,b=m%p;
        if(a<b)return 0;
        ans=ans*F[a]%p*Fn[b]%MOD*Fn[a-b]%p;
        n /= p;
        m /= p;
    }
    return ans;
}
inline LL cel(int a,int b){
    LL n=p[b].x-p[a].x,m=p[b].y-p[a].y;
    if((n+m)%3ll!=0) return 0ll;
    LL step=(m+n)/3ll;
    n-=step,m-=step;
    if(n<0||m<0) return 0ll;
    return Lucas(n+m,m,MOD);
}
inline bool check(int a,int b){
    if(p[a].x>p[b].x||p[a].y>p[b].y) return 0;
    return 1;
}
int main(){
    init(MOD);
    int cas=0;
        while(~scanf("%I64d%I64d%d",&n,&m,&r)){
            for(int i=1;i<=r;i++){
                scanf("%I64d%I64d",&p[i].x,&p[i].y);
            }
            p[0].x=1,p[0].y=1;
            p[++r].x=n,p[r].y=m;
            sort(p,p+r+1);
            for(int i=1;i<=r;i++) f[i]=0;
            f[0]=1;
            for(int i=1;i<=r;i++){
                f[i]=cel(0,i);
                for(int j=1;j<i;j++)
                    if(check(j,i)) f[i]=((f[i]-f[j]*cel(j,i)%MOD)%MOD+MOD)%MOD;
            }
            printf("Case #%d: %I64d\n",++cas,(f[r]%MOD+MOD)%MOD);
        }
        return 0;
}