oj.uz

Submission #992525

# Submission time Handle Problem Language Result Execution time Memory
992525 2024-06-04T14:34:51 Z PedroBigMan Islands (IOI08_islands) C++14
90 / 100
 468 ms 131072 KB 
#pragma GCC optimization ("O3")
#pragma GCC optimization ("unroll-loops")
#pragma GCC optimize("Ofast")
#include <iostream>
#include <vector>
#include <cmath>
#include <algorithm>
#include <string>
#include <map>
#include <unordered_map>
#include <set>
#include <unordered_set>
#include <queue>
#include <deque>
#include <list>
#include <iomanip>
#include <stdlib.h>
#include <time.h>
#include <cstring>
using namespace std;
typedef int ll;
typedef long long int lld;
typedef unsigned long long int ull;
typedef long double ld;
#define REP(i,a,b) for(ll i=(ll) a; i<(ll) b; i++)
#define pb push_back
#define mp make_pair
#define pl pair<ll,ll>
#define ff first
#define ss second
#define whole(x) x.begin(),x.end()
#define DEBUG(i) cout<<"Pedro "<<i<<endl
#define INF 1000000000000000000LL
#define EPS ((ld)0.00000000001)
#define pi ((ld)3.141592653589793)
#define VV(vvvv,NNNN,xxxx); REP(iiiii,0,NNNN) {vvvv.pb(xxxx);}
ll mod=1000000007;

template<class A=ll> 
void Out(vector<A> a) {REP(i,0,a.size()) {cout<<a[i]<<" ";} cout<<endl;}

template<class A=ll>
void In(vector<A> &a, ll N) {A cur; REP(i,0,N) {cin>>cur; a.pb(cur);}} 

vector<ll> f, w; vector<bool> visited;
vector<vector<ll> > adj; ll nei;
lld bestpath;
vector<lld> rm; vector<lld> ps;

vector<ll> Floyd(ll s) //Floyd's Algorithm, O(N) time, O(1) memory, return <element of cycle,length od cycle>
{
    ll a=f[s]; ll b=f[f[s]];
    while(a!=b) {a=f[a]; b=f[f[b]];}
    vector<ll> ans; ans.pb(a); b=f[a];
    while(b!=a) {ans.pb(b); b=f[b];}
    return ans;
}

lld DFS(ll s)
{
    visited[s]=true; lld ans = 0; lld thisval;
    REP(i,0,adj[s].size())
    {
        nei = adj[s][i];
        if(visited[nei]) {continue;}
        thisval = (lld)(w[nei]+DFS(nei));
        bestpath = max(bestpath,ans+thisval);
        ans = max(ans,thisval);
    }
    return ans;
}

void Rolling_Max(ll M, ll L) //ans[i]=j means aj<ai, j<i and j is the largest that fullfills these conditions. If no j exists then a[i]=-1. 
{
	deque<ll> chain; 
    REP(i,0,L)
    {
        while(chain.size()>0 && ps[i]>=ps[chain[chain.size()-1]]) {chain.pop_back();}
        chain.pb(i);
        if(chain.front()==i-M) {chain.pop_front();}
        rm[i]=ps[chain.front()];
    }	
} 


int main()
{
    ios_base::sync_with_stdio(0);
    cin.tie(0); cout.tie(0);
	cout.precision(20);
	ll N; cin>>N; f = vector<ll>(N,-1); w = vector<ll>(N,-1); visited = vector<bool>(N,false); adj=vector<vector<ll> >(N,vector<ll>());
    REP(i,0,N) {cin>>f[i]>>w[i]; f[i]--; adj[i].pb(f[i]); adj[f[i]].pb(i);}
    lld ans = 0LL; ll node; ll C; rm = vector<lld>(2*N,0); ps = vector<lld>(2*N,0);
    REP(s,0,N)
    {
        if(visited[s]) {continue;}
        lld bestval = 0; lld thisval;
        vector<ll> cycle = Floyd(s); C = cycle.size();
        REP(i,0,cycle.size()) {node=cycle[i]; visited[node]=true;}
        vector<lld> val(C,0);
        REP(i,0,cycle.size()) 
        {
            bestpath = 0;
            val[i]=DFS(cycle[i]);
            bestval = max(bestval,bestpath);
        }
        lld cursum=0;
        REP(i,0,2*C) {node = cycle[i%C]; ps[i]=(cursum+val[i%C]); cursum+=w[node];}
        Rolling_Max(C-1,2*C);
        REP(i,0,C)
        {
            thisval=rm[i+C-1] -(ps[i]-val[i])+val[i];
            bestval = max(thisval,bestval);
        }
        ans+=bestval;
    }
    cout<<ans<<endl;
    return 0;
}

Compilation message

islands.cpp:1: warning: ignoring '#pragma GCC optimization' [-Wunknown-pragmas]
    1 | #pragma GCC optimization ("O3")
      | 
islands.cpp:2: warning: ignoring '#pragma GCC optimization' [-Wunknown-pragmas]
    2 | #pragma GCC optimization ("unroll-loops")
      |

Subtask #1
33.0 / 33.0

# Verdict Execution time Memory Grader output
1 Correct 0 ms 348 KB Output is correct
2 Correct 0 ms 348 KB Output is correct
3 Correct 0 ms 348 KB Output is correct
4 Correct 0 ms 348 KB Output is correct
5 Correct 0 ms 348 KB Output is correct
6 Correct 0 ms 344 KB Output is correct
7 Correct 0 ms 348 KB Output is correct
8 Correct 0 ms 348 KB Output is correct
9 Correct 0 ms 348 KB Output is correct
10 Correct 0 ms 348 KB Output is correct
11 Correct 0 ms 348 KB Output is correct

Subtask #2
3.0 / 3.0

# Verdict Execution time Memory Grader output
1 Correct 1 ms 348 KB Output is correct
2 Correct 1 ms 348 KB Output is correct

Subtask #3
4.0 / 4.0

# Verdict Execution time Memory Grader output
1 Correct 1 ms 604 KB Output is correct
2 Correct 1 ms 860 KB Output is correct

Subtask #4
10.0 / 10.0

# Verdict Execution time Memory Grader output
1 Correct 4 ms 1884 KB Output is correct
2 Correct 10 ms 3916 KB Output is correct
3 Correct 6 ms 2396 KB Output is correct
4 Correct 3 ms 1368 KB Output is correct

Subtask #5
10.0 / 10.0

# Verdict Execution time Memory Grader output
1 Correct 13 ms 5468 KB Output is correct
2 Correct 20 ms 10076 KB Output is correct

Subtask #6
10.0 / 10.0

# Verdict Execution time Memory Grader output
1 Correct 41 ms 19804 KB Output is correct
2 Correct 45 ms 23284 KB Output is correct
3 Correct 54 ms 25704 KB Output is correct

Subtask #7
10.0 / 10.0

# Verdict Execution time Memory Grader output
1 Correct 77 ms 38388 KB Output is correct
2 Correct 98 ms 49860 KB Output is correct
3 Correct 118 ms 60784 KB Output is correct
4 Correct 132 ms 64704 KB Output is correct

Subtask #8
10.0 / 10.0

# Verdict Execution time Memory Grader output
1 Correct 164 ms 98248 KB Output is correct
2 Correct 427 ms 101148 KB Output is correct
3 Correct 159 ms 76192 KB Output is correct
4 Correct 206 ms 98504 KB Output is correct
5 Correct 218 ms 98244 KB Output is correct
6 Correct 468 ms 91220 KB Output is correct
7 Correct 212 ms 102964 KB Output is correct

Subtask #9
0 / 10.0

# Verdict Execution time Memory Grader output
1 Runtime error 189 ms 131072 KB Execution killed with signal 9
2 Halted 0 ms 0 KB -