結果

問題 No.1600 Many Shortest Path Problems
ユーザー chocoruskchocorusk
提出日時 2021-07-15 01:33:18
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 1,211 ms / 4,000 ms
コード長 5,962 bytes
コンパイル時間 3,826 ms
コンパイル使用メモリ 207,688 KB
実行使用メモリ 89,784 KB
最終ジャッジ日時 2024-07-04 01:48:15
合計ジャッジ時間 34,410 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge4
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 9 ms
19,328 KB
testcase_01 AC 9 ms
19,200 KB
testcase_02 AC 10 ms
19,200 KB
testcase_03 AC 10 ms
19,328 KB
testcase_04 AC 1,211 ms
83,976 KB
testcase_05 AC 1,188 ms
84,180 KB
testcase_06 AC 10 ms
19,200 KB
testcase_07 AC 10 ms
19,200 KB
testcase_08 AC 10 ms
19,328 KB
testcase_09 AC 10 ms
19,200 KB
testcase_10 AC 528 ms
26,596 KB
testcase_11 AC 615 ms
29,408 KB
testcase_12 AC 950 ms
44,404 KB
testcase_13 AC 1,064 ms
69,992 KB
testcase_14 AC 1,172 ms
83,976 KB
testcase_15 AC 9 ms
19,200 KB
testcase_16 AC 9 ms
19,328 KB
testcase_17 AC 939 ms
56,972 KB
testcase_18 AC 1,161 ms
84,104 KB
testcase_19 AC 9 ms
19,200 KB
testcase_20 AC 9 ms
19,328 KB
testcase_21 AC 896 ms
56,964 KB
testcase_22 AC 9 ms
19,200 KB
testcase_23 AC 9 ms
19,200 KB
testcase_24 AC 1,158 ms
83,980 KB
testcase_25 AC 9 ms
19,200 KB
testcase_26 AC 9 ms
19,328 KB
testcase_27 AC 9 ms
19,328 KB
testcase_28 AC 9 ms
19,200 KB
testcase_29 AC 802 ms
60,456 KB
testcase_30 AC 836 ms
59,692 KB
testcase_31 AC 896 ms
56,804 KB
testcase_32 AC 861 ms
56,828 KB
testcase_33 AC 9 ms
19,200 KB
testcase_34 AC 10 ms
19,200 KB
testcase_35 AC 846 ms
84,424 KB
testcase_36 AC 803 ms
89,784 KB
testcase_37 AC 10 ms
19,200 KB
testcase_38 AC 895 ms
59,560 KB
testcase_39 AC 9 ms
19,200 KB
testcase_40 AC 890 ms
59,948 KB
testcase_41 AC 709 ms
60,460 KB
testcase_42 AC 690 ms
59,564 KB
testcase_43 AC 778 ms
59,564 KB
testcase_44 AC 793 ms
58,672 KB
testcase_45 AC 805 ms
60,456 KB
testcase_46 AC 856 ms
59,820 KB
testcase_47 AC 891 ms
59,436 KB
testcase_48 AC 837 ms
59,568 KB
testcase_49 AC 9 ms
19,200 KB
testcase_50 AC 9 ms
19,200 KB
testcase_51 AC 9 ms
19,200 KB
testcase_52 AC 9 ms
19,200 KB
testcase_53 AC 9 ms
19,200 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <cstdio>
#include <cstring>
#include <iostream>
#include <string>
#include <cmath>
#include <bitset>
#include <vector>
#include <map>
#include <set>
#include <queue>
#include <deque>
#include <algorithm>
#include <complex>
#include <unordered_map>
#include <unordered_set>
#include <random>
#include <cassert>
#include <fstream>
#include <utility>
#include <functional>
#include <time.h>
#include <stack>
#include <array>
#include <list>
#include <atcoder/all>
#define popcount __builtin_popcount
using namespace std;
using namespace atcoder;
typedef long long ll;
typedef pair<int, int> P;
struct LCA{
    vector<vector<int>> g;
    vector<int> d;
    vector<vector<int>> p;
    int log;
    int n;
    LCA(const vector<vector<int>> &g):n(g.size()), g(g), d(g.size()){
        log=0;
        while(1<<log<=n) log++;
        p.resize(log, vector<int>(n));
    }
    void dfs(int x, int prev){
        for(auto y:g[x]){
            if(y==prev) continue;
            d[y]=d[x]+1;
            p[0][y]=x;
            dfs(y, x);
        }
    }
    void build(){
        dfs(0, -1);
        for(int i=1; i<log; i++){
            for(int j=0; j<n; j++){
                p[i][j]=p[i-1][p[i-1][j]];
            }
        }
    }
    int lca(int a, int b){
        if(d[a]>d[b]) swap(a, b);
        int dd=d[b]-d[a], i=0;
        int a1=a, b1=b;
        while(dd){
            if(dd&1) b1=p[i][b1];
            dd>>=1;
            i++;
        }
        if(a1==b1) return a1;
        for(int j=log-1; j>=0; j--){
            if(p[j][a1]!=p[j][b1]){
                a1=p[j][a1], b1=p[j][b1];
            }
        }
        return p[0][a1];
    }
    int dist(int a, int b){
        return d[a]+d[b]-2*d[lca(a, b)];
    }
};
struct unionfind{
    vector<int> par, sz;
    unionfind() {}
    unionfind(int n):par(n), sz(n, 1){
        for(int i=0; i<n; i++) par[i]=i;
    }
    int find(int x){
        if(par[x]==x) return x;
        return par[x]=find(par[x]);
    }
    void unite(int x, int y){
        x=find(x); y=find(y);
        if(x==y) return;
        if(sz[x]>sz[y]) swap(x, y);
        par[x]=y;
        sz[y]+=sz[x];
    }
    bool same(int x, int y){
        return find(x)==find(y);
    }
    int size(int x){
        return sz[find(x)];
    }
};
template<typename Monoid>
struct SegmentTree{
	using F=function<Monoid(Monoid, Monoid)>;
	int sz;
	vector<Monoid> seg;
	const F f;
	const Monoid e;

	SegmentTree(int n, const F f, const Monoid &e): f(f), e(e){
		sz=1;
		while(sz<n) sz<<=1;
		seg.resize(2*sz, e);
	}

	SegmentTree(int n, const F f, const Monoid &e, vector<Monoid> v): f(f), e(e){
		sz=1;
		while(sz<n) sz<<=1;
		seg.resize(2*sz, e);
		for(int i=0; i<n; i++) seg[i+sz]=v[i];
		for(int i=sz-1; i>=1; i--){
			seg[i]=f(seg[2*i], seg[2*i+1]);
		}
	}

	void update(int k, const Monoid &x){
		k+=sz;
		seg[k]=x;
		while(k>1){
			k>>=1;
			seg[k]=f(seg[2*k], seg[2*k+1]);
		}
	}

	Monoid query(int a, int b){
		a+=sz, b+=sz;
		Monoid ret=e;
		for(;a<b; a>>=1, b>>=1){
			if(b&1) ret=f(ret, seg[--b]);
			if(a&1) ret=f(ret, seg[a++]);
		}
		return ret;
	}

	/*Monoid query(int a, int b){ //演算が非可換のとき
		a+=sz, b+=sz;
		Monoid retl=e, retr=e;
		for(;a<b; a>>=1, b>>=1){
			if(b&1) retr=f(seg[--b], retr);
			if(a&1) retl=f(retl, seg[a++]);
		}
		return f(retl, retr);
	}*/

	Monoid operator[](const int &k) const{
		return seg[k+sz];
	}
};
int n, m;
int a[200020], b[200020];
int q;
int x[200020], y[200020], z[200020];
vector<vector<int>> g;
vector<vector<int>> ge;
vector<int> e;
using mint=modint1000000007;
mint d[200020];
int d0[200020];
mint p2[200020];
//mint ans[200020];
int l[200020], r[200020], rev[200020];
vector<int> v1[200020], v2[200020], vr[200020];
int mnl[200020], mnr[200020];
bool t[200020];
int main()
{
	cin>>n>>m;
	g.resize(n);
	ge.resize(n);
	for(int i=0; i<m; i++){
		cin>>a[i]>>b[i];
		a[i]--; b[i]--;
	}
	unionfind uf(n);
	for(int i=0; i<m; i++){
		if(uf.same(a[i], b[i])){
			e.push_back(i);
		}else{
			g[a[i]].push_back(b[i]);
			g[b[i]].push_back(a[i]);
			ge[a[i]].push_back(i);
			ge[b[i]].push_back(i);
			uf.unite(a[i], b[i]);
			t[i]=1;
		}
	}
	p2[0]=1;
	for(int i=0; i<m; i++) p2[i+1]=p2[i]*2;
	LCA lca(g);
	lca.build();
	int ord=0;
	auto dfs=[&](auto dfs, int x, int p)->void{
		l[x]=ord++;
		rev[l[x]]=x;
		for(int i=0; i<g[x].size(); i++){
			int y=g[x][i];
			if(y==p) continue;
			int j=ge[x][i];
			if(y==a[j]) swap(a[j], b[j]);
			d[y]=d[x]+p2[j+1];
			//d0[y]=d0[x]+1;
			dfs(dfs, y, x);
		}
		r[x]=ord;
		vr[r[x]].push_back(x);
	};
	dfs(dfs, 0, -1);
	const int INF=1e9;
	SegmentTree<int> seg(n, [](int x, int y){ return min(x, y);}, INF);
	for(auto i:e){
		if(l[a[i]]>l[b[i]]) swap(a[i], b[i]);
		v1[l[a[i]]].push_back(i);
		v2[l[b[i]]].push_back(i);
	}
	int mn1[200020];
	fill(mn1, mn1+n, INF);
	for(int i=0; i<n; i++){
		int r1=r[rev[i]];
		mnl[rev[i]]=seg.query(i, r1);
		for(auto j:v1[i]){
			int l1=l[b[j]];
			mn1[l1]=min(mn1[l1], j);
			seg.update(l1, mn1[l1]);
		}
	}
	SegmentTree<int> seg2(n, [](int x, int y){ return min(x, y);}, INF);
	fill(mn1, mn1+n, INF);
	for(int i=n; i>=1; i--){
		for(auto j:v2[i]){
			int l1=l[a[j]];
			mn1[l1]=min(mn1[l1], j);
			seg2.update(l1, mn1[l1]);
		}
		for(auto j:vr[i]){
			mnr[j]=seg2.query(l[j], r[j]);
		}
	}
	auto calc=[&](int x, int y){
		return d[x]+d[y]-2*d[lca.lca(x, y)];
	};
	cin>>q;
	for(int i=0; i<q; i++){
		cin>>x[i]>>y[i]>>z[i];
		x[i]--; y[i]--; z[i]--;
		if(!t[z[i]] || lca.dist(x[i], y[i])!=lca.dist(x[i], a[z[i]])+lca.dist(y[i], a[z[i]]) || lca.dist(x[i], y[i])!=lca.dist(x[i], b[z[i]])+lca.dist(y[i], b[z[i]])){
			mint ans=calc(x[i], y[i]);
			cout<<ans.val()<<endl;
		}else{
			int bz=b[z[i]];
			int k=min(mnl[bz], mnr[bz]);
			if(k==INF){
				cout<<-1<<endl;
				continue;
			}
			int a1=a[k], b1=b[k];
			if(l[bz]<=l[y[i]] && l[y[i]]<r[bz]) swap(x[i], y[i]);
			if(l[bz]<=l[b1] && l[b1]<r[bz]) swap(a1, b1);
			mint ans=calc(x[i], a1)+calc(y[i], b1)+p2[k+1];
			cout<<ans.val()<<endl;
		}
	}

    return 0;
}
0