HDU 5725/2016多校1C Game

题意：给一个图，其中包含一些点G表示障碍，并且两个G不会共行共列或有任何相交的点或边。求任意两点(不包含G)间的最短距离的期望

题解：

所求的期望即最短距离和除以点对数，点对数即非G点的数量的平方。

首先G的分布很有特点，也就较稀疏，只有很少的路径会被影响，也就是官方题解所说的两个点之间的路径最多被一个G所影响，并且被影响也只是会多走两条边，

那么我们可以先计算不考虑G时两点之间最短路径总和，再加上G的影响

对于任意两点之间的最短路径和 ，我最初的想法是求不考虑G时的任意两点之间的路径和 再减去每一个G到任意一点的，再加上两个G之间的，不过这些维护起来较复杂=。= 

所以改为用dp维护，某一个点到其左上方每一个非G的距离和，再处理右上方的，左下和右下是对称的关系所以乘2即可，但这样会把正上方距离的重复算两遍，所以干脆将行内和列内的单独处理，dp只维护左且上/有且上的值。 

而对于G值影响的路径数，开始以为只有行内列内会影响

# # # # # # 

#G # # # #

# # # # # #

# #G # # #

# # # # # #

# # #G # #

# # # # # # 

而像如上例子中 最左边的G的正上和最右边的G的正下直接的最短路径 也是会受到影响。

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
#include<iostream>
#include<algorithm>
#include<vector>
#include<list>
#include<set>
#include<map>
#include<stack>
#include<queue>
#include<deque>
#define mem(x,y) memset(x,y,sizeof(x))
#define pb push_back
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
typedef pair<int,int> pii;
#define bug puts("===========");
#define zjc puts("");
const double pi=(acos(-1.0));
const double eps=1e-8;
const ll INF=1e18+10;
const ll inf=1e9+10;
const int mod=1e9+7;
const int maxn=1000+10;
const int maxm=1e6+5;
char a[maxn][maxn];//原矩阵
ll hang[maxn],lie[maxn];///lie[i]表示i列出现的G的位置
ll flag[maxn];///
ll dp[maxn][maxn];///唯一其左且上的点 到这个点的距离
int n,m;
ll go(){
    ll ans=0;
    mem(flag,0);
    mem(dp,0);
    for(int i=1;i<=n;i++){
        ll sum=0;
        for(int j=1;j<=m;j++){
            dp[i][j]=dp[i][j-1]+(i-1)*(j-1)-sum;
            if(j!=1) dp[i][j] += i*(i-1)/2;
            if(flag[j-1]&&flag[j-1]!=i) {
                dp[i][j]-=i-flag[j-1]+1;
                sum++;
            }
            if(a[i][j]!='G') ans+=dp[i][j];
            else flag[j]=i;
        }
    }
    return ans;
}
ll qsum(ll x){
    return x*(x+1)/2;
}
ll fuck()
{
    ll ans=0;
    ans+=go();
    for(int i=1;i<=n;i++)
        for(int j=1;j<=m/2;j++) swap(a[i][j],a[i][m-j+1]);
    ans+=go();
    return ans;
}
ll z[maxn];
int main()
{
    int T_T;
    for(int i=1;i<maxn;i++) z[i]=z[i-1]+qsum(i-1);
    scanf("%d",&T_T);
    while(T_T--){
        scanf("%d%d",&n,&m);
        ll ans=0,cnt=n*m;
        mem(hang,0); mem(lie,0);
        for(int i=1;i<=n;i++){
            scanf("%s",a[i]+1);
            for(int j=1;j<=m;j++)
                if(a[i][j]=='G'){
                    cnt--;
                    ans += z[n-1]+3*(n-i)*(i-1);
                    ans += z[m-1]+3*(m-j)*(j-1);
                    hang[i]=j;
                    int now=i;
                    if(hang[i-1]){
                        if(j>hang[i-1]){
                            while(hang[now-1]&&hang[now]>hang[now-1]) ans+=2*(m-j)*(hang[--now]-1);
                        }else{
                            while(hang[now-1]&&hang[now]<hang[now-1]) ans+=2*(m-hang[--now])*(j-1);
                        }
                    }
                    break;
                }
                else if(j==m)  ans+=z[m];
        }
        for(int j=1;j<=m;j++)
            for(int i=1;i<=n;i++)
                if(a[i][j]=='G') {
                    lie[j]=i;
                    int now=j;
                    if(lie[j-1]){
                        if(i>lie[j-1]){
                            while(lie[now-1]&&lie[now]>lie[now-1]) ans+=2*(n-i)*(lie[--now]-1);
                        }else{
                            while(lie[now-1]&&lie[now]<lie[now-1]) ans+=2*(n-lie[--now])*(i-1);
                        }
                    }
                    break;
                }
                else if(i==n) ans += z[n];
        cnt=cnt*cnt;
        ans+=fuck();
        printf("%.4f\n",(double)ans*2/cnt);
    }
    return 0;
}