#include <bits/stdc++.h>
#include <ext/pb_ds/assoc_container.hpp>
#include <ext/pb_ds/tree_policy.hpp>
using namespace std;
using namespace __gnu_pbds;
typedef tree<int, null_type, less<int>, rb_tree_tag,tree_order_statistics_node_update> ordered_set;
#define scan(x) do{while((x=getchar())<'0'); for(x-='0'; '0'<=(_=getchar()); x=(x<<3)+(x<<1)+_-'0');}while(0)
char _;
#define complete_unique(a) a.erase(unique(a.begin(),a.end()),a.end())
#define all(a) a.begin(),a.end()
#define println printf("\n");
#define readln(x) getline(cin,x);
#define pb push_back
#define endl "\n"
#define INT_INF 0x3f3f3f3f
#define LL_INF 0x3f3f3f3f3f3f3f3f
#define MOD 1000000007
#define mp make_pair
#define fastio cin.tie(0); cin.sync_with_stdio(0);
#define MAXN 100005
#define MAXK 205
typedef unsigned long long ull;
typedef long long ll;
typedef long double ld;
typedef unordered_map<int,int> umii;
typedef pair<int,int> pii;
typedef pair<double,double> pdd;
typedef pair<ll,ll> pll;
typedef pair<int,pii> triple;
typedef int8_t byte;
mt19937 g1(time(0));
int randint(int a, int b){return uniform_int_distribution<int>(a, b)(g1);}
ll randlong(ll a,ll b){return uniform_int_distribution<long long>(a, b)(g1);}
ll gcd(ll a, ll b){return b == 0 ? a : gcd(b, a % b);}
ll lcm(ll a, ll b){return a*b/gcd(a,b);}
ll fpow(ll b, ll exp, ll mod){if(exp == 0) return 1;ll t = fpow(b,exp/2,mod);if(exp&1) return t*t%mod*b%mod;return t*t%mod;}
ll divmod(ll i, ll j, ll mod){i%=mod,j%=mod;return i*fpow(j,mod-2,mod)%mod;}
int len,K,arr[MAXN],par[MAXK][MAXN],dq[MAXN],l,r,lst;
ll dp[2][MAXN],psa[MAXN];
inline pll inter(int a, int b, int k){
ll m1 = psa[a], b1 = -(psa[a]*psa[a])+dp[!k][a];
ll m2 = psa[b], b2 = -(psa[b]*psa[b])+dp[!k][b];
return mp(b2-b1,m1-m2);
}
int main(){
scanf("%d %d",&len,&K); K++;
for(int i=1; i<=len; i++){
scanf(" %d",&arr[i]);
psa[i] = psa[i-1]+arr[i];
}
for(int k=2; k<=K; k++){
lst = !lst;
l = r = dq[0] = 0;
for(int i=1; i<=len; i++){
while(r-l+1 >= 2 && dp[!lst][dq[l+1]]+(psa[i]-psa[dq[l+1]])*psa[dq[l+1]] >= dp[!lst][dq[l]]+(psa[i]-psa[dq[l]])*psa[dq[l]]) ++l;
dp[lst][i] = dp[!lst][dq[l]]+(psa[i]-psa[dq[l]])*psa[dq[l]];
par[k][i] = dq[l];
while(r-l+1 >= 2){
pll f = inter(dq[r-1],i,lst), s = inter(dq[r-1],dq[r],lst);
ll fx = f.first*s.second;
ll sx = s.first*f.second;
if(fx <= sx) --r;
else break;
}
dq[++r] = i;
}
}
printf("%lld\n",dp[lst][len]);
vector<int> v;
int curr = len;
for(int i=K; i>=1; i--){
if(i < K) v.pb(curr);
curr = par[i][curr];
}
for(int i=(int)v.size()-1; i>=0; i--)
printf("%d ",v[i]);
}
/*
4 1
0 0 1 1
ans=1
*/
