钉子与小球

Posted 【對策局】

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了钉子与小球相关的知识,希望对你有一定的参考价值。

【题目描述】

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

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

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

技术分享
【输入描述】

第1行为整数n(2<=n<=50)。

以下n行依次为木板上从上至下n行钉子的信息,每行中‘*’表示钉子还在,‘.’ 表示钉子被拔去(最下面一排的钉子不会拔掉),注意在这n行中空格符可能出现在任何位置。

【输出描述】

共n+1行,每一行是一个既约分数(0写成0/1),为小球落在编号为0到编号为n这n+1个的格子中的概率m。

既约分数的定义:A/B是既约分数,当且仅当A、B为正整数且A和B没有大于1的公因子。

【样例输入】

4

*

* *

* . *

* * * *

【样例输出】

1/16

1/8

5/8

1/8

1/16

【数据范围及提示】

(2<=n<=50)

源代码:

#include<iostream>
using namespace std;
long long n,i1[61][61]={0},i2[61][61]={0}; //Long Long甚至也会爆掉!还得用cin、cout!
bool f[61][61]={0};
void x1(long long &x1,long long &y1,long long x2,long long y2) //辗转相除法进行坑爹的处理。
{
    if (!y2)
      return;
    if (x2&1)
      y2=y2<<1;
    else
      x2=x2>>1;
    if (!y1)
    {
        x1=x2;
        y1=y2;
    }
    else
    {
        long long t1=y1,t2=y2;
        while (t1%t2)
        {
            long long t=t1;
            t1=t2;
            t2=t%t2;
        }
        long long t=(y1*y2)/t2;
        t1=x1*(t/y1)+x2*(t/y2);
        t2=t;
        x1=t1;
        y1=t2;
        while (t1%t2)
        {
            long long t=t1;
            t1=t2;
            t2=t%t2;
        }
        if (!(x1%t2))
        {
            x1/=t2;
            y1/=t2;
        }
    }
}
int main()
{
    cin>>n;
    for (int a=1;a<=n;a++)
      for (int b=1;b<=a;b++)
      {
          char t;
          cin>>t;
          while (t!=*&&t!=.)
            cin>>t;
          if (t==*)
            f[a][b]=true;
      }
    for (int a=1;a<=n+1;a++)
      f[n+1][a]=true; 
    if (f[1][1])
      i1[1][1]=i2[1][1]=1;
    else
    {
        int t1=3,t2=2;
        while (!f[t1][t2]&&(t1+2)<=n+1)
        {
            t1+=2;
            t2++;
        }
        if (f[t1][t2])
          x1(i1[t1][t2],i2[t1][t2],2,1);
    }
    for (int a=1;a<n;a++)
      for (int b=1;b<=a;b++)
      {
        if (f[a+1][b])
          x1(i1[a+1][b],i2[a+1][b],i1[a][b],i2[a][b]);
        else
        {
            int t1=a+3,t2=b+1;
            while (!f[t1][t2]&&(t1+2)<=n+1)
            {
                t1+=2;
                t2++;
            }
            if (f[t1][t2])
              x1(i1[t1][t2],i2[t1][t2],i1[a][b],i2[a][b]);
        }
        if (f[a+1][b+1])
          x1(i1[a+1][b+1],i2[a+1][b+1],i1[a][b],i2[a][b]);
        else
        {
            int t1=a+3,t2=b+2;
            while (!f[t1][t2]&&(t1+2)<=n+1)
            {
                t1+=2;
                t2++;
            }
            if (f[t1][t2])
              x1(i1[t1][t2],i2[t1][t2],i1[a][b],i2[a][b]);
        }
      }
    for (int a=1;a<=n;a++)
    {
        x1(i1[n+1][a],i2[n+1][a],i1[n][a],i2[n][a]);
        x1(i1[n+1][a+1],i2[n+1][a+1],i1[n][a],i2[n][a]);
    }
    for (int a=1;a<=n+1;a++)
      if (i1[n+1][a])
        cout<<i1[n+1][a]<<"/"<<i2[n+1][a]<<endl;
      else
        cout<<"0/1"<<endl;
    return 0;
} //水完这道题整个人都不好了。

以上是关于钉子与小球的主要内容,如果未能解决你的问题,请参考以下文章

[POJ1189][BZOJ1867][CODEVS1709]钉子和小球

[bzoj1867][Noi1999][钉子和小球] (动态规划)

bzoj1867钉子和小球

bzoj1867: [Noi1999]钉子和小球(DP)

bzoj千题计划189:bzoj1867: [Noi1999]钉子和小球

bzoj 1867: [Noi1999]钉子和小球dp