from math import sqrt,sin,cos,tan,ceil,radians,floor,gcd,exp,log,log10,log2,factorial,fsum
from heapq import heapify, heappop, heappush
from bisect import bisect_left, bisect_right
from copy import deepcopy
import copy
import random
from collections import deque,Counter,defaultdict
from itertools import permutations,combinations
from decimal import Decimal,ROUND_HALF_UP
#tmp = Decimal(mid).quantize(Decimal('0'), rounding=ROUND_HALF_UP)
from functools import lru_cache, reduce
#@lru_cache(maxsize=None)
from operator import add,sub,mul,xor,and_,or_,itemgetter
import sys
input = sys.stdin.readline
# .rstrip()
INF = 10**18
mod1 = 10**9+7
mod2 = 998244353

#DecimalならPython
#再帰ならPython!!!!!!!!!!!!!!!!!!!!!!!!!!


'''
dfsじゃないといけないかも
頂点Gまでの経路
まあ普通にP回まで探索すればいいけど
それまでに通った辺をどうやって管理するか
重複していいならリストでいいのか
逆に、P回移動したときにたどりつく頂点を列挙して
そこからbfsするとか
ちょっと厳しいか?

なるほど
前後からか
ちょっと答え見るの早いかもな
まあ今回はしょうがないとして
'''

N,M,P = map(int, input().split())
S,Go = map(int, input().split())

G = [[] for _ in range(N+1)]

for _ in range(M):
    u,v = map(int, input().split())
    G[u].append(v)
    G[v].append(u)

Q = deque()
Q.append((S,0))

seen = [INF]*(N+1)

while len(Q) > 0:
    pos,cnt = Q.popleft()
    
    seen[pos] = cnt
    
    for nx in G[pos]:
        if seen[nx] == INF:
            Q.append((nx,cnt+1))
            seen[nx] = cnt+1

Q = deque()
Q.append((Go,0))

seen2 = [INF]*(N+1)

while len(Q) > 0:
    pos,cnt = Q.popleft()
    
    seen2[pos] = cnt
    
    for nx in G[pos]:
        if seen2[nx] == INF:
            Q.append((nx,cnt+1))
            seen2[nx] = cnt+1

ans = []
for i in range(1,N+1):
    if seen[i]+seen2[i] > P:
        continue
    if (seen[i]+seen2[i])%2 == P%2:
        ans.append(i)

ans.sort()

if len(ans) == 0:
    print(-1)

print(len(ans))
for a in ans:
    print(a)