AtCoder の問題を解いたのですがsegmentation faultの原因がわかりません。

https://atcoder.jp/contests/agc033/tasks/agc033_a
にある問題を解こうとしたのですが以下のコードを実行するとsegmentation fault
となりました。原因が分かれば教えていただきたいです。

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef long double ld;
typedef pair<int, int> P;
const ll mod = 1000000007;
const ll INF = 1e+14;
#define rep(i,n) for(int i=0;i<n;i++)
#define per(i,n) for(int i=n-1;i>=0;i--)
#define Rep(i,sta,n) for(int i=sta;i<n;i++)
#define rep1(i,n) for(int i=1;i<=n;i++)
#define per1(i,n) for(int i=n;i>=1;i--)
#define Rep1(i,sta,n) for(int i=sta;i<=n;i++)
#define _GLIBCXX_DEBUG

int main(){
    int h,w; cin>>h>>w;
    vector<string> field(h);
    rep(i,h){
        cin>>field[i];
    }

    ll dist1[1010][1010];
    rep(i,1010){
        rep(j,1010){
            dist1[i][j]=INF;
        }
    }

    rep(i,h){
        rep(j,w){
            if(field[i][j]=='.') continue;
            else {
                ll dist[1010][1010];
                rep(k,1010){
                    rep(l,1010){
                        dist[k][l]=INF;
                    }
                }
                queue<P> que;
                //初期条件
                dist[i][j]=0;
                que.push({i,j});
                int dx[4]={-1,1,0,0};
                int dy[4]={0,0,-1,1};
                while (que.empty()){
                    P v=que.front();
                    que.pop();
                    rep(k,4){
                        int nh=i+dx[k];
                        int nw=j+dy[k];
                        if(nh<0||nh>=h||nw<0||nw>=w) continue;
                        if(dist[nh][nw]!=INF) continue;
                        if(field[nh][nw]=='#') continue;
                        dist[nh][nw]=dist[v.first][v.second]+1;
                        dist1[nh][nw]=min(dist1[nh][nw],dist[nh][nw]);
                        que.push({nh,nw});
                    }

                }

            }
        }
    }

    ll ans=0;
    rep(i,h){
        rep(j,w){
            if(dist1[i][j]==INF) continue;
            else{
                ans=max(ans,dist1[i][j]);
            }
        }
    }
    cout<<ans<<endl;



}