#include <bits/stdc++.h>
#define int long long
using namespace std;
int n, m;
int s,d;
int xx,yy;
int exgcd(int a, int b, int &x, int &y)// 拓欧
{
    if (b == 0)
    {
        x = 1;
        y = 0;
        return a;
    }
    int d = exgcd(b, a % b, y, x);
    y -= (a / b) * x;
    return d;
}
int e;
void extend_gcd(int a, int b, int& d, int& x, int& y)
{
    if(!b){ d = a; x = 1; y = 0; }
    else { extend_gcd(b, a%b,d, y, x); y -= x*(a/b);}
}


void solve () {
    cin>>n>>m;
    int sum=0;
    for(int i=1;i<=n;i++){
        int x;
        cin>>x;
        sum+=x;
    }
    int a=n,b=(n+1)*n/2;
    int k=gcd(a,b);
    exgcd(a,b,s,d);
    s=(s%m+m)%m;
    d=(d%m+m)%m;
    // extend_gcd(a,b,e,s,d);
    // int b1=b/k;
    // s=(s+b1)*(e/k);
    // s=(s % b1 + b1) % b1;
    // d=(e - a*s) / b;
    int res=exgcd(k,m,xx,yy);
    xx=(((-sum%m)/res*xx)%m+m)%m;
    cout<<(((sum)%m+a*s*xx%m+b*d*xx%m)%m)<<endl;
    cout<<(s*xx%m)<<" "<<(d*xx%m)<<endl;
    // // work(a,b,(a*s+b*d)%m,s,d);
    // cerr<<s<< " "<<d<<endl;
    // int mi1=sum%m;
    // if(mi1==0){
    //     cout<<0<<endl<<0<<" "<<0<<endl;
    //     return ;
    // }
    // int mi2=m-sum%m;
    // if(mi2%k==0){
    //     cout<<0<<endl<<mi2/k*s<<" "<<mi2/k*d<<endl;
    //     return ;
    // }
    // else{
    //     mi2/=k;
    //     mi2=(sum+(mi2+1)*k)%m;
    //     if(mi1<=mi2){
    //         cout<<mi1<<endl<<0<<" "<<0<<endl;
    //     }
    //     else{
    //         cout<<mi2<<endl<< (mi2+1)*s<<" "<<(mi2+1)*d <<endl;
    //     }
    // }

}
signed main () {
    int T = 1; 
    std::ios::sync_with_stdio(false), cin.tie(0), cout.tie(0);
    //cin>>T;
    while (T --) solve ();
    return 0;
}