钉子和小球
Time Limit: 1000MS | Memory Limit: 10000K | |
Total Submissions: 8073 | Accepted: 2516 |
Description
有一个三角形木板,竖直立放,上面钉着n(n+1)/2颗钉子,还有(n+1)个格子(当n=5时如图1)。每颗钉子和周围的钉子的距离都等于d,每个格子的宽度也都等于d,且除了最左端和最右端的格子外每个格子都正对着最下面一排钉子的间隙。
让一个直径略小于d的小球中心正对着最上面的钉子在板上自由滚落,小球每碰到一个钉子都可能落向左边或右边(概率各1/2),且球的中心还会正对着下一颗将要碰上的钉子。例如图2就是小球一条可能的路径。
我们知道小球落在第i个格子中的概率pi=pi= ,其中i为格子的编号,从左至右依次为0,1,...,n。
现在的问题是计算拔掉某些钉子后,小球落在编号为m的格子中的概率pm。假定最下面一排钉子不会被拔掉。例如图3是某些钉子被拔掉后小球一条可能的路径。
让一个直径略小于d的小球中心正对着最上面的钉子在板上自由滚落,小球每碰到一个钉子都可能落向左边或右边(概率各1/2),且球的中心还会正对着下一颗将要碰上的钉子。例如图2就是小球一条可能的路径。
我们知道小球落在第i个格子中的概率pi=pi= ,其中i为格子的编号,从左至右依次为0,1,...,n。
现在的问题是计算拔掉某些钉子后,小球落在编号为m的格子中的概率pm。假定最下面一排钉子不会被拔掉。例如图3是某些钉子被拔掉后小球一条可能的路径。
Input
第1行为整数n(2 <= n <= 50)和m(0 <= m <= n)。以下n行依次为木板上从上至下n行钉子的信息,每行中'*'表示钉子还在,'.'表示钉子被拔去,注意在这n行中空格符可能出现在任何位置。
Output
仅一行,是一个既约分数(0写成0/1),为小球落在编号为m的格子中的概pm。既约分数的定义:A/B是既约分数,当且仅当A、B为正整数且A和B没有大于1的公因子。
Sample Input
5 2 * * . * * * * . * * * * * * *
Sample Output
7/16
Source
原题链接:http://poj.org/problem?id=1189
题目是中文就不多啰嗦了,但是有一点,如果某个点没有钉子,那么小球会落到下面第二层的位置。
理解以后就和POJ1136 The Triangle差不多了。
由于概率每次都要处于2,并且分数加法有点麻烦,所以,开始时一个数的值设为2^n,
这样便可简化运算.
运算方法有两种,
一.dp[i][j]+=(dp[i-1][j-1]+dp[i-1][j])/2;
二.dp[i+1][j]+=dp[i][j]/2,dp[i+1][j+1]+=dp[i][j]/2;
相对来说第二种方法更简单一点.
样例过程:
32
16 0
8 8 0
4 0 20 0
2 2 10 10 0
1 2 14 10 5 0
AC代码:
/**
* 行有余力,则来刷题!
* 博客链接:http://blog.csdn.net/hurmishine
*
*/
#include <iostream>
#include <cstdio>
#include <cstring>
using namespace std;
typedef long long LL;
const int maxn=50+5;
char a[maxn][maxn];
LL dp[maxn][maxn];
int n,m;
LL GCD(LL x,LL y)
{
if(y==0)
return x;
return GCD(y,x%y);
}
int main()
{
//freopen("C:\\Documents and Settings\\Administrator\\桌面\\data.txt","r",stdin);
while(cin>>n>>m)
{
for(int i=0;i<n;i++)
{
for(int j=0;j<=i;j++)
cin>>a[i][j];
}
memset(dp,0,sizeof(dp));
LL maxx=(LL)1<<n;
dp[0][0]=maxx;
for(int i=0;i<n;i++)
{
for(int j=0;j<=i;j++)
{
if(a[i][j]=='*')
{
dp[i+1][j]+=dp[i][j]/2;
dp[i+1][j+1]+=dp[i][j]/2;
}
else
{
dp[i+2][j+1]+=dp[i][j];
//dp[i][j]=0;
}
}
}
/**
for(int i=0;i<=n;i++)
{
for(int j=0;j<=i;j++)
cout<<dp[i][j]<<" ";
cout<<endl;
}
*/
//cout<<dp[n][m]<<endl;
LL gcd=GCD(maxx,dp[n][m]);
cout<<dp[n][m]/gcd<<"/"<<maxx/gcd<<endl;
}
return 0;
}