#include <stdio.h>
#include <complex.h>
#include <math.h>

void matwrite(int n, double complex A[n][n]){
	int i,j;
	for(i=0;i<n;i++) {
		for(j=0;j<n;j++) {
			printf("(%g %g) ",
				creal(A[i][j]),cimag(A[i][j]));
		}
		printf("\n");
	}
}

double complex dotprod(int n, double complex u[n],
	double complex v[n]){
	double complex s = 0;
	int i;
	for(i=0;i<n;i++){
		s+=u[i]*conj(v[i]);
	}
	return s;
}

void multmatvec(int n,double complex v[n],
	double complex A[n][n],
	double complex u[n]){
	int i,j;
	for(i=0;i<n;i++) v[i]=0;
	for(i=0;i<n;i++) for(j=0;j<n;j++) {
		v[i]+=A[i][j]*u[j];
	}
}

int main(){
	int n;
	scanf("%d",&n);
	printf("n=%d\n",n);
	double complex A[n][n];
	int i,j;
	for(i=0;i<n;i++) for(j=0;j<n;j++){
		double x,y;
		scanf("%lg %lg",&x,&y);
		A[i][j]=x+I*y;
	}
	matwrite(n,A);
	double complex X[n],Y[n];
	int k;
	for(i=0;i<n;i++) X[i]=2*i+I*(2*i+1);
	for(k=0;k<10;k++) {
		double complex sqxdotx=sqrt(creal(dotprod(n,X,X)));
		for(i=0;i<n;i++) X[i]/=sqxdotx;
		multmatvec(n,Y,A,X);
		double complex lambda=dotprod(n,Y,X);
		printf("lambda=%g %g\n",
			creal(lambda),cimag(lambda));
		for(i=0;i<n;i++) X[i]=Y[i];
	}
	return 0;
}