lightoj 1265 概率dp

題意：英文不難（因爲我都能看懂。。），大意就是在一個荒島上只有你和其他的一種吃人的動物和被吃的動物，然後題中說明了相遇時發生的事件，食肉動物數目爲0的時候可以逃生，問最後成功逃生的概率

思路：推推推的就寫出了轉移方程，因爲C（n，2）忘了除以2錯了一天。。。

代碼：

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <iostream>
#include <vector>
#include <cmath>
using namespace std;
#define ll long long
#define mem(a,b) memset(a,b,sizeof(a));
const int maxn = 1010;
double dp[maxn][maxn];
void init(){
	mem(dp,0);
	for(int i = 0;i < maxn;i ++)dp[0][i] = 1;
	for(int i = 2;i < maxn;i ++){
		dp[i][0] = dp[i-2][0]*(i-1)/(double)(i+1);
		for(int j = 1;j < maxn;j ++){
			double pa,pb,ma,mb;
			if(j >= 2){
				ma = (j+i+1)*(i+j) - j*(j-1);
				mb = (j+i+1)*(i+j) - j*(j-1) - j;
				ma /= 2.0;
				mb /= 2.0;
			}
			else {
				ma = (j+i+1)*(i+j);
				mb = (j+i+1)*(i+j) - j;
				ma /= 2.0;
				mb /= 2.0;
			}
			pa = j*dp[i][j-1] + dp[i][j-1]*i*j + i*(i-1)*dp[i-2][j]/2.0;
			pb = dp[i][j-1]*i*j + i*(i-1)*dp[i-2][j]/2.0;
			dp[i][j] = max(pa/ma,pb/mb);
		}
	}
}
int main(){
	int T;//printf("%d\n",'M');
	cin>>T;
	int cas = 1;
	init();
	while(T --){
		int a,b;
		cin>>a>>b;
		printf("Case %d: %.12f\n",cas ++,dp[a][b]);
	}
    return 0;
}