While Grisha was celebrating New Year with Ded Moroz, Misha gifted Sasha a small rectangular pond of size n × m, divided into cells of size 1 × 1, inhabited by tiny evil fishes (no more than one fish per cell, otherwise they'll strife!).

The gift bundle also includes a square scoop of size r × r, designed for fishing. If the lower-left corner of the scoop-net is located at cell (x, y), all fishes inside the square (x, y)...(x + r - 1, y + r - 1) get caught. Note that the scoop-net should lie completely inside the pond when used.

Unfortunately, Sasha is not that skilled in fishing and hence throws the scoop randomly. In order to not frustrate Sasha, Misha decided to release k fishes into the empty pond in such a way that the expected value of the number of caught fishes is as high as possible. Help Misha! In other words, put k fishes in the pond into distinct cells in such a way that when the scoop-net is placed into a random position among (n - r + 1)·(m - r + 1) possible positions, the average number of caught fishes is as high as possible.


-----Input-----

The only line contains four integers n, m, r, k (1 ≤ n, m ≤ 10^5, 1 ≤ r ≤ min(n, m), 1 ≤ k ≤ min(n·m, 10^5)).


-----Output-----

Print a single number — the maximum possible expected number of caught fishes.

You answer is considered correct, is its absolute or relative error does not exceed 10^{ - 9}. Namely, let your answer be a, and the jury's answer be b. Your answer is considered correct, if $\frac{|a - b|}{\operatorname{max}(1,|b|)} \leq 10^{-9}$.


-----Examples-----
Input
3 3 2 3

Output
2.0000000000

Input
12 17 9 40

Output
32.8333333333



-----Note-----

In the first example you can put the fishes in cells (2, 1), (2, 2), (2, 3). In this case, for any of four possible positions of the scoop-net (highlighted with light green), the number of fishes inside is equal to two, and so is the expected value. [Image]

Difficulty: interview

Starter code:

Solutions:
import heapq as hq
from queue import PriorityQueue

import math
n,m,r, k= input().split()
N = int(n)
M = int(m)
R = int(r)
K = int(k)

q = PriorityQueue()



for i in range(1,math.floor((N+1)/2) + 1):
    maxi = min(min(i,N-i+1),min(R,N-R+1)) * min(min(R,M-R+1),math.ceil(M/2))
    num = M - (2 * min(min(R,M-R+1),math.ceil(M/2))-2)
    mult = 2
    if(i > math.floor(N/2)):
        mult = 1
    q.put((-maxi,num * mult,i))
    #print(str(maxi) + " " + str(num) + " " + str(mult))


ans = 0
while(K > 0):
    pop = q.get()
    #print(pop)
    a = -1 * pop[0]
    b = pop[1]
    c = pop[2]
    d = min(min(c,N-c+1),min(R,N-R+1))
    if(d != a):
       # if(q.)
       # if(q.get(-(a - d)) != )
        mult = 2
        if (c > N / 2):
            mult = 1
        q.put((-(a - d),2*mult,c))
    ans += a * min(b,K)
    K -= b;

tot = (N-R+1) * (M-R+1)
#print("ANS = " + str(ans))
#print("FINANS = " + str(ans/tot))
print(str(ans/tot))


'''

d = []
for i in range(0,N):
    d.append([])
    for j in range(0,M):
        d[i].append(0)

tot = 0
for i in range(0,N-R+1):
    for j in range(0,M-R+1):
        for k in range(i,i+R):
            for l in range(j,j+R):
                d[k][l] += 1
                tot += 1


print(N-R+1)*(M-R+1) * (R*R)
print(tot)
print()
for i in d:
    print(i)

'''


class TupleHeap(object):
    def __init__(self):
        self.heap = [0]
    def mult(self, i):
        return self.heap[i][0] * self.heap[i][1]
    def __str__(self):
        return " ".join(map(str, self.heap))
    def add(self, var):
        self.heap.append(var)
        if len(self.heap) > 2:
            # надо сортировать снизу
            el = len(self.heap) - 1
            par = el // 2
            while (el != 1) and (self.mult(el) > self.mult(par)):
                self.heap[el], self.heap[par] = self.heap[par], self.heap[el]
                el = par
                par = el // 2
    def pop(self):
        if len(self.heap) == 1:
            return None  # Нечего брать из кучи
        ans = self.heap[1]
        if len(self.heap) == 2:
            self.heap.pop()
            return ans
        self.heap[1] = self.heap.pop()
        # делаем сортировку вниз
        el = 1
        var_el = self.mult(el)
        child1 = el * 2
        if child1 < len(self.heap):
            var_child1 = self.mult(child1)
        else:
            var_child1 = -1
        child2 = el * 2 + 1
        if child2 < len(self.heap):
            var_child2 = self.mult(child2)
        else:
            var_child2 = -1
        while (var_child2 > var_el) or (var_child1 > var_el):
            if var_child1 > var_child2:
                self.heap[el], self.heap[child1] = self.heap[child1], self.heap[el]
                el = child1
            else:
                self.heap[el], self.heap[child2] = self.heap[child2], self.heap[el]
                el = child2
            var_el = self.mult(el)
            child1 = el * 2
            if child1 < len(self.heap):
                var_child1 = self.mult(child1)
            else:
                var_child1 = -1
            child2 = el * 2 + 1
            if child2 < len(self.heap):
                var_child2 = self.mult(child2)
            else:
                var_child2 = -1
        return ans

n, m, r , k = map(int, input().split())

max_x = min(r, n - r + 1)
max_y = min(r, m - r + 1)
col_max_x = n - 2 * (max_x - 1)
col_max_y = m - 2 * (max_y - 1)
Nr = (n - r + 1) * (m - r + 1)


tuple_heap = TupleHeap()
tuple_heap.add((max_x, max_y))
used = {(max_x, max_y)}
SUMM = 0

curr = tuple_heap.pop()
while (curr is not None) and (curr[0] > 0) and (curr[1] > 0) and (k > 0):
    # кладём в кучу следующие элементы (если их там нет)
    new_tuple1 = (curr[0] - 1, curr[1])
    if new_tuple1 not in used:
        used.add(new_tuple1)
        tuple_heap.add(new_tuple1)
    new_tuple2 = (curr[0], curr[1] - 1)
    if new_tuple2 not in used:
        used.add(new_tuple2)
        tuple_heap.add(new_tuple2)
    # вычисляем кол-во элементов поля для такого curr
    x = 2 if curr[0] != max_x else col_max_x
    y = 2 if curr[1] != max_y else col_max_y
    col = min(x * y, k)                          # размещаем рыбок в квадраты, но не более 1 рыбки на квадрат
    k = k - col                           # сколько ещё рыбок надо разместить
    SUMM = SUMM + col * curr[0] * curr[1]
    curr = tuple_heap.pop()

#print(tuple_heap)
#print(used)
#print(max_x, max_y)
#print(col_max_x, col_max_y)
#print("SUMM =", SUMM)
#print("Nr =", Nr)
print(SUMM/Nr)

from heapq import *
n, m, r, k = map(int, input().split())
u, v = n // 2, m // 2
h = []
g = lambda z, l: min(z + 1, l - z, l - r + 1, r)
def f(x, y):
    if 0 <= x < n and 0 <= y < m:
        s = g(x, n) * g(y, m)
        heappush(h, (-s, x, y))
f(u, v)
t = 0
for i in range(k):
    s, x, y = heappop(h)
    t -= s
    if x <= u: f(x - 1, y)
    if x == u and y <= v: f(x, y - 1)
    if x >= u: f(x + 1, y)
    if x == u and y >= v: f(x, y + 1)
print(t / (n - r + 1) / (m - r + 1))

import queue
q = queue.PriorityQueue()
n, m, r, k = map(int, input().split())
a = min(r, n - r + 1)
b = min(r, m - r + 1)
u = n - 2 * a + 2
v = m - 2 * b + 2
for i in range(1, a + 1): q.put((-b * i, i))
t = 0
while k:
    s, i = q.get()
    q.put((s + i, i))
    d = min((u if i == a else 2) * (v if s == -b * i else 2), k)
    t -= s * d
    k -= d
print(t / (n - r + 1) / (m - r + 1))

# https://codeforces.com/problemset/problem/912/D
import heapq
from heapq import heappush as push
from heapq import heappop  as pop

base = 100009
Q=[]
used={}
ans=0

dx = [0,0,-1,1]
dy = [-1,1,0,0]

def f(x, y):
    x1, x2 = min(x, r-1),   max(0, x-n+r)
    y1, y2 = min(y, r-1),   max(0, y-m+r)
    
    return -(x1-x2+1)*(y1-y2+1)
    
def cal(i, j):
    return i*base+j

def get(ind):
    return ind//base, ind%base

def is_ok(x, y, x_, y_):
    if x+x_>=0 and x+x_<n and y+y_>=0 and y+y_<m:
        return True
    return False
   
n, m, r, k = map(int, input().split())    
x0, y0 = n//2, m//2 
Q=[(f(x0, y0), cal(x0, y0))]

while k > 0:
    val, ind = pop(Q)
    
    if ind in used:
        continue
        
    x, y = get(ind)
    used[ind]=True
    val  = -val
    #print((x, y), val)
    ans += val
    for x_, y_ in zip(dx, dy):
        if is_ok(x, y, x_, y_):
            push(Q, (f(x+x_, y+y_), cal(x+x_, y+y_)))
    k-=1
    
print(ans/((n-r+1)*(m-r+1)))