#include<bits/stdc++.h>
using namespace std;
struct Edge
{
	int end,rev;
	bool del;
};
int n,m,S,T;
vector<Edge> t[1100];
bool usd[1100];
vector<int> g[1100];
int val[1100];
bool vis[1100];
vector<vector<int>> opt;
void addEdge(int u,int v)
{
	// cout<<"AddEdge "<<u<<" "<<v<<endl;
	t[u].push_back((Edge){v,t[v].size(),0});
	t[v].push_back((Edge){u,t[u].size()-1,0});
}
namespace bipolarOrientation
{
	int dfn[1100],clk;
	vector<int> nxt[1100];
	bool vis[1100];
	bool shrink(int u,int f,vector<int> &chn)
	{
		// cout<<"Shrink u="<<u<<" f="<<f<<endl;
		dfn[u]=++clk;
		bool rchS=(u==S);
		int aid=-1,fid=-1;
		for(int i=0;i<t[u].size();i++)
		{
			int v=t[u][i].end;
			// cout<<"u="<<u<<" v="<<v<<" del="<<t[u][i].del<<endl;
			if(v==f)
				fid=i;
			if(v==f||t[u][i].del)
				continue;
			if(!dfn[v])
				rchS|=shrink(v,u,chn);
			else if(aid==-1||dfn[v]<dfn[t[u][aid].end])
				aid=i;
		}
		if(!rchS)
		{
			// cout<<"Shrinking "<<u<<" fid="<<fid<<" aid="<<aid<<endl;
			// if(aid!=-1)
			// 	cout<<"aid="<<aid<<" end="<<t[u][aid].end<<" dfn="<<dfn[t[u][aid].end]<<endl;
			if(!(fid!=-1&&aid!=-1&&dfn[t[u][aid].end]<dfn[u]))
			{
				// cout<<"Err at "<<u<<endl;
				cout<<-1<<endl;
				exit(0);
			}
			nxt[t[u][fid].end].push_back(u);
			nxt[t[u][aid].end].push_back(u);
			addEdge(t[u][fid].end,t[u][aid].end);
			for(int i=0;i<t[u].size();i++)
				t[u][i].del=t[t[u][i].end][t[u][i].rev].del=1;
		}
		else
			chn.push_back(u);
		return rchS;
	}
	void expand(int u,vector<int> &res)
	{
		vis[u]=1,res.push_back(u);
		for(int v:nxt[u])
		{
			if(!vis[v])
				expand(v,res);
		}
	}
	vector<int> solve()
	{
		vector<int> chn(0),res(0);
		shrink(T,0,chn);
		for(int u:chn)
			expand(u,res);
		return res;
	}
}
bool findPathOut(int u,int T,vector<int> &res)
{
	if(u==T)
		return 1;
	vis[u]=1;
	for(int i=0;i<t[u].size();i++)
	{
		int v=t[u][i].end;
		if(!vis[v]&&!usd[v])
		{
			if(findPathOut(v,T,res))
			{
				res.push_back(u);
				return 1;
			}
		}
	}
	return 0;
}
void findPathIn(int u,vector<int> &res)
{
	res.push_back(u);
	int mnf=0;
	for(int i=0;i<g[u].size();i++)
	{
		int v=g[u][i];
		if(usd[v]&&(mnf==0||val[v]<val[mnf]))
			mnf=v;
	}
	if(mnf==0||val[mnf]>val[u])
		return;
	swap(val[u],val[mnf]);
	findPathIn(mnf,res);
}
void appendPoint(int u)
{
	memset(vis,0,sizeof(vis));
	vector<int> pth(0);
	findPathOut(S,u,pth);
	reverse(pth.begin(),pth.end());
	for(int i=0;i<pth.size();i++)
	{
		swap(val[pth[i]],val[i==pth.size()-1?u:pth[i+1]]);
	}
	for(int i=0;i<t[u].size();i++)
	{
		int v=t[u][i].end;
		if(usd[v])
			g[u].push_back(v);
	}
	usd[u]=1;
	findPathIn(u,pth);
	opt.push_back(pth);
}
int main()
{
	cin>>n>>m>>S>>T;
	for(int i=1;i<=n;i++)
		cin>>val[i];
	for(int i=1;i<=m;i++)
	{
		int a,b;
		cin>>a>>b;
		addEdge(a,b);
	}
	vector<int> bpo=bipolarOrientation::solve();
	// cout<<"> ";
	// for(int i=0;i<bpo.size();i++)
	// 	cout<<bpo[i]<<" ";
	// cout<<endl;
	usd[bpo.back()]=1;
	for(int i=bpo.size()-2;i>=0;i--)
		appendPoint(bpo[i]);
	cout<<opt.size()<<endl;
	for(int i=0;i<opt.size();i++)
	{
		cout<<opt[i].size()<<" ";
		for(int j=0;j<opt[i].size();j++)
			cout<<opt[i][j]<<" ";
		cout<<endl;
	}
	// for(int i=1;i<=n;i++)
	// 	cout<<val[i]<<" ";
	// cout<<endl;
}