import sys
input = sys.stdin.readline

H,W,K,P=map(int,input().split())
mod=10**9+7

F=[input().split() for i in range(K)]
for i in range(K):
    F[i][0]=int(F[i][0])
    F[i][1]=int(F[i][1])

MAX=0
ind=0

for i in range(1<<K):
    DP=[[0]*(W+1) for i in range(H+1)]
    DP[0][0]=1

    p=0

    for j in range(K):
        if i & (1<<j) == 0:
            x,y,_ = F[j]
            DP[x][y]=-1
        else:
            p+=1

    if p>P:
        continue

    for j in range(H+1):
        for k in range(W+1):
            if DP[j][k]!=-1:
                if j-1>=0 and DP[j-1][k]!=-1:
                    DP[j][k]+=DP[j-1][k]
                if k-1>=0 and DP[j][k-1]!=-1:
                    DP[j][k]+=DP[j][k-1]

    if DP[-1][-1]>MAX:
        MAX=DP[-1][-1]
        ind=i
            
print(MAX%mod)

for j in range(K):
    if ind & (1<<j) != 0:
        print(F[j][2])