#include<cstdio>
#include<algorithm>
#include<vector>
#define mid ((l+r)>>1)
#define ll long long
using namespace std;
const int N=5e5+10;
const ll inf=3e17;
int n,k;
ll s[N],L,a[N],ans;
ll calc(int l,int r)
{
	if(l>=r) return 0;
	return a[mid]*(mid-l)-(s[mid-1]-s[l-1])+(s[r-1]-s[mid])-a[mid]*(r-mid-1);
}
namespace subtask1
{
	ll dp[N];
	int g[N];
	ll calc2(int i,int j){return dp[i]+2*calc(i,j);}
	int head,tail,lst[N];
	struct node{
		int l,r,k;
		ll res;
	}t[N];
	node create(int l,int r,int k){return node{l,r,k,calc2(k,l)};}
	void work(ll k)
	{
		int hd,tl;
		dp[1]=g[1]=0;
		t[hd=tl=1]=create(2,n+1,1);
		for(int i=2;;i++)
		{
			int w=t[hd].k;
			dp[i]=calc2(w,i)+k;
			g[i]=g[w]+1,lst[i]=w;
			if(i==n+1) break;
			if(t[hd].l==t[hd].r) hd++;
			else t[hd].l++,t[hd].res=calc2(t[hd].k,i+1);
			while(hd<=tl&&calc2(i,t[tl].l)<=t[tl].res) tl--;
			if(tl<hd) t[++tl]=create(i+1,n+1,i);
			else
			{
				int l=t[tl].l,r=t[tl].r+1,w=t[tl].k;
				while(l<r)
				{
					if(calc2(w,mid)>=calc2(i,mid)) r=mid;
					else l=mid+1;
				}
				t[tl].r=l-1;
				t[++tl]=create(l,n+1,i);
			}
		}
	}
	vector<int> work2()
	{
		int w=n+1;
		vector<int> v;
		while(w)
		{
			v.push_back(w);
			w=lst[w];
		}
		reverse(v.begin(),v.end());
		return v;
	}
	vector<int> solve()
	{
		ll l=-1,r=inf;
		while(l<r)
		{
			work(mid*2+1);
			if(g[n+1]>k) l=mid+1;
			else r=mid;
		}
		l=l*2+1;
		work(l);
		vector<int> v1=work2();
		if(g[n+1]==k) return v1;
		work(l-2);
		vector<int> v2=work2();
		if(g[n+1]==k) return v2;
		int m=v2.size();
		for(int i=1,j=0;i<m;i++)
		{
			while(v1[j]<v2[i]) j++;
			if(v2[i-1]>=v1[j-1])
			{
				if(m-i+j==k+1)
				{
					vector<int> v;
					for(int w=0;w<j;w++) v.push_back(v1[w]);
					for(int w=i;w<m;w++) v.push_back(v2[w]);
					return v;
				}
			}
		}
		return vector<int>();
	}
}
vector<int> ANS;
namespace subtask2
{
	ll f[N][2];
	int lst[N][2];
	void work(int l,int r,int L,int R,int p)
	{
		if(l>r) return;
		ll s=inf;
		for(int i=L;i<=R&&i<=mid;i++)
		{
			ll x=f[i][p^1]+calc(i,mid);
			if(x<s) s=x,lst[mid][p]=i;
		}
		f[mid][p]=s;
		work(l,mid-1,L,lst[mid][p],p);
		work(mid+1,r,lst[mid][p],R,p);
	}
	void solve(int l,int r,vector<pair<int,int>> v)
	{
		if(l>r) return;
		int m=(int)v.size();
		for(int i=v[0].first;i<=v[0].second;i++) f[i][0]=calc(mid,i);
		for(int i=1;i<m;i++)
			work(v[i].first,v[i].second,v[i-1].first,v[i-1].second,i&1);
		auto k=v.back();
		int w=k.first,p=(m-1)&1;
		vector<int> s;
		for(int i=k.first;i<=k.second;i++) f[i][p]+=calc(i,mid+n);
		for(int i=k.first+1;i<=k.second;i++) if(f[i][p]<f[w][p]) w=i;
		int flag=0;
		if(f[w][p]<ans) ans=f[w][p],flag=1;
		for(int i=m-1;i>=0;i--) s.push_back(w),w=lst[w][i&1];
		reverse(s.begin(),s.end());
		if(flag)
		{
			s.push_back(mid);
			ANS=s;
			s.pop_back();
		}
		if(l==r) return;
		vector<pair<int,int>> v1(m),v2(m);
		for(int i=0;i<m;i++)
		{
			v1[i]=make_pair(v[i].first,s[i]);
			v2[i]=make_pair(s[i],v[i].second);
		}
		solve(l,mid-1,v1);
		solve(mid+1,r,v2);
	}
	void solve(vector<int> v)
	{
		for(int i=0;i<=k;i++) v.push_back(v[i]+n);
		vector<pair<int,int>> V(k-1);
		int w=0;
		for(int i=0;i<k;i++)
			if(v[i+1]-v[i]<=n/k)
			{
				w=i;
				break;
			}
		for(int i=0;i<k-1;i++)
			V[i]=make_pair(v[w+i+1],v[w+i+2]);
		solve(v[w],v[w+1],V);
	}
}
void work(vector<int> &v)
{
	for(int i=0;i<(int)v.size();i++)
		if(v[i]>n) v[i]-=n;
	sort(v.begin(),v.end());
}
ll b[N];
int main()
{
	scanf("%d%d%lld",&n,&k,&L);
	for(int i=1;i<=n;i++) scanf("%lld",&a[i]);
	for(int i=n+1;i<=n*2;i++) a[i]=a[i-n]+L;
	for(int i=1;i<=2*n;i++) s[i]=s[i-1]+a[i];
	if(k==1)
	{
		ll ans=inf;
		int w=0;
		for(int i=1;i<=n;i++)
		{
			ll x=calc(i,i+n);
			if(x<ans) ans=x,w=i;
		}
		printf("%lld\n%d\n",ans,w);
		return 0;
	}
	vector<int> v=subtask1::solve();
	ans=0;
	for(int i=1;i<(int)v.size();i++) ans+=calc(v[i-1],v[i]);
	v.pop_back();
	ANS=v;
	subtask2::solve(v);
	work(ANS);
	printf("%lld\n",ans);
	ANS.push_back(ANS[0]+n);
	for(int i=0;i<k;i++)
	{
		int x=(ANS[i]+ANS[i+1])>>1;
		if(x>n) x-=n;
		b[i]=a[x];
	}
	sort(b,b+k);
	for(int i=0;i<k;i++) printf("%lld ",b[i]);
	printf("\n");
}