POJ1860_钉子小球_DP

题目描述

钉子和小球

Time Limit:1000MS     Memory Limit:10000KB     64bit IO Format:%I64d & %I64u

Submit  Status

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

解题报告

满眼钉子很瞎眼，我们可以把小球可放的位置理解成小洞。

小洞对应的位置有钉子便直接落到下一层，否则就左右两边放。

二维数组多好，偏要用一位数组模拟,实在脑残...

代码

#include <iostream>
#include <stdio.h>
#include <cstring>
using namespace std;

bool inp[13000];
__int64 trangle[14000];
__int64 gcd(__int64 a,__int64 b)                                                                                   //计算a,b的最大公约数
{
        if(b==0)
           return a;
        else
           return gcd(b,a%b);
}
/*************************/
int main()
{
    int n,m,num,now;//n为钉子行，m为第m个坑
    char temp;
    memset(trangle,0,sizeof(trangle));
        scanf("%d%d",&n,&m);
    if(n==0){printf("0/1");return 0;}
    num=n*(n+1)/2;
    for(int i=1;i<=num;i++){
            scanf("%c",&temp);
        while(temp!='*'&&temp!='.')
            scanf("%c",&temp);
        inp[i]=(temp=='*'?true:false);
    }
    trangle[1]=((__int64) 1)<<n;
    if(inp[1]){
                trangle[2]+=trangle[1]>>1;
                trangle[3]+=trangle[1]>>1;
              }
        else
        {
            trangle[5]=trangle[1];
        }
    //trangle[分子分母][第几个]
    for(int hang=1;hang<=n;hang++)
    for(int lie=1;lie<=hang+1;lie++){
        now=(hang+1)*(hang)/2+lie;
        if(trangle[now]==0)
                continue;
        if(inp[now]){
            trangle[now+hang+1]+=trangle[now]>>1;
            trangle[now+hang+2]+=trangle[now]>>1;

        }
        else{
            trangle[now+hang*2+4]+=trangle[now];
        }

    }
    __int64 all=((__int64) 1)<<n;
    if(trangle[num+m+1]==0)
        printf("0/1");
    else{
    __int64 t=gcd(all,trangle[num+m+1]);
    printf("%I64d/",trangle[num+m+1]/t);
    printf("%I64d\n",all/t);
    }
    return 0;


}