hdu 1267 下沙的沙子有几粒？

   摘要: f(m,n) =( ((m+n)!)/((m!)*n!) )*(1-n/(m+1))

 

 

#include<iostream>
using namespace std;
int main()
{
    int i,j,len=1,r=0,temp=0;
    int a[201][400]={0};
    a[1][0]=1;
    for(i=2;i<=200;i++)//求1……200的阶乘
    {
        for(j=0;j<len;j++)
            a[i][j] = a[i-1][j] * i;
        for(j=0,r=0;j<len;j++)
        {
            temp = a[i][j] + r;
            a[i][j] = temp % 10;
            r = temp /10;
        }
        while(r)
        {
            a[i][len++] = r%10;
            r /= 10;
        }
    }
    int m,n;
    while(cin>>m>>n)
    {
        r = 0;
        len = 400;//改成300就错了
        int b[500]={0};
        for(j=0;j<len;j++)
        {
            temp = a[m+n][j]*(m+1-n) + r;
            b[j] = temp % 10;
            r = temp /10;
        }
        while(r)
        {
            b[len++] = r%10;
            r /= 10;
        }

        for(j=len-1,r=0;j>=0;j--)
        {            //除法从高位到低位
            temp=r*10+b[j];
            b[j]=temp/(m+1);
            r=temp%(m+1);
        }
        while(!b[len])        //处理高位的零位
            len--;
            
        for(i=m;i>1;i--)//除m!
        {    
            for(j=len,r=0;j>=0;j--)  
                {            //除法从高位到低位
                temp=r*10+b[j];
                b[j]=temp/i;
                r=temp%i;
                }
            while(!b[len])        //处理高位的零位
            len--;
        }
        
        for(i=n;i>1;i--)//除n!
        {    
            for(j=len,r=0;j>=0;j--)  
                {            //除法从高位到低位
                temp=r*10+b[j];
                b[j]=temp/i;
                r=temp%i;
                }
            while(!b[len])        //处理高位的零位
            len--;
        }
                
        for(j=len;j>=0;j--)
            cout<<b[j];
        cout<<endl;
    }
    return 0;
}