#include<bits/stdc++.h>
using namespace std;
using LL=long long;
using u64=unsigned long long;
using LF=long double;
#define F(i,l,r) for(int i=l;i<=r;++i)
#define G(i,r,l) for(int i=r;i>=l;--i)
#define CT const&
#define IL inline
#define PB push_back
template<class T>
IL void tomn(T &x,T CT y){y<x?x=y,0:0;}
template<class T>
IL void tomx(T &x,T CT y){x<y?x=y,0:0;}
#define DEBUG(x) cerr<<"line:"<<__LINE__<<" "#x"="<<x<<endl
#define CUT cerr<<"**********\n"
const LF EPS=1e-10L;
IL int dcmp(LF x){return fabs(x)<=EPS?0:(x<0?-1:1);}
const LL INF=4e18L;

const int N=85;
void exgcd(int a,int b,int &x,int &y,int &g){
	if(!b)return x=1,y=0,g=a,void();
	exgcd(b,a%b,y,x,g),y-=(a/b)*x;
}
IL int mod(int x,int y){return x%=y,x+(x<0)*y;}
IL int dv(int x,int y){return (x-mod(x,y))/y;}
struct D{int x,m;}f[N][N][N];
bool operator<(D x,D y){return x.x<y.x;}
int n,g[N],ans[N];
void dfs(int l,int r,int x){
	if(l==r)return;
	D d=f[l][r][x%(r-l)];
	F(i,l,r-1)ans[i]+=(d.x-x)/(r-l);
	dfs(l,d.m,d.x),dfs(d.m+1,r,d.x);
}
IL bool ck(int l,int r,int x){return x<=f[l][r][x%(r-l)].x;}
IL void to(int l,int r,D x){tomx(f[l][r][x.x%(r-l)],x);}
int main(){
#ifdef LOCAL
	freopen("B.in","r",stdin);
//	freopen(".out","w",stdout);
#endif
	scanf("%d",&n);
	if(n==1)puts("Yes"),exit(0);
	F(i,1,n)scanf("%d",g+i);
	F(i,0,N-1)F(j,0,N-1)F(k,0,N-1)f[i][j][k]=D{-1};
	G(l,n,1)F(r,l,n){
		if(l==r){
			f[l][l][0]=D{g[l]};
			continue;
		}if(l==r-1){
			if(g[l]==g[r])tomx(f[l][r][0],D{g[l],l});
			continue;
		}
		F(m,l,r-1){
			if(m==l){
				if(ck(l+1,r,g[l]))to(l,r,D{g[l],m});
				continue;
			}if(m==r-1){
				if(ck(l,r-1,g[r]))to(l,r,D{g[r],m});
				continue;
			}int p=m-l,q=r-m-1,x,y,g,L,M=r-l;
			exgcd(p,q,x,y,g),L=p/g*q,M/=__gcd(L,M);
//			int mn1=max(dv(-x*g+q-1,q),dv(y*g+p-1,p));
//			int mn2=max(dv(x*g+q-1,q),dv(-y*g+p-1,p));
			F(i,0,L-1){
				int j=i%p,k=i%q,a=f[l][m][j].x,b=f[m+1][r][k].x;
				if(a<0||b<0)continue;
				int x2=(a-b)/g*x,y2=(b-a)/g*y,X=max(dv(-x2*g+q-1,q),dv(-y2*g+p-1,p));
				int c=a-p*(x2+q/g*X);
				{
//					cerr<<x2<<" "<<y2<<endl;
//					int x2=x,y2=y;
//					x2*=(a-b)/g,y2*=-(a-b)/g;
//					cerr<<a-x2*p<<" "<<b-y2*q<<endl;
//					x2+=q/g*X,y2+=p/g*X;
//					cerr<<x2<<" "<<y2<<endl;
//					cerr<<a-x2*p<<" "<<b-y2*q<<endl;
				}
//				int c=a-p/g*(a-b)*(x+q/g*(a>=b?mn1:-mn2));
//				cerr<<a<<" "<<p<<" "<<b<<" "<<q<<" "<<c<<endl;
				F(j,0,M-1){
					if(c<0)break;
					to(l,r,D{c,m}),c-=L;
				}
			}
		}
	}if(!ck(1,n,0))puts("No"),exit(0);
	puts("Yes"),dfs(1,n,0);
	F(i,1,n-1)printf("%d ",ans[i]);
}