Submission #894990

#	Time	Username	Problem	Language	Result	Execution time	Memory
894990		vjudge1	Split the sequence (APIO14_sequence)	C++17	50 / 100	2056 ms	131072 KiB

sequence

#include <bits/stdc++.h>
using namespace std;
#include <ext/pb_ds/assoc_container.hpp>
#include <ext/pb_ds/tree_policy.hpp>
using namespace __gnu_pbds;
  
#define ordered_set tree<int, null_type,less<int>, rb_tree_tag,tree_order_statistics_node_update>
 
// #pragma GCC optimize("Ofast,unroll-loops,fast-math")
// #pragma GCC target("popcnt")


typedef long long ll;
typedef unsigned long long ull;
typedef long double ld;
typedef pair<ll,ll> pll;
 
#define sz size()
#define ff first
#define ss second
#define all(a) a.begin(),a.end()
#define pb push_back

const int mod = ll(1e9)+7; //(b + (a%b)) % b (to mod -1%(10^9+7) correctly in c++ its -1 but its suppose to be 10^9+6
const ll MOD = 998244353;  // (a%mod)*(binpow(b,mod-2,mod) = (a/b)%mod
const int N = ll(1e5)+100;
const int M = ll(2e5) + 100;
const long long inf = 5e18;
const long double eps = 1e-15L;
 
ll lcm(ll a, ll b) { return (a / __gcd(a,b))*b; }

ll binpow(ll a, ll b, ll m) { ll res=1; a %= m; while(b>0){ if(b&1)res=(res * a) % m; a=(a * a) % m; b/=2; } return res%m;}

void Freopen(string Key){ freopen((Key+".in").c_str(), "r", stdin); freopen((Key+".out").c_str(), "w", stdout); }

#define int ll
int n, k, a[N];

ll pr[N], dp[N][205], last[N][205];

struct line {
	ld m, b;
	int id;
};

vector<line> vec;
vector<ld> posx;

ld interx(line a, line b) {
	return (ld)(b.b - a.b) / (a.m - b.m);
}

void add(line a) {
	while(vec.sz >= 2 && (vec.back().m == a.m || interx(vec[vec.sz - 2], vec.back()) > interx(vec[vec.sz - 2], a))) {
		vec.pop_back(), posx.pop_back();
	}
	if(vec.empty()) posx.pb(-inf);
	else posx.pb(interx(vec.back(), a));
	vec.pb(a);
}

pair<ll, int> get(ll x) {
	ll l = 0, r = posx.sz - 1;
	while(l <= r) {
		int mid = (l + r) / 2;
		if(posx[mid] <= x) l = mid + 1;
		else r = mid - 1;
	}
	l --;
	return {x * vec[l].m + vec[l].b, vec[l].id};
}

void Baizho() {
	cin>>n>>k;
	for(int i = 1; i <= n; i ++) {
		cin>>a[i];
		pr[i] = pr[i - 1] + a[i];
	}
	for(int i = 0; i <= n; i ++) for(int j = 0; j <= k; j ++) dp[i][j] = -inf;
	dp[0][0] = 0;
	add({0, 0, 0});
	for(int kk = 1; kk <= k; kk ++) {
		for(int i = kk; i <= n; i ++) {
			for(int j = 0; j < i; j ++) {
				if(dp[j][kk - 1] == -inf) continue;
//				+ pr[n] * pr[i] - pr[n] * pr[j]
				if(dp[j][kk - 1] - pr[i] * pr[i] + pr[i] * pr[j] > dp[i][kk]) {
					dp[i][kk] = dp[j][kk - 1] - pr[i] * pr[i] + pr[i] * pr[j];
					last[i][kk] = j;
				}
			}
		}
	}
	ll pos = 0;
	for(int i = k; i <= n; i ++) {
		if(dp[i][k] + pr[n] * pr[i] > dp[pos][k] + pr[n] * pr[pos]) {
			pos = i;
		}
	}
	cout<<dp[pos][k] + pr[n] * pr[pos]<<"\n";
	vector<int> ans;
	while(pos) {
		ans.pb(pos);
		pos = last[pos][k --];
	}
	reverse(all(ans));
	for(auto x : ans) cout<<x<<" ";
}
 
signed main() {		
// 	Freopen("nondec");
    ios_base::sync_with_stdio(false);   
    cin.tie(0);cout.tie(0); 
//   	precalc();
   	
    int ttt = 1;
//    cin>>ttt;

    for(int i=1; i<=ttt; i++) {Baizho(); }
}

Compilation message (stderr)

sequence.cpp: In function 'void Freopen(std::string)':
sequence.cpp:35:34: warning: ignoring return value of 'FILE* freopen(const char*, const char*, FILE*)' declared with attribute 'warn_unused_result' [-Wunused-result]
   35 | void Freopen(string Key){ freopen((Key+".in").c_str(), "r", stdin); freopen((Key+".out").c_str(), "w", stdout); }
      |                           ~~~~~~~^~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
sequence.cpp:35:76: warning: ignoring return value of 'FILE* freopen(const char*, const char*, FILE*)' declared with attribute 'warn_unused_result' [-Wunused-result]
   35 | void Freopen(string Key){ freopen((Key+".in").c_str(), "r", stdin); freopen((Key+".out").c_str(), "w", stdout); }
      |                                                                     ~~~~~~~^~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

• Subtask 1: 11 / 11
• Subtask 2: 11 / 11
• Subtask 3: 11 / 11
• Subtask 4: 17 / 17
• Subtask 5: 0 / 21
• Subtask 6: 0 / 29