CodeForces 398B

http://acm.hdu.edu.cn/showproblem.php?pid=6415

ZOJ 3822

这3题的模型好像都一样，但做法都不一样。。。

二维n*n,随机填格子，期望多少次每行每列都至少有一个格子被填。一开始有些格子已经被填

又纠结了半天，最后尝试了期望DP万能方法倒着DP想了想终于做出来了。

倒着dp又是解方程。。。

设f[i][j]为占领i行j列，到占领n行n列期望要多少次。

f[i][j]=(f[i][j]*i*j+f[i+1][j]*(n-i)*j+f[i][j+1]*i*(n-j)+f[i+1][j+1]*(n-i)*(n-j))/n/n

然后倒着推，解这个方程，就能该解出f[i][j]了。

由于一开始已经放了一些格子，那么只要到那时的占领情况就是答案了。

#include<bits/stdc++.h>
#define maxl 2010

int n,m,rcnt,ccnt;
int rin[maxl],cin[maxl];
double f[maxl][maxl];

inline void prework()
{
	scanf("%d%d",&n,&m);
	int r,c;
	for(int i=1;i<=m;i++)
	{
		scanf("%d%d",&r,&c);
		if(!rin[r]) rin[r]=true,rcnt++;
		if(!cin[c]) cin[c]=true,ccnt++;
	}
}

inline void mainwork()
{
	double tmp;
	for(int i=n;i>=0;i--)
		for(int j=n;j>=0;j--)
		{
			if(i==n && j==n) continue;
			tmp=1.0*i*j/n/n;
			if(i<n)
				tmp+=(f[i+1][j]+1)*1.0*(n-i)*j/n/n;
			if(j<n)
				tmp+=(f[i][j+1]+1)*1.0*i*(n-j)/n/n;
			if(i<n && j<n)
				tmp+=(f[i+1][j+1]+1)*1.0*(n-i)*(n-j)/n/n;
			f[i][j]=tmp/(1.0-1.0*i*j/n/n);
		}
}

inline void print()
{
	printf("%.4f",f[rcnt][ccnt]);
}

int main()
{
	prework();
	mainwork();
	print();
	return 0;
}