hdu 5725(计数)

最新推荐文章于 2018-07-14 16:05:00 发布

top628LJ

最新推荐文章于 2018-07-14 16:05:00 发布

阅读量660

点赞数

分类专栏： 2016多校

本文链接：https://blog.csdn.net/u013849646/article/details/52004910

版权

2016多校专栏收录该内容

27 篇文章 0 订阅

订阅专栏

题目要求出所有合法点对间的最短路径的平均值，因此我们应当求出所有合法最短点对的最短路径之和，再除以合法点对个数。题目中Guard之间有着很不自然的制约关系，每个Guard的周围和同行、列都不能有其余的 Guard,但不可路过的格子却只有本身一格。由此发现，任意两点间的最短路径只会至多被一个Guard所干扰。因此可以想到，我们可以先算出没有Guard制约的所有最短路径长度之和，再分别对于每个 Guard求出它对最短路径长度的贡献（即原本必须经过它的最短路要多绕两格）。

#include <bits/stdc++.h>

#pragma comment(linker, "/STACK:102400000,102400000")
using namespace std;

#define LL long long
#define pii pair<int,int>
#define MP make_pair
#define ls i << 1
#define rs ls | 1
#define md (ll + rr >> 1)
#define lson ll, md, ls
#define rson md + 1, rr, rs
#define Pi acos(-1.0)
#define mod 1000000007
#define eps 1e-10
#define inf 0x3f3f3f3f
#define N 1010
#define M 800020

LL ans;
void f(int n, int m, int g[]){
    LL sum = 0;
    for(int i = 1; i <= n; ++i){
        if(g[i] > g[i-1]) sum += g[i] - 1;
        else if(g[i] == 0) sum = 0;
        else sum = g[i] - 1;
        ans += 1LL * sum * (m - g[i]) * 4;
    }
}
void g(int n, int m, int g[]){
    LL sum = 0;
    for(int i = n; i >= 1; --i){
        if(g[i] < g[i+1]) sum = 0;
        ans += 1LL * sum * (m - g[i]) * 4;
        if(g[i] > g[i+1]) sum += g[i] - 1;
        else if(g[i] == 0) sum = 0;
        else sum = g[i] - 1;
    }
}
int gx[N], gy[N], x[N], y[N];
char s[N][N];
int main(){
    int cas;
    scanf("%d", &cas);
    while(cas--){
        int n, m;
        scanf("%d%d", &n, &m);
        memset(x, 0, sizeof x);
        memset(y, 0, sizeof y);
        memset(gx, 0, sizeof gx);
        memset(gy, 0, sizeof gy);
        for(int i = 1; i <= n; ++i)
            scanf("%s", s[i] + 1);
        int tot = 0;
        for(int i = 1; i <= n; ++i)
            for(int j = 1; j <= m; ++j){
                if(s[i][j] == 'G')
                    gx[i] = j, gy[j] = i;
                else{
                    tot++;
                    x[i]++, y[j]++;
                }
            }
        ans = 0;
        for(int i = 1; i <= n; ++i)
            for(int j = 1; j <= n; ++j)
                ans += 1LL * abs(i - j) * x[i] * x[j];
        for(int i = 1; i <= m; ++i)
            for(int j = 1; j <= m; ++j)
                ans += 1LL * abs(i - j) * y[i] * y[j];
        f(n, m, gx);
        f(m, n, gy);
        g(n, m, gx);
        g(m, n, gy);
        printf("%.4f\n", 1.0 * ans / tot / tot);
    }
    return 0;
}