hdu 5194 DZY Loves Balls(组合数学)

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

题目大意：给你n个黑球和m个白球，每次取一个球，黑球代表1，白球代表0，求p/q，q表示不同字符串的个数，p表示01子串的总个数。q比较好求，q=（m+n）的全排列除以（m的全排列乘n的全排列），直接算（m+n）的全排列会超出longlong 的范围，m+n最大是24，所以需要一点小技巧，具体看代码就懂，（m+n）的全排列中包含着m或n的全排列，所以不用乘那么多次。求p的话，先把所有的1排成一排，然后依次枚举有i（1<=i<=min(m,n)）个01子串的字符串，i个01串只需要i个0，如果0还有剩余，还需将这些0插入到串中，但是前面没有0的1不能插0，所以实在i+1地方插入m-i个0，需要采用dp的方法预处理一下，dp[i][j]就表示在i个位置中插入j个0有多少种可能，初始化dp数组为0，dp[1][i]和dp[i][1]=1。然后就能求得答案。

这题用dfs是会超时的，应为大概有150个例子，也没有什么好的剪枝方法，但是可以先写个dfs作为对拍程序知道一些数据的正确答案，比如12 12输出6/1。

AC代码：

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<queue>
#include<set>
#include<stack>
using namespace std;
long long c(long long x,long long y)
{
    long long i,r,j;
    x=min(x,y-x);
    if(x==0)
    {
        return 1;
    }
    r=1;
    for(i=1,j=y;i<=x;i++,j--)
    {
        r*=j;
    }
    for(i=x;i>=1;i--)
    {
        r/=i;
    }
    return r;
}
long long f(long long da,long long xiao)
{
    long long t;
    while(da%xiao!=0)
    {
        t=da;
        da=xiao;
        xiao=t%xiao;
    }
    return xiao;
}
long long dp[15][15];
int main()
{
	long long n,m,z,az,an,am,q,i,p,sum,y,j,sw,x,k;
	memset(dp,0,sizeof(dp));
	for(i=0;i<=15;i++)
	{
	    dp[1][i]=1;
	    dp[i][0]=1;
	}
	for(i=2;i<=13;i++)
	{
	    for(j=1;j<=12;j++)
	    {
	        for(k=0;k<=j;k++)
	        {
	            dp[i][j]+=dp[i-1][j-k];
	        }
	    }
	}
	while(cin>>n>>m)
	{
	    z=m+n;
	    az=an=am=1;
	    for(i=z;i>n;i--)
	    {
	        az*=i;
	    }
	    for(i=1;i<=m;i++)
	    {
	        am*=i;
	    }
	    q=az/(am);
	    //cout<<"q:"<<q<<endl;
	    p=0;
	    for(i=1;i<=n&&i<=m;i++)
	    {
	        sw=m-i;
	        sum=dp[i+1][sw];
	        p+=sum*c(i,n)*i;
	        //cout<<"p:"<<p<<endl;
	    }
	    //cout<<"p:"<<p<<endl;
	    y=f(max(p,q),min(p,q));
	    cout<<p/y<<"/"<<q/y<<endl;
	}
	return 0;
}