import itertools as iter
import collections as coll
import heapq as hq
import bisect as bis
from decimal import Decimal as dec
from functools import cmp_to_key
import math
import sys
#import pypyjit
#pypyjit.set_param('max_unroll_recursion=-1')
sys.setrecursionlimit(10 ** 6)
inp = sys.stdin.readline
input = lambda : inp()[:-1]
getN = lambda : int(inp())
getNs = lambda : map(int, inp().split())
getList = lambda : list(map(int, inp().split()))
getStrs = lambda n : [input() for _ in [0] * n]
getEdges = lambda n : [[x - 1 for x in getNs()] for _ in [0] * n]
getWEdges = lambda n : [[x - (i < 2) for i, x in enumerate(getNs())] for _ in [0] * n]
def yexit(): print("Yes"); exit(0)
def nexit(): print("No"); exit(0)
pi = 3.141592653589793
mod = 1000000007
MOD = 998244353
INF = 4611686018427387903
dxs = [1, 0, -1, 0];  dys = [0, 1, 0, -1]
#di = coll.defaultdict(int)


"""
Main Code
"""

def solve(N, K, edges, a):
    graph = [[] for _ in [0] * N]
    for x, y in edges:
        graph[x].append(y)
        graph[y].append(x)
    depth = [-1] * N
    depth[0] = 0
    stk = [0]
    while(stk):
        v = stk.pop()
        nd = depth[v] + 1
        for nv in graph[v]:
            if(depth[nv] != -1):
                continue
            depth[nv] = nd
            stk.append(nv)
    lis = [i for i in range(N) if(depth[i] & 1)]
    g = 0
    for k in lis:
        g ^= a[k] % (K + 1)
    return 'K' if(g) else 'P'

for _ in [0] * getN():
    N, K = getNs()
    edges = getEdges(N - 1)
    a = getList()
    print(solve(N, K, edges, a))