Submission #842424

#	Time	Username	Problem	Language	Result	Execution time	Memory
842424		arnold518	Overtaking (IOI23_overtaking)	C++17	19 / 100	3 ms	15708 KiB

overtaking

#include "overtaking.h"
#include <bits/stdc++.h>
using namespace std;

typedef long long ll;
typedef pair<int, int> pii;
typedef pair<ll, ll> pll;

const int MAXN = 1000;
const ll INF = 1e18;

int L, N;
pll P[MAXN+10];
int M, S[MAXN+10];
ll X;

int C[MAXN+10];
ll A[MAXN+10], B[MAXN+10], A2[MAXN+10][MAXN+10], ans[MAXN+10][MAXN+10];

struct Frac
{
    ll a, b;
    Frac(ll a, ll b) : a(a), b(b) {}
};
bool operator <= (Frac p, Frac q)
{
    if(p.a<0) p.a*=-1, p.b*=-1;
    if(q.a<0) q.a*=-1, q.b*=-1;
    return (__int128)p.b*q.a<=(__int128)q.b*p.a;
}

struct CHT
{
    vector<pair<pll, int>> V;
    CHT() {}
    Frac cross(pll p, pll q) { return {p.first-q.first, q.second-p.second}; }
    Frac f(pll p, ll y) { return {p.first, y-p.second}; }
    ll f2(pll p, ll y) { return (y-p.second+p.first-1)/p.first; }
    void push(pll p, int q)
    {
        while(!V.empty() && (V.back().first.first==p.first || cross(V[V.size()-1].first, p)<=Frac(1, 0))) V.pop_back();
        while(V.size()>1 && (V.back().first.first==p.first || cross(V[V.size()-2].first, V[V.size()-1].first)<=cross(V[V.size()-1].first, p))) V.pop_back();
        V.push_back({p, q});
    }
    pll query(ll y)
    {
        //while(V.size()>1 && f(V[V.size()-2].first, y)<=f(V[V.size()-1].first, y)) V.pop_back();
        Frac t={1, INF}; int ret=0;
        for(int i=0; i<V.size(); i++) if(f(V[i].first, y)<=t) t=f(V[i].first, y), ret=i;
        if(V.empty()) return {0, 0};
        else
        {
            return {V[ret].second, f2(V[ret].first, y)};
        }
    }
};
vector<ll> V1;
vector<CHT> V2;

void init(int _L, int _N, vector<ll> _T, vector<int> _W, int _X, int _M, vector<int> _S)
{
    L=_L; N=_N;
    for(int i=1; i<=N; i++) P[i]=pll(_W[i-1], _T[i-1]);
    M=_M; X=_X;
    for(int i=0; i<M; i++) S[i]=_S[i];
    M--;

    sort(P+1, P+N+1);
    reverse(P+1, P+N+1);
    while(N && P[N].first<=X) N--;

    //for(int i=1; i<=N; i++) printf("%lld %lld\n", P[i].first, P[i].second);

    for(int i=1; i<=N; i++) C[i]=i;
    for(int i=1; i<=N; i++) A[i]=P[i].second;
    for(int i=1; i<=N; i++) A2[i][0]=A[i];
    for(int i=1; i<=M; i++)
    {
        for(int j=1; j<=N; j++) B[j]=A[j]+(ll)(S[i]-S[i-1])*P[j].first;
        sort(C+1, C+N+1, [&](const int &p, const int &q) { return pll(A[p], B[p])<pll(A[q], B[q]); });
        ll t=0;
        for(int j=1; j<=N; j++) A[C[j]]=t=max(t, B[C[j]]);
        for(int j=1; j<=N; j++) A2[j][i]=A[j];
    }

    for(int i=M; i>=0; i--)
    {
        sort(C+1, C+N+1, [&](const int &p, const int &q) { return A2[p][i]<A2[q][i]; });
        CHT cht;
        for(int j=1; j<=N;)
        {
            auto p=cht.query(A2[C[j]][i]); p.second+=S[i];
            //printf("%d %lld : %lld %lld\n", C[j], A2[C[j]][i], p.first, p.second);
            ll val;
            if(p.first==0) val=A2[C[j]][i]+X*(L-S[i]);
            else
            {
                if(p.second>L) val=A2[C[j]][i]+X*(L-S[i]);
                else
                {
                    int it=lower_bound(S, S+M+1, p.second)-S;
                    val=ans[p.first][it];
                }
            }
            ll tt=A2[C[j]][i];
            for(; j<=N && A2[C[j]][i]==tt; j++)
            {
                ans[C[j]][i]=val;
                //printf("!%d : %lld %lld\n", C[j], P[C[j]].first-X, A2[C[j]][i]);
                cht.push({P[C[j]].first-X, A2[C[j]][i]}, C[j]);
            }
            if(i==0) V1.push_back(tt), V2.push_back(cht);
        }
        //printf("=======\n");
    }
    /*
    for(int i=0; i<=M; i++)
    {
        for(int j=1; j<=N; j++) printf("%lld ", ans[j][i]);
        printf("\n");
    }
    for(int i=0; i<=M; i++)
    {
        for(int j=1; j<=N; j++) printf("%lld ", A2[j][i]);
        printf("\n");
    }
    */
}

ll arrival_time(ll Y)
{
    auto it=lower_bound(V1.begin(), V1.end(), Y);
    if(it==V1.begin()) return Y+X*L;
    it--;
    auto &cht = V2[it-V1.begin()];
    /*
    int lo=0, hi=cht.V.size();
    while(lo+1<hi)
    {
        int mid=lo+hi>>1;
        if(cht.f(cht.V[mid].first, Y)<=cht.f(cht.V[mid-1].first, Y)) lo=mid;
        else hi=mid;
    }
    pll p;
    if(cht.V.empty()) p={0, 0};
    else p=pll(cht.V[lo].second, cht.f2(cht.V[lo].first, Y));
*/
    //for(auto it : cht.V) printf("%d %lld\n", it.second, cht.f({P[it.second].first-X, A2[it.second][0]}, Y));
    //for(int i=1; i<=N; i++) if(A2[i][0]<Y) printf("%d %lld %.10lf\n", i, cht.f({P[i].first-X, A2[i][0]}, Y), (double)(Y-A2[i][0])/(P[i].first-X));
    //for(auto it : cht.V) printf("!%d %lld %lld\n", it.second, it.first.first, it.first.second);
    pll p=cht.query(Y);
    //printf("%lld %lld\n", p.first, p.second);
    if(p.first==0) return Y+X*L;
    else
    {
        if(p.second>L) return Y+X*L;
        else
        {
            int it=lower_bound(S, S+M+1, p.second)-S;
            //printf("?%d %d\n", p.first ,it);
            return ans[p.first][it];
        }
    }
}

Compilation message (stderr)

overtaking.cpp: In member function 'pll CHT::query(ll)':
overtaking.cpp:49:23: warning: comparison of integer expressions of different signedness: 'int' and 'std::vector<std::pair<std::pair<long long int, long long int>, int> >::size_type' {aka 'long unsigned int'} [-Wsign-compare]
   49 |         for(int i=0; i<V.size(); i++) if(f(V[i].first, y)<=t) t=f(V[i].first, y), ret=i;
      |                      ~^~~~~~~~~

• Subtask 0: 0 / 0
• Subtask 1: 9 / 9
• Subtask 2: 10 / 10
• Subtask 3: 0 / 20
• Subtask 4: 0 / 26
• Subtask 5: 0 / 35