AOJ B 钉子小球模拟题

B - 钉子和小球
Time Limit: 1000MS
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

以前用的记忆化搜索，这次用的简单的DP，KD的是数据错了！！这题没有A过，悲催= =

YXG你小心点给我！

#include<iostream>
#define ull unsigned long long
using namespace std;

ull Inf;

ull gcd( ull a,ull b )
{
        if( b==0 )return a;
        return gcd( b,a%b );
}

int main()
{
    int n,m;
    char map[160][160];
    int map1[160][160];
    ull sum[160][160];
    while( scanf( "%d %d",&n,&m )!=EOF )
    {
           int i,j;
           for( i=1;i<=n;i++ )
           for( j=1;j<=i;j++ )
                scanf( " \n%c",&map[i][j] );
           
           for( i=1;i<=n;i++ )
           for( j=1;j<=i;j++ )
                if( map[i][j]=='*')
                    map1[i][j]=1;
                else
                    map1[i][j]=0;
           
           memset( sum,0,sizeof(sum) );
           sum[1][1]=1;
           for( i=1;i<=n;i++ )
           for( j=1;j<=i;j++ )
                if( map1[i][j]==1 )
                {
                    sum[i+1][j]+=sum[i][j];
                    sum[i+1][j+1]+=sum[i][j];
                }
                else
                    sum[i+2][j+1]+=4*sum[i][j];
           Inf=0;
           for( i=1;i<=n+1;i++ )
                Inf+=sum[n+1][i];
           ull ans=gcd( Inf,sum[n+1][m+1] );
           
           printf( "%llu/%llu\n",sum[n+1][m+1]/ans,Inf/ans );
    }
    return 0;
}