###スニペット始まり###
import sys, re
from copy import copy, deepcopy
from math import ceil, floor, sqrt,factorial, gcd, pi, degrees, radians, sin, asin, cos, acos, tan, atan2
from statistics import mean, median
from collections import Counter, deque, defaultdict
from heapq import heapify, heappop, heappush
from itertools import permutations, accumulate, product, combinations, combinations_with_replacement
from bisect import bisect, bisect_left, bisect_right
from functools import reduce, lru_cache
from string import ascii_uppercase, ascii_lowercase
from decimal import Decimal, ROUND_HALF_UP #四捨五入用

def input(): return sys.stdin.readline().rstrip('\n')
#easy-testのpypyでは再帰数を下げる。
if __file__=='prog.py':
    sys.setrecursionlimit(10**5)
else:
    sys.setrecursionlimit(10**6)

def lcm(a, b): return a * b // gcd(a, b)

#3つ以上の最大公約数/最小公倍数。Nを要素数、Mを数値の大きさとして、O(NlogM)
def gcd_v2(l: list): return reduce(gcd, l)
def lcm_v2(l: list): return reduce(lcm, l)

#nPk
def nPk(n, k): return factorial(n) // factorial(n - k)

#逆元
def modinv(a, mod=10**9+7): return pow(a, mod-2, mod)
INF = float('inf')
MOD = 10 ** 9 + 7
###スニペット終わり###

N, M=map(int, input().split())
A=list(map(int, input().split()))

es=[]
for _ in range(M):
    a, b = map(int, input().split())
    es.append([a,b])

K=int(input())
B=list(map(int, input().split()))

lamps=[0]*N #0がoff、1がon
for b in B:
    lamps[b-1]=1

graph = [[] for _ in range(N)]
input_edge=[0]*N
for i in range(M):
    a,b=es[i]
    if A[a-1]<A[b-1]:
        graph[a-1].append(b-1)
        if lamps[a-1]:
            input_edge[b-1]+=1
    elif A[a-1]>A[b-1]:
        graph[b-1].append(a-1)  
        if lamps[b-1]:
            input_edge[b-1]+=1

ans = [i for i in range(N) if input_edge[i] == 0 and lamps[i]==1]
que=deque(ans)
while que:
    q = que.popleft()
    lamps[q]=0 #消す。
    for e in graph[q]:
        input_edge[e] -= 1
        lamps[e]=0 if lamps[e]==1 else 1
        if input_edge[e] == 0 and lamps[e]==1:
            que.append(e)
            ans.append(e)
print(len(ans))
for a in ans:
    print(a+1)