Test #1:
6
Test #2:
18
Test #3:
180

C++ 代码

#include <iostream>
#include <stdio.h>
#include <string.h>
using namespace std;
#define MAX 100
#define BASE 10000
void multiply(int a[],int Max,int b)  //大数乘小数
{
    int i,array=0;
    for (i=Max-1; i>=0; i--)
    {
        array+=b*a[i];
        a[i] = array%BASE;
        array /= BASE;
    }
}

void divide(int a[], int Max, int b)  //大数除小数
{
    int i,div=0;
    for (i=0;i<Max; i++)
    {
        div = div*BASE + a[i];
        a[i] = div / b;
        div %= b;
    }
}
int fact[205][MAX];
void setFact ()
{
     fact[0][MAX-1] = fact[1][MAX-1] = 1;
     for ( int i = 2; i <= 200; ++ i )
     {
           memcpy ( fact[i] , fact[i-1] , MAX * sizeof ( int ) );
           multiply ( fact[i] , MAX , i );
     }
}
void outPut ( int ctl[MAX] )
{
     int i = 0;
     while ( i < MAX && ctl[i] == 0 )
     {
             i ++ ; //去前导0
     }
     printf ( "%d", ctl[i++] );
     while ( i < MAX )
     {
             printf ( "%04d", ctl[i++] );
     }
     putchar ( '\n' );
}
int res[MAX];
int main ()
{
     int M,N;
     int ca = 1;
     setFact();
     while ( cin >> M >> N , M + N )
     {
             printf ( "Test #%d:\n",ca++ );
             if ( N > M )
             {
                  puts ( "0" );
                  continue;
             }
             memcpy ( res , fact[M+N] , MAX * sizeof ( int ) );     // 阶乘 ( m + n )!
             multiply ( res, MAX, M - N + 1 );                             //  ( m + n )! * ( m-n+1 )
             divide ( res, MAX, M + 1 );                                   //     ( m + n )! * ( m-n+1 )  / ( m+ 1 )
             outPut ( res );
     }
     return 0;
}

本文转自NewPanderKing51CTO博客，原文链接：http://www.cnblogs.com/newpanderking/archive/2011/07/31/2122530.html ，如需转载请自行联系原作者