Submission #53783

# Submission time Handle Problem Language Result Execution time Memory
53783 2018-07-01T07:27:10 Z 노영훈(#1436) Islands (IOI08_islands) C++11
18 / 100
1380 ms 190752 KB
#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef pair<int, ll> pil;
const int MX=1000010;
const ll inf=2e17;

int n, A[MX], C[MX];
vector<pil> G[MX];
bool done[MX];

ll D[MX];

bool vis[MX];

// list all vertices
void dfs1(int v, vector<int> &V){
    done[v]=true;
    V.push_back(v);
    for(pil &e:G[v]){
        int x=e.first;
        if(done[x]) continue;
        dfs1(x,V);
    }
}

// find cycle 
int found=-1;
int dep[MX], now;
void dfs2(int v, int p, vector<int> &V){
    dep[v]=++now;
    for(pil &e:G[v]){
        int x; ll c; tie(x,c)=e;
        if(p==x) continue;
        if(dep[x]>=0)
            found=x;
        else
            dfs2(x,v,V);
        if(found>=0){
            if(dep[v]>=dep[found]) V.push_back(v);
            return;
        }
    }
    --now;
}

// farthest point
ll dist[MX]; 
int dfs3(int v){
    int res=v;
    for(pil &e:G[v]){
        int x; ll c; tie(x,c)=e;
        if(dist[x]>=0) continue;
        dist[x]=dist[v]+c;
        int y=dfs3(x);
        if(dist[y]>dist[res]) res=y;
    }
    return res;
}
ll solve_tree(vector<int> &V){
    int s=V[0];

    for(int v:V) dist[v]=-1;
    dist[s]=0;
    int u=dfs3(s);

    for(int v:V) dist[v]=-1;
    dist[u]=0;
    int v=dfs3(u);

    return dist[v];
}

bool incyc[MX];
ll arm[MX];
// longest path from v in cyc
void dfs4(int v){
    arm[v]=0;
    for(pil &e:G[v]){
        int x; ll c; tie(x,c)=e;
        if(incyc[x] || arm[x]>=0) continue;
        dfs4(x);
        arm[v]=max(arm[v], arm[x]+c);
    }
}

ll Val[MX], Len[MX], Sum[MX];

vector<int> V, cyc;
ll solve(int s){
    // cout<<"\nSOLVE ON "<<s<<'\n';
    V.clear(), cyc.clear();
    dfs1(s, V);

    found=-1, now=0;
    for(int v:V) dep[v]=-1;
    dfs2(s, -1, cyc);

    if(cyc.empty())
        return solve_tree(V);

    for(int v:V) arm[v]=-1;
    for(int v:cyc) incyc[v]=true, arm[v]=0;
    for(int v:cyc) dfs4(v);

    int sz=cyc.size();
    Sum[0]=0;
    for(int i=0; i<sz; i++){
        int v=cyc[i];
        Val[i]=arm[v];
        int nxt=cyc[(i+1)%sz];
        for(pil &e:G[v])
            if(e.first==nxt) Len[i]=e.second;
        Sum[i+1]=Sum[i]+Len[i];
    }
    // ll res=0;
    deque<int> Q;
    for(int i=0; i<sz-1; i++){
        while(!Q.empty()){
            int x=Q.back();
            ll a=Val[x]-Sum[x];
            ll b=Val[i]-Sum[i];
            if(a<=b) Q.pop_back();
            else break;
        }
        Q.push_back(i);
    }
    ll res=Val[Q.front()]-Sum[Q.front()]+Val[sz-1]-Sum[sz-1];
    while(!Q.empty()){
        int x=Q.back();
        ll a=Val[x]-Sum[x];
        ll b=Val[sz-1]-Sum[sz-1];
        if(a<=b) Q.pop_back();
        else break;
    }
    Q.push_back(sz-1);
    for(int i=sz; i<=2*sz-2; i++){
        while(!Q.empty() && Q.front()<=i-sz) Q.pop_front();
        int x=Q.front();
        ll now=Val[x]-Sum[x]+ Sum[i%sz]+Sum[sz]+Val[i%sz];
        res=max(res, now);
    }
    return res;
    for(int i=0; i<sz; i++)
        for(int j=i+1; j<sz; j++){
            ll dist=max(Sum[j]-Sum[i], Sum[sz]+Sum[i]-Sum[j]);
            ll now=dist+Val[i]+Val[j];

            // printf("%d - %d : %lld (%lld)\n", cyc[i], cyc[j], now, dist);
            res=max(res, now);
        }
    return res;
}

int main(){
    ios::sync_with_stdio(0); cin.tie(0);
    cin>>n;
    for(int u=1; u<=n; u++){
        int v, c; cin>>v>>c;
        A[u]=v, C[u]=c;
        if(v<u && A[v]==u){
            C[u]=max(C[u], C[v]);
            C[v]=-1;
        }
    }
    for(int i=1; i<=n; i++){
        if(C[i]<0) continue;
        G[i].push_back({A[i], C[i]});
        G[A[i]].push_back({i, C[i]});
    }
    ll ans=0;
    for(int i=1; i<=n; i++){
        if(done[i]) continue;
        ans+=solve(i);
    }
    cout<<ans;
    return 0;
}
# Verdict Execution time Memory Grader output
1 Incorrect 22 ms 23932 KB Output isn't correct
2 Incorrect 23 ms 23932 KB Output isn't correct
3 Correct 21 ms 24076 KB Output is correct
4 Correct 21 ms 24076 KB Output is correct
5 Correct 21 ms 24100 KB Output is correct
6 Incorrect 21 ms 24100 KB Output isn't correct
7 Incorrect 24 ms 24124 KB Output isn't correct
8 Incorrect 22 ms 24128 KB Output isn't correct
9 Incorrect 23 ms 24252 KB Output isn't correct
10 Correct 25 ms 24252 KB Output is correct
11 Correct 25 ms 24252 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 21 ms 24300 KB Output is correct
2 Correct 21 ms 24300 KB Output is correct
# Verdict Execution time Memory Grader output
1 Incorrect 22 ms 24316 KB Output isn't correct
2 Halted 0 ms 0 KB -
# Verdict Execution time Memory Grader output
1 Correct 32 ms 25540 KB Output is correct
2 Incorrect 45 ms 28408 KB Output isn't correct
3 Halted 0 ms 0 KB -
# Verdict Execution time Memory Grader output
1 Incorrect 51 ms 29944 KB Output isn't correct
2 Halted 0 ms 0 KB -
# Verdict Execution time Memory Grader output
1 Incorrect 144 ms 41252 KB Output isn't correct
2 Halted 0 ms 0 KB -
# Verdict Execution time Memory Grader output
1 Correct 191 ms 54652 KB Output is correct
2 Incorrect 289 ms 81980 KB Output isn't correct
3 Halted 0 ms 0 KB -
# Verdict Execution time Memory Grader output
1 Correct 393 ms 99768 KB Output is correct
2 Correct 1380 ms 152600 KB Output is correct
3 Incorrect 360 ms 152600 KB Output isn't correct
4 Halted 0 ms 0 KB -
# Verdict Execution time Memory Grader output
1 Correct 585 ms 157136 KB Output is correct
2 Correct 609 ms 157136 KB Output is correct
3 Correct 574 ms 190752 KB Output is correct
4 Correct 373 ms 190752 KB Output is correct
5 Incorrect 499 ms 190752 KB Output isn't correct
6 Halted 0 ms 0 KB -