【NOIP2012模拟10.25】旅行

题目

给定一个n行m列的字符矩阵，’.’代表空地，’X’代表障碍。移动的规则是：每秒钟以上下左右四个方向之一移动一格，不能进入障碍。
计算：在空地中随机选择起点和终点（可以重合，此时最短耗时为0），从起点移动到终点最短耗时的平均值。
每一行每一列至多有1个障碍，并且障碍不在对角线方向相邻。以下矩阵是不合法的：
.X

X.

分析

50%的数据，全都是空地，答案就是所有两点间的曼哈顿距离和。
很容易求，({第i行空地的数量}*{第j行空地的数量}*|i-j|*2)
因为每一行每一列至多有1个障碍，并且障碍不在对角线方向相邻。得出从A走到B的耗时，要么等于AB的曼哈顿距离，要么等于AB的曼哈顿距离+2。
现在考虑什么情况下距离加二。
观察下图，规律显而易见，我就不多说了，

#include <cmath>
#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <algorithm>
#include <queue>
const int maxlongint=2147483647;
const int mo=1000000007;
const int N=1005;
using namespace std;
int a[N][N],n,m,sum,sx[N],sy[N],xx[N],xy[N];
double ans;
int main()
{
	scanf("%d%d",&n,&m);
	for(int i=1;i<=n;i++)
		for(int j=1;j<=m;j++)
		{
			char c=getchar();
			while(c!='.' && c!='X') c=getchar();
			if(c=='.')
			{
				sx[i]++;
				sy[j]++;
				sum++;
			}
			else
			{
				xx[i]=j;
				xy[j]=i;	
			}
		}
	for(int i=1;i<=n;i++)
		for(int j=1;j<=n;j++)
			ans+=abs(i-j)*sx[i]*sx[j]*1.0/sum/sum;
	for(int i=1;i<=m;i++)
		for(int j=1;j<=m;j++)
			ans+=abs(i-j)*sy[i]*sy[j]*1.0/sum/sum;
	for(int i=1;i<=n;i++)
	if(xx[i])
	{
		int mn=maxlongint,mx=-1;
		for(int j=i;j>=1;j--)
		{
			if(xx[j]>mx && xx[j])
			{
				ans+=1.0*(xx[i]-1)*(m-xx[j])*2/sum/sum;
				mx=xx[j];
			}
			else
			break;
		}
		for(int j=i;j>=1;j--)
		{
			if(xx[j]<mn && xx[j])
			{
				ans+=1.0*(xx[j]-1)*(m-xx[i])*2/sum/sum;
				mn=xx[j];
			}
			else
			break;
		}
		mn=xx[i],mx=xx[i];
		for(int j=i+1;j<=n;j++)
		{
			if(xx[j]>mx && xx[j])
			{
				ans+=1.0*(xx[i]-1)*(m-xx[j])*2/sum/sum;
				mx=xx[j];
			}
			else
			break;
		}
		for(int j=i+1;j<=n;j++)
		{
			if(xx[j]<mn && xx[j])
			{
				ans+=1.0*(xx[j]-1)*(m-xx[i])*2/sum/sum;
				mn=xx[j];
			}
			else
			break;
		}
	}
	for(int i=1;i<=m;i++)
	if(xy[i])
	{
		int mn=maxlongint,mx=-1;
		for(int j=i;j>=1;j--)
		{
			if(xy[j]>mx && xy[j])
			{
				ans+=1.0*(xy[i]-1)*(n-xy[j])*2/sum/sum;
				mx=xy[j];
			}
			else
			break;
		}
		for(int j=i;j>=1;j--)
		{
			if(xy[j]<mn && xy[j])
			{
				ans+=1.0*(xy[j]-1)*(n-xy[i])*2/sum/sum;
				mn=xy[j];
			}
			else
			break;
		}
		mn=xy[i],mx=xy[i];
		for(int j=i+1;j<=m;j++)
		{
			if(xy[j]>mx && xy[j])
			{
				ans+=1.0*(xy[i]-1)*(n-xy[j])*2/sum/sum;
				mx=xy[j];
			}
			else
			break;
		}
		for(int j=i+1;j<=m;j++)
		{
			if(xy[j]<mn && xy[j])
			{
				ans+=1.0*(xy[j]-1)*(n-xy[i])*2/sum/sum;
				mn=xy[j];
			}
			else
			break;
		}
	}
	printf("%.4lf",ans);
}