【GDOI2014模拟】网格

题目

某城市的街道呈网格状，左下角坐标为A(0, 0)，右上角坐标为B(n, m)，其中n >= m。现在从A(0, 0)点出发，只能沿着街道向正右方或者正上方行走，且不能经过图示中直线左上方的点，即任何途径的点(x, y)都要满足x >= y，请问在这些前提下，到达B(n, m)有多少种走法。

分析

首先，我们知道：如果现在从(0, 0)点出发，只能沿着街道向正右方或者正上方行走时，到（n,m）点(n>=m)的方案数是\(C^{m}_{n+m}\)。
发现，任何途径的点(x, y)都要满足x>=y就是不经过粉色线x=y+1

有一条违法的路径：

我们把路径按照粉色线做个对称，
我们发现，到橙色点的路径都一一对应每一条违法路径，也就是说，违法路径的方案数就是到橙色点的方案数。
所以，合法方案数=随便走的方案数-违法方案数=\(C^{m}_{n+m}-C^{m-1}_{n+m}\)
注意：如果直接算可能会超时，要把式子简化，而且高精度要压位。

#include <cmath>
#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <algorithm>
#include <queue>
const int maxlongint=2147483647;
const int mo=10000;
using namespace std;
int f[20000],n,m,g[20000];
inline int times(int x)
{
    memcpy(g,f,sizeof(g));
    f[1]=0;
    f[0]=g[0]+4;
    for(int i=1;i<=f[0];i++)
    {
        f[i]+=g[i]*x;
        f[i+1]=f[i]/mo;
        f[i]-=f[i+1]*mo;
    }
    while(!f[f[0]] && f[0]>1)
        f[0]--;
}
inline int div(int x)
{
    memcpy(g,f,sizeof(g));
    int k=0;
    for(int i=g[0];i>=1;i--)
    {
        f[i]=(k*mo+g[i])/x;
        k=(k*mo+g[i])-x*f[i];
    }
    f[0]=g[0];
    while(!f[f[0]] && f[0]>1)
        f[0]--;
}
int main()
{
    scanf("%d%d",&n,&m);
    f[0]=1;
    f[1]=1;
    int j=2;
    for(int i=n+2;i<=n+m;i++)
    {
        if(j<=m)
        {
            if(i%j==0)
            {
                times(i/j);
                j++;
            }
                else
                    times(i);
        }
        
    }
    times(n+1-m);
    for(int i=j;i<=m;i++)
    {
        div(i);
    }
    for(int i=f[0];i>=1;i--)
    {
        if(f[0]!=i)
        printf("%04d",f[i]);
            else
                printf("%d",f[i]);
    }
}