Language: C++

// This amazing code is by Eric Sunli Chen.
#include<bits/stdc++.h>
using namespace std;
template<typename T> bool get_int(T &x)
{
	char t=getchar();
	bool neg=false;
	x=0;
	for(; (t>'9'||t<'0')&&t!='-'&&t!=EOF; t=getchar());
	if(t=='-')neg=true,t=getchar();if(t==EOF)return false;
	for(; t<='9'&&t>='0'; t=getchar())x=x*10+t-'0';
	if(neg)x=-x;return true;
}
template<typename T> void print_int(T x)
{
	if(x<0)putchar('-'),x=-x;
	short a[20]= {},sz=0;
	while(x>0)a[sz++]=x%10,x/=10;
	if(sz==0)putchar('0');
	for(int i=sz-1; i>=0; i--)putchar('0'+a[i]);
}
#define ff first
#define ss second
#define pb push_back
#define mp make_pair
#define get1(a) get_int(a)
#define get2(a,b) (get1(a)&&get1(b))
#define get3(a,b,c) (get1(a)&&get2(b,c))
#define printendl(a) print_int(a),puts("")
typedef long long LL;
typedef unsigned long long uLL;
typedef pair<int,int> pii;
const int inf=0x3f3f3f3f;
const LL Linf=1ll<<61;
const double pi=acos(-1.0);

const int maxn=200111;
const int maxm=400111;

struct edge{int to,cap;};
namespace flow
{
	int head[maxn],nxt[maxm<<1],iter[maxn],tot;
	edge e[maxm<<1];
	void clear()
	{
		memset(head,0,sizeof(head));
		tot=1;
	}
	max_flow(){clear();}
	void addedge(int u,int v,int cap)
	{
		nxt[++tot]=head[u];
		head[u]=tot;
		e[tot]=edge{v,cap};
		nxt[++tot]=head[v];
		head[v]=tot;
		e[tot]=edge{u,0};
	}
	 
	int lvl[maxn],q[maxn];
	bool bfs(int s,int t)
	{
		memset(lvl,-1,sizeof(lvl));
		lvl[s]=0;int rr=0;q[rr++]=s;
		for(int fr=0;fr<rr;fr++)
		{
			int x=q[fr];if(x==t)return true;
			for(int i=head[x];i;i=nxt[i])if(e[i].cap&&lvl[e[i].to]==-1)
			{
				lvl[e[i].to]=lvl[x]+1;
				q[rr++]=e[i].to;
			}
		}
		return false;
	}
	int dfs(int x,int t,int f)
	{
		if(x==t||f==0)return f;
		int ret=0;
		for(int&i=iter[x];i;i=nxt[i])if(e[i].cap&&lvl[e[i].to]==lvl[x]+1)
		{
			int d=dfs(e[i].to,t,min(f,e[i].cap));
			e[i].cap-=d;e[i^1].cap+=d;
			f-=d;ret+=d;
			if(!f)break;
		}
		return ret;
	}
	int maxflow(int s,int t)
	{
		int flow=0;
		while(bfs(s,t))
		{
			memcpy(iter,head,sizeof(head));
			flow+=dfs(s,t,inf);
		}
		return flow;
	}
};

vector<int> e[maxn],g[maxn];
int n,piv[maxn],fa[maxn];

bool use[maxn];
void dfs(int x)
{
	use[x]=1;
	for(auto&v:g[x])if(!use[v])
	{
		fa[v]=x;
		dfs(v);
	}
}
int main()
{
	get1(n);
	int S=1,T=n+n;
	for(int i=2;i<=n;i++)flow::addedge(S,i,1);
	for(int i=1;i<n;i++)
	{
		int cnt,x;get1(cnt);
		flow::addedge(i+n,T,1);
		while(cnt--)
		{
			get1(x);
			e[i].pb(x);
			if(x>1)flow::addedge(x,i+n,1);
		}
	}
	if(flow::maxflow(S,T)<n-1)
	{
		puts("-1");
		return 0;
	}
	
	for(int i=2;i<=n;i++)for(int j=flow::head[i];j;j=flow::nxt[j])if(flow::e[j].to!=S&&!flow::e[j].cap)
		piv[flow::e[j].to-n]=i;
	
	for(int i=1;i<n;i++)for(int j=0;j<(int)e[i].size();j++)if(e[i][j]!=piv[i])g[e[i][j]].pb(piv[i]);
	
	dfs(1);
	for(int i=2;i<=n;i++)if(!fa[i])
	{
		puts("-1");
		return 0;
	}
	for(int i=1;i<n;i++)printf("%d %d\n",piv[i],fa[piv[i]]);
	return 0;
}

Compiler output:
a.cc:50:9: error: ISO C++ forbids declaration of 'max_flow' with no type [-fpermissive]
   50 |         max_flow(){clear();}
      |         ^~~~~~~~
a.cc: In function 'int flow::max_flow()':
a.cc:50:28: warning: no return statement in function returning non-void [-Wreturn-type]
   50 |         max_flow(){clear();}
      |                            ^