#include #include #include #include #include using namespace atcoder; using namespace std; using ll = long long; using ull = unsigned long long; const ull MOD = 1000000007; using mll = static_modint<1000000007>; #define rep(i,n) for(int i=0; i<(n); i++) struct DSU{ vector V; DSU(int n){ V.resize(n); rep(i,n) V[i] = i; } int leader(int a){ if(V[a] == a) return a; return V[a] = leader(V[a]); } void merge(int r,int c){ V[leader(c)] = leader(r); } }; struct Edge{ int u,v,i; }; int N, M; vector J; vector> E; vector flows; vector cost; vector P; vector> P2; vector Eidx; vector depi; vector dep; vector iEidx; vector flowIdx; void read_graph(){ cin >> N >> M; E.resize(N); cost.resize(M); cost[0] = 2; for(int i=1; i> u >> v; u--; v--; J.push_back({u,v,i}); if(G1.leader(u) == G1.leader(v)){ flows.push_back({u,v,i}); } else{ G1.merge(u,v); E[u].push_back({u,v,i}); E[v].push_back({v,u,i}); } } } void initLCA(){ vector I; P.assign(N,-1); I.push_back(0); depi.assign(N,0); dep.assign(N,0); Eidx.assign(N,-1); rep(i,I.size()){ int p = I[i]; for(Edge e : E[p]) if(P[p] != e.v){ P[e.v] = p; dep[e.v] = dep[p] + cost[e.i]; depi[e.v] = depi[p] + 1; Eidx[e.v] = e.i; I.push_back(e.v); } } P2.assign(20,vector(N)); rep(i,N) P2[0][i] = P[i]; P2[0][0] = 0; rep(d,19) rep(i,N) P2[d+1][i] = P2[d][P2[d][i]]; } int LCA(int u,int v){ if(depi[u] < depi[v]) swap(u,v); for(int i=19; i>=0; i--) if(depi[u]-depi[v] >= (1<=0; i--) if(P2[i][u] != P2[i][v]){ u = P2[i][u]; v = P2[i][v]; } return P2[0][u]; } int dist0(int u,int v){ int g = LCA(u,v); return depi[u] + depi[v] - depi[g] * 2; } mll dist1(int u,int v){ int g = LCA(u,v); return dep[u] + dep[v] - dep[g] * 2; } void build_flows(){ iEidx.assign(M,-1); rep(i,N) if(Eidx[i] != -1) iEidx[Eidx[i]] = i; flowIdx.assign(N,-1); DSU G2(N); for(Edge e : flows){ int g = G2.leader(LCA(e.u,e.v)); for(int s : {e.u,e.v}){ int p = G2.leader(s); while(p != g){ flowIdx[p] = e.i; G2.merge(P[p],p); p = G2.leader(p); } } } } int path_includes(int u,int v,int z){ if(iEidx[z] == -1) return -1; if((LCA(u,v)^LCA(u,J[z].u)^LCA(v,J[z].u)) != J[z].u) return -1; if((LCA(u,v)^LCA(u,J[z].v)^LCA(v,J[z].v)) != J[z].v) return -1; int res = flowIdx[iEidx[z]]; if(res == -1) return -2; return res; } mll shortest_using(int u,int v,int z){ mll ans = cost[z]; //cout << "u = " << u << ", v = " << v << ", z = " << z << endl; //cout << "e = " << J[z].u << ", " << J[z].v << endl; //cout << "d1 = " << (dist0(u,J[z].u) + dist0(v,J[z].v)) << ", d2 = " << (dist0(u,J[z].v) + dist0(v,J[z].u)) << endl; if(dist0(u,J[z].u) + dist0(v,J[z].v) < dist0(u,J[z].v) + dist0(v,J[z].u)) ans += dist1(u,J[z].u) + dist1(v,J[z].v); else ans += dist1(u,J[z].v) + dist1(v,J[z].u); return ans; } int main(){ read_graph(); initLCA(); build_flows(); int Q; cin >> Q; rep(q,Q){ int u,v,z; cin >> u >> v >> z; u--; v--; z--; int inc = path_includes(u,v,z); if(inc == -2){ cout << "-1\n"; continue; } if(inc == -1){ cout << dist1(u,v).val() << "\n"; continue; } auto e = J[inc]; mll ans = shortest_using(u,v,e.i); cout << ans.val() << "\n"; } return 0; } struct ios_do_not_sync{ ios_do_not_sync(){ ios::sync_with_stdio(false); cin.tie(nullptr); } } ios_do_not_sync_inst;