#include<iostream.h>
#include<math.h>
#include<iomanip.h>
void main()
{
int i=0;
int k=0;
int n=0;
float x[4]={1,1,2,3};
float y[4]={1,14,23,33};
float a[4]={1,1,2,3};
float u[4]={1,0,2,3};
float v[4]={1,1,2,3};
float t[4]={1,1,2,3};
float sy=0;
float u1=0;
float s1=0;
float f=0;
float q=0;
float p=0;
float s0=0;
cin>>n;
cout<<n<<endl;
/*for(i=1;i<=4;i++)
{
cin>>x[i];
}
for(i=1;i<=4;i++)
{
cin>>y[i];
} */
for(i=1;i<=4;i++)
{
sy+=y[i];
u1+=x[i];
}
a[0]=sy/4;
q=1;
u[0]=u1/4;
v[0]=0;
s0=4;
for(k=0;k<=n-1;k++)
{
s1=0;
sy=0;
for(i=1;i<=4;i++)
{
f=(x[i]-u[k])*q-v[0]*p;
sy+=f*y[i];
s1+=f*f;
u[k+1]+=x[i]*f*f;
p=q;
q=f;
}
a[k+1]=(sy)/(s1);
u[k+1]=u[k+1]/(s1);
v[k+1]=(s1)/(s0);
(s0)=(s1);
cout<<a[k+1]<<" ";
}
cout<<endl;
for(i=1;i<=4;i++)
{
cin>>t[i];
s0=1;
s1=t[i]-u[0];
f=a[0]*(s0)+a[1]*(s1);
for(k=1;k<=n-1;k++)
{
sy=(t[i]-u[k])*(s1)-v[k]*(s0);
f=f+a[k+1]*(sy);
s0=s1;
s1=f;
}
cout<<"f="<<f<<endl;
}
}
