#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(2e5)+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];
int last[N][205];
struct line {
ll m, b;
int id;
};
vector<line> vec;
ll val(ll x, line a) {
return x * a.m + a.b;
}
bool bad(line a, line b, line c){
return (b.b - a.b) * (a.m - c.m) >= (c.b - a.b) * (a.m - b.m);
}
void add(line a) {
while(vec.sz >= 2 && bad(vec[vec.sz - 2], vec.back(), a)) {
vec.pop_back();
}
vec.pb(a);
}
pair<ll, int> get(ll x) {
if(vec.size() == 1) return {val(x, vec[0]), vec[0].id};
int l = 0, r = vec.sz - 2;
while(l <= r) {
int mid = (l + r) / 2;
if(val(x, vec[mid]) <= val(x, vec[mid + 1])) l = mid + 1;
else r = mid - 1;
}
return {val(x, vec[l]), 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 ++) dp[i] = -inf;
dp[0] = 0;
for(int kk = 1; kk <= k; kk ++) {
add({pr[kk - 1], dp[kk - 1], kk - 1});
for(int i = kk; i <= n; i ++) {
// need to find j
auto got = get(pr[i]);
if(dp[i] != -inf) add({pr[i], dp[i], i});
if(got.ff - pr[i] * pr[i] > dp[i]) {
dp[i] = got.ff - pr[i] * pr[i];
last[i][kk] = got.ss;
}
// 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;
// }
// }
}
vec.clear();
}
ll pos = 0;
for(int i = k; i <= n; i ++) {
if(dp[i] + pr[n] * pr[i] > dp[pos] + pr[n] * pr[pos]) {
pos = i;
}
}
cout<<dp[pos] + 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<<" ";
cout<<"\n";
}
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(); }
}
