[ZOJ3781]Paint the Grid Reloaded

思路：

先用DFS缩点，然后BFS找出每个点出发能到达的最长路，取$min$。
注意多组数据，初始化一定要仔细，刚开始存边的$e$忘记初始化，一直WA，调了半个晚上。
一开始和网上的题解对拍$T=1$一直对不出错，后来测了大数据才发现是初始化的问题。

#include<queue>
#include<vector>
#include<cstring>
#include<iostream>
const int inf=0x7fffffff;
const int dx[]={-1,1,0,0},dy[]={0,0,-1,1};
int n,m;
const int N=42,V=1601;
bool a[N][N],b[N][N];
int g[N][N];
bool map[V][V];
int cnt,ans;
std::vector<int> e[V];
inline void init() {
    cnt=0,ans=inf;
    memset(a,0,sizeof a);
    memset(g,0,sizeof g);
    memset(map,0,sizeof map);
    memset(b,0,sizeof b);
    for(int i=1;i<=n;i++) {
        b[i][0]=b[i][m+1]=true;
    }
    for(int j=1;j<=m;j++) {
        b[0][j]=b[n+1][j]=true;
    }
    for(int i=0;i<V;i++) e[i].clear();
}
inline void add_edge(const int u,const int v) {
    if(!map[u][v]) {
        e[u].push_back(v);
        map[u][v]=true;
    }
}
void dfs(const int x,const int y) {
    b[x][y]=true;
    for(int i=0;i<4;i++) {
        int xx=x+dx[i],yy=y+dy[i];
        if(g[xx][yy]&&g[xx][yy]!=g[x][y]) {
            add_edge(g[x][y],g[xx][yy]);
            add_edge(g[xx][yy],g[x][y]);
        }
        if(b[xx][yy]) continue;
        if(a[x][y]==a[xx][yy]) {
            g[xx][yy]=g[x][y];
            dfs(xx,yy);
        }
    }
}
bool v[V];
void bfs(const int s) {
    memset(v,0,sizeof v);
    std::queue<std::pair<int,int> > q;
    q.push(std::make_pair(s,0));
    v[s]=true;
    int max=0;
    while(!q.empty()) {
        int x=q.front().first,d=q.front().second;
        q.pop();
        for(unsigned i=0;i<e[x].size();i++) {
            int &y=e[x][i];
            if(v[y]) continue;
            q.push(std::make_pair(y,d+1));
            v[y]=true;
        }
        max=std::max(max,d);
    }
    ans=std::min(ans,max);
}
int main() {
    std::ios_base::sync_with_stdio(false);
    std::cin.tie(NULL);
    int T;
    std::cin>>T;
    while(T--) {
        std::cin>>n>>m;
        init();
        for(int i=1;i<=n;i++) {
            for(int j=1;j<=m;j++) {
                char ch;
                std::cin>>ch;
                a[i][j]=ch=='O';
            }
        }
        for(int i=1;i<=n;i++) {
            for(int j=1;j<=m;j++) {
                if(b[i][j]) continue;
                g[i][j]=++cnt;
                dfs(i,j);
            }
        }
        for(int i=1;i<=cnt;i++) bfs(i);
        std::cout<<ans<<std::endl;
    }
    return 0;
}