Submission #248760

#	Time	Username	Problem	Language	Result	Execution time	Memory
248760		MarcoMeijer	Wiring (IOI17_wiring)	C++14	100 / 100	843 ms	41440 KiB

wiring

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

//macros
typedef long long ll;
typedef pair<ll, ll> ii;
typedef pair<ll, ll> lll;
typedef tuple<ll, ll, ll> iii;
typedef vector<ll> vi;
typedef vector<ii> vii;
typedef vector<iii> viii;
typedef vector<ll> vll;
typedef vector<lll> vlll;
#define REP(a,b,c) for(ll a=ll(b); a<ll(c); a++)
#define RE(a,c) REP(a,0,c)
#define RE1(a,c) REP(a,1,c+1)
#define REI(a,b,c) REP(a,b,c+1)
#define REV(a,b,c) for(ll a=ll(c-1); a>=ll(b); a--)
#define INF 1e18
#define pb push_back
#define fi first
#define se second

const ll MX = 1e6;

vector<int> f[2];
ll N[2], n = 0;
ll a[MX], b[MX], sm[MX], ch[MX], CENT[MX];
ll dp[MX];

struct Seg {
    void create(int size) {
        sz = size;
        SEG.assign(sz*4, 0);
        LAZY.assign(sz*4, 0);
    }

    void add(int i, int j, ll v, ll lazy=0, int p=0, int l=0, int r=-1) {
        if(r == -1) r=sz-1;
        SEG[p] += lazy;
        LAZY[p] += lazy;
        if(j < l || i > r) return;
        if(i <= l && j >= r) {
            SEG[p] += v;
            LAZY[p] += v;
            return;
        }
        int m=(l+r)/2;
        add(i,j,v,LAZY[p],p*2+1,l,m);
        add(i,j,v,LAZY[p],p*2+2,m+1,r);
        SEG[p] = min(SEG[p*2+1], SEG[p*2+2]);
        LAZY[p] = 0;
    }
    ll get(int i, int j, ll lazy=0, int p=0, int l=0, int r=-1) {
        if(r == -1) r=sz-1;
        SEG[p] += lazy;
        LAZY[p] += lazy;
        if(j < l || i > r) return INF;
        if(i <= l && j >= r) return SEG[p];
        int m=(l+r)/2;
        ll a=get(i,j,LAZY[p],p*2+1,l,m);
        ll b=get(i,j,LAZY[p],p*2+2,m+1,r);
        LAZY[p] = 0;
        return min(a,b);
    }

    vi SEG, LAZY;
    int sz;
};

Seg seg[2];

ll min_total_length(std::vector<int> R, std::vector<int> B) {
    f[0] = R; f[1] = B;
    RE(i,2) N[i]=f[i].size();

    // fill a and b
    priority_queue<ii, vii, greater<ii>> pq;
    RE(i,2) RE(j,N[i]) pq.push({f[i][j], i});
    while(!pq.empty()) {
        ii p = pq.top(); pq.pop();
        a[n] = p.fi; b[n] = p.se;
        n++;
    }

    sm[0] = 0; ch[0] = 0; CENT[0] = 0;
    REP(i,1,n+1) {
        ch[i] = ch[i-1];
        sm[i] = sm[i-1]+a[i-1];
        CENT[i] = CENT[i-1];
        if(b[i] != b[i-1]) ch[i]++, CENT[i]=i;
    }

    RE(i,2) seg[i].create(n);
    REP(i,1,n+1) {
        dp[i] = INF;
        ll ed = i-1;
        ll cent = CENT[ed];
        ll beg = CENT[max(0ll, cent-1)];
        if(cent == ed) {
            REP(j,beg,cent) RE(k,2) seg[k].add(beg,j,-a[j]);
            REP(j,beg,cent) seg[0].add(j,j,-a[cent  ]*j);
            REP(j,beg,cent) seg[1].add(j,j,-a[cent-1]*j);
        }
        ll switchJ = 2*cent-ed-1;
        dp[i] = min(dp[i], seg[0].get(beg,switchJ-1) + sm[ed+1]-sm[cent] + a[cent  ]*switchJ);
        dp[i] = min(dp[i], seg[1].get(switchJ,cent-1) + sm[ed+1]-sm[cent] + a[cent-1]*switchJ);
        RE(j,2) seg[j].add(i,i,dp[i]);
        if(dp[i] < dp[i-1]) RE(j,2) seg[j].add(i-1,i-1,dp[i] - dp[i-1]);
    }

    return dp[n];
}

• Subtask 1: 7 / 7
• Subtask 2: 13 / 13
• Subtask 3: 10 / 10
• Subtask 4: 25 / 25
• Subtask 5: 45 / 45