#543422 (C++14) No.1215 都市消滅ビーム

提出ソース

結果

問題	No.1215 都市消滅ビーム
ユーザー	autumn-eel
提出日時	2020-08-30 08:28:28
言語	C++14 (gcc 12.3.0 + boost 1.83.0)
結果	AC
実行時間	3,217 ms / 6,000 ms
コード長	3,909 bytes
コンパイル時間	2,565 ms
コンパイル使用メモリ	194,684 KB
実行使用メモリ	29,084 KB
最終ジャッジ日時	2024-11-15 06:12:27
合計ジャッジ時間	35,773 ms
ジャッジサーバーID （参考情報）	judge5 / judge2

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

テストケース

テストケース表示

入力	結果	実行時間実行使用メモリ
testcase_00	AC	2 ms 5,248 KB
testcase_01	AC	2 ms 5,248 KB
testcase_02	AC	2 ms 5,248 KB
testcase_03	AC	2 ms 5,248 KB
testcase_04	AC	2 ms 5,248 KB
testcase_05	AC	2 ms 5,248 KB
testcase_06	AC	2 ms 5,248 KB
testcase_07	AC	3 ms 5,248 KB
testcase_08	AC	2 ms 5,248 KB
testcase_09	AC	2 ms 5,248 KB
testcase_10	AC	2 ms 5,248 KB
testcase_11	AC	3 ms 5,248 KB
testcase_12	AC	2 ms 5,248 KB
testcase_13	AC	571 ms 11,620 KB
testcase_14	AC	722 ms 10,680 KB
testcase_15	AC	1,774 ms 20,164 KB
testcase_16	AC	511 ms 9,324 KB
testcase_17	AC	740 ms 19,160 KB
testcase_18	AC	792 ms 15,388 KB
testcase_19	AC	76 ms 5,248 KB
testcase_20	AC	1,459 ms 17,708 KB
testcase_21	AC	22 ms 5,248 KB
testcase_22	AC	1,156 ms 14,308 KB
testcase_23	AC	1,072 ms 17,888 KB
testcase_24	AC	901 ms 12,968 KB
testcase_25	AC	568 ms 12,408 KB
testcase_26	AC	1,116 ms 18,592 KB
testcase_27	AC	1,225 ms 23,168 KB
testcase_28	AC	746 ms 10,564 KB
testcase_29	AC	656 ms 10,928 KB
testcase_30	AC	353 ms 19,816 KB
testcase_31	AC	278 ms 8,856 KB
testcase_32	AC	66 ms 18,032 KB
testcase_33	AC	74 ms 15,428 KB
testcase_34	AC	3,141 ms 28,932 KB
testcase_35	AC	3,182 ms 28,932 KB
testcase_36	AC	3,118 ms 28,932 KB
testcase_37	AC	2,595 ms 29,084 KB
testcase_38	AC	1,770 ms 28,956 KB
testcase_39	AC	3,217 ms 28,892 KB
testcase_40	AC	2 ms 5,248 KB
testcase_41	AC	2 ms 5,248 KB

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

ソースコード

raw source code

#include <bits/stdc++.h>
#define rep(i,n)for(int i=0;i<(n);i++)
using namespace std;
typedef long long ll;
typedef pair<ll,ll>P;

class LCA{
public:
	vector<vector<pair<int,int>>>E;
	vector<vector<int>>par;
	vector<int>d1;
	vector<long long>d2;
	int MAX_LOG=22;
	int n;
	LCA(){}
	LCA(int N){
		n=N;
		E=vector<vector<pair<int,int>>>(n);
		par=vector<vector<int>>(MAX_LOG,vector<int>(n));
		d1=vector<int>(n);
		d2=vector<long long>(n);
	}
	void add_edge(int a,int b,int c){
		E[a].push_back(make_pair(c,b));
		E[b].push_back(make_pair(c,a));
	}
	void add_edge(int a,int b){
		E[a].push_back(make_pair(1,b));
		E[b].push_back(make_pair(1,a));
	}
private:
	void dfs(int u,int p,int a,long long b){
		par[0][u]=p;d1[u]=a;d2[u]=b;
		for(auto&v:E[u]){
			if(v.second!=p)dfs(v.second,u,a+1,b+v.first);
		}
	}
	bool flag=false;
	void init(){
		dfs(0,-1,0,0);
		for(int i=1;i<MAX_LOG;i++)for(int j=0;j<n;j++){
			if(par[i-1][j]==-1)par[i][j]=-1;
			else par[i][j]=par[i-1][par[i-1][j]];
		}
		flag=true;
	}
public:
	int lca(int u,int v){
		if(!flag)init();
		if(d1[u]>d1[v])swap(u,v);
		for(int i=0;i<MAX_LOG;i++){
			if((d1[v]-d1[u])>>i&1)v=par[i][v];
		}
		if(u==v)return u;
		for(int i=MAX_LOG-1;i>=0;i--){
			if(par[i][u]!=par[i][v]){
				u=par[i][u];
				v=par[i][v];
			}
		}
		return par[0][u];
	}
	long long dist(int u,int v){
		int l=lca(u,v);
		return d2[u]+d2[v]-d2[l]*2;
	}
};

template<class T>
class BIT{
public:
	vector<T>bit;
	T cnt;
	BIT(int n){
		bit=vector<T>(n+10);
		cnt=0;
	}
	void add(int k,T x){
		k++;
		while(k<(int)bit.size()){
			bit[k]+=x;
			k+=k&-k;
		}
		cnt+=x;
	}
	T sum(int k){
		k++;
		T res=0;
		while(k){
			res+=bit[k];
			k-=k&-k;
		}
		return res;
	}
	T sum2(int k){
		return cnt-sum(k);
	}
};

LCA lca;

int n,K;
int c[200000];
ll d[200000];

int L[200000],R[200000];//lca from left,right
int b[200000];

ll calc1(ll t){//score only
	ll ans=0;
	BIT<int>bit(K);
	vector<P>query;
	ll S=0;
	rep(i,K-2){
		int dep=lca.d1[lca.lca(L[i],c[K-1])];
		S+=d[i];
		query.push_back(P(t-dep-S,b[i]));
	}
	S=0;
	rep(i,K-1){
		S+=d[i];
		if(S+lca.d1[L[i]]<=t)ans++;
	}
	sort(query.begin(),query.end());
	vector<P>sum;
	ll s=0;
	for(int i=K-1;i>=0;i--){
		s+=d[i];
		sum.push_back(P(s,i));
	}
	sort(sum.begin(),sum.end());
	int g=0;
	for(auto p:query){
		while(g<sum.size()&&sum[g].first<=p.first){
			bit.add(sum[g].second,1);
			g++;
		}
		ans+=bit.sum2(p.second-1);
	}
	return ans;
}

struct st{
	ll a,b,c;
};

ll calc2(ll t){//score+depth
	ll ans=0;
	BIT<int>bit(K);
	vector<st>query;
	query.push_back({t,0,K});
	ll S=0;
	rep(i,K-2){
		S+=d[i];
		query.push_back({t-S,i+2,b[i]});
	}
	sort(query.begin(),query.end(),[](st a,st b){return a.a<b.a;});
	vector<P>sum;
	ll s=0;
	for(int i=K-1;i>=0;i--){
		s+=d[i];
		sum.push_back(P(s+lca.d1[R[i]],i));
	}
	sort(sum.begin(),sum.end());
	int g=0;
	for(auto p:query){
		while(g<sum.size()&&sum[g].first<=p.a){
			bit.add(sum[g].second,1);
			g++;
		}
		ans+=bit.sum(p.c-1)-bit.sum(p.b-1);
	}
	return ans;
}

ll C(ll t){
	ll cnt=0;
	cnt+=calc1(t);
	cnt+=calc2(t);
	//~ cout<<t<<' '<<calc1(t)<<' '<<calc2(t)<<endl;
	return cnt;
}

int main(){
	cin>>n>>K;
	lca=LCA(n);
	rep(i,K){
		scanf("%d",&c[i]);c[i]--;
	}
	rep(i,K)scanf("%lld",&d[i]);
	rep(i,n-1){
		int a,b;scanf("%d%d",&a,&b);a--;b--;
		lca.add_edge(a,b);
	}
	lca.lca(0,1);
	for(int i=K-1;i>=0;i--){
		if(i==K-1)R[i]=c[i];
		else R[i]=lca.lca(R[i+1],c[i]);
	}
	rep(i,K){
		if(i==0)L[i]=c[i];
		else L[i]=lca.lca(L[i-1],c[i]);
	}
	rep(i,K-2){
		int ok=K-1,ng=i+1;
		int node=lca.lca(c[K-1],L[i]);
		while(ok-ng>1){
			int t=(ok+ng)/2;
			if(lca.lca(L[i],R[t])==node)ok=t;
			else ng=t;
		}
		b[i]=ok;
	}
	ll E=K*(ll)(K+1)/2+1;
	//~ cout<<"E="<<E<<endl;
	E=(E+1)/2-1;
	//~ cout<<"E="<<E<<endl;
	ll ok=LLONG_MAX/3,ng=LLONG_MIN/3;
	while(ok-ng>1){
		ll t=(ok+ng)/2;
		if(C(t)>=E)ok=t;
		else ng=t;
	}
	cout<<ok<<endl;
}
/*
3 2
2 3
1 1
1 2
2 3
*/

/*
7 3
4 5 7
1 10 100
1 2
2 3
3 4
3 5
1 6
6 7
*/

yukicoder

結果

テストケース

ソースコード