h, w, k, p = map(int, input().split())
dat = []
for i in range(k):
    x, y, n = input().split()
    x = int(x)
    y = int(y)
    dat.append([x, y, n])
ret = 0
num = 0
for i in range(1 << k):
    if p < bin(i).count("1"):
        continue
    cnt = [[0] * (w + 1) for j in range(h + 1)]
    cnt[0][0] = 1
    for j in range(k):
        if not (i >> j) & 1:
            cnt[dat[j][0]][dat[j][1]] = -1
    for x in range(h + 1):
        for y in range(w + 1):
            if cnt[x][y] == -1:
                continue
            if x + 1 <= h and cnt[x + 1][y] != -1:
                cnt[x + 1][y] += cnt[x][y]
            if y + 1 <= w and cnt[x][y + 1] != -1:
                cnt[x][y + 1] += cnt[x][y]
    if cnt[h][w] > num:
        num = cnt[h][w]
        ret = i
print(num % (10**9 + 7))
for i in range(k):
    if (ret >> i) & 1:
        print(dat[i][2])