XVII Open Cup named after E.V. Pankratiev Stage 14, Grand Prix of Tatarstan, Sunday, April 2, 2017 Problem L. Canonical duel

题目：Problem L. Canonical duel
Input file: standard input
Output file: standard output
Time limit: 2 seconds
Memory limit: 256 megabytes
In the game «Canonical duel» board N × M is used. Some of the cells of the board contain turrets. A
turret is the unit with four cannons oriented to north, east, south and west. Each turret can shoot exactly
once. When turret is hit by the cannon of another turret, its activated. When turret is activated, all four
cannons shoot simultaneously, then self-destruction process causes the turret to disappear.
Given the board with some turrets. You may add exactly one turret on one of cells which does not
contains a turret and activate this new turret. Your goal is to destroy maximum number of turrets.
Input
First line of the input contains two positive integers N and M, does not exceeding 2000 — size of the
board.
Each of the next N lines contain exactly M chars: ‘+’ denotes that cell is occupied by a turret, and ‘.’
that cell is empty.
Output
In the first line print maximum number of existing turrets you may destroy, then in second line print two
space-separated integers — row and column of place where turret can be set. If it is impossible to destroy
ever one turret in such a way, print only one line containing a zero; if several solutions exist, print any of
them.
Examples

standard input	standard output
3 4 ++.. +... ..++	5 2 4
4 5 ++... ..+.. ....+ ...++	5 4 1
3 3 +++ +++ +++	0

思路：

　　用并查集维护每个点所能炸的点数量，然后处理出每个空位置上下左右最近的炮台的位置。之后枚举放置位置即可。

　　

#include <bits/stdc++.h>

using namespace std;

#define MP make_pair
#define PB push_back
typedef long long LL;
typedef pair<int,int> PII;
const double eps=1e-8;
const double pi=acos(-1.0);
const int K=1e6+7;
const int mod=1e9+7;

char ss[2005][2005];
int n,m,ans,ansx,ansy,f[4001000],cnt[4001000];
int dir[4001000][5],vis[4001000];
int idx(int i,int j)
{
    if(i<1||j<1||i>n||j>m) return 0;
    return (i-1)*m+j;
}
int fd(int x)
{
    return f[x]==x?x:f[x]=fd(f[x]);
}
void join(int y,int x)
{
    int fx=fd(x),fy=fd(y);
    if(fx!=fy)
        f[fy]=fx,cnt[fx]+=cnt[fy];
}
int main(void)
{
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++)
        scanf("%s",&ss[i][1]);
    for(int i=1;i<=n;i++)
    for(int j=1;j<=m;j++)
        f[idx(i,j)]=idx(i,j),cnt[idx(i,j)]=ss[i][j]=='+';
    for(int i=1;i<=n;i++)
    for(int j=1;j<=m;j++)
    {
        dir[idx(i,j)][1]=ss[i-1][j]=='+'?idx(i-1,j):dir[idx(i-1,j)][1];
        dir[idx(i,j)][3]=ss[i][j-1]=='+'?idx(i,j-1):dir[idx(i,j-1)][3];
    }
    for(int i=n;i;i--)
    for(int j=m;j;j--)
    {
        dir[idx(i,j)][2]=ss[i+1][j]=='+'?idx(i+1,j):dir[idx(i+1,j)][2];
        dir[idx(i,j)][4]=ss[i][j+1]=='+'?idx(i,j+1):dir[idx(i,j+1)][4];
    }
    for(int i=1;i<=n;i++)
    for(int j=1;j<=m;j++)
    if(ss[i][j]=='+')
    {
        if(dir[idx(i,j)][1]) join(idx(i,j),dir[idx(i,j)][1]);
        if(dir[idx(i,j)][3]) join(idx(i,j),dir[idx(i,j)][3]);
        //printf("f[%d][%d]:%d
",i,j,f[idx(i,j)]);
    }
//    for(int i=1;i<=n;i++)
//    for(int j=1;j<=m;j++)
//    {
//        printf("dir[%d][%d]:",i,j);
//        for(int k=1;k<=4;k++)
//        printf("%d ",dir[idx(i,j)][k]);
//        printf("
");
//    }
//    for(int i=1;i<=n;i++)
//    for(int j=1;j<=m;j++)
//    if(f[idx(i,j)]==idx(i,j))
//        printf("cnt[%d][%d]:%d
",i,j,cnt[idx(i,j)]);
    for(int i=1;i<=n;i++)
    for(int j=1;j<=m;j++)
    if(ss[i][j]=='.')
    {
        int sum=0;
        for(int k=1;k<=4;k++)
        if(!vis[fd(dir[idx(i,j)][k])])
            vis[fd(dir[idx(i,j)][k])]=1,sum+=cnt[fd(dir[idx(i,j)][k])];
        for(int k=1;k<=4;k++)
            vis[fd(dir[idx(i,j)][k])]=0;
        //printf("sum:%d %d %d
",i,j,sum);
        if(sum>ans)
            ans=sum,ansx=i,ansy=j;
    }
    if(ans) printf("%d
%d %d
",ans,ansx,ansy);
    else printf("0
");
    return 0;
}