#596716 (C++17) No.1326 ふたりのDominator

提出ソース

結果

問題	No.1326 ふたりのDominator
ユーザー	kmjp
提出日時	2020-12-23 02:01:00
言語	C++17 (gcc 12.3.0 + boost 1.83.0)
結果	AC
実行時間	808 ms / 2,000 ms
コード長	4,793 bytes
コンパイル時間	3,545 ms
コンパイル使用メモリ	247,848 KB
実行使用メモリ	117,488 KB
最終ジャッジ日時	2024-09-23 03:02:01
合計ジャッジ時間	13,165 ms
ジャッジサーバーID （参考情報）	judge5 / judge3

このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示

入力	結果	実行時間実行使用メモリ
testcase_00	AC	22 ms 63,452 KB
testcase_01	AC	20 ms 59,572 KB
testcase_02	AC	21 ms 63,528 KB
testcase_03	AC	21 ms 62,112 KB
testcase_04	AC	22 ms 61,860 KB
testcase_05	AC	20 ms 59,352 KB
testcase_06	AC	21 ms 63,448 KB
testcase_07	AC	25 ms 62,816 KB
testcase_08	AC	25 ms 64,100 KB
testcase_09	AC	25 ms 64,008 KB
testcase_10	AC	25 ms 60,292 KB
testcase_11	AC	25 ms 62,560 KB
testcase_12	AC	638 ms 93,088 KB
testcase_13	AC	634 ms 92,928 KB
testcase_14	AC	647 ms 93,380 KB
testcase_15	AC	639 ms 93,076 KB
testcase_16	AC	634 ms 92,868 KB
testcase_17	AC	605 ms 92,888 KB
testcase_18	AC	678 ms 101,428 KB
testcase_19	AC	416 ms 90,816 KB
testcase_20	AC	537 ms 101,180 KB
testcase_21	AC	599 ms 101,152 KB
testcase_22	AC	808 ms 117,488 KB
testcase_23	AC	522 ms 90,376 KB
testcase_24	AC	610 ms 93,360 KB

権限があれば一括ダウンロードができます

ソースコード

raw source code

#include <bits/stdc++.h>
using namespace std;
typedef signed long long ll;

#define _P(...) (void)printf(__VA_ARGS__)
#define FOR(x,to) for(x=0;x<(to);x++)
#define FORR(x,arr) for(auto& x:arr)
#define FORR2(x,y,arr) for(auto& [x,y]:arr)
#define ALL(a) (a.begin()),(a.end())
#define ZERO(a) memset(a,0,sizeof(a))
#define MINUS(a) memset(a,0xff,sizeof(a))
template<class T> bool chmax(T &a, const T &b) { if(a<b){a=b;return 1;}return 0;}
template<class T> bool chmin(T &a, const T &b) { if(a>b){a=b;return 1;}return 0;}
//-------------------------------------------------------


template<int um> class UF {
	public:
	vector<int> par,rank,cnt;
	UF() {par=rank=vector<int>(um,0); cnt=vector<int>(um,1); for(int i=0;i<um;i++) par[i]=i;}
	void reinit() {int i; FOR(i,um) rank[i]=0,cnt[i]=1,par[i]=i;}
	int operator[](int x) {return (par[x]==x)?(x):(par[x] = operator[](par[x]));}
	int count(int x) { return cnt[operator[](x)];}
	int operator()(int x,int y) {
		if((x=operator[](x))==(y=operator[](y))) return x;
		cnt[y]=cnt[x]=cnt[x]+cnt[y];
		if(rank[x]>rank[y]) return par[x]=y;
		rank[x]+=rank[x]==rank[y]; return par[y]=x;
	}
};

class SCC_ART {
public:
	static const int MV = 110000;
	int NV,time;
	vector<vector<int> > E;
	vector<int> low,num,isart,gr;
	vector<pair<int,int>> ES;
	vector<vector<int> > SC; // out
	vector<int> ART; // out
	vector<set<pair<int,int>> > BE; // out, Edgeset
	UF<MV> uf;
	
	void init(int NV=MV) { this->NV=NV; E.clear(); E.resize(NV);}
	void add_edge(int x,int y) { E[x].push_back(y); E[y].push_back(x);}
	void dfs(int cur,int pre) {
		low[cur]=num[cur]=++time;
		FORR(e,E[cur]) {
			if(e!=pre && num[e]<num[cur]) ES.push_back({min(e,cur),max(e,cur)});
			if(num[e]) low[cur]=min(low[cur],num[e]);
			else {
				dfs(e,cur);
				low[cur]=min(low[cur],low[e]);
				if((num[cur]==1&&num[e]>2) || (num[cur]!=1&&low[e]>=num[cur])) isart[cur]=1;
				if(low[e]>=num[cur]) {
					BE.push_back(set<pair<int,int>>());
					while(1) {
						pair<int,int> te=ES.back();
						ES.pop_back();
						BE.back().insert(te);
						if(min(cur,e)==te.first && max(cur,e)==te.second) break;
					}
					
				}
			}
		}
		if(isart[cur]) ART.push_back(cur);
	}
	void scc() {
		ART.clear();uf.reinit();SC.clear();BE.clear();ES.clear();
		low=num=isart=gr=vector<int>(NV,0);
		for(int i=0;i<NV;i++) if(!num[i]) time=0,dfs(i,-1);
		sort(ART.begin(),ART.end());
		for(int i=0;i<NV;i++) FORR(r,E[i]) if(isart[i]==0 && isart[r]==0) uf(i,r);
		for(int i=0;i<NV;i++) if(uf[i]==i) gr[i]=SC.size(), SC.push_back(vector<int>());
		for(int i=0;i<NV;i++) gr[i]=gr[uf[i]], SC[gr[uf[i]]].push_back(i);
	}
};

int N,M,Q;
SCC_ART scc;
map<pair<int,int>,int> ME;
int A[101010],B[101010];
vector<int> E[303030];
set<int> NE[303030];
UF<303030> uf;
int isart[301010];
int isarts[301010];
vector<int> V[303030];


int P[21][300005],D[300005];

void dfs(int cur,int ps) {
	ps+=isart[cur];
	isarts[cur]=ps;
	FORR(e,NE[cur]) if(e!=P[0][cur]) D[e]=D[cur]+1, P[0][e]=cur, dfs(e,ps);
}
int getpar(int cur,int up) {
	int i;
	FOR(i,20) if(up&(1<<i)) cur=P[i][cur];
	return cur;
}

int lca(int a,int b) {
	int ret=0,i,aa=a,bb=b;
	if(D[aa]>D[bb]) swap(aa,bb);
	for(i=19;i>=0;i--) if(D[bb]-D[aa]>=1<<i) bb=P[i][bb];
	for(i=19;i>=0;i--) if(P[i][aa]!=P[i][bb]) aa=P[i][aa], bb=P[i][bb];
	return (aa==bb)?aa:P[0][aa];               // vertex
}


void solve() {
	int i,j,k,l,r,x,y; string s;
	
	cin>>N>>M;
	scc.init(N);
	FOR(i,M) {
		cin>>A[i]>>B[i];
		A[i]--;
		B[i]--;
		ME[{A[i],B[i]}]=ME[{B[i],A[i]}]=i;
		scc.add_edge(A[i],B[i]);
		E[A[i]].push_back(N+2*i);
		E[N+2*i].push_back(A[i]);
		E[B[i]].push_back(N+2*i+1);
		E[N+2*i+1].push_back(B[i]);
		E[N+2*i].push_back(N+2*i+1);
		E[N+2*i+1].push_back(N+2*i);
	}
	scc.scc();
	FORR(a,scc.ART) isart[a]=1;
	FOR(i,N) if(isart[i]==0) {
		FORR(e,E[i]) uf(i,e);
	}
	FOR(i,N) if(isart[i]) assert(uf.count(i)==1);
	
	FORR(e,scc.BE) {
		int pre=-1;
		FORR(s,e) {
			x=s.first;
			y=s.second;
			i=ME[{x,y}];
			if(A[i]==y) swap(x,y);
			uf(N+i*2,N+i*2+1);
			if(pre>=0) uf(pre,N+i*2);
			pre=N+i*2;
		}
	}
	FOR(i,N+2*M) {
		FORR(e,E[i]) if(uf[i]!=uf[e]) {
			NE[uf[e]].insert(uf[i]);
			NE[uf[i]].insert(uf[e]);
			
		}
	}
	dfs(uf[0],0);
	FOR(i,19) FOR(x,N+M*2) P[i+1][x]=P[i][P[i][x]];
	cin>>Q;
	while(Q--) {
		cin>>x>>y;
		x=uf[x-1];
		y=uf[y-1];
		
		if(x==y) {
			cout<<0<<endl;
			continue;
		}
		int lc=lca(x,y);
		if(lc==x||lc==y) {
			if(lc==y) swap(x,y);
			int ret=isarts[y]-isarts[x]-isart[y];
			cout<<ret<<endl;
		}
		else {
			int ret=isarts[x]-isarts[lc]+isarts[y]-isarts[lc]-isart[x]-isart[y]+isart[lc];
			cout<<ret<<endl;
		}
		
	}
	
}


int main(int argc,char** argv){
	string s;int i;
	if(argc==1) ios::sync_with_stdio(false), cin.tie(0);
	FOR(i,argc-1) s+=argv[i+1],s+='\n'; FOR(i,s.size()) ungetc(s[s.size()-1-i],stdin);
	cout.tie(0); solve(); return 0;
}

yukicoder

結果

テストケース

ソースコード