#include<iostream>
#include<cstdio>
#include<cmath>
#define N 502
using namespace std;
const long double eps=1e-7;
int read()
{
	char c=0;
	int sum=0;
	while (c<'0'||c>'9') c=getchar();
	while ('0'<=c&&c<='9') sum=sum*10+c-'0',c=getchar();
	return sum;
}
struct reads
{
	long double a,b;
};
int d,n,sn;
long double A[N+1],B[N+1],dA[N+1],dB[N+1],ans[N+1],X[N+1][N+1],delta[N+1][N+1],G[N+1][N+1],res[N+1],dis[N+1];
void gauss()
{
	long double d;
	for (int i=1;i<=sn;++i)
		for (int j=i+1;j<=sn;++j)
		{
			d=G[j][i]/G[i][i];
			for (int k=i;k<=sn+2;++k) G[j][k]-=G[i][k]*d;
		}
	for (int i=sn;i>=1;--i)
	{
		A[i]=G[i][sn+1],B[i]=G[i][sn+2];
		for (int j=i+1;j<=sn;++j) A[i]-=G[i][j]*A[j],B[i]-=G[i][j]*B[j];
		A[i]/=G[i][i],B[i]/=G[i][i];
	}
	return;
}
int main()
{
	long double tA=0,tB=0,tC=0,dst,rest;
	bool op=1;
	d=read(),n=read(),sn=min(n,d);
	for (int i=1;i<=n;++i)
	{
		for (int j=1;j<=d;++j) cin>>X[i][j];
		for (int j=1;j<=d;++j) res[i]+=X[i][j]*X[i][j];
		cin>>dis[i],dis[i]=dis[i]*dis[i]-res[i];
	}
	for (int i=1;i<=sn;++i)
		for (int j=1;j<=sn;++j)
			for (int k=1;k<=d;++k)
				delta[i][j]+=X[i][k]*X[j][k];
	for (int i=1;i<=sn;++i)
	{
		for (int j=1;j<=sn;++j) G[i][j]=delta[i][j];
		G[i][sn+1]=0.5,G[i][sn+2]=-dis[i]/2;
	}
	gauss();
	for (int i=1;i<=d;++i)
	{
		for (int j=1;j<=sn;++j) dA[i]+=A[j]*X[j][i],dB[i]+=B[j]*X[j][i];
		tA+=dA[i]*dA[i],tB+=2*dA[i]*dB[i],tC+=dB[i]*dB[i];
	}
	tB--,dst=(-tB-sqrt(tB*tB-4*tA*tC))/(2*tA);
	for (int i=1;i<=d;++i) ans[i]=dA[i]*dst+dB[i];
	for (int i=sn+1;i<=n;++i)
	{
		rest=0;
		for (int j=1;j<=d;++j) rest+=ans[j]*X[i][j];
		op&=(abs(rest-(dst-dis[i])/2)/max(rest,(long double)(1))<eps);
	}
	if (!op)
	{
		dst=(-tB+sqrt(tB*tB-4*tA*tC))/(2*tA);
		for (int i=1;i<=d;++i) ans[i]=dA[i]*dst+dB[i];
	}
	for (int i=1;i<=d;++i) printf("%0.6Lf ",ans[i]);
	puts("");
	return 0;
}