Submission #555719

#	Time	Username	Problem	Language	Result	Execution time	Memory
555719		urosk	K blocks (IZhO14_blocks)	C++14	14 / 100	1 ms	340 KiB

blocks

// __builtin_popcount(x)
// __builtin_popcountll(x)
#define here cerr<<"===========================================\n"
#include <bits/stdc++.h>
#define ld double
#define ll long long
#define ull unsigned long long
#define llinf 100000000000000000LL // 10^17
#define iinf 2000000000 // 2*10^9
#define pb push_back
#define popb pop_back
#define fi first
#define sc second
#define endl '\n'
#define pii pair<int,int>
#define pll pair<ll,ll>
#define pld pair<ld,ld>
#define sz(a) int(a.size())
#define all(a) a.begin(),a.end()
#define ceri(a,l,r) {for(ll i_ = l;i_<=r;i_++) cerr<<a[i_]<< " ";cerr<<endl;}
using namespace std;

#define maxn 100005
#define maxk 105
#define lg 25
ll n,k;
ll a[maxn];
ll dp[maxn];
ll ndp[maxn];
ll st[maxn][lg];
ll powi[maxn];
bool test = 0;
ll get(ll l,ll r){
    ll len = (r-l+1);
    ll k = powi[len];
    ll ans = max(st[l][k],st[r-(1<<k)+1][k]);
    if(test){
        ll g =0 ;
        for(ll i = l;i<=r;i++) g = max(g,a[i]);
        if(g!=ans){
            cerr<<"MIN: "<<l<< " "<<r<<endl;
            exit(0);
        }
    }
    return ans;
}
void reshi(ll l,ll r,ll optl,ll optr){
    if(l>r) return;
    ll mid = (l+r)/2;
    ndp[mid] = dp[mid-1] + a[mid];
    ll opt = mid-1;
    for(ll i = optl;i<=min(mid-1,optr);i++){
        if(dp[i] + get(i+1,mid) < ndp[mid]){
            opt = i;
            ndp[mid] = dp[i] + get(i+1,mid);
        }
    }
    //cerr<<"i: "<<mid<< " "<<optl<< " "<<optr<<" "<<opt<<endl;
    //cerr<<l<< " "<<r<< " "<<optl<< " "<<optr<<" "<<mid<< " "<<opt<<endl;
    reshi(l,mid-1,opt,optr);
    reshi(mid+1,r,optl,opt);
}
ll brut(){
    vector<vector<ll> > v(n+1,vector<ll>(k+1,0));
    v[0][0] = 0;
    for(ll i = 1;i<=n;i++) v[i][0] = llinf;
    for(ll j = 1;j<=k;j++){
        v[0][j] = llinf;
        for(ll i = 1;i<=n;i++){
            v[i][j] = llinf;
            ll mx = a[i];
            for(ll f = i-1;f>=0;f--){
                v[i][j] = min(v[i][j],v[f][j-1] + mx);
                mx = max(mx,a[f]);
            }
        }
        for(ll i = 0;i<=n;i++) cerr<<v[i][j]<< " ";
        cerr<<endl;
    }
    return v[n][k];
}
mt19937_64 rng(chrono::steady_clock::now().time_since_epoch().count());
ll rnd(ll l,ll r){
    return uniform_int_distribution<ll>(l,r)(rng);
}
bool gen = 0;
int main(){
    test = 0;
	ios_base::sync_with_stdio(false);cerr.tie(0);cout.tie(0);cin.tie(0);
    if(gen==0){
        cin >> n >> k;
        for(ll i = 1;i<=n;i++) cin >> a[i];
    }else{
        n = 7;
        k = 4;
        for(ll i = 1;i<=n;i++) a[i] = rnd(1,30);
        cerr<<n<< " "<<k<<endl;
        for(ll i = 1;i<=n;i++) cerr<<a[i]<< " ";
        cerr<<endl;
    }
    for(ll i = 1;i<=n;i++) st[i][0] = a[i];
    for(ll j = 1;j<lg;j++){
        for(ll i = 1;i+(1<<j)-1<=n;i++){
            st[i][j] = max(st[i][j-1],st[i+(1<<(j-1))][j-1]);
        }
    }

    for(ll j = 0;j<lg;j++){
        if((1<<j)>n) break;
        powi[(1<<j)] = j;
    }
    for(ll i = 1;i<=n;i++){
        if(powi[i]==0) powi[i] = powi[i-1];
    }
    fill(dp,dp+n+1,llinf);
    dp[0] = 0;
    for(ll ii = 1;ii<=k;ii++){
        fill(ndp,ndp+n+1,llinf);
        //here;
        reshi(1,n,0,n);
        //cerr<<"dp: ";
        for(ll i = 1;i<=n;i++) dp[i] = ndp[i];
        for(ll i = 0;i<ii;i++) dp[i] = llinf;
        //ceri(ndp,0,n);
    }
    //here;
    if(test)cerr<<"brut: "<<brut()<<endl;
    cout<<dp[n]<<endl;
	return 0;
}
/*
5 2
1 2 3 4 5
out: 6
7 4
20 1 10 19 12 8 22
out: 53
*/

• Subtask 1: 14 / 14
• Subtask 2: 0 / 18
• Subtask 3: 0 / 21
• Subtask 4: 0 / 47