#include <bits/stdc++.h>
using namespace std;
int Qread()
{
	int x=0;bool f=false;char ch=getchar();
	while(ch<'0'||ch>'9') f|=(ch=='-'),ch=getchar();
	while(ch>='0'&&ch<='9') x=x*10+(ch^48),ch=getchar();
	return f?-x:x;
}
typedef pair<int,int> pr;
int n,bk1[200010],bk2[200010];
pr p1[200010],p2[200010];
int X[400010],Y[400010],totx,toty;
int fl1,fl2;
int pp[200010],p_[200010],fa[200010];

#define ls (pos<<1)
#define rs (pos<<1|1)
#define mid (l+r>>1)
int s[2000010],mx[200010],mn[200010];
void pushdown(int pos){if(s[pos]) mx[ls]=mn[ls]=s[ls]=mx[rs]=mn[rs]=s[rs]=s[pos],s[pos]=0;}
void update(int pos){mx[pos]=max(mx[ls],mx[rs]),mn[pos]=min(mn[ls],mn[rs]);}
int get_nd(int x,int pos=1,int l=1,int r=toty)
{
	if(l==r) return s[pos];
	if(x<=mid) return pushdown(pos),get_nd(x,ls,l,mid);
	else return pushdown(pos),get_nd(x,rs,mid+1,r);
}
void range_change(int L,int R,int x,int pos=1,int l=1,int r=toty)
{
	if(L<=l&&r<=R) return s[pos]=x,void();
	if(r<L||R<l) return;
	return pushdown(pos),range_change(L,R,x,ls,l,mid),range_change(L,R,x,rs,mid+1,r),update(pos);
}
int range_max(int L,int R,int pos=1,int l=1,int r=toty)
{
	if(L<=l&&r<=R) return mx[pos];
	if(r<L||R<l) return -20070610;
	else return pushdown(pos),max(range_max(L,R,ls,l,mid),range_max(L,R,rs,mid+1,r));
}
int range_min(int L,int R,int pos=1,int l=1,int r=toty)
{
	if(L<=l&&r<=R) return mn[pos];
	if(r<L||R<l) return 20070610;
	else return pushdown(pos),min(range_min(L,R,ls,l,mid),range_min(L,R,rs,mid+1,r));
}
void show(int pos=1,int l=1,int r=toty)
{
	if(l==r) printf("%d ",s[pos]);
	else pushdown(pos),show(ls,l,mid),show(rs,mid+1,r);
}
#undef ls
#undef rs
#undef mid

#define x first
#define y second
struct Node{pr loc;int id;};
bool operator<(Node A,Node B)
{
	if(A.loc.y!=B.loc.y) return A.loc.y<B.loc.y;
	if(A.loc.x!=B.loc.x) return A.loc.x<B.loc.x;
	return false;
}
set<Node> S;set<Node>::iterator it;
int main()
{
	freopen("T5.in","r",stdin);
    freopen("T5.out","w",stdout);

	n=Qread();
	for(int i=1;i<=n;i++) p1[i].x=Qread(),p1[i].y=Qread();
	for(int i=1;i<=n;i++) p2[i].x=Qread(),p2[i].y=Qread();

	for(int i=1;i<=n;i++) X[++totx]=p1[i].x,Y[++toty]=p1[i].y;
	for(int i=1;i<=n;i++) X[++totx]=p2[i].x,Y[++toty]=p2[i].y;

	sort(X+1,X+totx+1);totx=unique(X+1,X+totx+1)-X-1;
	sort(Y+1,Y+toty+1);toty=unique(Y+1,Y+toty+1)-Y-1;

	for(int i=1;i<=n;i++) p1[i].x=lower_bound(X+1,X+totx+1,p1[i].x)-X,p1[i].y=lower_bound(Y+1,Y+toty+1,p1[i].y)-Y;
	for(int i=1;i<=n;i++) p2[i].x=lower_bound(X+1,X+totx+1,p2[i].x)-X,p2[i].y=lower_bound(Y+1,Y+toty+1,p2[i].y)-Y;

	// for(int i=1;i<=n;i++) printf("(%d %d)",p1[i].x,p1[i].y);printf("\n");
	// for(int i=1;i<=n;i++) printf("(%d %d)",p2[i].x,p2[i].y);printf("\n");

	for(int i=1;i<=n;i++) bk2[i]=i;
	sort(bk2+1,bk2+n+1,[&](int a,int b){return p1[a]<p1[b];});
	for(int i=1;i<=n;i++) bk1[bk2[i]]=i;
	sort(bk2+1,bk2+n+1,[&](int a,int b){return p2[a]<p2[b];});
	sort(p1+1,p1+n+1),sort(p2+1,p2+n+1);

	// for(int i=1;i<=n;i++) printf("(%d %d)",p1[i].x,p1[i].y);printf("\n");
	// for(int i=1;i<=n;i++) printf("(%d %d)",p2[i].x,p2[i].y);printf("\n");

	fl1=fl2=1;
	while(fl1<=n&&fl2<=n)
	{
		if(p2[fl2].x<=p1[fl1].x)
		{
			it=S.lower_bound(Node{make_pair(0,p2[fl2].y),fl2});
			if(it==S.begin()) return printf("No\n"),0;
			it--;
			pp[it->id]=fl2,p_[fl2]=it->id;S.erase(it);
			fl2++;
		}
		else S.insert(Node{p1[fl1],fl1}),fl1++;
	}
	if(fl1<=n) return printf("No\n"),0;
	while(fl2<=n)
	{
		it=S.lower_bound(Node{make_pair(0,p2[fl2].y),fl2});
		if(it==S.begin()) return printf("No\n"),0;
		it--;
		pp[it->id]=fl2;S.erase(it);
		fl2++;
	}
	
	// for(int i=1;i<=n;i++) printf("[%d %d]*[%d %d]\n",p1[i].x,p2[pp[i]].x,p1[i].y,p2[pp[i]].y);

	if(n>=1000) printf("--------\n");

	fl1=fl2=1;
	int tl,tr;
	while(fl1<=n&&fl2<=n)
	{
		if(p2[fl2].x<=p1[fl1].x)
		{
			range_change(p1[p_[fl2]].y,p2[fl2].y-1,fa[p_[fl2]]);
			fl2++;
		}
		else
		{
			tl=get_nd(p1[fl1].y),tr=get_nd(p2[pp[fl1]].y-1);
			if(tl!=tr||range_max(p1[fl1].y,p2[pp[fl1]].y-1)!=range_min(p1[fl1].y,p2[pp[fl1]].y-1)) return printf("No\n"),0;
			if(tl&&p2[pp[fl1]].x>p2[pp[tl]].x) return printf("No\n"),0;
			range_change(p1[fl1].y,p2[pp[fl1]].y-1,fl1);
			fa[fl1]=tl;
			fl1++;
		}

		// printf("%d %d:\n",fl1,fl2);
		// show();printf("\n");
	}

	printf("Yes\n");
	for(int i=1;i<=n;i++) printf("%d\n",bk2[pp[bk1[i]]]);
	return 0;
}