NC16856. [NOI1999]钉子和小球
描述
有一个三角形木板,竖直立放,上面钉着n(n+1)/2颗钉子,还有(n+1)个格子(当n=5时如图1)。每颗钉子和周围的钉子的距离都等于d,每个格子的宽度也都等于d,且除了最左端和最右端的格子外每个格子都正对着最下面一排钉子的间隙。
让一个直径略小于d的小球中心正对着最上面的钉子在板上自由滚落,小球每碰到一个钉子都可能落向左边或右边(概率各1/2),且球的中心还会正对着下一颗将要碰上的钉子。例如图2就是小球一条可能的路径。
我们知道小球落在第i个格子中的概率pi= ,其中i为格子的编号,从左至右依次为0,1,...,n。
图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
C++(g++ 7.5.0) 解法, 执行用时: 3ms, 内存消耗: 464K, 提交时间: 2023-05-14 20:53:38
#include<bits/stdc++.h> #define ll long long using namespace std; ll f[60][60]; int main() { int n,m; cin>>n>>m; ll q=(ll)1<<n; f[0][0]=q; for(int i=0;i<n;i++) { for(int j=0;j<=i;j++) { char t; cin>>t; if(t=='*')f[i+1][j]+=f[i][j]/2,f[i+1][j+1]+=f[i][j]/2; else if(t=='.')f[i+2][j+1]+=f[i][j]; } } if(!f[n][m])cout<<"0/1"<<endl; else { ll w=__gcd(f[n][m],q); cout<<f[n][m]/w<<'/'<<q/w<<endl; } }
C++14(g++5.4) 解法, 执行用时: 4ms, 内存消耗: 480K, 提交时间: 2019-01-07 18:12:56
#include<iostream> using namespace std; long long n,m,f[55][55],tot; char map[55][55]; int main() { int i,j,k; cin>>n>>m;tot=1ll<<n; for(i=1;i<=n;i++)for(j=1;j<=i;j++)cin>>map[i][j]; f[1][1]=tot; for(i=1;i<=n;i++) for(j=1;j<=i;j++) if(map[i][j]=='*') f[i+1][j]+=f[i][j]/2,f[i+1][j+1]+=f[i][j]/2; else f[i+2][j+1]+=f[i][j]; if(f[n+1][m+1]==0){cout<<0<<'/'<<1;return 0;} while(f[n+1][m+1]%2==0){f[n+1][m+1]/=2;tot/=2;} cout<<f[n+1][m+1]<<'/'<<tot; return 0; }