Submission #354561

#	Time	Username	Problem	Language	Result	Execution time	Memory
354561		arnold518	Toll (APIO13_toll)	C++14	78 / 100	308 ms	16224 KiB

toll

#pragma GCC optimize ("Ofast")
#pragma GCC optimize ("O3")
#pragma GCC optimize ("unroll-loops")
 
#include <bits/stdc++.h>
using namespace std;
 
typedef long long ll;
typedef pair<int, int> pii;
typedef pair<ll, ll> pll;
 
const int MAXN = 1e5;
const int MAXM = 3e5;
const int MAXK = 20;
 
int N, M, K;
 
struct Edge
{
	int u, v, w;
};
Edge E[MAXM+10];
pii A[MAXK+10];
ll B[MAXN+10];
 
struct UF
{
	int par[MAXN+10], sz[MAXN+10];
	void init()
	{
		memset(par, -1, (N+1)<<2);
		memset(sz, 0, (N+1)<<2);
	}
	inline int Find(int x)
	{
		if(-1==par[x]) return x;
		return par[x]=Find(par[x]);
	}
	void Union(int x, int y)
	{
		x=Find(x); y=Find(y);
		if(sz[x]>sz[y]) swap(x, y);
		if(sz[x]==sz[y]) sz[y]++;
		par[x]=y;
	}
}uf, uf2;
 
ll ans, val=0;
vector<pii> adj2[MAXN+10];
vector<int> adj[MAXN+10];
int root[MAXN+10];
 
ll sz[MAXN+10];
vector<int> ST;
 
int col[MAXN+10];
ll D[MAXN+10];

int par2[MAXN+10];
int dep2[MAXN+10];
int Q[MAXN+10];
 
int main()
{
	scanf("%d%d%d", &N, &M, &K);
	for(register int i=1; i<=M; i++) scanf("%d%d%d", &E[i].u, &E[i].v, &E[i].w);
	for(register int i=0; i<K; i++) scanf("%d%d", &A[i].first, &A[i].second);
	for(register int i=1; i<=N; i++) scanf("%d", &B[i]);
 
	sort(E+1, E+M+1, [&](const Edge &p, const Edge &q)
	{
		return p.w<q.w;
	});
 
	uf.init();
	vector<Edge> EE;
	for(register int i=1; i<=M; i++)
	{
		if(uf.Find(E[i].u)==uf.Find(E[i].v)) continue;
		uf.Union(E[i].u, E[i].v);
		EE.push_back(E[i]);
	}
	for(register int i=0; i<EE.size(); i++) E[i+1]=EE[i];
	M=EE.size();
	EE.clear();
	assert(M==N-1);
 
	for(register int i=1; i<=M; i++)
	{
		adj2[E[i].u].push_back({E[i].v, E[i].w});
		adj2[E[i].v].push_back({E[i].u, E[i].w});
	}
 
	vector<int> VV;
	VV.push_back(1);
	for(register int i=0; i<K; i++)
	{
		VV.push_back(A[i].first);
		VV.push_back(A[i].second);
	}
	sort(VV.begin(), VV.end());
	VV.erase(unique(VV.begin(), VV.end()), VV.end());
 
	sort(VV.begin(), VV.end());
 
	uf.init();
	for(auto it : VV) col[it]=it;
 
	for(register int i=1; i<=N; i++) D[i]=B[i];
	for(register int i=1; i<=M; i++)
	{
		if(col[uf.Find(E[i].u)] && col[uf.Find(E[i].v)])
		{
			EE.push_back(E[i]);
		}
		else
		{
			int p=uf.Find(E[i].u), q=uf.Find(E[i].v);
			D[p]+=D[q];
			uf.par[q]=p;
			col[p]|=col[q];
		}
	}
 
	for(auto &it : EE)
	{
		it.u=lower_bound(VV.begin(), VV.end(), col[uf.Find(it.u)])-VV.begin()+1;
		it.v=lower_bound(VV.begin(), VV.end(), col[uf.Find(it.v)])-VV.begin()+1;
	}
	for(register int i=0; i<EE.size(); i++) E[i+1]=EE[i];
	M=EE.size();
 
	for(register int i=0; i<K; i++)
	{
		A[i].first=lower_bound(VV.begin(), VV.end(), A[i].first)-VV.begin()+1;
		A[i].second=lower_bound(VV.begin(), VV.end(), A[i].second)-VV.begin()+1;
	}
 
	for(register int i=0; i<VV.size(); i++)
	{
		B[i+1]=D[uf.Find(VV[i])];
	}
	N=VV.size();
/*
	printf("VV\n");
	for(auto it : VV) printf("%d ", it); printf("\n");
	printf("E\n");
	for(register int i=1; i<=M; i++) printf("%d %d %d\n", E[i].u, E[i].v, E[i].w);
	printf("B\n");
	for(register int i=1; i<=N; i++) printf("%lld ", B[i]); printf("\n");
	printf("A\n");
	for(register int i=0; i<K; i++) printf("%d %d\n", A[i].first, A[i].second);
*/
 
	assert(N<=100);
	for(register int i=1; i<(1<<K); i++)
	{
		uf.init();
		bool flag=true;
		for(register int it=0; it<K; it++)
		{
			if(!(i&(1<<it))) continue;
			if(uf.Find(A[it].first)==uf.Find(A[it].second))
			{
				flag=false;
				break;
			}
			uf.Union(A[it].first, A[it].second);
			adj[A[it].first].push_back(A[it].second);
			adj[A[it].second].push_back(A[it].first);
		}
 
		if(!flag)
		{
			for(register int j=1; j<=N; j++) adj[j].clear();
			continue;
		}
 
		uf2.init();
		for(register int j=1; j<=N; j++) root[j]=j;
 
		for(register int j=1; j<=M; j++)
		{
			if(uf.Find(E[j].u)==uf.Find(E[j].v))
			{
				Q[j]=1;
			}
			else
			{
				Q[j]=0;
				uf.Union(E[j].u, E[j].v);
				adj[E[j].u].push_back(E[j].v);
				adj[E[j].v].push_back(E[j].u);
				//printf("!%d %d %d\n", E[j].u, E[j].v, 0);
			}
		}
	 	
 
		vector<pii> ST1;
		vector<int> ST2;
		ST1.reserve(N);
		ST2.reserve(N);
		ST1.push_back({1, 1});
		ST2.push_back(0);
		for(int i=1; i<=N; i++) sz[i]=B[i];
		dep2[1]=1;
		while(!ST1.empty())
		{
			int now=ST1.back().first, bef=ST1.back().second;
			int j=ST2.back();
			par2[now]=bef;
 
			if(j>=adj[now].size())
			{
				for(auto nxt : adj[now])
				{
					if(nxt==bef) continue;
					sz[now]+=sz[nxt];
				}
				ST1.pop_back();
				ST2.pop_back();
				continue;
			}
			int nxt=adj[now][j];
			if(nxt==bef) ST2[ST2.size()-1]++;
			else
			{
				dep2[nxt]=dep2[now]+1;
				ST2[ST2.size()-1]++;
				ST1.push_back({nxt, now});
				ST2.push_back(0);
			}
		}
 
		val=0;
		for(register int j=1; j<=M; j++)
		{
			if(Q[j])
			{
				int u=E[j].u, v=E[j].v;
				u=root[uf2.Find(u)]; v=root[uf2.Find(v)];
				while(u!=v)
				{
					if(dep2[u]>dep2[v]) swap(u, v);
					val+=sz[v]*E[j].w;
					int t2=root[uf2.Find(par2[v])];
					uf2.Union(v, par2[v]);
					root[uf2.Find(v)]=t2;
					v=t2;
				}
			}
			else
			{
				int uu=root[uf2.Find(E[j].u)], vv=root[uf2.Find(E[j].v)];
				if(dep2[uu]>dep2[vv]) swap(uu, vv);
				uf2.Union(vv, uu);
				root[uf2.Find(vv)]=uu;
			}
		}
 
		ans=max(ans, val);
 
		for(register int j=1; j<=N; j++) vector<int>().swap(adj[j]);
		if(K>15) break;
	}
	printf("%lld\n", ans);
}

Compilation message (stderr)

toll.cpp: In function 'int main()':
toll.cpp:68:43: warning: format '%d' expects argument of type 'int*', but argument 2 has type 'll*' {aka 'long long int*'} [-Wformat=]
   68 |  for(register int i=1; i<=N; i++) scanf("%d", &B[i]);
      |                                          ~^   ~~~~~
      |                                           |   |
      |                                           |   ll* {aka long long int*}
      |                                           int*
      |                                          %lld
toll.cpp:83:25: warning: comparison of integer expressions of different signedness: 'int' and 'std::vector<Edge>::size_type' {aka 'long unsigned int'} [-Wsign-compare]
   83 |  for(register int i=0; i<EE.size(); i++) E[i+1]=EE[i];
      |                        ~^~~~~~~~~~
toll.cpp:130:25: warning: comparison of integer expressions of different signedness: 'int' and 'std::vector<Edge>::size_type' {aka 'long unsigned int'} [-Wsign-compare]
  130 |  for(register int i=0; i<EE.size(); i++) E[i+1]=EE[i];
      |                        ~^~~~~~~~~~
toll.cpp:139:25: warning: comparison of integer expressions of different signedness: 'int' and 'std::vector<int>::size_type' {aka 'long unsigned int'} [-Wsign-compare]
  139 |  for(register int i=0; i<VV.size(); i++)
      |                        ~^~~~~~~~~~
toll.cpp:213:8: warning: comparison of integer expressions of different signedness: 'int' and 'std::vector<int>::size_type' {aka 'long unsigned int'} [-Wsign-compare]
  213 |    if(j>=adj[now].size())
      |       ~^~~~~~~~~~~~~~~~~
toll.cpp:65:7: warning: ignoring return value of 'int scanf(const char*, ...)', declared with attribute warn_unused_result [-Wunused-result]
   65 |  scanf("%d%d%d", &N, &M, &K);
      |  ~~~~~^~~~~~~~~~~~~~~~~~~~~~
toll.cpp:66:40: warning: ignoring return value of 'int scanf(const char*, ...)', declared with attribute warn_unused_result [-Wunused-result]
   66 |  for(register int i=1; i<=M; i++) scanf("%d%d%d", &E[i].u, &E[i].v, &E[i].w);
      |                                   ~~~~~^~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
toll.cpp:67:39: warning: ignoring return value of 'int scanf(const char*, ...)', declared with attribute warn_unused_result [-Wunused-result]
   67 |  for(register int i=0; i<K; i++) scanf("%d%d", &A[i].first, &A[i].second);
      |                                  ~~~~~^~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
toll.cpp:68:40: warning: ignoring return value of 'int scanf(const char*, ...)', declared with attribute warn_unused_result [-Wunused-result]
   68 |  for(register int i=1; i<=N; i++) scanf("%d", &B[i]);
      |                                   ~~~~~^~~~~~~~~~~~~

• Subtask 1: 16 / 16
• Subtask 2: 18 / 18
• Subtask 3: 22 / 22
• Subtask 4: 22 / 22
• Subtask 5: 0 / 22