#include<iostream>
using namespace std;
const long mod=1e9+7;
long N,A,B,C;
long X[3][3],E[3][3],T[3][3];
main()
{
	cin>>N>>A>>B>>C;
	N--;
	X[0][0]=X[1][1]=X[2][2]=1;
	X[0][1]=X[1][2]=X[2][0]=mod-1;
	E[0][0]=E[1][1]=E[2][2]=1;
	while(N)
	{
		if(N&1)
		{
			for(int i=0;i<3;i++)for(int j=0;j<3;j++)
			{
				T[i][j]=E[i][j];
				E[i][j]=0;
			}
			for(int i=0;i<3;i++)for(int j=0;j<3;j++)
			{
				for(int k=0;k<3;k++)(E[i][j]+=T[i][k]*X[k][j])%=mod;
			}
		}
		N>>=1;
		for(int i=0;i<3;i++)for(int j=0;j<3;j++)
		{
			T[i][j]=X[i][j];
			X[i][j]=0;
		}
		for(int i=0;i<3;i++)for(int j=0;j<3;j++)
		{
			for(int k=0;k<3;k++)(X[i][j]+=T[i][k]*T[k][j])%=mod;
		}
	}
	for(int i=0;i<3;i++)
	{
		cout<<(E[i][0]*A+E[i][1]*B+E[i][2]*C)%mod<<(i==2?"\n":" ");
	}
}