#include<bits/stdc++.h>
#pragma GCC optimize ("O3")
#pragma GCC target ("sse4")
#include "simurgh.h"
#define MAXN 507
#define MAXM 300007
using namespace std;
int n,m,dep[MAXN],edges;
vector< pair<int,int> > v[MAXN];
bool intree[MAXM],vis[MAXN];
pair<int,int> edge[MAXM],parent[MAXN];
int value[MAXM],st[MAXM],sumtree;
vector<int> ans;
unordered_set<int> s;
int questions;
int calc(vector<int> rem,vector<int> add){
questions++;
for(int i:rem)s.erase(i);
for(int i:add)s.insert(i);
vector<int> w;
for(int i:s){
w.push_back(i-1);
}
for(int i:rem)s.insert(i);
for(int i:add)s.erase(i);
return count_common_roads(w);
}
void dfs(int x){
vis[x]=true;
for(int i=0;i<v[x].size();i++){
if(vis[v[x][i].first])continue;
dep[v[x][i].first]=dep[x]+1;
parent[v[x][i].first]={x,v[x][i].second};
dfs(v[x][i].first);
intree[v[x][i].second]=true;
s.insert(v[x][i].second);
}
}
vector<int> path,toask;
int total;
void dfs2(int x){
vis[x]=true;
for(int i=0;i<v[x].size();i++){
if(vis[v[x][i].first] or !intree[v[x][i].second]){
continue;
}
toask.push_back(v[x][i].second-1);
total-=(value[v[x][i].second]-1);
dfs2(v[x][i].first);
}
}
bool ok;
void up(int x,int y){
if(x==y)return;
if(value[parent[x].second]==0 or ok){
if(value[parent[x].second])ok=false;
path.push_back(parent[x].second);
}
up(parent[x].first,y);
}
int howmany(int root,int l,int r){
toask.clear();
for(int i=1;i<=n;i++)vis[i]=false;
vis[root]=true;
for(int i=l;i<=r;i++)vis[v[root][i].first]=true;
for(int i=l;i<=r;i++){
toask.push_back(v[root][i].second-1);
}
total=0;
dfs2(root);
for(int i=l;i<=r;i++){
dfs2(v[root][i].first);
}
total+=count_common_roads(toask);
return total;
}
void solve(int root,int l,int r,int s){
if(l==r){
value[v[root][l].second]=2;
return;
}
int mid=(l+r)/2;
int lt=howmany(root,l,mid);
if(lt>0){
solve(root,l,mid,lt);
}
if(s-lt){
solve(root,mid+1,r,s-lt);
}
}
vector<int> find_roads(int N,vector<int> from,vector<int> to){
n=N; m=int(from.size());
double start=clock();
for(int i=1;i<=m;i++){
from[i-1]++; to[i-1]++;
edge[i]={from[i-1],to[i-1]};
v[from[i-1]].push_back({to[i-1],i});
v[to[i-1]].push_back({from[i-1],i});
}
dfs(1);
sumtree=calc({},{});
for(int i=1;i<=m;i++){
if(intree[i])continue;
int u=from[i-1],w=to[i-1];
if(dep[u]>dep[w])swap(u,w);
path.clear(); ok=true;
up(w,u);
int mins=sumtree,maxs=sumtree,br=0,ask=0;
for(int f:path){
if(value[f]==0){
st[f]=calc({f},{i});
br++;
mins=min(mins,st[f]);
maxs=max(maxs,st[f]);
}
}
if(br==0)continue;
if(mins==maxs){
for(int f:path){
if(value[f]==0)continue;
st[f]=calc({f},{i});
ask=value[f];
mins=min(mins,st[f]);
maxs=max(maxs,st[f]);
break;
}
}
if(ask==0)ask=1;
for(int f:path){
if(value[f]!=0)continue;
if(mins==maxs)value[f]=ask;
else{
if(st[f]==mins)value[f]=2;
else value[f]=1;
}
}
}
for(int i=1;i<=m;i++){
if(intree[i] and value[i]==0)value[i]=2;
}
double fin=clock();
if((fin-start)/CLOCKS_PER_SEC>0.3)cout<<1/0;
for(int i=1;i<=n;i++){
solve(i,0,v[i].size()-1,howmany(i,0,v[i].size()-1));
double fn=clock();
if((fn-start)/CLOCKS_PER_SEC>2.5)cout<<1/0;
}
for(int i=1;i<=m;i++){
if(value[i]==2){
ans.push_back(i-1);
}
}
return ans;
}
Compilation message
simurgh.cpp: In function 'void dfs(int)':
simurgh.cpp:46:15: warning: comparison of integer expressions of different signedness: 'int' and 'std::vector<std::pair<int, int> >::size_type' {aka 'long unsigned int'} [-Wsign-compare]
46 | for(int i=0;i<v[x].size();i++){
| ~^~~~~~~~~~~~
simurgh.cpp: In function 'void dfs2(int)':
simurgh.cpp:65:15: warning: comparison of integer expressions of different signedness: 'int' and 'std::vector<std::pair<int, int> >::size_type' {aka 'long unsigned int'} [-Wsign-compare]
65 | for(int i=0;i<v[x].size();i++){
| ~^~~~~~~~~~~~
simurgh.cpp: In function 'std::vector<int> find_roads(int, std::vector<int>, std::vector<int>)':
simurgh.cpp:206:43: warning: division by zero [-Wdiv-by-zero]
206 | if((fin-start)/CLOCKS_PER_SEC>0.3)cout<<1/0;
| ~^~
simurgh.cpp:212:43: warning: division by zero [-Wdiv-by-zero]
212 | if((fn-start)/CLOCKS_PER_SEC>2.5)cout<<1/0;
| ~^~
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
1 ms |
2384 KB |
correct |
2 |
Correct |
1 ms |
2384 KB |
correct |
3 |
Correct |
1 ms |
2384 KB |
correct |
4 |
Correct |
2 ms |
2384 KB |
correct |
5 |
Correct |
1 ms |
2384 KB |
correct |
6 |
Correct |
1 ms |
2384 KB |
correct |
7 |
Correct |
1 ms |
2384 KB |
correct |
8 |
Correct |
1 ms |
2384 KB |
correct |
9 |
Correct |
1 ms |
2384 KB |
correct |
10 |
Correct |
1 ms |
2384 KB |
correct |
11 |
Correct |
1 ms |
2384 KB |
correct |
12 |
Correct |
1 ms |
2384 KB |
correct |
13 |
Correct |
1 ms |
2640 KB |
correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
1 ms |
2384 KB |
correct |
2 |
Correct |
1 ms |
2384 KB |
correct |
3 |
Correct |
1 ms |
2384 KB |
correct |
4 |
Correct |
2 ms |
2384 KB |
correct |
5 |
Correct |
1 ms |
2384 KB |
correct |
6 |
Correct |
1 ms |
2384 KB |
correct |
7 |
Correct |
1 ms |
2384 KB |
correct |
8 |
Correct |
1 ms |
2384 KB |
correct |
9 |
Correct |
1 ms |
2384 KB |
correct |
10 |
Correct |
1 ms |
2384 KB |
correct |
11 |
Correct |
1 ms |
2384 KB |
correct |
12 |
Correct |
1 ms |
2384 KB |
correct |
13 |
Correct |
1 ms |
2640 KB |
correct |
14 |
Correct |
6 ms |
2384 KB |
correct |
15 |
Correct |
6 ms |
2384 KB |
correct |
16 |
Correct |
5 ms |
2384 KB |
correct |
17 |
Correct |
5 ms |
2384 KB |
correct |
18 |
Correct |
2 ms |
2384 KB |
correct |
19 |
Correct |
5 ms |
2384 KB |
correct |
20 |
Correct |
5 ms |
2384 KB |
correct |
21 |
Correct |
4 ms |
2384 KB |
correct |
22 |
Correct |
3 ms |
2384 KB |
correct |
23 |
Correct |
3 ms |
2384 KB |
correct |
24 |
Correct |
2 ms |
2556 KB |
correct |
25 |
Correct |
1 ms |
2384 KB |
correct |
26 |
Correct |
2 ms |
2384 KB |
correct |
27 |
Correct |
3 ms |
2384 KB |
correct |
28 |
Correct |
2 ms |
2384 KB |
correct |
29 |
Correct |
1 ms |
2384 KB |
correct |
30 |
Correct |
3 ms |
2384 KB |
correct |
31 |
Correct |
3 ms |
2556 KB |
correct |
32 |
Correct |
3 ms |
2384 KB |
correct |
33 |
Correct |
3 ms |
2384 KB |
correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
1 ms |
2384 KB |
correct |
2 |
Correct |
1 ms |
2384 KB |
correct |
3 |
Correct |
1 ms |
2384 KB |
correct |
4 |
Correct |
2 ms |
2384 KB |
correct |
5 |
Correct |
1 ms |
2384 KB |
correct |
6 |
Correct |
1 ms |
2384 KB |
correct |
7 |
Correct |
1 ms |
2384 KB |
correct |
8 |
Correct |
1 ms |
2384 KB |
correct |
9 |
Correct |
1 ms |
2384 KB |
correct |
10 |
Correct |
1 ms |
2384 KB |
correct |
11 |
Correct |
1 ms |
2384 KB |
correct |
12 |
Correct |
1 ms |
2384 KB |
correct |
13 |
Correct |
1 ms |
2640 KB |
correct |
14 |
Correct |
6 ms |
2384 KB |
correct |
15 |
Correct |
6 ms |
2384 KB |
correct |
16 |
Correct |
5 ms |
2384 KB |
correct |
17 |
Correct |
5 ms |
2384 KB |
correct |
18 |
Correct |
2 ms |
2384 KB |
correct |
19 |
Correct |
5 ms |
2384 KB |
correct |
20 |
Correct |
5 ms |
2384 KB |
correct |
21 |
Correct |
4 ms |
2384 KB |
correct |
22 |
Correct |
3 ms |
2384 KB |
correct |
23 |
Correct |
3 ms |
2384 KB |
correct |
24 |
Correct |
2 ms |
2556 KB |
correct |
25 |
Correct |
1 ms |
2384 KB |
correct |
26 |
Correct |
2 ms |
2384 KB |
correct |
27 |
Correct |
3 ms |
2384 KB |
correct |
28 |
Correct |
2 ms |
2384 KB |
correct |
29 |
Correct |
1 ms |
2384 KB |
correct |
30 |
Correct |
3 ms |
2384 KB |
correct |
31 |
Correct |
3 ms |
2556 KB |
correct |
32 |
Correct |
3 ms |
2384 KB |
correct |
33 |
Correct |
3 ms |
2384 KB |
correct |
34 |
Correct |
704 ms |
3936 KB |
correct |
35 |
Correct |
679 ms |
3684 KB |
correct |
36 |
Correct |
478 ms |
3408 KB |
correct |
37 |
Correct |
20 ms |
2384 KB |
correct |
38 |
Correct |
714 ms |
3656 KB |
correct |
39 |
Correct |
519 ms |
3632 KB |
correct |
40 |
Correct |
380 ms |
3408 KB |
correct |
41 |
Correct |
455 ms |
3664 KB |
correct |
42 |
Correct |
453 ms |
3688 KB |
correct |
43 |
Correct |
177 ms |
3192 KB |
correct |
44 |
Correct |
127 ms |
2896 KB |
correct |
45 |
Correct |
160 ms |
3068 KB |
correct |
46 |
Correct |
114 ms |
2896 KB |
correct |
47 |
Correct |
34 ms |
2640 KB |
correct |
48 |
Correct |
6 ms |
2384 KB |
correct |
49 |
Correct |
12 ms |
2612 KB |
correct |
50 |
Correct |
35 ms |
2640 KB |
correct |
51 |
Correct |
150 ms |
2896 KB |
correct |
52 |
Correct |
122 ms |
2896 KB |
correct |
53 |
Correct |
102 ms |
2896 KB |
correct |
54 |
Correct |
172 ms |
3192 KB |
correct |
55 |
Correct |
200 ms |
2896 KB |
correct |
56 |
Correct |
199 ms |
3092 KB |
correct |
57 |
Correct |
219 ms |
3100 KB |
correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
1 ms |
2384 KB |
correct |
2 |
Correct |
1 ms |
2640 KB |
correct |
3 |
Runtime error |
2523 ms |
11624 KB |
Execution killed with signal 4 |
4 |
Halted |
0 ms |
0 KB |
- |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
1 ms |
2384 KB |
correct |
2 |
Correct |
1 ms |
2384 KB |
correct |
3 |
Correct |
1 ms |
2384 KB |
correct |
4 |
Correct |
2 ms |
2384 KB |
correct |
5 |
Correct |
1 ms |
2384 KB |
correct |
6 |
Correct |
1 ms |
2384 KB |
correct |
7 |
Correct |
1 ms |
2384 KB |
correct |
8 |
Correct |
1 ms |
2384 KB |
correct |
9 |
Correct |
1 ms |
2384 KB |
correct |
10 |
Correct |
1 ms |
2384 KB |
correct |
11 |
Correct |
1 ms |
2384 KB |
correct |
12 |
Correct |
1 ms |
2384 KB |
correct |
13 |
Correct |
1 ms |
2640 KB |
correct |
14 |
Correct |
6 ms |
2384 KB |
correct |
15 |
Correct |
6 ms |
2384 KB |
correct |
16 |
Correct |
5 ms |
2384 KB |
correct |
17 |
Correct |
5 ms |
2384 KB |
correct |
18 |
Correct |
2 ms |
2384 KB |
correct |
19 |
Correct |
5 ms |
2384 KB |
correct |
20 |
Correct |
5 ms |
2384 KB |
correct |
21 |
Correct |
4 ms |
2384 KB |
correct |
22 |
Correct |
3 ms |
2384 KB |
correct |
23 |
Correct |
3 ms |
2384 KB |
correct |
24 |
Correct |
2 ms |
2556 KB |
correct |
25 |
Correct |
1 ms |
2384 KB |
correct |
26 |
Correct |
2 ms |
2384 KB |
correct |
27 |
Correct |
3 ms |
2384 KB |
correct |
28 |
Correct |
2 ms |
2384 KB |
correct |
29 |
Correct |
1 ms |
2384 KB |
correct |
30 |
Correct |
3 ms |
2384 KB |
correct |
31 |
Correct |
3 ms |
2556 KB |
correct |
32 |
Correct |
3 ms |
2384 KB |
correct |
33 |
Correct |
3 ms |
2384 KB |
correct |
34 |
Correct |
704 ms |
3936 KB |
correct |
35 |
Correct |
679 ms |
3684 KB |
correct |
36 |
Correct |
478 ms |
3408 KB |
correct |
37 |
Correct |
20 ms |
2384 KB |
correct |
38 |
Correct |
714 ms |
3656 KB |
correct |
39 |
Correct |
519 ms |
3632 KB |
correct |
40 |
Correct |
380 ms |
3408 KB |
correct |
41 |
Correct |
455 ms |
3664 KB |
correct |
42 |
Correct |
453 ms |
3688 KB |
correct |
43 |
Correct |
177 ms |
3192 KB |
correct |
44 |
Correct |
127 ms |
2896 KB |
correct |
45 |
Correct |
160 ms |
3068 KB |
correct |
46 |
Correct |
114 ms |
2896 KB |
correct |
47 |
Correct |
34 ms |
2640 KB |
correct |
48 |
Correct |
6 ms |
2384 KB |
correct |
49 |
Correct |
12 ms |
2612 KB |
correct |
50 |
Correct |
35 ms |
2640 KB |
correct |
51 |
Correct |
150 ms |
2896 KB |
correct |
52 |
Correct |
122 ms |
2896 KB |
correct |
53 |
Correct |
102 ms |
2896 KB |
correct |
54 |
Correct |
172 ms |
3192 KB |
correct |
55 |
Correct |
200 ms |
2896 KB |
correct |
56 |
Correct |
199 ms |
3092 KB |
correct |
57 |
Correct |
219 ms |
3100 KB |
correct |
58 |
Correct |
1 ms |
2384 KB |
correct |
59 |
Correct |
1 ms |
2640 KB |
correct |
60 |
Runtime error |
2523 ms |
11624 KB |
Execution killed with signal 4 |
61 |
Halted |
0 ms |
0 KB |
- |