POJ 2096 Collecting Bugs G++ 概率dp 没掌握

#include <cstdio>
#include <iostream>
#include <cstring>
using namespace std;
//英语     看博友分析     抄博友好程序       POJac       百练wa       概率dp      没掌握 
double dp[1010][1010];
int main()
{
	int n,s;
	while(cin>>n>>s)
	{
		memset(dp,0,sizeof(dp));
		for(int i=n;i>=0;i--)
		{
			for(int j=s;j>=0;j--)
			{
				double sum=n*s-i*j;
				if(i==n && j==s)
				{
					continue;
				}
				dp[i][j]+=(dp[i][j]+1.0)*(double)i*(double)j/sum;
				dp[i][j]+=(dp[i+1][j]+1.0)*(double)(n-i)*(double)j/sum;
				dp[i][j]+=(dp[i+1][j+1]+1.0)*(double)(n-i)*(double)(s-j)/sum;
				dp[i][j]+=(dp[i][j+1]+1.0)*(double)i*(double)(s-j)/sum;
			}
		} 	
		printf("%.4f\n",dp[0][0]);	
	}	
	return 0;
}