oj.uz

Submission #554676

# Submission time Handle Problem Language Result Execution time Memory
554676 2022-04-29T06:57:11 Z kwongweng Split the sequence (APIO14_sequence) C++17
100 / 100
 1306 ms 86540 KB 
/*
Solution for APIO 2014 - Sequence
Tags : dp, Convex-hull Trick (CHT)
*/
 
#include <bits/stdc++.h>
using namespace std;
 
#pragma GCC target ("avx2")
#pragma GCC optimization ("Ofast")
#pragma GCC optimization ("unroll-loops")
 
typedef long long ll;
typedef vector<int> vi;
typedef pair<ll, ll> ii;
typedef vector<ii> vii;
typedef long double ld;
typedef pair<ll, ll> pll;
#define FOR(i, a, b) for(int i = a; i < b; i++)
#define ROF(i, a, b) for(int i = a; i >= b; i--)
#define ms memset
#define pb push_back
#define fi first
#define se second
 
ll MOD = 1000000007;
 
ll power(ll base, ll n){
	if (n == 0) return 1;
	if (n == 1) return base;
	ll halfn = power(base, n/2);
	if (n % 2 == 0) return (halfn * halfn) % MOD;
	return (((halfn * halfn) % MOD) * base) % MOD;
}
 
ll inverse(ll n){
	return power(n, MOD-2);
}
 
ll add(ll a, ll b){
	return (a+b) % MOD;
}
 
ll mul(ll a, ll b){
	a %= MOD;
	return (a*b) % MOD;
}
 
ll gcd(ll a, ll b){
    if (a == 1) return 1;
    if (a == 0) return b;
    return gcd(b%a, a);
}
 
 
const int N = 100001;
vector<ll> m(N), c(N);
 
ld g(ii l){
    ll l1 = l.fi; ll l2 = l.se;
    ld a = c[l1]-c[l2];
    ld b = m[l2]-m[l1];
    if (b==0) return -MOD;
    return (ld) a/b;
}
 
void solve(){
    // CHT to optimise O(n^2 * k) into O(n*k)
    int n, k; cin >> n >> k;
    vector<ll> a(n+1); FOR(i,1,n+1) cin >> a[i];
    vector<ll> s(n+1); FOR(i,1,n+1) s[i]=a[i]+s[i-1];
    ll dp[n+1][2]; // replace usual dp[n+1][k+1] to reduce memory usage for last subtask
    int pos[n+1][k+1];
    ms(dp,-1,sizeof(dp));
    ms(pos,-1,sizeof(pos));
    FOR(i,1,n+1){
        dp[i][0]=s[i]*(s[n]-s[i]);
        // only 1 component
    }
    FOR(j,1,k+1){
        // m_l = 2*s[l]+s[n], m_l non-decreasing
        // x = s[i]
        // c_l = dp[l][j-1] - s[l] * (s[n]+s[l])
        // f_l(x) = x^2 + m_l x + c_l
        // l1 < l2, g(l1, l2) = (c_l1-c_l2)/(m_l2-m_l1)
        // f_l1(x) <= f_l2(x) <==>  g(l1, l2) <= x
        // f_l1(x) >= f_l2(x) || f_l2(x) <= f_l3(x)
        // g(l1, l2) >= x || g(l2, l3) <= x
        // g(l1, l2) >= g(l2, l3) <==> l2 ignored
        
        FOR(i,1,n+1){
            dp[i][j%2]=-1; // dp values to compute later on
            m[i] = 2*s[i]+s[n];
            c[i] = dp[i][(j+1)%2] - s[i] * (s[n] + s[i]);
        }
        int best_l = 1;
        // edge case : i = 2, only 1 possible value
        ll val = dp[1][(j+1)%2] + (s[2]-s[1])*(s[n]-s[2]+s[1]);
        if (val > dp[2][j%2] && dp[1][(j+1)%2] != -1){
            dp[2][j%2]=val;
            pos[2][j]=best_l;
        }
        list <ii> li; // stores pairs with increasing g(l1, l2) from left to right.
        int r_l = 1;
        FOR(i,3,n+1){
            best_l = r_l;
            ld cur;
            if (dp[i-1][(j+1)%2] != -1){
                r_l = i-1;
                best_l = i-1;
                ii c = {i-2, i-1};
                cur = g(c);
                while (!li.empty()){
                    ii u = li.back();
                    ld val = g(u);
                    if (val >= cur){
                        c.fi = u.fi;
                        cur = g(c);
                        // merge u and c since u.se, or c.fi can be ignored
                    }else{
                        break;
                    }
                    li.pop_back();
                }
                li.pb(c);
            }
            best_l = r_l;
            cur = s[i];
            while (!li.empty()){
                if (cur < g(li.front())){
                    ii u = li.front();
                    best_l = u.fi;
                    break;
                }
                li.pop_front();
                // remove values from the left since they are now smaller than cur
            }
            ll val1 = dp[best_l][(j+1)%2] + (s[i]-s[best_l])*(s[n]-s[i]+s[best_l]);
            if (val1 > dp[i][j%2]){
                dp[i][j%2]=val1;
                //cout<<i<<" "<<j<<" "<<val<<" "<<val-dp[l][j-1]<<'\n';
                pos[i][j]=best_l;
            }
            /*
            FOR(l,1,i){
                //O(K*N^2) brute force
                ll val = dp[l][j-1] + (s[i]-s[l])*(s[n]-s[i]+s[l]);
                if (val > dp[i][j] && dp[l][j-1] != -1){
                    dp[i][j]=val;
                    //cout<<i<<" "<<j<<" "<<val<<" "<<val-dp[l][j-1]<<'\n';
                    pos[i][j]=l;
                }
                
            }
            */
        }
    }
    cout << dp[n][k%2]/2 << '\n';
    int cur = n;
    vi b;
    ROF(i,k,1){
        cur = pos[cur][i];
        b.pb(cur);
    }
    ROF(i,k-1,0){
        cout << b[i]<<" ";
    }
    cout << '\n';
    return;
}
 
int main() {
    ios::sync_with_stdio(false);
    int TC = 1;
    //cin >> TC;
    FOR(i, 1, TC+1){
        //cout << "Case #" << i << ": ";
        solve();
    }
}

Compilation message

sequence.cpp:10: warning: ignoring '#pragma GCC optimization' [-Wunknown-pragmas]
   10 | #pragma GCC optimization ("Ofast")
      | 
sequence.cpp:11: warning: ignoring '#pragma GCC optimization' [-Wunknown-pragmas]
   11 | #pragma GCC optimization ("unroll-loops")
      |

Subtask #1
11.0 / 11.0

# Verdict Execution time Memory Grader output
1 Correct 1 ms 1880 KB contestant found the optimal answer: 108 == 108
2 Correct 1 ms 1876 KB contestant found the optimal answer: 999 == 999
3 Correct 1 ms 1904 KB contestant found the optimal answer: 0 == 0
4 Correct 1 ms 1896 KB contestant found the optimal answer: 1542524 == 1542524
5 Correct 1 ms 1876 KB contestant found the optimal answer: 4500000000 == 4500000000
6 Correct 1 ms 1880 KB contestant found the optimal answer: 1 == 1
7 Correct 1 ms 1880 KB contestant found the optimal answer: 1 == 1
8 Correct 1 ms 1876 KB contestant found the optimal answer: 1 == 1
9 Correct 1 ms 1876 KB contestant found the optimal answer: 100400096 == 100400096
10 Correct 1 ms 1880 KB contestant found the optimal answer: 900320000 == 900320000
11 Correct 1 ms 1880 KB contestant found the optimal answer: 3698080248 == 3698080248
12 Correct 1 ms 1880 KB contestant found the optimal answer: 3200320000 == 3200320000
13 Correct 1 ms 1880 KB contestant found the optimal answer: 140072 == 140072
14 Correct 1 ms 1880 KB contestant found the optimal answer: 376041456 == 376041456
15 Correct 1 ms 1880 KB contestant found the optimal answer: 805 == 805
16 Correct 1 ms 1880 KB contestant found the optimal answer: 900189994 == 900189994
17 Correct 1 ms 1880 KB contestant found the optimal answer: 999919994 == 999919994

Subtask #2
11.0 / 11.0

# Verdict Execution time Memory Grader output
1 Correct 1 ms 1880 KB contestant found the optimal answer: 1093956 == 1093956
2 Correct 1 ms 1876 KB contestant found the optimal answer: 302460000 == 302460000
3 Correct 1 ms 1884 KB contestant found the optimal answer: 122453454361 == 122453454361
4 Correct 1 ms 1876 KB contestant found the optimal answer: 93663683509 == 93663683509
5 Correct 1 ms 1876 KB contestant found the optimal answer: 1005304678 == 1005304678
6 Correct 1 ms 1876 KB contestant found the optimal answer: 933702 == 933702
7 Correct 1 ms 1904 KB contestant found the optimal answer: 25082842857 == 25082842857
8 Correct 1 ms 1876 KB contestant found the optimal answer: 687136 == 687136
9 Correct 1 ms 1876 KB contestant found the optimal answer: 27295930079 == 27295930079
10 Correct 1 ms 1876 KB contestant found the optimal answer: 29000419931 == 29000419931

Subtask #3
11.0 / 11.0

# Verdict Execution time Memory Grader output
1 Correct 1 ms 1876 KB contestant found the optimal answer: 610590000 == 610590000
2 Correct 1 ms 1904 KB contestant found the optimal answer: 311760000 == 311760000
3 Correct 3 ms 2004 KB contestant found the optimal answer: 1989216017013 == 1989216017013
4 Correct 1 ms 1900 KB contestant found the optimal answer: 1499437552673 == 1499437552673
5 Correct 3 ms 2004 KB contestant found the optimal answer: 1019625819 == 1019625819
6 Correct 3 ms 2028 KB contestant found the optimal answer: 107630884 == 107630884
7 Correct 3 ms 2004 KB contestant found the optimal answer: 475357671774 == 475357671774
8 Correct 2 ms 1876 KB contestant found the optimal answer: 193556962 == 193556962
9 Correct 1 ms 1876 KB contestant found the optimal answer: 482389919803 == 482389919803
10 Correct 2 ms 1876 KB contestant found the optimal answer: 490686959791 == 490686959791

Subtask #4
17.0 / 17.0

# Verdict Execution time Memory Grader output
1 Correct 2 ms 1904 KB contestant found the optimal answer: 21503404 == 21503404
2 Correct 1 ms 1876 KB contestant found the optimal answer: 140412195 == 140412195
3 Correct 14 ms 2672 KB contestant found the optimal answer: 49729674225461 == 49729674225461
4 Correct 2 ms 1876 KB contestant found the optimal answer: 37485571387523 == 37485571387523
5 Correct 15 ms 2756 KB contestant found the optimal answer: 679388326 == 679388326
6 Correct 12 ms 2644 KB contestant found the optimal answer: 4699030287 == 4699030287
7 Correct 12 ms 2644 KB contestant found the optimal answer: 12418819758185 == 12418819758185
8 Correct 11 ms 2644 KB contestant found the optimal answer: 31093317350 == 31093317350
9 Correct 4 ms 2004 KB contestant found the optimal answer: 12194625429236 == 12194625429236
10 Correct 6 ms 2260 KB contestant found the optimal answer: 12345131038664 == 12345131038664

Subtask #5
21.0 / 21.0

# Verdict Execution time Memory Grader output
1 Correct 4 ms 2552 KB contestant found the optimal answer: 1818678304 == 1818678304
2 Correct 4 ms 2644 KB contestant found the optimal answer: 1326260195 == 1326260195
3 Correct 110 ms 10344 KB contestant found the optimal answer: 4973126687469639 == 4973126687469639
4 Correct 4 ms 2544 KB contestant found the optimal answer: 3748491676694116 == 3748491676694116
5 Correct 76 ms 7252 KB contestant found the optimal answer: 1085432199 == 1085432199
6 Correct 88 ms 7920 KB contestant found the optimal answer: 514790755404 == 514790755404
7 Correct 92 ms 8380 KB contestant found the optimal answer: 1256105310476641 == 1256105310476641
8 Correct 67 ms 7252 KB contestant found the optimal answer: 3099592898816 == 3099592898816
9 Correct 85 ms 7968 KB contestant found the optimal answer: 1241131419367412 == 1241131419367412
10 Correct 99 ms 9428 KB contestant found the optimal answer: 1243084101967798 == 1243084101967798

Subtask #6
29.0 / 29.0

# Verdict Execution time Memory Grader output
1 Correct 27 ms 8444 KB contestant found the optimal answer: 19795776960 == 19795776960
2 Correct 29 ms 8916 KB contestant found the optimal answer: 19874432173 == 19874432173
3 Correct 1014 ms 86540 KB contestant found the optimal answer: 497313449256899208 == 497313449256899208
4 Correct 28 ms 9432 KB contestant found the optimal answer: 374850090734572421 == 374850090734572421
5 Correct 1306 ms 86240 KB contestant found the optimal answer: 36183271951 == 36183271951
6 Correct 889 ms 62412 KB contestant found the optimal answer: 51629847150471 == 51629847150471
7 Correct 976 ms 66992 KB contestant found the optimal answer: 124074747024496432 == 124074747024496432
8 Correct 756 ms 55876 KB contestant found the optimal answer: 309959349080800 == 309959349080800
9 Correct 856 ms 62696 KB contestant found the optimal answer: 124113525649823701 == 124113525649823701
10 Correct 1016 ms 78540 KB contestant found the optimal answer: 124309619349406845 == 124309619349406845