QOJ.ac

QOJ

ID题目提交者结果用时内存语言文件大小提交时间测评时间
#60174#4809. Maximum RangexiaoyaowudiWA 77ms89992kbC++179.1kb2022-11-03 10:31:272022-11-03 10:31:28

Judging History

你现在查看的是最新测评结果

  • [2023-08-10 23:21:45]
  • System Update: QOJ starts to keep a history of the judgings of all the submissions.
  • [2022-11-03 10:31:28]
  • 评测
  • 测评结果:WA
  • 用时:77ms
  • 内存:89992kb
  • [2022-11-03 10:31:27]
  • 提交

answer

#include <iostream>
#include <vector>
#include <tuple>
#include <algorithm>
#include <functional>
constexpr int N(400010),inf(1.01e9);
using t3i=std::tuple<int,int,int,int>;
using veci=std::vector<int>;
int bcnt,n,m;
t3i mxm[N],mnm[N];
veci pts[N],ts[N<<1],gs[N];
// #define DEBUG {std::cout<<__FUNCTION__<<std::endl;}
#define DEBUG
std::vector<std::pair<int,int>> eds[N],es[N];
void tarjan(int u)
{
	static std::tuple<int,int,int> stk[N];
	static int top,dfn[N],low[N],dcnt;
	dfn[u]=low[u]=++dcnt;stk[++top]={u,-1,-1};
	for(auto [v,w]:es[u])
	{
		if(!dfn[v])
		{
			tarjan(v);low[u]=std::min(low[u],low[v]);
			if(low[v]==dfn[u])
			{
				int k=++bcnt;
				pts[k].emplace_back(u);
				mxm[k]={-inf,0,0,0},mnm[k]={inf,0,0,0};
				while(1)
				{
					auto [x,y,w]=stk[top--];
					// std::cerr<<x<<" "<<y<<" "<<w<<std::endl;
					if(~y) eds[k].emplace_back(x,y),mxm[k]=std::max(mxm[k],{w,x,y,k+n}),mnm[k]=std::min(mnm[k],{w,x,y,k+n});
					else
					{
						pts[k].emplace_back(x);
						if(x==v) break;
					}
				}
				if(eds[k].size()>1)
				{
					for(int v:pts[k]) ts[k+n].emplace_back(v),ts[v].emplace_back(k+n);
				}
			}
		}
		else if(dfn[v]<dfn[u]) stk[++top]={u,v,w},low[u]=std::min(low[u],dfn[v]);
	}
}
// std::pair<veci,veci> calc(const veci &vs,const veci &vall)
// {
// 	DEBUG;
// 	// if(m==98516)
// 	// {
// 	// 	for(int u:vs) if(std::find(vall.begin(),vall.end(),u)==vall.end())
// 	// 	{
// 	// 		std::cout<<"can't find "<<u<<std::endl;
// 	// 	}
// 	// }
// 	int l=vs.size(),u=vs[0],v=vs[l-1];
// 	static int idx[N],fv[N];
// 	// static bool mk[N];
// 	// for(int i:vall) fu[i]=fv[i]=0,mk[i]=false;
// 	for(int i:vall) fv[i]=0,idx[i]=-100;
// 	for(int i(0);i<l;++i) idx[vs[i]]=i;
// 	std::function<void(int,int,int*)> dfs=[&dfs](int u,int fa,int *ff){ff[u]=fa;if(mk[u]) return;for(int v:gs[u]) if(!ff[v]) dfs(v,u,ff);};
// 	for(int k:gs[u]) if(k==v)
// 	{
// 		return {vs,{u,v}};
// 	}
// 	fv[v]=-1;for(int k:gs[v]) if(!mk[k] && !fv[k]) dfs(k,v,fv);
// 	// if(m==98516) std::cout<<"❤"<<std::endl;
// 	if(fv[u])
// 	{
// 		// if(n==99996 && m==100000) std::cout<<"❤gyc"<<std::endl;
// 		veci res;
// 		for(int i(u);~i;i=fv[i]) res.emplace_back(i);
// 		// if(m==98516)
// 		// {
// 		// 	// std::cout<<"end"<<std::endl;
// 		// 	for(int i(0);i<std::min<int>(5,vs.size());++i) std::cout<<vs[vs.size()-i-1]<<" ";std::cout<<std::endl;
// 		// 	for(int i(0);i<std::min<int>(5,res.size());++i) std::cout<<res[res.size()-i-1]<<" ";std::cout<<std::endl;
// 		// }
// 		return {vs,res};
// 	}
// 	// if(m==98516) std::cout<<"❤"<<std::endl;
// 	int p1(l-2);while(!fv[vs[p1]]) ++p1;
// 	// if(m==98516) std::cout<<p1<<" "<<l<<std::endl;
// 	// if(m==98516) std::cout<<"❤"<<std::endl;
// 	int p2(p1+1);while(!fu[vs[p2]]) ++p2;
// 	// if(m==98516) std::cout<<"❤"<<std::endl;
// 	// if(m==98516) std::cout<<p1<<" "<<p2<<" "<<l<<std::endl;
// 	veci r1;for(int i(vs[p2]);~i;i=fu[i]) r1.emplace_back(i);std::reverse(r1.begin(),r1.end());r1.insert(r1.end(),vs.begin()+p2+1,vs.end());
// 	veci r2;r2.insert(r2.begin(),vs.begin(),vs.begin()+p1);for(int i(vs[p1]);~i;i=fv[i]) r2.emplace_back(i);
// 	// if(m==98516) for(int i(0);i<std::min<int>(5,r1.size());++i) std::cout<<r1[i]<<" ";std::cout<<std::endl;
// 	// if(m==98516) for(int i(0);i<std::min<int>(5,r2.size());++i) std::cout<<r2[i]<<" ";std::cout<<std::endl;
// 	return {r1,r2};
// }
std::pair<veci,veci> calc(const veci &vs,const veci &vall)
{
	DEBUG;
	int l(vs.size()),u(vs[0]),v(vs[l-1]),k(vs[1]);
	static int idx[N];
	for(int i:vall) idx[i]=-1;
	for(int i(0);i<l;++i) idx[vs[i]]=i;
	veci r1({u}),r2({u,k});int p1(u),p2(k);
	static int fa1[N],fa2[N];
	static bool vis1[N],vis2[N];
	for(int i:vall) vis1[i]=vis2[i]=false,fa1[i]=fa2[i]=0;
	int *f1(fa1),*f2(fa2);bool *v1(vis1),*v2(vis2);
	// std::cout<<u<<" "<<v<<" "<<k<<std::endl;
	// for(int i:vs) std::cout<<i<<" ";std::cout<<std::endl;
	while(1)
	{
		// if(n==99983 && m==99998) std::cout<<p1<<" "<<p2<<std::endl;
		if(idx[p1]>idx[p2]) std::swap(p1,p2),std::swap(f1,f2),std::swap(v1,v2),r1.swap(r2);
		if(p2==v)
		{
			for(int i(idx[p1]+1);i<l;++i) r1.emplace_back(vs[i]);
			break;
		}
		int mxi(-1);
		std::function<void(int,int)> dfs=[&](int u,int fa)->void{v1[u]=true;f1[u]=fa;mxi=std::max(mxi,idx[u]);if(idx[u]>=idx[p2]) return;
			for(int t:gs[u]) if(!v1[t]) dfs(t,u);};
		dfs(p1,0);
		int len(r1.size());
		for(int t(vs[mxi]);t!=p1;t=f1[t]) r1.emplace_back(t);
		std::reverse(r1.begin()+len,r1.end());
		p1=vs[mxi];
	}
	for(int k:r1)
	{
		if(k!=u && k!=v && std::find(r2.begin(),r2.end(),k)!=r2.end())
		{
			if(idx[k]!=-1) std::cout<<"G!"<<std::endl;
			std::cout<<"multiple point found"<<std::endl;
		}
	}
	return {r1,r2};
}
veci solve_path(veci pts,int u,int v,int frb,int pre)
{
	DEBUG;
	static int fs[N];
	static bool vis[N];
	fs[pre]=0;
	std::function<void(int,int)> dfs=[&](int u,int fa)->void{fs[u]=fa;vis[u]=true;for(int v:gs[u]) if(!vis[v] && v!=frb) dfs(v,u);};
	dfs(v,pre);
	veci ret;
	for(int k:pts) vis[k]=false;
	for(int i(u);i;i=fs[i]) ret.emplace_back(i);
	// if(n==99996 && m==100000){std::cout<<*ret.rbegin()<<" "<<v<<" "<<pre<<std::endl;}
	return ret;
}
std::pair<veci,veci> solve1(int idx,int u,int v,int s)
{
	if(n==99983 && m==99998) std::cout<<"solve 3"<<std::endl;
	DEBUG;
	idx-=n;
	for(auto [u,v]:eds[idx])
	{
		gs[u].emplace_back(v),gs[v].emplace_back(u);
		if(n==99983 && m==99998) std::cout<<"e "<<u<<" "<<v<<std::endl;
	}
	std::pair<veci,veci> res;
	if(u==s)
	{
		res=calc({s,v},pts[idx]);
	}
	else
	{
		res=calc(solve_path(pts[idx],s,v,u,u),pts[idx]);
	}
	for(int u:pts[idx]) gs[u].clear();
	return res;
}
std::pair<veci,veci> solve2(int idx,int u,int v)
{
	DEBUG;
	// if(m==98516) std::cout<<u<<" "<<v<<std::endl;
	idx-=n;
	for(auto [u,v]:eds[idx]) gs[u].emplace_back(v),gs[v].emplace_back(u);
	auto res=calc(solve_path(pts[idx],u,v,0,0),pts[idx]);
	for(int u:pts[idx]) gs[u].clear();
	return res;
}
veci solve3(int idx,int u,int v,int s,int t)
{
	if(n==99983 && m==99998) std::cout<<"solve 3"<<std::endl;
	DEBUG;
	idx-=n;
	for(auto [u,v]:eds[idx])
	{
		gs[u].emplace_back(v),gs[v].emplace_back(u);
		// if(n==99983 && m==99998) std::cout<<"e "<<u<<" "<<v<<std::endl;
	}
	auto hp=[](std::pair<veci,veci> pr,veci vall)->veci
	{
		auto [p1,p2]=pr;
		std::reverse(p2.begin(),p2.end());
		p1.pop_back();p2.pop_back();
		p1.insert(p1.end(),p2.begin(),p2.end());
		for(int v:vall) gs[v].clear();
		return p1;
	};
	if(v==s || v==t) return hp(calc({u,v,(s^t^v)},pts[idx]),pts[idx]);
	else if(u==s || u==t) return hp(calc({v,u,(s^t^u)},pts[idx]),pts[idx]);
	auto pt=solve_path(pts[idx],v,s,t,t);
	int hu=-1,len=pt.size();
	for(int k(0);k<len;++k) if(pt[k]==u){hu=k;break;}
	// if(n==99997) std::cout<<hu<<std::endl;
	if(~hu)
	{
		pt.erase(pt.begin()+1,pt.begin()+hu);
	}
	else pt.insert(pt.begin(),u);
	// for(int i=0;i<std::min<int>(pt.size(),10);++i) std::cout<<pt[i]<<" ";std::cout<<std::endl;
	return hp(calc(pt,pts[idx]),pts[idx]);
}
int main()
{
	DEBUG;
	std::ios::sync_with_stdio(false);
	std::cin.tie(0);std::cin>>n>>m;
	for(int i(1),u,v,w;i<=m;++i) std::cin>>u>>v>>w,es[u].emplace_back(v,w),es[v].emplace_back(u,w);
	tarjan(1);
	// if(m==98516) std::cout<<"❤"<<std::endl;
	static bool vis[N<<1];
	std::function<void(int,t3i&,t3i&)> dfs1=[&dfs1](int u,t3i &mx,t3i &mn)->void
	{
		if(u>n) mx=std::max(mxm[u-n],mx),mn=std::min(mnm[u-n],mn);vis[u]=true;
		for(int v:ts[u]) if(!vis[v]) dfs1(v,mx,mn);
	};
	int ans=-1,su,sv,sk,tu,tv,tk;
	for(int i(1);i<=n;++i) if(!vis[i] && ts[i].size())
	{
		t3i mn{inf,0,0,0},mx{-inf,0,0,0};
		dfs1(i,mx,mn);
		auto [w1,e1u,e1v,v1]=mn;auto [w2,e2u,e2v,v2]=mx;
		if(w2-w1>ans)
		{
			ans=w2-w1;
			su=e1u;sv=e1v;sk=v1;
			tu=e2u;tv=e2v;tk=v2;
		}
	}
	std::cout<<ans<<std::endl;
	if(sk!=tk)
	{
		static int ff[N<<1];
		std::function<void(int,int)> dfs2=[&dfs2](int u,int fa)->void{ff[u]=fa;for(int v:ts[u])if(fa!=v) dfs2(v,u);};
		dfs2(tk,0);
		static int vs[N<<1],vc;
		for(int i(sk);i;i=ff[i]) vs[++vc]=i;
		auto [r1,r2]=solve1(sk,su,sv,vs[2]);
		std::reverse(r1.begin(),r1.end());std::reverse(r2.begin(),r2.end());
		for(int i=3;i<vc;i+=2)
		{
			auto [c1,c2]=solve2(vs[i],vs[i-1],vs[i+1]);
			c1.erase(c1.begin());c2.erase(c2.begin());
			r1.insert(r1.end(),c1.begin(),c1.end());
			r2.insert(r2.end(),c2.begin(),c2.end());
		}
		auto [c1,c2]=solve1(tk,tu,tv,vs[vc-1]);
		// std::reverse(c1.begin(),c1.end());std::reverse(c2.begin(),c2.end());
		c1.erase(c1.begin());c2.erase(c2.begin());
		DEBUG;
		r1.insert(r1.end(),c1.begin(),c1.end());
		r2.insert(r2.end(),c2.begin(),c2.end());
		std::reverse(r2.begin(),r2.end());r1.erase(--r1.end());r2.erase(--r2.end());
		// if(n==99996 && m==100000){
		// 	for(int i(0);i<std::min<int>(5,r1.size());++i) std::cout<<r1[i]<<" ";std::cout<<std::endl;
		// 	for(int i(0);i<std::min<int>(5,r2.size());++i) std::cout<<r2[i]<<" ";std::cout<<std::endl;
		// }
		r1.insert(r1.end(),r2.begin(),r2.end());
		std::cout<<r1.size()<<std::endl;
		for(int v: r1) std::cout<<v<<" ";std::cout<<std::endl;
	}
	else
	{
		auto r1=solve3(sk,su,sv,tu,tv);
		std::cout<<r1.size()<<std::endl;
		for(int v:r1) std::cout<<v<<" ";std::cout<<std::endl;
	}
	return 0;
}

详细

Test #1:

score: 100
Accepted
time: 8ms
memory: 70652kb

input:

5 7
1 2 1
1 3 -2
2 3 1
3 4 3
4 5 1
1 5 -1
2 5 2

output:

5
4
1 3 4 5 

result:

ok ok

Test #2:

score: 0
Accepted
time: 68ms
memory: 84284kb

input:

99997 100000
12238 99016 352755196
99016 25485 -412473602
25485 2440 991507552
2440 31171 -181894654
36970 2440 -800167579
2440 41865 -148191946
96629 31171 847888506
36970 95740 395546542
27992 2440 647886610
99016 29557 369124914
80795 27992 -673871966
36970 3509 573208857
57672 29557 874406776
41...

output:

1959330954
56
31171 2440 27992 90367 94522 24616 56899 57732 80256 964 80875 37318 41262 21467 79804 22393 7822 43423 93639 10503 73944 25485 99016 29557 57672 69259 68883 44442 4697 96048 94991 86274 74677 92617 86665 76022 72089 22074 96230 87712 51491 72825 3463 84407 67966 89628 84997 54073 6852...

result:

ok ok

Test #3:

score: 0
Accepted
time: 53ms
memory: 87564kb

input:

99997 100000
41884 21178 -431811360
41884 42699 -450057006
36523 21178 582079730
21178 96679 615552614
63637 21178 498974417
96679 5108 235820276
75058 41884 220112636
35148 42699 589595309
36523 18002 -637739861
65854 5108 -312755792
45137 41884 -511118771
5108 31311 554050951
25335 35148 -28341059...

output:

1968439328
40
70298 33264 26071 90144 25926 52252 51434 69337 7577 5108 50088 6204 28694 41126 87303 83047 26981 54901 59612 14678 35287 78274 18331 89860 71024 99686 98098 23692 87673 42699 41884 45137 85192 38202 83711 83919 55330 71151 98733 99716 

result:

ok ok

Test #4:

score: 0
Accepted
time: 76ms
memory: 89992kb

input:

99984 99999
33974 29867 335681778
33974 87468 348956829
83048 87468 320849805
29867 69456 -424530698
72457 69456 -950650074
53838 83048 755969166
85914 69456 569454441
51728 87468 -202158773
15970 29867 -865071002
15970 94894 697607001
94894 74694 616318126
33974 11496 -89287579
53838 34365 -6577379...

output:

1985932414
194
11395 92249 27866 23026 83298 1652 76013 55585 53162 41141 5525 41602 63721 42073 54251 11958 88466 28353 49296 99532 17325 75547 55727 29864 31939 52174 85999 70541 55140 80600 7311 63247 51807 97023 29296 25952 23014 84486 96848 6053 59406 91605 30363 85067 58162 86558 90720 49758 5...

result:

ok ok

Test #5:

score: 0
Accepted
time: 67ms
memory: 87504kb

input:

99988 99992
8584 11873 -811540160
68064 11873 -930246087
11873 60056 916668870
68064 82193 -859062523
60056 75072 790866030
27767 75072 357619485
75072 78221 411650300
39636 82193 264106928
6675 60056 933851261
71747 78221 -508471038
11873 92771 -665232168
34402 27767 -906494982
11873 42714 63734230...

output:

1932268861
30
93559 75593 83884 20678 16385 62720 25891 75176 34179 64174 38274 8778 99809 28745 6675 60056 75072 91957 3186 29873 5185 88236 50628 2518 83283 23798 89787 8975 26922 21107 

result:

ok ok

Test #6:

score: 0
Accepted
time: 69ms
memory: 87844kb

input:

99996 99996
58191 98120 261718607
91298 98120 471683748
58191 68921 217652908
67441 91298 -731916804
78177 68921 810185021
98120 54747 -35446486
78177 2822 -409569426
91298 68058 -897038977
68921 39067 892161204
30165 78177 379543758
32418 98120 -139944101
11281 68921 422411872
37751 32418 331606200...

output:

1752928792
25
16812 65221 81854 64478 74019 46983 4369 77560 30739 1917 58464 3308 3783 48440 58191 68921 11281 34081 78216 29915 18913 85525 42296 25934 67249 

result:

ok ok

Test #7:

score: 0
Accepted
time: 51ms
memory: 88676kb

input:

99996 100000
39127 4358 657531703
4358 66528 484843263
47215 4358 -856669390
47215 26179 -147254695
24822 39127 -635228854
81984 26179 600617794
24822 60559 327733708
39127 23879 286268283
95563 81984 -766366787
96587 24822 723252700
23879 13711 -303309809
60559 38379 992907085
60559 6012 -15086498
...

output:

1948904917
51
3750 42006 67183 31390 97465 24485 4358 61378 82877 24890 38051 91154 63827 13465 90743 30145 41482 70789 91374 66385 37095 8030 96587 24822 60559 38379 67362 37704 53517 23254 1066 28267 79904 54151 24450 79459 52647 10570 24822 39127 23879 13711 16150 60168 12024 27325 73608 58919 19...

result:

ok ok

Test #8:

score: -100
Wrong Answer
time: 77ms
memory: 83244kb

input:

99983 99998
360 38113 273639182
29807 360 -492749399
360 45494 960572841
67090 45494 -168787586
38113 61765 -90469418
71988 360 -556152065
67090 77653 704061103
30847 38113 542389160
84363 30847 295740326
30847 62591 -916431414
86104 77653 878763485
45494 11422 -795069866
86104 64096 714130240
61765...

output:

1972142685
solve 3
multiple point found
multiple point found
multiple point found
multiple point found
multiple point found
multiple point found
multiple point found
multiple point found
multiple point found
multiple point found
multiple point found
multiple point found
multiple point found
multiple...

result:

wrong output format Expected integer, but "solve" found