// Arayi
#include <bits/stdc++.h>
#define fr first
#define sc second
#define MP make_pair
#define ad push_back
#define PB push_back
#define fastio \
ios_base::sync_with_stdio(false); \
cin.tie(0); \
cout.tie(0);
#define lli long long int
#define y1 arayikhalatyan
#define j1 jigglypuff
#define ld long double
#define itn int
#define pir pair<int, int>
#define all(x) (x).begin(), (x).end()
#define str string
#define enl endl
#define en endl
#define cd complex<long double>
#define vcd vector<cd>
#define vii vector<int>
#define vlli vector<lli>
using namespace std;
lli gcd(lli a, lli b) { return (b == 0LL ? a : gcd(b, a % b)); }
ld dist(ld x, ld y1, ld x2, ld y2)
{
return sqrt((x - x2) * (x - x2) + (y1 - y2) * (y1 - y2));
}
lli S(lli a)
{
return (a * (a + 1LL)) / 2;
}
mt19937 rnd(363542);
char vow[] = {'a', 'e', 'i', 'o', 'u'};
int dx[] = {0, -1, 0, 1, -1, -1, 1, 1, 0};
int dy[] = {-1, 0, 1, 0, -1, 1, -1, 1, 0};
char dc[] = {'R', 'D', 'L', 'U'};
const int N = 1e5;
const lli mod = 1e9 + 7;
const ld pi = acos(-1);
const int T = 200;
lli bp(lli a, lli b = mod - 2LL)
{
lli ret = 1;
while (b)
{
if (b & 1)
ret *= a, ret %= mod;
a *= a;
a %= mod;
b >>= 1;
}
return ret;
}
ostream &operator<<(ostream &c, pir a)
{
c << a.fr << " " << a.sc;
return c;
}
template <class T>
void maxi(T &a, T b)
{
a = max(a, b);
}
template <class T>
void mini(T &a, T b)
{
a = min(a, b);
}
int n, k;
int p[N][201];
lli nx[N], nw[N];
lli a[N], pr[N];
ld f(int i, int i1)
{
ld x1;
ld a = pr[i1] - pr[i];
if (a != 0)
x1 = ((ld)(nx[i1] - nx[i])) / a;
else if (nx[i1] > nx[i])
x1 = 1e10;
else
x1 = -1;
return x1;
}
bool col[N];
int main()
{
fastio;
cin >> n >> k;
for (int i = 0; i < n; i++)
{
cin >> a[i];
pr[i] = a[i];
if (i)
pr[i] += pr[i - 1];
}
for (int i = 0; i < n; i++)
nx[i] = pr[i] * (pr[n - 1] - pr[i]); //, cout << dp[i][1] << " ";
// cout << endl;
int ii = 0;
if (k == 1)
{
for (int i = 0; i < n; i++)
if (nx[i] > nx[ii])
ii = i;
cout << nx[ii] << endl;
cout << ii + 1 << endl;
return 0;
}
for (int x = 2; x <= k; x++)
{
vector<pair<ld, int>> fp;
fp.ad(MP(1e10, 0));
for (int i = 1; i < n; i++)
{
while ((int)fp.size() > 1 && f(fp.back().sc, i) > fp.back().fr)
fp.pop_back();
fp.ad(MP(f(fp.back().sc, i), i));
}
int i1 = 0;
for (int i = 0; i < n - 1; i++)
{
nw[i] = pr[i] * (pr[n - 1] - pr[i]);
ld x1 = pr[n - 1] - pr[i];
while (i1 + 1 < (int)fp.size() && x1 <= fp[i1 + 1].fr && fp[i1 + 1].sc < i)
i1++;
// cout << fp[i1].sc << "-";
p[i][x] = fp[i1].sc;
nw[i] += nx[fp[i1].sc];
nw[i] -= pr[fp[i1].sc] * x1;
// cout << dp[i][x] << " ";
if (x == k && nw[i] >= nw[ii])
ii = i;
}
swap(nx, nw);
}
cout << nx[ii] << endl;
for (int i = 0; i < k; i++)
{
col[ii] = 1;
cout << ii + 1 << " ";
ii = p[ii][k - i];
while (ii < n - 2 && a[ii + 1] == 0 && col[ii + 1] == 0)
ii++;
}
return 0;
}
/*
__
*(><)*
\/ /
||/
--||
||
/\
/ \
/ \
*/
