容易将题意转化成，每次选最大的 $a$ 个加 $1$，然后再将最大的 $b$ 个减 $1$，重复操作 $d$ 次。

显然只有最开始的 $a$ 个是会被操作的，其他的贡献就是初始值的平方。只保留最大的 $a$ 个，每次操作变成选最小的 $b-a$ 个加 $1$。

然后就变成了下面这个问题：

$n$ 个数，每次给最小的 $m$ 个数加 $1$，求操作 $k$ 次后数列。

$m<n\le 10^5,k\le 10^9$。

要总共给 $n$ 个数加 $mk$，并且每个数最多加 $k$。容易发现，如果存在 $x<y$ 且 $x$ 加的数没有到 $k$，那么 $y$ 不可能被加任何数，否则不满足每次给最小的加的条件。

考虑二分答案，判断是否所有小于 $mid$ 的数都变成 $a_i+k$ 和 $mid$ 的最小值的代价是否 $\le mk$。二分要求出最大的满足条件的 $mid$。最后二分出来花了 $sum$ 的代价，给 $mk-sum$ 个等于 $mid$ 的数加 $1$ 就行了。

::::info[code]

#include <bits/stdc++.h>
using namespace std;

namespace z {

#define int long long
const int N = 5e5 + 5, mod = 1e9 + 7;
int n, m, x, y, a[N], b[N];
void main() {

    ios::sync_with_stdio(false);
    cin.tie(nullptr);cout.tie(nullptr);
    cin >> n >> m >> x >> y;
    int ans = 0, k = x - y;
    for(int i = 1; i <= n; i++) cin >> a[i];
    sort(a + 1, a + n + 1, greater<int>());
    auto chk = [&](int t) -> bool {
        int sum = 0;
        for(int i = 1; i <= x; i++) {
            sum += max(0ll, min(t - a[i], m));
        }
        return sum <= m * k;
    };
    int l = 0, r = 2e9, t = -1;
    while(l <= r) {
        int mid = l + r >> 1;
        if(chk(mid)) l = mid + 1, t = mid;
        else r = mid - 1;
    }  
    int sum = 0;
    for(int i = 1; i <= x; i++) {
        if(a[i] < t) sum += min(m, t - a[i]), a[i] = min(t, a[i] + m); 
    }
    int rem = m * k - sum;
    for(int i = 1; i <= x && rem; i++) if(a[i] == t) a[i]++, rem--;
    for(int i = 1; i <= n; i++) ans += a[i] * a[i] % mod;
    cout << ans % mod << '\n';

}

#undef int

}


int main()
{
    z::main();
    return 0;
}

::::