#include <algorithm>
#include <cstdio>

using namespace std ;

const int N=2e5+10 ;

struct Edge { int to , nxt , v ; };
Edge e[N<<2] ;
int lst[N] ;
inline void add ( int u , int v )
{
	static int cnt=1 ;
	e[++cnt]=(Edge){v,lst[u]},lst[u]=cnt,
	e[++cnt]=(Edge){u,lst[v]},lst[v]=cnt;
}

int sz[N] ;

int fa[N] , sa , sb ;
int n ;

inline void dfs1 ( int x , int f )
{
	int i , y ;
	for ( i=lst[x],y=e[i].to ; i ; i=e[i].nxt,y=e[i].to )
		if ( y^f ) dfs1(y,x),sz[x]^=1;
	return ;
}
inline int dfs2 ( int x , int f )
{
	int i , y , ret=1 , s=0 ;
	for ( i=lst[x],y=e[i].to ; i ; i=e[i].nxt,y=e[i].to )
		if ( y^f ) ret&=dfs2(y,x),s^=1;
	if ( s^sz[x-n] ) return 0 ;
	return ret ;
}

int ans[N] ;
inline void Euler ( int x )
{
	for ( int &i=lst[x] ; i ; i=e[i].nxt ) 
		if ( !e[i].v ) 
		{
			// printf("From %d to %d\n",x,e[i].to);
			if ( e[i].to==x+n ) ans[x]=1;
			if ( e[i].to==x-n ) ans[x-n]=-1;
			e[i].v=e[i^1].v=1,Euler(e[i].to);
		}
	return ;
}

int main ()
{
	int i ;
	scanf("%d",&n);
	for ( i=1 ; i<=n ; i++ ) 
	{
		scanf("%d",fa+i);
		if ( ~fa[i] ) add(i,fa[i]);
		sa=(~fa[i]?sa:i);
	}
	for ( i=1 ; i<=n ; i++ ) 
	{
		scanf("%d",fa+n+i);
		if ( ~fa[i+n] ) add(i+n,fa[i+n]+n);
		sb=(~fa[i+n]?sb:i+n);
	}
	dfs1(sa,sa);
	if ( !dfs2(sb,sb) ) return puts("IMPOSSIBLE")*0;
	// for ( i=1 ; i<=n ; i++ ) 
	// 	printf("%d ",v[i]); puts("");
	for ( i=1 ; i<=n ; i++ ) 
		if ( !sz[i] ) add(i,i+n);
	int S=2*n+1 ;
	add(S,sa),add(S,sb);
	Euler(S);
	puts("POSSIBLE");
	for ( i=1 ; i<=n ; i++ ) 
		printf("%d ",ans[i]);
	return 0 ;
}