POJ1189:钉子和小球(dp)

钉子和小球
Time Limit: 1000MS Memory Limit: 10000K
Total Submissions: 7903 Accepted: 2420

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是某些钉子被拔掉后小球一条可能的路径。 

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

思路:用2^n(路径总数)个球模拟一下就行了。

# include <iostream>
# include <cstdio>
# include <cstring>
# define LL long long

using namespace std;
LL gcd(LL a, LL b)
{
    return b==0?a:gcd(b, a%b);
}

int main()
{
    LL dp[53][53];
    bool map[53][53];
    int n, m;
    while(~scanf("%d%d",&n,&m))
    {
        getchar();
        memset(dp, 0, sizeof(dp));
        dp[1][1] = (LL)1<<n;//不能写成1<<n。
        for(int i=1; i<=n; ++i)
        {
            int k=1;
            string s;
            getline(cin, s);
            for(int j=0; j<s.size(); ++j)
            {
                if(s[j]=='*')
                    map[i][k++] = true;
                else if(s[j] == '.')
                    map[i][k++] = false;
                if(k == i+1)
                    break;
            }
        }
        for(int i=1; i<=n; ++i)
            for(int j=1; j<=i; ++j)
            {
                if(map[i][j])//如果有钉,分两半。
                {
                    dp[i+1][j] += (dp[i][j]>>1);
                    dp[i+1][j+1] += (dp[i][j]>>1);
                }
                else//如果没钉,直接掉下去。
                    dp[i+2][j+1] += dp[i][j];
            }
        LL g = gcd(dp[n+1][m+1], dp[1][1]);
        printf("%lld/%lld\n",dp[n+1][m+1]/g, dp[1][1]/g);
    }
    return 0;
}


转载于:https://www.cnblogs.com/junior19/p/6729987.html

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值