kmjp's blog

競技プログラミング参加記です

Codeforces #949 : Div2 E. Turtle and Intersected Segments

またなんかややこし気な設定。
https://codeforces.com/contest/1981/problem/E

問題

N個の区間がある。
各区間にはパラメータA[i]が設定されている。

ここでN頂点のグラフを考える。
もし区間iと区間jが共通部分を持つとき、点iと点jの間に、コストの|A[i]-A[j]|の辺があるとする。

このグラフのMSTのコストを求めよ。

解法

区間を端から走査し、区間の集合を保持しておく。
区間の追加・削除のタイミングで、追加・削除対象の区間iと、A[j]の値が近い区間jの間にある辺を追加しよう。
こうするとO(N)本の辺の候補が出るので、あとはMSTを求めるだけ。

int T,N;
int L[505050],R[505050],A[505050];
template<int um> class UF {
	public:
	vector<int> par,rank,cnt,G[um];
	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 num=um) {int i; FOR(i,num) 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;
	}
};
UF<505050> uf;

void solve() {
	int i,j,k,l,r,x,y; string s;
	
	cin>>T;
	while(T--) {
		cin>>N;
		uf.reinit(N);
		map<int,vector<int>> ev;
		FOR(i,N) {
			cin>>L[i]>>R[i]>>A[i];
			L[i]=L[i]*2;
			R[i]=R[i]*2+1;
			ev[L[i]].push_back(i);
			ev[R[i]].push_back(i);
		}
		vector<pair<ll,pair<int,int>>> Es;
		set<pair<ll,int>> S;
		FORR2(a,b,ev) {
			if(a%2==0) {
				FORR(v,b) S.insert({A[v],v});
				FORR(v,b) {
					auto it=S.lower_bound({A[v],v});
					if(it!=S.begin()) {
						auto p=*prev(it);
						Es.push_back({A[v]-p.first,{v,p.second}});
					}
					it++;
					if(it!=S.end()) {
						auto p=*it;
						Es.push_back({p.first-A[v],{v,p.second}});
					}
				}
			}
			else {
				FORR(v,b) {
					auto it=S.lower_bound({A[v],v});
					if(it!=S.begin()) {
						auto p=*prev(it);
						Es.push_back({A[v]-p.first,{v,p.second}});
					}
					it++;
					if(it!=S.end()) {
						auto p=*it;
						Es.push_back({p.first-A[v],{v,p.second}});
					}
				}
				FORR(v,b) S.erase({A[v],v});
			}
		}
		ll ret=0;
		sort(ALL(Es));
		FORR2(a,v,Es) {
			if(uf[v.first]!=uf[v.second]) {
				uf(v.first,v.second);
				ret+=a;
			}
		}
		if(uf.count(0)!=N) ret=-1;
		cout<<ret<<endl;
	}
		
		
}

まとめ

コードは若干ややこしいけど、考え方はそうでもない。