CodeForces 57 D.Journey（组合数学）

Description

给出一$n\times m$的地图，$X$表示该位置被占，否则是空，问任意两个空位置的最短距离的期望值，两个位置是可达的当且仅当两个位置均为空且共边，保证任意两个$X$不会同行同列也不会对角相邻

Input

第一行两个整数$n,m$表示地图行列数，之后输入一$n\times m$的地图$(2\le n,m\le 1000)$

Output

输出任意两个空位置最短距离的期望值，误差不超过$10^{-6}$

Sample Input

2 2
$..$
$.X$

Sample Output

0.888888888889

Solution

首先不考虑$X$对最短路的影响，那么两个位置$(x_1,y_1)$和$(x_2,y_2)$的最短距离即为$|x_2-x_1|+|y_2-y_1|$，对于位置$(x,y)$，其到这$nm$个点的最短距离之和为:

$\sum\limits_{i=1}^n\sum\limits_{j=1}^m(|x-i|+|y-j|)=m\sum\limits_{i=1}^n|x-i|+n\sum\limits_{j=1}^m|y-j|=m\cdot (C_x^2+C_{n-x+1}^2)+n\cdot (C_y^2+C_{m-y+1}^2)$

根据该公式可以求出任意两点最短距离和$S_1$，所有$X$到其他点的最短距离和$S_2$，任意两个$X$的最短距离和$S_3$

注意其中$S_2,S_3$在考虑$X$对最短路的影响时也不变，因为任意两个$X$不同行不同列不对角相邻，这样就不会出现连续的（此处连续指的是把共边或对角相邻的点视为相邻）的$X$把这个$X$和其他点完全隔开，如此求可以求出不考虑$X$对最短路影响时任意两个空位置的最短距离和$S_1-2(S_2-S_3)-S_3=S_1-2S_2+S_3$

但是注意到$X$会对两个空位置之间的最短距离造成影响，对于$(x,y)$处的$X$，第$x$行和第$y$列只有这一个$X$，那么第$x$行前$y-1$个空位置到后$m-y$个位置的最短路被$X$隔离开了，必须绕路才可达，绕路距离至少多$2$，且多$2$方案存在，因为$X$一圈都没有其他$X$，所以只需多走两步绕过$X$即可，这样对每个点对都会多$2$步距离，总共就会多$4\cdot (y-1)\cdot (m-y)​$步距离，同理考虑第$y$列会多出$4\cdot (x-1)\cdot (n-x)​$步距离

上面只考虑了$X$对其所在行和所在列的影响，但是注意到如果一些相邻行都有$X$且这些$X$的列数单调，或者一些相邻列都有$X$且这些$X$的行数单调，例如第$i$行到第$j$行都有$X$，且这些$X$的列数为$a_i<a_{i+1}<...<a_j$，那么第$i$行前$a_i-1$个元素到第$k$行后$n-a_k$个元素的所有最短路径均会被这些$X$阻挡，必须要绕过一个$X$才可达，故最短距离步数增加$2$，总的多出$4\cdot (a_i-1)\cdot (m-a_k)$步，$i<k\le j$，以此累计多出的步数加到答案里即可

现在求出了所有空位置最短距离之和，除以$(nm-num)^2$即为答案，其中$num$为 X X <script type="math/tex" id="MathJax-Element-69">X</script>的个数

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;
typedef pair<int,int>P;
const int INF=0x3f3f3f3f,maxn=1005;
int n,m,x[maxn],y[maxn],row[maxn],col[maxn];
char s[maxn][maxn]; 
int C(int n)
{
    return n*(n-1)/2;
}
ll Solve(int x,int y)
{
    return (ll)m*(C(x)+C(n-x+1))+(ll)n*(C(y)+C(m-y+1));
}
int main()
{
    while(~scanf("%d%d",&n,&m))
    {
        int res=0;
        memset(row,0,sizeof(row));
        memset(col,0,sizeof(col));
        for(int i=1;i<=n;i++)
        {
            scanf("%s",s[i]+1);
            for(int j=1;j<=m;j++)
                if(s[i][j]=='X')
                    row[i]=j,col[j]=i,x[res]=i,y[res++]=j;
        }
        ll sum0=0,sum1=0,sum2=0;
        for(int i=1;i<=n;i++)
            for(int j=1;j<=m;j++)
                sum0+=Solve(i,j);
        for(int i=0;i<res;i++)
            for(int j=0;j<res;j++)
                sum1+=abs(x[j]-x[i])+abs(y[j]-y[i]);
        for(int i=0;i<res;i++)sum2+=Solve(x[i],y[i]);
        sum2=2ll*sum2-sum1;
        sum0-=sum2;
        for(int i=1;i<=n;i++)
            if(row[i])
            {
                sum0+=4ll*(row[i]-1)*(m-row[i]);
                for(int j=i-1;j>=1;j--)
                    if(!row[j]||row[j]<row[j+1])break;
                    else sum0+=4ll*(row[i]-1)*(m-row[j]);
                for(int j=i+1;j<=n;j++)
                    if(!row[j]||row[j]<row[j-1])break;
                    else sum0+=4ll*(row[i]-1)*(m-row[j]);
            }
        for(int i=1;i<=m;i++)
            if(col[i])
            {
                sum0+=4ll*(col[i]-1)*(n-col[i]);
                for(int j=i-1;j>=1;j--)
                    if(!col[j]||col[j]<col[j+1])break;
                    else sum0+=4ll*(col[i]-1)*(n-col[j]);
                for(int j=i+1;j<=m;j++)
                    if(!col[j]||col[j]<col[j-1])break;
                    else sum0+=4ll*(col[i]-1)*(n-col[j]);
            }       
        double temp=1.0*(n*m-res)*(n*m-res);
        printf("%.10f\n",1.0*sum0/temp);
    }
    return 0;
}