oj.uz

Submission #491320

# Submission time Handle Problem Language Result Execution time Memory
491320 2021-12-01T13:30:15 Z blue Treatment Project (JOI20_treatment) C++17
35 / 100
 519 ms 41852 KB 
#include <iostream>
#include <algorithm>
#include <vector>
#include <deque>
#include <set>
#include <queue>
#include <cassert>
using namespace std;
 
using ll = long long;
using pll = pair<ll, ll>;
using vpll = vector<pll>;
using vi = vector<int>;
using vll = vector<ll>;
using pii = pair<int, int>;
 
struct project
{
    int T;
    int L;
    int R;
    int C;
    int I;
};
 
 
 
vector<project> P;
 
const int maxM = 100'000;
vi query_ans;
 
vi visited(maxM, 0);
 
const ll INF = 1'000'000'000'000'000'000LL;
 
vll dp(maxM, INF);
 
struct pos
{
    int i;
};
 
bool operator < (pos A, pos B)
{
    if(dp[A.i] == dp[B.i]) return A.i < B.i;
    else return dp[A.i] < dp[B.i];
}
 
 
 
struct segtree
{
    // deque<int> v[220'000];
    vector<int> v[1 << 18];
 
    segtree()
    {
        ;
    }
 
    segtree(int L, int R)
    {
        ;
    }
 
    void add(int i, int l, int r, int I)
    {
        v[i].push_back(I);
        if(l != r)
        {
            if(I <= (l+r)/2) add(2*i, l, (l+r)/2, I);
            else add(2*i + 1, (l+r)/2+1, r, I);
        }
    }
 
    void getU(int i, int l, int r, int L, int R, ll mxv)
    {
        if(R < l || r < L || R < L) return;
        else if(L <= l && r <= R)
        {
            while(1)
            {
                if(v[i].empty()) break;
                if(P[v[i].back()].L + P[v[i].back()].T > mxv) break;
                if(!visited[v[i].back()])
                {
                    query_ans.push_back(v[i].back());
                    visited[v[i].back()] = 1;
                }
                v[i].pop_back();
            }
        }
        else
        {
            getU(2*i, l, (l+r)/2, L, R, mxv);
            getU(2*i+1, (l+r)/2+1, r, L, R, mxv);
        }
    }
 
    void getD(int i, int l, int r, int L, int R, ll mxv)
    {
        if(R < l || r < L || R < L) return;
        else if(L <= l && r <= R)
        {
            while(1)
            {
                if(v[i].empty()) break;
                if(P[v[i].back()].L - P[v[i].back()].T > mxv) break;
                if(!visited[v[i].back()])
                {
                    query_ans.push_back(v[i].back());
                    visited[v[i].back()] = 1;
                }
                v[i].pop_back();
            }
        }
        else
        {
            getD(2*i, l, (l+r)/2, L, R, mxv);
            getD(2*i+1, (l+r)/2+1, r, L, R, mxv);
        }
    }
};
 
segtree U, D;
 
int main()
{
    ios_base::sync_with_stdio(false);
    cin.tie(NULL);
 
    int N, M;
    cin >> N >> M;
 
    P = vector<project>(M);
    for(int i = 0; i < M; i++) cin >> P[i].T >> P[i].L >> P[i].R >> P[i].C;
 
    sort(P.begin(), P.end(), [] (project x, project y)
    {
        return x.T < y.T;
    });
 
    for(int i = 0; i < M; i++) P[i].I = i;
 
 
    sort(P.begin(), P.end(), [] (project x, project y)
    {
        return x.L + x.T <= y.L + y.T;
    });
 
    for(int i = 0; i < M; i++)
        U.add(1, 0, M-1, P[i].I);
 
    sort(P.begin(), P.end(), [] (project x, project y)
    {
        return x.L - x.T < y.L - y.T;
    });
 
    for(int i = 0; i < M; i++)
        D.add(1, 0, M-1, P[i].I);
 
    for(int f = 0; f < (1 << 18); f++)
    {
        reverse(U.v[f].begin(), U.v[f].end());
        reverse(D.v[f].begin(), D.v[f].end());
    }
 
 
 
    sort(P.begin(), P.end(), [] (project x, project y)
    {
        return x.I < y.I;
    });
    set<pos> tbv;
 
    for(int i = 0; i < M; i++)
    {
        if(P[i].L == 1)
        {
            dp[i] = P[i].C;
            visited[i] = 1;
            tbv.insert(pos{i});
        }
    }
 
    while(!tbv.empty())
    {
        int u = tbv.begin()->i;
        tbv.erase(tbv.begin());
        // cerr << "visiting " << u << ' ' << P[u].L << ' ' << P[u].R <<  '\n';
 
        query_ans.clear();
        U.getU(1, 0, M-1, u+1, M-1, P[u].R + P[u].T + 1);
        D.getD(1, 0, M-1, 0, u-1, P[u].R - P[u].T + 1);
 
        for(int v: query_ans)
        {
            // cerr << u << " -> " << v << '\n';
            assert(P[v].L <= P[u].R + 1 - abs(P[u].T - P[v].T));
            assert(dp[v] == INF);
            dp[v] = dp[u] + P[v].C;
            tbv.insert(pos{v});
        }
    }
 
    ll final_ans = INF;
    for(int i = 0; i < M; i++)
        if(P[i].R == N)
            final_ans = min(final_ans, dp[i]);
 
    if(final_ans >= INF) final_ans = -1;
 
    cout << final_ans << '\n';
}

Subtask #1
0 / 4.0

# Verdict Execution time Memory Grader output
1 Correct 519 ms 41852 KB Output is correct
2 Correct 326 ms 41824 KB Output is correct
3 Runtime error 65 ms 31760 KB Execution killed with signal 11
4 Halted 0 ms 0 KB -

Subtask #2
5.0 / 5.0

# Verdict Execution time Memory Grader output
1 Correct 9 ms 13736 KB Output is correct
2 Correct 8 ms 13772 KB Output is correct
3 Correct 11 ms 13696 KB Output is correct
4 Correct 9 ms 13772 KB Output is correct
5 Correct 11 ms 13808 KB Output is correct
6 Correct 9 ms 13772 KB Output is correct
7 Correct 9 ms 13804 KB Output is correct
8 Correct 10 ms 13772 KB Output is correct
9 Correct 9 ms 13772 KB Output is correct
10 Correct 8 ms 13772 KB Output is correct
11 Correct 8 ms 13772 KB Output is correct
12 Correct 8 ms 13772 KB Output is correct
13 Correct 9 ms 13772 KB Output is correct
14 Correct 8 ms 13744 KB Output is correct
15 Correct 8 ms 13772 KB Output is correct
16 Correct 9 ms 13772 KB Output is correct
17 Correct 8 ms 13772 KB Output is correct
18 Correct 8 ms 13772 KB Output is correct
19 Correct 9 ms 13772 KB Output is correct

Subtask #3
30.0 / 30.0

# Verdict Execution time Memory Grader output
1 Correct 9 ms 13736 KB Output is correct
2 Correct 8 ms 13772 KB Output is correct
3 Correct 11 ms 13696 KB Output is correct
4 Correct 9 ms 13772 KB Output is correct
5 Correct 11 ms 13808 KB Output is correct
6 Correct 9 ms 13772 KB Output is correct
7 Correct 9 ms 13804 KB Output is correct
8 Correct 10 ms 13772 KB Output is correct
9 Correct 9 ms 13772 KB Output is correct
10 Correct 8 ms 13772 KB Output is correct
11 Correct 8 ms 13772 KB Output is correct
12 Correct 8 ms 13772 KB Output is correct
13 Correct 9 ms 13772 KB Output is correct
14 Correct 8 ms 13744 KB Output is correct
15 Correct 8 ms 13772 KB Output is correct
16 Correct 9 ms 13772 KB Output is correct
17 Correct 8 ms 13772 KB Output is correct
18 Correct 8 ms 13772 KB Output is correct
19 Correct 9 ms 13772 KB Output is correct
20 Correct 20 ms 15188 KB Output is correct
21 Correct 18 ms 15112 KB Output is correct
22 Correct 20 ms 15052 KB Output is correct
23 Correct 20 ms 15060 KB Output is correct
24 Correct 24 ms 15308 KB Output is correct
25 Correct 22 ms 15060 KB Output is correct
26 Correct 19 ms 15052 KB Output is correct
27 Correct 19 ms 15128 KB Output is correct
28 Correct 21 ms 15328 KB Output is correct
29 Correct 22 ms 15052 KB Output is correct
30 Correct 21 ms 15180 KB Output is correct
31 Correct 15 ms 15052 KB Output is correct
32 Correct 20 ms 15308 KB Output is correct
33 Correct 21 ms 15284 KB Output is correct
34 Correct 23 ms 15160 KB Output is correct
35 Correct 22 ms 15248 KB Output is correct
36 Correct 20 ms 15328 KB Output is correct
37 Correct 21 ms 15160 KB Output is correct

Subtask #4
0 / 61.0

# Verdict Execution time Memory Grader output
1 Correct 519 ms 41852 KB Output is correct
2 Correct 326 ms 41824 KB Output is correct
3 Runtime error 65 ms 31760 KB Execution killed with signal 11
4 Halted 0 ms 0 KB -