LightOJ 1265 Island of Survival（概率DP）

Description
一座孤岛上有一个人n只老虎m只鹿，两只老虎遇见都死掉，两只鹿遇见什么也不会发生，老虎遇见人或鹿都会把他们吃掉，人遇见鹿可以选择杀死鹿或不杀死鹿，任意两只动物相遇的概率相同，问人活下来的最大几率
Input
第一行一整数T表示用例组数，每组用例输入两个整数n和m分别表示老虎的数量和鹿的数量(T<=200,0<=t,d<=1000)
Output
对于每组用例，输出人活下来的最大几率
Sample Input
4
0 0
1 7
2 0
0 10
Sample Output
Case 1: 1
Case 2: 0
Case 3: 0.3333333333
Case 4: 1
Solution
令dp[i][j]表示i只老虎j只鹿时人活下来的最大几率，那么有
dp[0][i]=1,0<=i<=1000（人杀掉所有鹿）
dp[1][0]=0（老虎吃掉人）
dp[i][j]=(C(i,2)*dp[i-2][j]+C(j,2)*dp[i][j]++i*j*dp[i][j-1]+i*0+j*dp[i][j])/C(i+j+1,2)
（第一项表示两只老虎相遇会死掉两只老虎，第二项表示两只鹿相遇什么都没发生，第三项表示一只老虎遇见一只鹿会死掉一只鹿，第四项表示人遇到老虎会死掉，第五项表示人遇见鹿不会杀掉鹿（这样可以减小人遇见老虎的几率））
解该方程即得到dp[i][j]
Code

#include<cstdio>
#include<iostream>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<vector>
#include<queue>
#include<map>
#include<set>
#include<ctime>
using namespace std;
typedef long long ll;
#define INF 0x3f3f3f3f
#define maxn 1111
int T,n,m;
double dp[maxn][maxn];
int main()
{
    memset(dp,0,sizeof(dp));
    for(int i=0;i<=1000;i++)dp[0][i]=1;
    dp[1][0]=0;
    for(int i=1;i<=1000;i++)
        for(int j=0;j<=1000;j++)
        { 
            double a=(i+j+1)*(i+j)/2,b=0;
            if(i>=2)b+=1.0*i*(i-1)/2*dp[i-2][j];//T&T 
            if(j>=2)a-=1.0*j*(j-1)/2;//D&D 
            if(j>0)
            {
                b+=1.0*i*j*dp[i][j-1];//T&D 
                a-=j;//D&P P not kill D 
            }
            dp[i][j]=b/a; 
        }
    int res=1;
    scanf("%d",&T);
    while(T--)
    {
        scanf("%d%d",&n,&m);
        printf("Case %d: %.10f\n",res++,dp[n][m]);
    }
    return 0;
}