Submission #987097

#	Submission time	Handle	Problem	Language	Result	Execution time	Memory
987097	2024-05-21T21:01:28 Z	Pyqe	Thousands Islands (IOI22_islands)	C++17	100 / 100	113 ms	43764 KB

islands

#include <bits/stdc++.h>
#include "islands.h"

using namespace std;

#define mp make_pair
#define fr first
#define sc second

const long long inf=1e18;
long long n,m,nn=0,nn2[2],nn3[2],cnn,dh[100069],an[100069],cyc[100069],ex[100069],ex2[2][100069],ex3[2][100069],cex[100069];
vector<pair<long long,long long>> al[100069],nx[100069],nx2[100069];
vector<long long> vl[100069];
bitset<100069> vtd,vtd2;
pair<long long,long long> pr[100069];
bitset<200069> spc[2];

void dfs(long long x)
{
	long long i,j,sz=al[x].size(),sz2,k,l,p;
	
	vtd[x]=1;
	vtd2[x]=1;
	for(i=0;i<sz;i++)
	{
		l=al[x][i].fr;
		p=al[x][i].sc;
		if(!vtd[l])
		{
			dh[l]=dh[x]+1;
			pr[l]={x,p};
			dfs(l);
		}
		if(vtd2[l])
		{
			if(dh[l]>dh[an[x]]&&!cyc[x])
			{
				an[x]=l;
				nx[x]={{l,p}};
			}
		}
		else
		{
			if(cyc[l])
			{
				if(!cyc[x])
				{
					nx[x].clear();
				}
				nx2[x].push_back({l,p});
			}
			cyc[x]=min(cyc[x]+cyc[l],2ll);
			if(vl[x].empty())
			{
				if(dh[an[l]]>dh[an[x]]&&!cyc[x])
				{
					an[x]=an[l];
					nx[x]={{l,p}};
				}
			}
			else
			{
				if(!cyc[l])
				{
					if(!cyc[x])
					{
						nx[x].clear();
					}
					cyc[x]=min(cyc[x]+1,2ll);
					nx[x].push_back({l,p});
				}
				sz2=vl[x].size();
				for(j=0;j<sz2;j++)
				{
					k=vl[x][j];
					cyc[k]=max(cyc[k],1ll);
				}
				vl[x].clear();
			}
		}
	}
	vtd2[x]=0;
	if(vtd2[an[x]]&&!cyc[x])
	{
		vl[an[x]].push_back(x);
	}
}

variant<bool,vector<int>> find_journey(int on,int om,vector<int> ka,vector<int> la)
{
	long long i,ii,k,l,p,sz,e,kk,ll,rt=0;
	pair<long long,long long> mne={inf,-1};
	vector<int> sq;
	
	n=on;
	m=om;
	for(i=1;i<=m;i++)
	{
		k=ka[i-1]+1;
		l=la[i-1]+1;
		al[k].push_back({l,i});
	}
	dh[0]=-1;
	dfs(1);
	if(cyc[1]!=2)
	{
		return false;
	}
	for(i=1;i<=n;i++)
	{
		if(cyc[i]==2&&dh[i]>dh[rt])
		{
			rt=i;
		}
	}
	for(k=rt;k!=1;k=pr[k].fr)
	{
		p=pr[k].sc;
		nn++;
		ex[nn]=p;
	}
	reverse(ex+1,ex+nn+1);
	for(ii=0;ii<2;ii++)
	{
		for(k=rt;!nx2[k].empty();k=l)
		{
			l=nx2[k].back().fr;
			p=nx2[k].back().sc;
			nn2[ii]++;
			ex2[ii][nn2[ii]]=p;
		}
		kk=k;
		for(;1;k=l)
		{
			l=nx[k][0].fr;
			p=nx[k][0].sc;
			nn3[ii]++;
			ex3[ii][nn3[ii]]=p;
			if(l==an[k])
			{
				break;
			}
		}
		ll=l;
		cnn=0;
		for(k=kk;k!=ll;k=pr[k].fr)
		{
			p=pr[k].sc;
			cnn++;
			cex[cnn]=p;
		}
		reverse(cex+1,cex+cnn+1);
		for(i=1;i<=cnn;i++)
		{
			nn3[ii]++;
			ex3[ii][nn3[ii]]=cex[i];
		}
		for(i=1;i<=n;i++)
		{
			reverse(nx[i].begin(),nx[i].end());
		}
		if(!nx2[rt].empty())
		{
			nx2[rt].pop_back();
		}
		mne=min(mne,{dh[ll],ii});
	}
	for(ii=0;ii<2;ii++)
	{
		for(i=1;i<=nn3[ii];i++)
		{
			k=ex3[ii][i];
			spc[ii][k]=1;
		}
	}
	for(i=1;i<=nn3[1];i++)
	{
		k=ex3[1][i];
		if(!spc[0][k])
		{
			break;
		}
	}
	if(nn3[0]==nn3[1]&&i>nn3[0])
	{
		for(i=1;i<=nn;i++)
		{
			sq.push_back(ex[i]);
		}
		for(i=1;i<=nn2[0];i++)
		{
			sq.push_back(ex2[0][i]);
		}
		for(i=1;i<=nn3[0];i++)
		{
			sq.push_back(ex3[0][i]);
		}
		for(i=nn2[0];i;i--)
		{
			sq.push_back(ex2[0][i]);
		}
		for(i=1;i<=nn2[1];i++)
		{
			sq.push_back(ex2[1][i]);
		}
		for(i=nn3[1];i;i--)
		{
			sq.push_back(ex3[1][i]);
		}
		for(i=nn2[1];i;i--)
		{
			sq.push_back(ex2[1][i]);
		}
		for(i=nn;i;i--)
		{
			sq.push_back(ex[i]);
		}
	}
	else
	{
		for(i=1;i<=nn3[1];i++)
		{
			k=ex3[1][i];
			if(spc[0][k])
			{
				break;
			}
		}
		if(i<=nn3[1])
		{
			e=mne.sc;
			for(i=1;i<=nn;i++)
			{
				sq.push_back(ex[i]);
			}
			for(i=1;i<=nn2[!e];i++)
			{
				sq.push_back(ex2[!e][i]);
			}
			for(i=1;i<=nn3[!e];i++)
			{
				sq.push_back(ex3[!e][i]);
			}
			for(i=nn2[!e];i;i--)
			{
				sq.push_back(ex2[!e][i]);
			}
			for(i=1;i<=nn2[e];i++)
			{
				sq.push_back(ex2[e][i]);
			}
			for(i=1;i<=nn3[e];i++)
			{
				k=ex3[e][i];
				if(spc[!e][k])
				{
					break;
				}
				sq.push_back(k);
			}
			kk=i;
			for(i=1;ex3[!e][i]!=k;i++);
			for(i=(i+nn3[!e]-2)%nn3[!e]+1;1;i=(i+nn3[!e]-2)%nn3[!e]+1)
			{
				k=ex3[!e][i];
				if(spc[e][k])
				{
					break;
				}
				sq.push_back(k);
			}
			for(i=kk;i<=nn3[e];i++)
			{
				k=ex3[e][i];
				if(!spc[!e][k])
				{
					sq.push_back(k);
				}
			}
			for(i=nn2[e];i;i--)
			{
				sq.push_back(ex2[e][i]);
			}
			for(i=1;i<=nn2[!e];i++)
			{
				k=ex2[!e][i];
				if(spc[e][k])
				{
					break;
				}
				sq.push_back(k);
			}
			if(i<=nn2[!e])
			{
				kk=k;
			}
			else
			{
				kk=ex3[!e][1];
			}
			for(i=1;ex3[e][i]!=kk;i++);
			for(i=(i+nn3[e]-2)%nn3[e]+1;1;i=(i+nn3[e]-2)%nn3[e]+1)
			{
				k=ex3[e][i];
				sq.push_back(k);
				if(k==kk)
				{
					break;
				}
			}
			for(i=nn2[!e];i;i--)
			{
				k=ex2[!e][i];
				if(!spc[e][k])
				{
					sq.push_back(k);
				}
			}
			for(i=nn;i;i--)
			{
				sq.push_back(ex[i]);
			}
		}
		else
		{
			for(i=1;i<=nn;i++)
			{
				sq.push_back(ex[i]);
			}
			for(i=1;i<=nn2[0];i++)
			{
				sq.push_back(ex2[0][i]);
			}
			for(i=1;i<=nn3[0];i++)
			{
				sq.push_back(ex3[0][i]);
			}
			for(i=nn2[0];i;i--)
			{
				sq.push_back(ex2[0][i]);
			}
			for(i=1;i<=nn2[1];i++)
			{
				sq.push_back(ex2[1][i]);
			}
			for(i=1;i<=nn3[1];i++)
			{
				sq.push_back(ex3[1][i]);
			}
			for(i=nn2[1];i;i--)
			{
				sq.push_back(ex2[1][i]);
			}
			for(i=1;i<=nn2[0];i++)
			{
				sq.push_back(ex2[0][i]);
			}
			for(i=nn3[0];i;i--)
			{
				sq.push_back(ex3[0][i]);
			}
			for(i=nn2[0];i;i--)
			{
				sq.push_back(ex2[0][i]);
			}
			for(i=1;i<=nn2[1];i++)
			{
				sq.push_back(ex2[1][i]);
			}
			for(i=nn3[1];i;i--)
			{
				sq.push_back(ex3[1][i]);
			}
			for(i=nn2[1];i;i--)
			{
				sq.push_back(ex2[1][i]);
			}
			for(i=nn;i;i--)
			{
				sq.push_back(ex[i]);
			}
		}
	}
	sz=sq.size();
	for(i=0;i<sz;i++)
	{
		sq[i]--;
	}
	return sq;
}

Subtask #1

5.0 / 5.0

#	Verdict	Execution time	Memory	Grader output
1	Correct	4 ms	12636 KB	Output is correct
2	Correct	3 ms	14684 KB	Output is correct
3	Correct	4 ms	12636 KB	Output is correct
4	Correct	3 ms	12636 KB	Output is correct
5	Correct	4 ms	12636 KB	Output is correct
6	Correct	3 ms	14680 KB	Output is correct
7	Correct	33 ms	22580 KB	Output is correct

Subtask #2

5.0 / 5.0

#	Verdict	Execution time	Memory	Grader output
1	Correct	3 ms	14684 KB	Output is correct
2	Correct	4 ms	12636 KB	Output is correct
3	Correct	4 ms	14684 KB	Output is correct
4	Correct	4 ms	14684 KB	Output is correct
5	Correct	4 ms	14684 KB	Output is correct
6	Correct	31 ms	22456 KB	Output is correct

Subtask #3

21.0 / 21.0

#	Verdict	Execution time	Memory	Grader output
1	Correct	5 ms	15196 KB	Output is correct
2	Correct	3 ms	12636 KB	Output is correct
3	Correct	3 ms	12636 KB	Output is correct
4	Correct	6 ms	14788 KB	Output is correct
5	Correct	4 ms	12892 KB	Output is correct
6	Correct	4 ms	14684 KB	Output is correct
7	Correct	4 ms	12636 KB	Output is correct
8	Correct	3 ms	14684 KB	Output is correct
9	Correct	5 ms	14684 KB	Output is correct
10	Correct	4 ms	14940 KB	Output is correct
11	Correct	4 ms	14780 KB	Output is correct
12	Correct	4 ms	14940 KB	Output is correct
13	Correct	3 ms	12636 KB	Output is correct
14	Correct	3 ms	12888 KB	Output is correct
15	Correct	3 ms	14684 KB	Output is correct
16	Correct	3 ms	12636 KB	Output is correct
17	Correct	23 ms	19696 KB	Output is correct
18	Correct	18 ms	19800 KB	Output is correct
19	Correct	4 ms	12636 KB	Output is correct
20	Correct	4 ms	14680 KB	Output is correct
21	Correct	5 ms	14680 KB	Output is correct
22	Correct	4 ms	14684 KB	Output is correct
23	Correct	31 ms	23636 KB	Output is correct

Subtask #4

24.0 / 24.0

#	Verdict	Execution time	Memory	Grader output
1	Correct	4 ms	14684 KB	Output is correct
2	Correct	5 ms	14936 KB	Output is correct
3	Correct	30 ms	21296 KB	Output is correct
4	Correct	33 ms	21836 KB	Output is correct
5	Correct	4 ms	14940 KB	Output is correct
6	Correct	5 ms	14940 KB	Output is correct
7	Correct	5 ms	14684 KB	Output is correct
8	Correct	4 ms	14684 KB	Output is correct
9	Correct	4 ms	14684 KB	Output is correct
10	Correct	5 ms	15068 KB	Output is correct
11	Correct	4 ms	12892 KB	Output is correct
12	Correct	5 ms	14936 KB	Output is correct
13	Correct	5 ms	15196 KB	Output is correct
14	Correct	5 ms	15096 KB	Output is correct
15	Correct	5 ms	15044 KB	Output is correct
16	Correct	4 ms	14940 KB	Output is correct
17	Correct	3 ms	12636 KB	Output is correct
18	Correct	5 ms	14940 KB	Output is correct
19	Correct	4 ms	14940 KB	Output is correct
20	Correct	36 ms	20056 KB	Output is correct
21	Correct	31 ms	21580 KB	Output is correct
22	Correct	3 ms	12892 KB	Output is correct
23	Correct	3 ms	12636 KB	Output is correct
24	Correct	4 ms	12636 KB	Output is correct
25	Correct	3 ms	14680 KB	Output is correct
26	Correct	5 ms	15172 KB	Output is correct
27	Correct	32 ms	22620 KB	Output is correct
28	Correct	35 ms	22696 KB	Output is correct
29	Correct	3 ms	12632 KB	Output is correct
30	Correct	39 ms	21060 KB	Output is correct
31	Correct	5 ms	14840 KB	Output is correct
32	Correct	33 ms	22604 KB	Output is correct
33	Correct	33 ms	21844 KB	Output is correct
34	Correct	19 ms	19804 KB	Output is correct
35	Correct	4 ms	14940 KB	Output is correct
36	Correct	30 ms	21328 KB	Output is correct
37	Correct	35 ms	24404 KB	Output is correct
38	Correct	4 ms	12632 KB	Output is correct
39	Correct	30 ms	19036 KB	Output is correct
40	Correct	5 ms	14936 KB	Output is correct
41	Correct	35 ms	20984 KB	Output is correct
42	Correct	34 ms	22612 KB	Output is correct
43	Correct	3 ms	12632 KB	Output is correct
44	Correct	3 ms	12892 KB	Output is correct
45	Correct	5 ms	14940 KB	Output is correct
46	Correct	18 ms	19032 KB	Output is correct

Subtask #5

45.0 / 45.0

#	Verdict	Execution time	Memory	Grader output
1	Correct	4 ms	12636 KB	Output is correct
2	Correct	3 ms	14684 KB	Output is correct
3	Correct	4 ms	12636 KB	Output is correct
4	Correct	3 ms	12636 KB	Output is correct
5	Correct	4 ms	12636 KB	Output is correct
6	Correct	3 ms	14680 KB	Output is correct
7	Correct	33 ms	22580 KB	Output is correct
8	Correct	3 ms	14684 KB	Output is correct
9	Correct	4 ms	12636 KB	Output is correct
10	Correct	4 ms	14684 KB	Output is correct
11	Correct	4 ms	14684 KB	Output is correct
12	Correct	4 ms	14684 KB	Output is correct
13	Correct	31 ms	22456 KB	Output is correct
14	Correct	5 ms	15196 KB	Output is correct
15	Correct	3 ms	12636 KB	Output is correct
16	Correct	3 ms	12636 KB	Output is correct
17	Correct	6 ms	14788 KB	Output is correct
18	Correct	4 ms	12892 KB	Output is correct
19	Correct	4 ms	14684 KB	Output is correct
20	Correct	4 ms	12636 KB	Output is correct
21	Correct	3 ms	14684 KB	Output is correct
22	Correct	5 ms	14684 KB	Output is correct
23	Correct	4 ms	14940 KB	Output is correct
24	Correct	4 ms	14780 KB	Output is correct
25	Correct	4 ms	14940 KB	Output is correct
26	Correct	3 ms	12636 KB	Output is correct
27	Correct	3 ms	12888 KB	Output is correct
28	Correct	3 ms	14684 KB	Output is correct
29	Correct	3 ms	12636 KB	Output is correct
30	Correct	23 ms	19696 KB	Output is correct
31	Correct	18 ms	19800 KB	Output is correct
32	Correct	4 ms	12636 KB	Output is correct
33	Correct	4 ms	14680 KB	Output is correct
34	Correct	5 ms	14680 KB	Output is correct
35	Correct	4 ms	14684 KB	Output is correct
36	Correct	31 ms	23636 KB	Output is correct
37	Correct	4 ms	14684 KB	Output is correct
38	Correct	3 ms	12636 KB	Output is correct
39	Correct	3 ms	12636 KB	Output is correct
40	Correct	3 ms	12636 KB	Output is correct
41	Correct	17 ms	16476 KB	Output is correct
42	Correct	5 ms	12892 KB	Output is correct
43	Correct	34 ms	27968 KB	Output is correct
44	Correct	38 ms	25032 KB	Output is correct
45	Correct	40 ms	30148 KB	Output is correct
46	Correct	4 ms	14684 KB	Output is correct
47	Correct	3 ms	14684 KB	Output is correct
48	Correct	3 ms	14684 KB	Output is correct
49	Correct	4 ms	14940 KB	Output is correct
50	Correct	92 ms	43764 KB	Output is correct
51	Correct	76 ms	33984 KB	Output is correct
52	Correct	64 ms	36084 KB	Output is correct
53	Correct	95 ms	40640 KB	Output is correct
54	Correct	94 ms	37556 KB	Output is correct
55	Correct	77 ms	39620 KB	Output is correct
56	Correct	103 ms	40388 KB	Output is correct
57	Correct	58 ms	29524 KB	Output is correct
58	Correct	110 ms	35292 KB	Output is correct
59	Correct	71 ms	25040 KB	Output is correct
60	Correct	51 ms	26448 KB	Output is correct
61	Correct	49 ms	23376 KB	Output is correct
62	Correct	10 ms	14428 KB	Output is correct
63	Correct	43 ms	22752 KB	Output is correct
64	Correct	22 ms	20056 KB	Output is correct
65	Correct	3 ms	14684 KB	Output is correct
66	Correct	3 ms	14680 KB	Output is correct
67	Correct	113 ms	35412 KB	Output is correct
68	Correct	43 ms	22612 KB	Output is correct
69	Correct	35 ms	24404 KB	Output is correct
70	Correct	4 ms	14936 KB	Output is correct
71	Correct	45 ms	21332 KB	Output is correct
72	Correct	4 ms	14680 KB	Output is correct
73	Correct	59 ms	24296 KB	Output is correct
74	Correct	105 ms	37944 KB	Output is correct
75	Correct	6 ms	15452 KB	Output is correct
76	Correct	41 ms	25928 KB	Output is correct
77	Correct	58 ms	40392 KB	Output is correct
78	Correct	48 ms	26448 KB	Output is correct
79	Correct	6 ms	14684 KB	Output is correct
80	Correct	90 ms	35048 KB	Output is correct
81	Correct	4 ms	14940 KB	Output is correct
82	Correct	33 ms	20812 KB	Output is correct
83	Correct	3 ms	14940 KB	Output is correct
84	Correct	4 ms	14680 KB	Output is correct
85	Correct	51 ms	25648 KB	Output is correct
86	Correct	4 ms	14940 KB	Output is correct
87	Correct	4 ms	12636 KB	Output is correct
88	Correct	6 ms	15196 KB	Output is correct
89	Correct	49 ms	26452 KB	Output is correct
90	Correct	50 ms	23552 KB	Output is correct
91	Correct	92 ms	37168 KB	Output is correct
92	Correct	3 ms	14684 KB	Output is correct