HDOJ 5794 (2016多校联合训练 Training Contest 6) A Simple Chess

题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=5794

第六场多校唯一"做"出来的一题(小编觉得，像 A Boring Question(01)和 A Simple Nim(03)这种打个表莫名其妙就A的题没啥好讲的)。

题意：给我们一个棋盘，棋盘上有r个障碍，我们能够从(1,1)出发，每一次能够向右下走日字，问我们走到终点的方案数。

很简单的一个容斥原理的运用，举个例子，我们从点i出发，途中有一个障碍点j，从点i到终点的方案数，就等于总方案数减去从j到终点的方案数*从点i到点j的方案数。

#include<cstdio>
#include<algorithm>
using namespace std;
typedef long long LL;
const LL mod=110119;
const LL maxn=100;
LL dp[maxn*10+20];
LL f[mod*10];
struct Point{LL x,y;}p[maxn+20];
bool cmp(Point a,Point b)
{
  if(a.x<b.x)return 1;
  if(a.x>b.x)return 0;
  return a.y<b.y;
}
LL Power_mod(LL a, LL b, LL p)
{
    LL res = 1;
    while(b)
	{
        if(b&1) res = (res * a) % p;
        a = (a*a) % p;
        b >>= 1;
    }
    return res;
}
LL Comb(LL a, LL b, LL p)
{
    if(a < b)   return 0;
    if(a == b)  return 1;
    LL ans=f[a]*Power_mod(f[b]*f[a-b]%p,p-2,p)%p;
    return ans;
}
LL Lucas(LL n, LL m, LL p)
{
    LL ans = 1;
    while(n && m && ans)
	{
    	ans = (ans*Comb(n%p, m%p, p)) % p;
        n /= p;
        m /= p;
    }
    return ans;
}
LL get(LL x,LL y)
{
	if((x+y)%3 || max(x,y)>2*min(x,y))return 0;
	LL a=(x+y)/3,b=(2*x-y)/3;
	if(a<0 || b<0)return 0;
	return Lucas(a,b,mod);
}
int main()
{    
	LL i,j,k,h,w,n,cas=0;
	f[0]=f[1]=1;
	for(i=2;i<=mod+20;i++)f[i]=f[i-1]*i%mod; 
	while(~scanf("%I64d%I64d%I64d",&h,&w,&n))
    {
	    for(i=0; i<n; i++) scanf("%I64d%I64d",&p[i].x,&p[i].y);
	    if((h+w-2)%3 || max(h-1,w-1)>2*min(h-1,w-1))
        {
	        printf("Case #%I64d: 0\n",++cas);
	        continue;
        }
	    sort(p,p+n,cmp),p[n].x=h,p[n].y=w;
	    for(i=0; i<=n; i++)
        {
	        dp[i] = get(p[i].x-1,p[i].y-1);
	        for(j=0; j<i; j++)
				if(p[j].x <= p[i].x && p[j].y <= p[i].y)
	            {
		            dp[i] -= get(p[i].x-p[j].x,p[i].y-p[j].y)*dp[j];
		            dp[i] = (dp[i]%mod+mod)%mod;    
	            }  
        }
    	printf("Case #%I64d: %I64d\n",++cas,dp[n]);
    }
	return 0;
}