#pragma GCC target("avx") #pragma GCC optimize("O3") #pragma GCC optimize("unroll-loops") #include // #include // #include // #include // using namespace __gnu_pbds; // #include // namespace multiprecisioninteger = boost::multiprecision; // using cint=multiprecisioninteger::cpp_int; using namespace std; using ll=long long; using datas=pair; using ddatas=pair; using tdata=pair; using vec=vector; using mat=vector; using pvec=vector; using pmat=vector; // using llset=tree,rb_tree_tag,tree_order_statistics_node_update>; #define For(i,a,b) for(i=a;i<(ll)b;++i) #define bFor(i,b,a) for(i=b,--i;i>=(ll)a;--i) #define rep(i,N) For(i,0,N) #define rep1(i,N) For(i,1,N) #define brep(i,N) bFor(i,N,0) #define brep1(i,N) bFor(i,N,1) #define all(v) (v).begin(),(v).end() #define allr(v) (v).rbegin(),(v).rend() #define vsort(v) sort(all(v)) #define vrsort(v) sort(allr(v)) #define uniq(v) vsort(v);(v).erase(unique(all(v)),(v).end()) #define endl "\n" #define popcount __builtin_popcountll #define eb emplace_back #define print(x) cout< ostream& operator<<(ostream& os,const pair& p){return os<<"("< ostream& operator<<(ostream& os,const vector& v){ os<<"{";bool f=false; for(auto& x:v){if(f)os<<",";os< ostream& operator<<(ostream& os,const set& v){ os<<"{";bool f=false; for(auto& x:v){if(f)os<<",";os< ostream& operator<<(ostream& os,const multiset& v){ os<<"{";bool f=false; for(auto& x:v){if(f)os<<",";os< ostream& operator<<(ostream& os,const map& v){ os<<"{";bool f=false; for(auto& x:v){if(f)os<<",";os< inline bool chmax(T& a,const T b){bool x=a inline bool chmin(T& a,const T b){bool x=a>b;if(x)a=b;return x;} #ifdef DEBUG void debugg(){cout<void debugg(const T& x,const Args&... args){cout<<" "<size;--i)modncrlistm[i-1]=modncrlistm[i]*i%mod; } return modncrlistp[n]*modncrlistm[r]%mod*modncrlistm[n-r]%mod; } ll modpow(ll a,ll n,const ll m=mod){ if(n<0)return 0; ll res=1; while(n>0){ if(n&1)res=res*a%m; a=a*a%m; n>>=1; } return res; } constexpr ll gcd(const ll a,const ll b) noexcept{return (!b)?abs(a):(a%b==0)?abs(b):gcd(b,a%b);} constexpr ll lcm(const ll a,const ll b) noexcept{return a/gcd(a,b)*b;} class modmatrix{ mat a; ll H,W; public: modmatrix(mat& g):a(g){ H=g.size(); W=g[0].size(); } modmatrix(ll i,ll j):a(i,vec(j,0)){H=i;W=j;} modmatrix(ll n):a(n,vec(n,0)){H=W=n;} inline vec& operator [](int k){ return a[k]; } auto begin(){return a.begin();} auto end(){return a.end();} modmatrix operator +=(modmatrix b){ ll i,j; rep(i,this->H)rep(j,this->W){ (*this)[i][j]+=b[i][j]; if((*this)[i][j]>=mod)(*this)[i][j]-=mod; } return (*this); } modmatrix operator -=(modmatrix b){ ll i,j; rep(i,this->H)rep(j,this->W){ (*this)[i][j]-=b[i][j]; if((*this)[i][j]<0)(*this)[i][j]+=mod; } return (*this); } modmatrix operator *=(modmatrix b){ ll i,j,k; assert(this->W==b.H); modmatrix c(this->H,b.W); rep(i,this->H){ rep(k,this->W)rep(j,b.W){ c[i][j]+=(*this)[i][k]*b[k][j]%mod; } rep(j,b.W)c[i][j]%=mod; } (*this)=c; return (*this); } modmatrix operator ^=(ll K){ assert(this->H==this->W); modmatrix c(this->H); ll i; rep(i,this->H)c[i][i]=1; if(K&1)c*=(*this); while(K){ K>>=1; (*this)*=(*this); if(K&1)c*=(*this); } this->a.swap(c.a); return (*this); } modmatrix operator +(modmatrix c){ return modmatrix(*this)+=c; } modmatrix operator -(modmatrix c){ return modmatrix(*this)-=c; } modmatrix operator *(modmatrix c){ return modmatrix(*this)*=c; } modmatrix operator ^(ll K){ return modmatrix(*this)^=K; } modmatrix del(ll eh,ll ew){ ll i,j; mat res; rep(i,H){ if(i==eh)continue; res.resize(res.size()+1); rep(j,W){ if(j==ew)continue; res.back().eb(a[i][j]); } } return res; } ll determinant(){ assert(H==W); ll i,j,k,ans=1; auto b(a); rep(i,H){ if(!b[i][i]){ For(j,i+1,H)if(b[j][i]){ swap(b[i],b[j]); ans*=-1; break; } if(j==H)return 0; } (ans*=b[i][i])%=mod; For(j,i+1,W)(b[i][j]*=modinv(b[i][i]))%=mod; For(j,i+1,H)if(b[j][i]){ ll x=mod-b[j][i]; b[j][i]=0; For(k,i+1,W){ b[j][k]+=x*b[i][k]%mod; if(b[j][k]>=mod)b[j][k]-=mod; } } } if(ans<0)ans+=mod; return ans; } modmatrix inv(){ assert(H==W); ll i,j,k; modmatrix b(a),c(H,W); rep(i,H)c[i][i]=1; rep(i,H){ if(!b[i][i]){ For(j,i+1,H)if(b[j][i]){ swap(b[i],b[j]); swap(c[i],c[j]); break; } if(j==H)assert(false); } ll x=modinv(b[i][i]); rep(j,W){ (b[i][j]*=x)%=mod; (c[i][j]*=x)%=mod; } rep(j,H){ if(i==j)continue; if(b[j][i]){ x=mod-b[j][i]; rep(k,W){ b[j][k]+=x*b[i][k]%mod; if(b[j][k]>=mod)b[j][k]-=mod; c[j][k]+=x*c[i][k]%mod; if(c[j][k]>=mod)c[j][k]-=mod; } } } } return c; } void out(){ for(auto x:a)output(x); } }; ll N,M,K,H,W,L; int main(){ startupcpp(); // int codeforces;cin>>codeforces;while(codeforces--){ ll i,j; cin>>N>>M; mat dp(2,vec(M+1,-inf)),nx(dp); vec a(N); rep(i,N)cin>>a[i]; vrsort(a); dp[0][0]=0; rep(i,N-1){ ll x=a[i]-a[i+1]; rep(j,M+1){ nx[0][j]=max(dp[0][j],dp[1][j]); if(j+x<=M)nx[1][j+x]=dp[0][j]+a[i]; } swap(dp,nx); } ll ans=-inf; rep(i,2)rep(j,M+1)chmax(ans,dp[i][j]); print(ans); // } }