Problem:
Score : 700 points

Problem Statement

N problems are proposed for an upcoming contest. Problem i has an initial integer score of A_i points.

M judges are about to vote for problems they like. Each judge will choose exactly V problems, independently from the other judges,
and increase the score of each chosen problem by 1.

After all M judges cast their vote, the problems will be sorted in non-increasing order of score, and the first P problems will be chosen for the problemset.
Problems with the same score can be ordered arbitrarily, this order is decided by the chief judge.

How many problems out of the given N have a chance to be chosen for the problemset?

Constraints

2 \le N \le 10^5

1 \le M \le 10^9

1 \le V \le N - 1

1 \le P \le N - 1

0 \le A_i \le 10^9

Input

Input is given from Standard Input in the following format:

N M V P
A_1 A_2 ... A_N

Output

Print the number of problems that have a chance to be chosen for the problemset.

Sample Input 1

6 1 2 2
2 1 1 3 0 2

Sample Output 1

5

If the only judge votes for problems 2 and 5, the scores will be 2 2 1 3 1 2.
The problemset will consist of problem 4 and one of problems 1, 2, or 6.

If the only judge votes for problems 3 and 4, the scores will be 2 1 2 4 0 2.
The problemset will consist of problem 4 and one of problems 1, 3, or 6.

Thus, problems 1, 2, 3, 4, and 6 have a chance to be chosen for the problemset. On the contrary, there is no way for problem 5 to be chosen.

Sample Input 2

6 1 5 2
2 1 1 3 0 2

Sample Output 2

3

Only problems 1, 4, and 6 have a chance to be chosen.

Sample Input 3

10 4 8 5
7 2 3 6 1 6 5 4 6 5

Sample Output 3

8

Solution in PyPy3 (2.4.0) (judge verdict: Accepted)

from heapq import heappush, heappop, heapify
from collections import deque, defaultdict, Counter
import itertools
from itertools import permutations, combinations
import sys
import bisect
import string
import math
import time
from fractions import gcd
#import random


def I(): return int(input())


def MI(): return map(int, input().split())


def LI(): return [int(i) for i in input().split()]


def LI_(): return [int(i)-1 for i in input().split()]


def StoI(): return [ord(i)-97 for i in input()]


def ItoS(nn): return chr(nn+97)


def input(): return sys.stdin.readline().rstrip()


def show(*inp, end='\n'):
    if show_flg:
        print(*inp, end=end)


YN = ['Yes', 'No']
MOD = 10**9+7
inf = float('inf')
IINF = 10**19
l_alp = string.ascii_lowercase
u_alp = string.ascii_uppercase
ts = time.time()
sys.setrecursionlimit(10**6)


# show_flg = False
show_flg = True


def main():
    N, M, V, P = MI()
    A = MI()
    A = sorted(A, reverse=True)
    minp = A[P-1]
    result = 0
    td = 0
    count = 0
    tpcount = 1
    tptd = 0
    pre = minp
    for i in range(P, N):
        if pre != A[i]:
            td += tptd
            count += tpcount
            tptd = 0
            tpcount = 0
        if M * min(N-V, count) >= (minp - A[i]) * count - td and minp <= A[i] + M:
            result += 1
            tpcount += 1
            tptd += minp - A[i]
            pre = A[i]
        else:
            break
    print(result + P)
    # N, M, V, P = 10**5, 10*19, 10**4, 10**4
    # A = [i for i in range(N)]
    # MAP = {}
    # A = sorted(A, reverse=True)
    # AH = A[:P]
    # AH_min = AH[-1]
    # AT = []
    # cur = 0
    # sent_c = 0
    # for i in range(P, N):
    #     if AH_min == A[i]:
    #         sent_c += 1
    #     # else:
    #     if A[i] not in MAP:
    #         MAP[A[i]] = cur
    #     cur += 1
    #     AT.append(A[i])
    # print(AH, AT)
    # # print('M=', M, 'V=', V)
    # ans = len(AH)
    # AT_len = len(AT)
    # # print(AT, MAP)
    # # show(len(AT), V - len(AH) + 1)

    # _V = V - len(AH) + 1
    # # print('_V', _V, V, AH_min)
    # for i in range(AT_len):
    #     if AH_min == AT[i]:
    #         ans += 1
    #     elif AH_min <= (AT[i] + M):
    #         remaining = AT_len - (MAP[AT[i]] + 1)
    #         print('AT[i]', AT[i], 'ans=', ans, 'i=', i, 'remain=', remaining +
    #               1, 'V=', _V, 'AH_MIN=', AH_min, 'MAP=', MAP[AT[i]])
    #         if remaining + 1 >= _V:
    #             ans += 1

    # print(ans)


if __name__ == '__main__':
    main()