寒假学习记录D24

题目描述

有一个三角形木板,竖直立放，上面钉着n(n+1)/2颗钉子，还有(n+1)个格子（当n=5时如图1）。每颗钉子和周围的钉子的距离都等于d，每个格子的宽度也都等于d，且除了最左端和最右端的格子外每个格子都正对着最下面一排钉子的间隙。

让一个直径略小于d的小球中心正对着最上面的钉子在板上自由滚落，小球每碰到一个钉子都可能落向左边或右边（概率各1/2），且球的中心还会正对着下一颗将要碰上的钉子。例如图2就是小球一条可能的路径。

我们知道小球落在第i个格子中的概率pi=

 ，其中i为格子的编号，从左至右依次为0,1,...,n。

现在的问题是计算拔掉某些钉子后，小球落在编号为m的格子中的概率pm。假定最下面一排钉子不会被拔掉。例如图3是某些钉子被拔掉后小球一条可能的路径。

              图1                                  图2                              图3

输入描述:

第1行为整数n（2 ≤ n ≤ 50）和m（0 ≤ m ≤ n）。以下n行依次为木板上从上至下n行钉子的信息，每行中‘*’表示钉子还在，‘.’表示钉子被拔去，注意在这n行中空格符可能出现在任何位置。

输出描述:

仅一行，是一个既约分数(0写成0/1)，为小球落在编号为m的格子中的概pm。既约分数的定义：A/B是既约分数，当且仅当A、B为正整数且A和B没有大于1的公因子。

示例1

输入

5 2
    *
   * .
  * * *
 * . * *
* * * * *

输出

7/16

#include <iostream>
   #include <cstdio>
   #include <cstring>
   
   using namespace std;
  
   char s[55][55];
   
   long long gcd(long long x,long long y){
      if(y==0) return x;
      else return gcd(y,x%y);
  }
  
  long long lc(long long x,long long y){
      return x/gcd(x,y)*y;
  }
  
  int main(){
      int n,m;
      while(scanf("%d %d",&n,&m)!=EOF){
          pair<long long,long long>  dp[2][150];
          for(int i=0;i<n;i++)
    for(int j=0;j<=i;j++){
              char c=' ';
              for(;c!='.'&&c!='*';scanf("%c",&c));
              s[i][j]=c;
          }
          int e=0;
          dp[e][n].first=dp[e][n].second=1;
          for(int i=0;i<n;i++){
              memset(dp[e^1],0,sizeof(long long)*150*2);
             for(int j=0;j<=i;j++){
                  int k=n-i+2*j;
                  if(s[i][j]=='*'&&dp[e][k].first!=0){
                      if(dp[e^1][k-1].first==0){
                          dp[e^1][k-1].first=dp[e][k].first;
                          dp[e^1][k-1].second=dp[e][k].second*2;
                      }
                      else{
                         long long t=lc(dp[e][k].second*2,dp[e^1][k-1].second);
                         dp[e^1][k-1].first=dp[e^1][k-1].first*(t/dp[e^1][k-1].second)+dp[e][k].first*(t/(dp[e][k].second*2));
                          dp[e^1][k-1].second=t;
                      }
                      if(dp[e^1][k+1].first==0){
                          dp[e^1][k+1].first=dp[e][k].first;
                          dp[e^1][k+1].second=dp[e][k].second*2;
                      }
                      else{
                          long long t=lc(dp[e][k].second*2,dp[e^1][k+1].second);
                          dp[e^1][k+1].first=dp[e^1][k+1].first*(t/dp[e^1][k+1].second)+dp[e][k].first*(t/(dp[e][k].second*2));
                          dp[e^1][k+1].second=t;
                     }
                      long long t=gcd(dp[e^1][k-1].first,dp[e^1][k-1].second);
                      if(t!=0&&t!=1){
                          dp[e^1][k-1].first/=t;
                          dp[e^1][k-1].second/=t;
                      }
                      t=gcd(dp[e^1][k+1].first,dp[e^1][k+1].second);
                     if(t!=0&&t!=1){
                         dp[e^1][k+1].first/=t;
                          dp[e^1][k+1].second/=t;
                      }
                  }
                  else if(s[i][j]=='.'&&dp[e][k].first!=0){
                      if(dp[e^1][k].first==0){
                          dp[e^1][k].first=dp[e][k].first;
                          dp[e^1][k].second=dp[e][k].second;
                      }
                      else{
                          long long t=lc(dp[e][k].second,dp[e^1][k].second);
                          dp[e^1][k].first=dp[e^1][k].first*(t/dp[e^1][k].second)+dp[e][k].first*(t/dp[e][k].second);
                          dp[e^1][k].second=t;
                      }
                      long long t=gcd(dp[e^1][k].first,dp[e^1][k].second);
                     if(t!=0&&t!=1){
                          dp[e^1][k].first/=t;
                          dp[e^1][k].second/=t;
                      }
                  }
              }
              for(int j=0,k=n-i+1;j<i;k+=2,j++){
                  if(dp[e][k].first!=0){
                      if(dp[e^1][k].first==0){
                          dp[e^1][k].first=dp[e][k].first;
                          dp[e^1][k].second=dp[e][k].second;
                      }
                      else{
                          long long t=lc(dp[e][k].second,dp[e^1][k].second);
                          dp[e^1][k].first=dp[e^1][k].first*(t/dp[e^1][k].second)+dp[e][k].first*(t/dp[e][k].second);
                          dp[e^1][k].second=t;
                      }
                      long long t=gcd(dp[e^1][k].first,dp[e^1][k].second);
                      if(t!=0&&t!=1){
                          dp[e^1][k].first/=t;
                         dp[e^1][k].second/=t;
                      }
               }
              }
             e^=1;
         }
        m*=2;
         if(dp[e][m].first==0)
             puts("0/1");
         else printf("%lld/%lld\n",dp[e][m].first,dp[e][m].second);
     }
     return 0;
 }