Submission #282944

# Submission time Handle Problem Language Result Execution time Memory
282944 2020-08-25T07:38:48 Z 최은수(#5745) Iqea (innopolis2018_final_C) C++17
36 / 100
1866 ms 262148 KB
#include<iostream>
#include<vector>
#include<algorithm>
#define ep emplace
#define eb emplace_back
#define fi first
#define se second
#define all(x) (x).begin(),(x).end()
using namespace std;
typedef long long ll;
typedef pair<int,int>pi;
typedef pair<ll,ll>pl;
const int inf=1e9+7;
const ll INF=1e18;
struct dsu
{
    int pa[100010];
    inline void init()
    {
        fill(pa,pa+100010,-1);
        return;
    }
    dsu(){init();}
    int par(int x)
    {
        if(pa[x]==-1)
            return x;
        return pa[x]=par(pa[x]);
    }
    inline void uni(int x,int y)
    {
        x=par(x);
        y=par(y);
        if(x==y)
            return;
        pa[x]=y;
        return;
    }
}uf1,uf2;
vector<int>adj[200010];
int sz[200010];
void sdfs(int x,int p)
{
    sz[x]=1;
    for(int&t:adj[x])
        if(t!=p)
            sdfs(t,x),sz[x]+=sz[t];
    return;
}
bool chk[200010];
int dis[200010][18];
void dfs(int x,int p,int d,int cd)
{
    dis[x][cd]=d;
    sz[x]=1;
    for(int&t:adj[x])
        if(!chk[t]&&t!=p)
            dfs(t,x,d+1,cd),sz[x]+=sz[t];
    return;
}
int get(int x,int p,int z)
{
    int ms=-1,mi=-1;
    for(int&t:adj[x])
        if(!chk[t]&&t!=p&&sz[t]>ms)
            ms=sz[t],mi=t;
    if(ms<=sz[z]/2)
        return x;
    return get(mi,x,z);
}
int cdep[200010];
int cpa[200010];
void dnc(int x)
{
    chk[x]=1;
    for(int&t:adj[x])
    {
        if(chk[t])
            continue;
        dfs(t,x,1,cdep[x]);
        int nc=get(t,x,t);
        cpa[nc]=x;
        cdep[nc]=cdep[x]+1;
        dnc(nc);
    }
    return;
}
vector<pair<unsigned int,unsigned char> >sv1[100010];
vector<pair<unsigned int,unsigned char> >sv2[100010];
vector<vector<int> >ind[100010];
int cnt1[100010],cnt2[100010];
int sz1[100010],sz2[100010];
int main()
{
    ios_base::sync_with_stdio(false);
    cin.tie(nullptr);
    int n;
    cin>>n;
    vector<pi>p(n);
    for(pi&t:p)
        cin>>t.fi>>t.se;
    vector<pi>sxv=p,syv=p;
    auto cmp1=[](const pi&x,const pi&y){return x<y;};
    auto cmp2=[](const pi&x,const pi&y){return pi(x.se,x.fi)<pi(y.se,y.fi);};
    sort(all(sxv),cmp1);
    sort(all(syv),cmp2);
    for(int i=1;i<n;i++)
        if(sxv[i].fi==sxv[i-1].fi&&sxv[i].se==sxv[i-1].se+1)
            uf1.uni(i,i-1);
    for(int i=1;i<n;i++)
        if(syv[i].fi==syv[i-1].fi+1&&syv[i].se==syv[i-1].se)
            uf2.uni(i,i-1);
    for(int i=1;i<n;i++)
    {
        if(sxv[i].fi==sxv[i-1].fi&&sxv[i].se==sxv[i-1].se+1)
        {
            int p1=uf2.par(lower_bound(all(syv),sxv[i-1],cmp2)-syv.begin())+n;
            int p2=uf2.par(lower_bound(all(syv),sxv[i],cmp2)-syv.begin())+n;
            adj[p1].eb(p2);
            adj[p2].eb(p1);
        }
    }
    for(int i=1;i<n;i++)
    {
        if(syv[i].fi==syv[i-1].fi+1&&syv[i].se==syv[i-1].se)
        {
            int p1=uf1.par(lower_bound(all(sxv),syv[i-1],cmp1)-sxv.begin());
            int p2=uf1.par(lower_bound(all(sxv),syv[i],cmp1)-sxv.begin());
            adj[p1].eb(p2);
            adj[p2].eb(p1);
        }
    }
    for(int i=0;i<n*2;i++)
    {
        sort(all(adj[i]));
        adj[i].erase(unique(all(adj[i])),adj[i].end());
    }
    sdfs(uf1.par(0),-1);
    int c1=get(uf1.par(0),-1,uf1.par(0));
    cpa[c1]=-1;
    dnc(c1);
    sdfs(uf2.par(0)+n,-1);
    int c2=get(uf2.par(0)+n,-1,uf2.par(0)+n);
    cpa[c2]=-1;
    dnc(c2);
    for(int i=0;i<n;i++)
    {
        auto&t=sxv[i];
        int p1=uf1.par(i);
        int p2=uf2.par(lower_bound(all(syv),t,cmp2)-syv.begin())+n;
        vector<int>cp1,cp2;
        while(p1!=-1)
            cp1.eb(p1),p1=cpa[p1];
        while(p2!=-1)
            cp2.eb(p2),p2=cpa[p2];
        for(int&u:cp1)
            cnt1[u]+=(int)cp2.size();
        for(int&v:cp2)
            cnt2[v-n]+=(int)cp1.size();
    }
    for(int i=0;i<n;i++)
        sv1[i].reserve(cnt1[i]);
    for(int i=0;i<n;i++)
    {
        auto&t=sxv[i];
        int p1=uf1.par(i);
        int p2=uf2.par(lower_bound(all(syv),t,cmp2)-syv.begin())+n;
        vector<int>cp1,cp2;
        while(p1!=-1)
            cp1.eb(p1),p1=cpa[p1];
        while(p2!=-1)
            cp2.eb(p2),p2=cpa[p2];
        for(int&u:cp1)
            for(int&v:cp2)
                sv1[u].eb(v-n+i%32768*131072,i/32768);
        sz1[i]=cp1.size();
        sz2[i]=cp2.size();
    }
    for(int i=0;i<n;i++)
    {
        for(auto&t:sv1[i])
        {
            int fi=t.fi%131072;
            int se=t.se*32768+t.fi/131072;
            if(sv2[fi].empty())
                sv2[fi].reserve(cnt2[fi]);
            sv2[fi].eb(i+se%32768*131072,se/32768);
        }
        sv1[i].clear(),sv1[i].shrink_to_fit();
    }
    int ct=0;
    for(int i=0;i<n;i++)
    {
        for(int j=0;j<(int)sv2[i].size();j++)
        {
            auto&t=sv2[i][j];
            int fi=t.fi%131072;
            int cid=t.se*32768+t.fi/131072;
            if(ind[cid].empty())
                ind[cid].resize(sz1[cid],vector<int>(sz2[cid],0));
            ind[cid][cdep[fi]][cdep[i+n]]=ct;
            if(j==(int)sv2[i].size()-1||fi!=(int)(sv2[i][j+1].fi%131072))
                ct++;
        }
        sv2[i].clear(),sv2[i].shrink_to_fit();
    }
    vector<int>arr(ct,0);
    int q;
    cin>>q;
    for(int qi=0;qi<q;qi++)
    {
        int t,x,y;
        cin>>t>>x>>y;
        int id=lower_bound(all(sxv),pi(x,y),cmp1)-sxv.begin();
        int p1=uf1.par(id);
        int p2=uf2.par(lower_bound(all(syv),pi(x,y),cmp2)-syv.begin())+n;
        if(t==1)
            for(int i=cdep[p1];i>=0;i--)
                for(int j=cdep[p2];j>=0;j--)
                    arr[ind[id][i][j]]=min(arr[ind[id][i][j]],dis[p1][i]+dis[p2][j]-inf);
        else
        {
            int ans=0;
            for(int i=cdep[p1];i>=0;i--)
                for(int j=cdep[p2];j>=0;j--)
                    ans=min(ans,arr[ind[id][i][j]]+dis[p1][i]+dis[p2][j]);
            if(ans==0)
                ans=-1;
            else
                ans+=inf;
            cout<<ans<<'\n';
        }
    }
    return 0;
}
# Verdict Execution time Memory Grader output
1 Correct 9 ms 12928 KB Output is correct
2 Correct 10 ms 12928 KB Output is correct
3 Correct 10 ms 12928 KB Output is correct
4 Correct 30 ms 16300 KB Output is correct
5 Correct 37 ms 16160 KB Output is correct
6 Correct 32 ms 17024 KB Output is correct
7 Correct 43 ms 17904 KB Output is correct
8 Correct 41 ms 18132 KB Output is correct
9 Correct 42 ms 18208 KB Output is correct
10 Correct 37 ms 17876 KB Output is correct
11 Correct 38 ms 17972 KB Output is correct
12 Correct 40 ms 18080 KB Output is correct
13 Correct 41 ms 18368 KB Output is correct
14 Correct 39 ms 17884 KB Output is correct
15 Correct 43 ms 17916 KB Output is correct
16 Correct 39 ms 17828 KB Output is correct
17 Correct 40 ms 18072 KB Output is correct
18 Correct 38 ms 17736 KB Output is correct
19 Correct 39 ms 17848 KB Output is correct
20 Correct 44 ms 19072 KB Output is correct
21 Correct 42 ms 18124 KB Output is correct
22 Correct 49 ms 20468 KB Output is correct
23 Correct 46 ms 19352 KB Output is correct
24 Correct 50 ms 20552 KB Output is correct
25 Correct 29 ms 17664 KB Output is correct
26 Correct 30 ms 15492 KB Output is correct
27 Correct 23 ms 14976 KB Output is correct
28 Correct 50 ms 20424 KB Output is correct
29 Correct 49 ms 20504 KB Output is correct
30 Correct 54 ms 20788 KB Output is correct
31 Correct 51 ms 20536 KB Output is correct
32 Correct 51 ms 20232 KB Output is correct
33 Correct 49 ms 20432 KB Output is correct
34 Correct 51 ms 20312 KB Output is correct
35 Correct 51 ms 20696 KB Output is correct
36 Correct 28 ms 17536 KB Output is correct
37 Correct 49 ms 20232 KB Output is correct
38 Correct 50 ms 20172 KB Output is correct
39 Correct 48 ms 20484 KB Output is correct
40 Correct 46 ms 20468 KB Output is correct
41 Correct 45 ms 19704 KB Output is correct
42 Correct 46 ms 19824 KB Output is correct
43 Correct 27 ms 17528 KB Output is correct
44 Correct 46 ms 19524 KB Output is correct
45 Correct 45 ms 19760 KB Output is correct
46 Correct 23 ms 15104 KB Output is correct
47 Correct 42 ms 18492 KB Output is correct
48 Correct 37 ms 17256 KB Output is correct
49 Correct 36 ms 17612 KB Output is correct
50 Correct 35 ms 17136 KB Output is correct
51 Correct 34 ms 17396 KB Output is correct
52 Correct 40 ms 18080 KB Output is correct
53 Correct 41 ms 18616 KB Output is correct
54 Correct 41 ms 18748 KB Output is correct
55 Correct 44 ms 18804 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 9 ms 12928 KB Output is correct
2 Correct 10 ms 12928 KB Output is correct
3 Correct 10 ms 12928 KB Output is correct
4 Correct 30 ms 16300 KB Output is correct
5 Correct 37 ms 16160 KB Output is correct
6 Correct 32 ms 17024 KB Output is correct
7 Correct 43 ms 17904 KB Output is correct
8 Correct 41 ms 18132 KB Output is correct
9 Correct 42 ms 18208 KB Output is correct
10 Correct 37 ms 17876 KB Output is correct
11 Correct 38 ms 17972 KB Output is correct
12 Correct 40 ms 18080 KB Output is correct
13 Correct 41 ms 18368 KB Output is correct
14 Correct 39 ms 17884 KB Output is correct
15 Correct 43 ms 17916 KB Output is correct
16 Correct 39 ms 17828 KB Output is correct
17 Correct 40 ms 18072 KB Output is correct
18 Correct 38 ms 17736 KB Output is correct
19 Correct 39 ms 17848 KB Output is correct
20 Correct 44 ms 19072 KB Output is correct
21 Correct 42 ms 18124 KB Output is correct
22 Correct 49 ms 20468 KB Output is correct
23 Correct 46 ms 19352 KB Output is correct
24 Correct 50 ms 20552 KB Output is correct
25 Correct 29 ms 17664 KB Output is correct
26 Correct 30 ms 15492 KB Output is correct
27 Correct 23 ms 14976 KB Output is correct
28 Correct 50 ms 20424 KB Output is correct
29 Correct 49 ms 20504 KB Output is correct
30 Correct 54 ms 20788 KB Output is correct
31 Correct 51 ms 20536 KB Output is correct
32 Correct 51 ms 20232 KB Output is correct
33 Correct 49 ms 20432 KB Output is correct
34 Correct 51 ms 20312 KB Output is correct
35 Correct 51 ms 20696 KB Output is correct
36 Correct 28 ms 17536 KB Output is correct
37 Correct 49 ms 20232 KB Output is correct
38 Correct 50 ms 20172 KB Output is correct
39 Correct 48 ms 20484 KB Output is correct
40 Correct 46 ms 20468 KB Output is correct
41 Correct 45 ms 19704 KB Output is correct
42 Correct 46 ms 19824 KB Output is correct
43 Correct 27 ms 17528 KB Output is correct
44 Correct 46 ms 19524 KB Output is correct
45 Correct 45 ms 19760 KB Output is correct
46 Correct 23 ms 15104 KB Output is correct
47 Correct 42 ms 18492 KB Output is correct
48 Correct 37 ms 17256 KB Output is correct
49 Correct 36 ms 17612 KB Output is correct
50 Correct 35 ms 17136 KB Output is correct
51 Correct 34 ms 17396 KB Output is correct
52 Correct 40 ms 18080 KB Output is correct
53 Correct 41 ms 18616 KB Output is correct
54 Correct 41 ms 18748 KB Output is correct
55 Correct 44 ms 18804 KB Output is correct
56 Correct 121 ms 16760 KB Output is correct
57 Correct 96 ms 17352 KB Output is correct
58 Correct 100 ms 16648 KB Output is correct
59 Correct 144 ms 17768 KB Output is correct
60 Correct 123 ms 18060 KB Output is correct
61 Correct 126 ms 18764 KB Output is correct
62 Correct 122 ms 17764 KB Output is correct
63 Correct 116 ms 18204 KB Output is correct
64 Correct 126 ms 18312 KB Output is correct
65 Correct 117 ms 18192 KB Output is correct
66 Correct 117 ms 18972 KB Output is correct
67 Correct 113 ms 17964 KB Output is correct
68 Correct 118 ms 18732 KB Output is correct
69 Correct 114 ms 18236 KB Output is correct
70 Correct 114 ms 18404 KB Output is correct
71 Correct 107 ms 17268 KB Output is correct
72 Correct 120 ms 18252 KB Output is correct
73 Correct 117 ms 18212 KB Output is correct
74 Correct 94 ms 17728 KB Output is correct
75 Correct 89 ms 15744 KB Output is correct
76 Correct 81 ms 15284 KB Output is correct
77 Correct 136 ms 20384 KB Output is correct
78 Correct 146 ms 20412 KB Output is correct
79 Correct 105 ms 17712 KB Output is correct
80 Correct 139 ms 20368 KB Output is correct
81 Correct 157 ms 20176 KB Output is correct
82 Correct 138 ms 19724 KB Output is correct
83 Correct 139 ms 19940 KB Output is correct
84 Correct 95 ms 17484 KB Output is correct
85 Correct 137 ms 19668 KB Output is correct
86 Correct 133 ms 19660 KB Output is correct
87 Correct 92 ms 15268 KB Output is correct
88 Correct 138 ms 18008 KB Output is correct
89 Correct 123 ms 17876 KB Output is correct
90 Correct 136 ms 17952 KB Output is correct
91 Correct 131 ms 17508 KB Output is correct
92 Correct 98 ms 17424 KB Output is correct
93 Correct 124 ms 18316 KB Output is correct
94 Correct 128 ms 19056 KB Output is correct
95 Correct 130 ms 18780 KB Output is correct
96 Correct 133 ms 18664 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 463 ms 78328 KB Output is correct
2 Correct 727 ms 144500 KB Output is correct
3 Runtime error 1133 ms 262148 KB Execution killed with signal 9
4 Halted 0 ms 0 KB -
# Verdict Execution time Memory Grader output
1 Correct 1400 ms 158320 KB Output is correct
2 Correct 1487 ms 167976 KB Output is correct
3 Correct 1427 ms 180332 KB Output is correct
4 Correct 1434 ms 190864 KB Output is correct
5 Correct 1403 ms 191368 KB Output is correct
6 Correct 1445 ms 174764 KB Output is correct
7 Correct 1388 ms 162928 KB Output is correct
8 Correct 1381 ms 196136 KB Output is correct
9 Correct 1479 ms 189208 KB Output is correct
10 Correct 1866 ms 239952 KB Output is correct
11 Runtime error 822 ms 262148 KB Execution killed with signal 9
12 Halted 0 ms 0 KB -
# Verdict Execution time Memory Grader output
1 Correct 9 ms 12928 KB Output is correct
2 Correct 10 ms 12928 KB Output is correct
3 Correct 10 ms 12928 KB Output is correct
4 Correct 30 ms 16300 KB Output is correct
5 Correct 37 ms 16160 KB Output is correct
6 Correct 32 ms 17024 KB Output is correct
7 Correct 43 ms 17904 KB Output is correct
8 Correct 41 ms 18132 KB Output is correct
9 Correct 42 ms 18208 KB Output is correct
10 Correct 37 ms 17876 KB Output is correct
11 Correct 38 ms 17972 KB Output is correct
12 Correct 40 ms 18080 KB Output is correct
13 Correct 41 ms 18368 KB Output is correct
14 Correct 39 ms 17884 KB Output is correct
15 Correct 43 ms 17916 KB Output is correct
16 Correct 39 ms 17828 KB Output is correct
17 Correct 40 ms 18072 KB Output is correct
18 Correct 38 ms 17736 KB Output is correct
19 Correct 39 ms 17848 KB Output is correct
20 Correct 44 ms 19072 KB Output is correct
21 Correct 42 ms 18124 KB Output is correct
22 Correct 49 ms 20468 KB Output is correct
23 Correct 46 ms 19352 KB Output is correct
24 Correct 50 ms 20552 KB Output is correct
25 Correct 29 ms 17664 KB Output is correct
26 Correct 30 ms 15492 KB Output is correct
27 Correct 23 ms 14976 KB Output is correct
28 Correct 50 ms 20424 KB Output is correct
29 Correct 49 ms 20504 KB Output is correct
30 Correct 54 ms 20788 KB Output is correct
31 Correct 51 ms 20536 KB Output is correct
32 Correct 51 ms 20232 KB Output is correct
33 Correct 49 ms 20432 KB Output is correct
34 Correct 51 ms 20312 KB Output is correct
35 Correct 51 ms 20696 KB Output is correct
36 Correct 28 ms 17536 KB Output is correct
37 Correct 49 ms 20232 KB Output is correct
38 Correct 50 ms 20172 KB Output is correct
39 Correct 48 ms 20484 KB Output is correct
40 Correct 46 ms 20468 KB Output is correct
41 Correct 45 ms 19704 KB Output is correct
42 Correct 46 ms 19824 KB Output is correct
43 Correct 27 ms 17528 KB Output is correct
44 Correct 46 ms 19524 KB Output is correct
45 Correct 45 ms 19760 KB Output is correct
46 Correct 23 ms 15104 KB Output is correct
47 Correct 42 ms 18492 KB Output is correct
48 Correct 37 ms 17256 KB Output is correct
49 Correct 36 ms 17612 KB Output is correct
50 Correct 35 ms 17136 KB Output is correct
51 Correct 34 ms 17396 KB Output is correct
52 Correct 40 ms 18080 KB Output is correct
53 Correct 41 ms 18616 KB Output is correct
54 Correct 41 ms 18748 KB Output is correct
55 Correct 44 ms 18804 KB Output is correct
56 Correct 121 ms 16760 KB Output is correct
57 Correct 96 ms 17352 KB Output is correct
58 Correct 100 ms 16648 KB Output is correct
59 Correct 144 ms 17768 KB Output is correct
60 Correct 123 ms 18060 KB Output is correct
61 Correct 126 ms 18764 KB Output is correct
62 Correct 122 ms 17764 KB Output is correct
63 Correct 116 ms 18204 KB Output is correct
64 Correct 126 ms 18312 KB Output is correct
65 Correct 117 ms 18192 KB Output is correct
66 Correct 117 ms 18972 KB Output is correct
67 Correct 113 ms 17964 KB Output is correct
68 Correct 118 ms 18732 KB Output is correct
69 Correct 114 ms 18236 KB Output is correct
70 Correct 114 ms 18404 KB Output is correct
71 Correct 107 ms 17268 KB Output is correct
72 Correct 120 ms 18252 KB Output is correct
73 Correct 117 ms 18212 KB Output is correct
74 Correct 94 ms 17728 KB Output is correct
75 Correct 89 ms 15744 KB Output is correct
76 Correct 81 ms 15284 KB Output is correct
77 Correct 136 ms 20384 KB Output is correct
78 Correct 146 ms 20412 KB Output is correct
79 Correct 105 ms 17712 KB Output is correct
80 Correct 139 ms 20368 KB Output is correct
81 Correct 157 ms 20176 KB Output is correct
82 Correct 138 ms 19724 KB Output is correct
83 Correct 139 ms 19940 KB Output is correct
84 Correct 95 ms 17484 KB Output is correct
85 Correct 137 ms 19668 KB Output is correct
86 Correct 133 ms 19660 KB Output is correct
87 Correct 92 ms 15268 KB Output is correct
88 Correct 138 ms 18008 KB Output is correct
89 Correct 123 ms 17876 KB Output is correct
90 Correct 136 ms 17952 KB Output is correct
91 Correct 131 ms 17508 KB Output is correct
92 Correct 98 ms 17424 KB Output is correct
93 Correct 124 ms 18316 KB Output is correct
94 Correct 128 ms 19056 KB Output is correct
95 Correct 130 ms 18780 KB Output is correct
96 Correct 133 ms 18664 KB Output is correct
97 Correct 463 ms 78328 KB Output is correct
98 Correct 727 ms 144500 KB Output is correct
99 Runtime error 1133 ms 262148 KB Execution killed with signal 9
100 Halted 0 ms 0 KB -