【NOIP 模拟题】掷骰子（dp）
掷骰子（dice.cpp）
【问题描述】
太郎和一只兔子正在玩一个掷骰子的游戏。有一个N个格子的长条棋盘，太郎和兔子轮流掷一个M面的骰子，骰子M面分别是1到M的数字，且掷到每一面的概率是相同的，掷到几，就往前走几步，当谁走到第N格时，谁就获胜了。游戏中还有一个规则“反弹”，就是当一位选手走到N格外时，他就会后退（就像飞行棋进营一样）。
假设现在一位选手在A格，当他掷出B时：1）若A+B<N,走到A+B格；2）若A+B=N,走到第N格，获胜；3）若A+B>N,走到N-(A+B-N)格。
假设现在太郎和兔子分别在x格、y格，接下来是太郎掷骰子，太郎想知道他赢得比赛的概率是多少。
【问题输入】
一行4格整数N、M、x、y。
【问题输出】
一行一个小数，表示太郎获胜的概率（保留六位小数）
【样例输入】
10 6 1 1
【样例输出】
0.541725
【数据范围】
30%的数据：10<=N<=100；
100%的数据：10<=N<=2000,M，x,y<=n-1。
————————————————————————————————
【题解】【dp】
【f[i][j]表示A在第i格，B在第j格时，A获胜的概率。特别的：若i=n或j=n时，概率是1或0】
【分四种情况考虑：
[1)若i+m<=n、j+m<=n，即i、j都不能一步获胜]
（A从走到i+1-i+m这之间的位置的可能性都是均等的，都是1/m，同理，B从走到j+1-j+m这之间的位置的可能性都是1/m,所以每一个f[i'][j']出现的概率都是1/（m²））
【2）i+m<=n、j+m>n [A不能一步到终点，B可能一步到终点]】
j只有1中可能上次从n转移过来，j共有m中取值；i也有m种取值，所以共有m²种可能，再加上i这一步可以到终点的一种可能性
【3）i+m>n,j+m<=n[A可能一步到终点，B不可能一步到终点]】
（这种情况下，与第2种情况相似）
【4）i+m>n,j+m>n[A、B都有可能一步到终点]】
在这种情况下，因为有反弹，soA、B的位置只可能在[n-m+1,n]中，两人一步到终点的几率都是1/m,
那么，A先到n的可能就是
等比数列求和：
化简得：
【5）最后，还有一点，若i=n-m,当i'=j=n的概率为1】
。
[转自： http://blog.csdn.net/lyd_7_29/article/details/52245309]
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
int n,m,x,y;
double f[2010][2010];//f[i][j]表示A在第i格,B在第j格时，A获胜的概率
inline double SUM(int x1,int y1,int x2,int y2)
{
return (f[x1][y1]-f[x1][y2]-f[x2][y1]+f[x2][y2]);
}
int main()
{
freopen("dice.in","r",stdin);
freopen("dice.out","w",stdout);
int i,j;
scanf("%d%d%d%d",&n,&m,&x,&y);
for(i=n;i>0;--i)
for(j=n;j>0;--j)
{
f[i][j]=f[i+1][j]+f[i][j+1]-f[i+1][j+1];
if(i==n||j==n) {f[i][j]+=(i==n); continue;}//有一个已到终点
if(i<=n-m&&j<=n-m) {f[i][j]+=SUM(i+1,j+1,i+m+1,j+m+1)/m/m; continue; }//A、B都不能一步到终点
if(i<=n-m&&j>n-m) {f[i][j]+=((m-1)*SUM(i+1,j,i+m+1,j+1)+(i+m==n))/m/m; continue;}//A不能一步到终点，B可能一步到终点
if(i>n-m&&j<=n-m) {f[i][j]+=((m-1)*SUM(i,j+1,i+1,j+m+1)+m)/m/m; continue; }//A可能一步到终点，B不能一步到终点
if(i>n-m&&j>n-m) {f[i][j]+=1.0*m/(2*m-1); continue;}//A、B都可能一步到终点
}
printf("%.6lf\n",SUM(x,y,x+1,y+1));
return 0;
}
