best coder #35-01<组合数学 || 概率数学>

问题描述
一个盒子里有n个黑球和m个白球。现在DZY每次随机从盒子里取走一个球,取了n+m次后,刚好取完。DZY用这种奇怪的方法生成了一个随机的01串S[1(n+m)]。如果DZY第i次取出的球是黑色的,那么S[i]=1,如果是白色的,那么S[i]=0。
DZY现在想知道,'01'在S串中出现的期望次数。
输入描述
输入有多组测试数据。 (TestCase150)
每行两个整数, n, m(1n,m12)
输出描述
对于每个测试数据,输出一行答案,格式为p/q(p,q互质)。
输入样例
1 1
2 3
输出样例
1/2
6/5
Hint
Case 1: S='01' or S='10', 所以期望次数 = 1/2
Case 2: S='00011' or S='00101' or S='00110' or S='01001' or S='01010' 
or S='01100' or S='10001' or S='10010' or S='10100' or S='11000',
所以期望次数 = (1+2+1+2+2+1+1+1+1+0)/10 = 12/10 = 6/5
*************************************************************************************************
本题的官方解释是数学:
第i位上出现0,第i+1位出现1,这种概率是:(m/(n+m))*(n/(n+m-1)),0的位置可以出现在第一个一直到倒数第二个;所以累加(n+m-1),相承,公式是n*m/(n+m);

我的做法不是这样,组合数学:
刚开始看作所有的1发个在最左边,0放在最右边,一共最多有n或m个01串,(个数是其中较小的那个),枚举串的个数从1到最大,当有1
个01串的时候,从n个1中取一个空当,用来插入0,从m个0中可以选1个0放到空档里,也可以选两个一直到选m个,所以对于每一个1都有m种插法,c(n,1)*c(m,1)*1;
当有两个01串的时候,取1,c(n,2),取0,第一个1有1到m-1种插法,第二个1有剩余的插法,组合数公式,c(m,2);c(n,2)*c(m,2)*2;
累加下去
总共的排列数就是c(n+m,n);
相比一下,约分,就得到答案了;


代码如下:
 1 #include<iostream>
 2 #include<cstdio>
 3 #include<algorithm>
 4 #include<string>
 5 #include<cstring>
 6 #define ll long long
 7 using namespace std;
 8 int N,M;
 9 long long k,t;
10 long long zuhe(ll n,ll x)
11 {
12     ll di=1,gao=1,ans=1;
13     for(ll i=1;i<=x;i++)
14     { 
15         di*=n;
16         n--;
17         gao*=i;
18     }
19     return di/gao;
20 }
21 void yuefen(ll a,ll b)
22 {
23     for(ll i=b;i>=1;i--)
24     {
25         if(a%i==0&&b%i==0)
26         {
27             k=a/i;
28             t=b/i;
29             return ;
30         }
31     }
32 }
33 int main()
34 {
35 //    freopen("in","r",stdin);
36 //    freopen("out.out","w",stdout);
37     while(scanf("%d%d",&N,&M)!=EOF)
38     {
39         k=0;
40         t=0;
41         for(ll i=1;i<=N;i++)
42         {
43             ll temp=0;
44             if(i>M)
45                 break;
46             t+=i*zuhe(N,i)*zuhe(M,i);
47         }
48         k=zuhe(N+M,M);
49         yuefen(k,t);
50         printf("%I64d/%I64d\n",t,k);
51     }
52     return 0;
53 }
View Code

 

posted on 2015-04-10 03:19  bai_yan  阅读(204)  评论(0编辑  收藏  举报

导航