###スニペット始まり### 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[b-1].append(a-1) input_edge[a-1]+=1 ans = [] que=deque([i for i in range(N) if input_edge[i] == 0]) while que: q = que.popleft() if not lamps[q]: for e in graph[q]: input_edge[e] -= 1 if input_edge[e] == 0: que.append(e) else: #切り替えする。 lamps[q]=0 #消す。 ans.append(q) for e in graph[q]: input_edge[e] -= 1 lamps[e]=0 if lamps[e]==1 else 1 if input_edge[e] == 0: que.append(e) print(len(ans)) for a in ans: print(a+1)