#include <bits/stdc++.h>
using namespace std;
#define N 105
#define MOD 998244353
mt19937 rand1(0);
int n=100,nw,w[N],rt[N],st[N],a[N][N],z[N][N];
int f(int x,int y) {return (x-1)*(x-2)/2+y;}
int rdm(int l,int r) {return uniform_int_distribution<int>(l,r)(rand1);}
int qry(int x1,int x2,int x3)
{printf("? %d %d %d\n",x1,x2,x3);fflush(stdout);int t;scanf("%d",&t);return t;}
int findRt(int u)
{if(u==rt[u]) return u;int t=findRt(rt[u]);w[u]^=w[rt[u]];return rt[u]=t;}
void W(int &x,int y) {x+=y;if(x>=MOD) x-=MOD;}
int qPow(int x,int y)
{int res=1;for(;y;y/=2,x=1ll*x*x%MOD) if(y&1) res=1ll*res*x%MOD;return res;}
void Gauss(int n)
{
	for(int i=1,t,t1;i<=n;++i)
	{
		if(!a[i][i]) for(int j=i+1;j<=n;++j) if(a[j][i])
		{for(int k=i;k<=n+1;++k) swap(a[i][k],a[j][k]);break;}
		t=qPow(MOD-a[i][i],MOD-2);
		for(int j=i+1;j<=n;++j)
		{
			t1=1ll*t*a[j][i]%MOD;
			for(int k=i;k<=n+1;++k) a[j][k]=(a[j][k]+1ll*t1*a[i][k])%MOD;
		}
	}
	for(int i=n;i;--i)
	{
		for(int j=i+1;j<=n;++j) W(a[i][n+1],MOD-1ll*a[i][j]*a[j][n+1]%MOD);
		a[i][n+1]=1ll*a[i][n+1]*qPow(a[i][i],MOD-2)%MOD;
	}
}
int main()
{
	for(int i=1;i<=5;++i) for(int j=1;j<i;++j) for(int k=1;k<j;++k)
		++nw,a[nw][f(i,j)]=a[nw][f(j,k)]=a[nw][f(i,k)]=1,a[nw][11]=qry(i,j,k);
	Gauss(10);for(int i=1;i<=5;++i) for(int j=1;j<i;++j) z[i][j]=a[f(i,j)][11];
	for(int i=6,t,t1,t2;i<=n;++i)
	{
		iota(rt+1,rt+i,1);fill(w+1,w+i,0);fill(z[i]+1,z[i]+i,-1);
		while(1)
		{
			st[0]=0;for(int j=1;j<i;++j) if(findRt(j)==j && z[i][j]==-1) st[++st[0]]=j;
			if(!st[0]) break;if(st[0]==1) {z[i][st[1]]=qry(i,i,st[1])/2;break;}
			t1=t2=st[rdm(1,st[0])];while(t1==t2) t2=st[rdm(1,st[0])];
			t=qry(i,t1,t2)-z[max(t1,t2)][min(t1,t2)];
			if(t==1) rt[t1]=t2,w[t1]=1;else z[i][t1]=z[i][t2]=t/2;
		}for(int j=1;j<i;++j) if(z[i][j]==-1) t=findRt(j),z[i][j]=z[i][t]^w[j];
	}for(int i=1;i<=n;++i) for(int j=1;j<i;++j) z[j][i]=z[i][j];printf("!\n");
	for(int i=1;i<=n;++i,putchar('\n')) for(int j=1;j<=n;++j) putchar(z[i][j]+48);
	fflush(stdout);return 0;
}