1164_上班路线

上班路线

• Description

  一位秘书在距离家以东n个街区，以北m个街区的一座大楼里工作。每天他要步行5个街区上班。（右图为n=3，m=2时的情形）。对于他来说可能有多少种不同的上班路线f(n,m)（上班是件敢时间的事情，他的行进速度保持不变，他需要尽可能早到达公司）

• Input

输入包括一行，两个正整数n，m ( 0 < n , m <= 50 )。

• Output

输出包括一行，一个数字f(n,m) mod 9901 。

• Sample Input

3 2

• Sample Output

10

• Hint

f(n,m) mod 9901 代表的意思是 f(n,m) 除以 9901 所得的余数

#include<iostream>
using namespace std;
#define N 100
#define M 9901
int prime[N];
int tab[N];
int total;
int num[N];

void prepare()
{
    int i,j;
    total=0;
    memset(tab,0,sizeof(tab));
    for(i=2;i<N;i++)
    {
        if(tab[i]==0)
        {
            prime[total++]=i;
            for(j=2*i;j<N;j+=i)
                tab[j]=1;
        }
    }
}

void fast(int n,int flag)
{
    int i;
    for(i=0;i<total&&prime[i]<=n;i++)
    {
        int t=n;
        while(t)
        {
            num[i]+=flag*t/prime[i];
            t/=prime[i];
        }
    }
}

int mod(int a,int b)
{
    int t=a,res=1;
    while(b)
    {
        if(b&1)res=res*t%M;
        t=t*t%M;
        b>>=1;//等同于b=b>>1,即b右移一位,最前面补0
    }
    return res;
}
int solve(int n,int m)
{
    int res=1,i;
    memset(num,0,sizeof(num));
    fast(n+m,1);
    fast(m,-1);
    fast(n,-1);
    for(i=0;i<total&&prime[i]<=n+m;i++)
    {
        res=res*mod(prime[i],num[i]);
        res=res%M;
        if(res==0)
            return 0;
    }
    return res;
}

int main()
{
    int n,m;
    prepare();
    while(scanf("%d%d",&n,&m)!=EOF)
        printf("%d\n",solve(n,m));
    return 0;
}