#include <bits/stdc++.h>
using namespace std;
typedef pair <int,int> pii;

const int RLEN=1<<18|1;
inline char nc() {
	static char ibuf[RLEN],*ib,*ob;
	(ib==ob) && (ob=(ib=ibuf)+fread(ibuf,1,RLEN,stdin));
	return (ib==ob) ? -1 : *ib++;
}
inline int rd() {
	char ch=nc(); int i=0,f=1;
	while(!isdigit(ch)) {if(ch=='-')f=-1; ch=nc();}
	while(isdigit(ch)) {i=(i<<1)+(i<<3)+ch-'0'; ch=nc();}
	return i*f;
}

const int N=1e5+50;

int n,ans[N];
vector <int> g[N];
struct Tree {
	int rt,col[N];
	vector <int> edge[N];
	inline int dfs(int x,int f) {
		vector <int> vec;
		int sum=1;
		for(auto v:edge[x]) {
			int t=dfs(v,x);
			sum^=1;
			if(t) vec.push_back(t);
		}
		if(sum) vec.push_back(x);
		while(vec.size()>2) {
			int u=vec.back(); vec.pop_back();
			int v=vec.back(); vec.pop_back();
			g[u].push_back(v); g[v].push_back(u);
		} return col[x]=sum,vec.size()?vec[0]:0;
	}
	inline void init() {
		for(int i=1;i<=n;i++) {
			int f=rd();
			if(~f) edge[f].push_back(i);
			else rt=i; 
		} dfs(rt,0);
	}
} T[2];
inline void dfs(int x,int c) {
	ans[x]=c;
	for(auto v:g[x])
		if(!~ans[v]) dfs(v,c^1);
}
int main() {
	n=rd(); 
	T[0].init(); 
	T[1].init();
	for(int i=1;i<=n;i++) 
		if(T[0].col[i]^T[1].col[i]) return puts("IMPOSSIBLE"),0;
	for(int i=1;i<=n;i++) ans[i]=-1;
	for(int i=1;i<=n;i++) if(!~ans[i]) dfs(i,0);
	puts("POSSIBLE");
	for(int i=1;i<=n;i++)
		printf("%d ",T[0].col[i] ? (ans[i]?1:-1) :0);
}