#include <bits/stdc++.h>
                            #include<iostream>
                            #include<cstdio>
                            #include<vector>
                            #include<queue>
                            #include<map>
                            #include<cstring>
                            #include<string>
                            #include <math.h>
                            #include<algorithm>
                        //    #include <boost/multiprecision/cpp_int.hpp>
                            #include<functional>
                     #define int long long
                            #define inf  1000000007
                            #define pa pair<int,int>
                            #define ll long long
                            #define pal pair<double,int>
                            #define ppa pair<string,pa>
                            #define ssa pair<string,int>
                            #define  mp make_pair
                            #define  pb push_back
                            #define EPS (1e-10)
                            #define equals(a,b) (fabs((a)-(b))<EPS)
                     
                            using namespace std;
                     //priority_queue<int, vector<int>, greater<int> > que;
                            class Point{
                            	public:
                            	double x,y;
                            	Point(double x=0,double y=0):x(x),y(y) {}
                            	Point operator + (Point p) {return Point(x+p.x,y+p.y);}
                            	Point operator - (Point p) {return Point(x-p.x,y-p.y);}
                            	Point operator * (double a) {return Point(x*a,y*a);}
                            	Point operator / (double a) {return Point(x/a,y/a);}
                            	double absv() {return sqrt(norm());}
                            	double norm() {return x*x+y*y;}
                            	bool operator < (const Point &p) const{
                            		return x != p.x ? x<p.x: y<p.y;
                            	}
                            	bool operator == (const Point &p) const{
                            		return fabs(x-p.x)<EPS && fabs(y-p.y)<EPS;
                            	}
                            };
                            typedef Point Vector;
                     
                            struct Segment{
                            Point p1,p2;
                            };
                     
                        double hen(Vector a){
                        if(fabs(a.x)<EPS && a.y>0) return acos(0);
                        else if(fabs(a.x)<EPS && a.y<0) return 3*acos(0);
                        else if(fabs(a.y)<EPS && a.x<0) return 2*acos(0);
                        else if(fabs(a.y)<EPS && a.x>0) return 0.0;
                        else if(a.y>0) return acos(a.x/a.absv());
                        else return 2*acos(0)+acos(-a.x/a.absv());
                     
                        }
                     
                string itos( int i ) {
                ostringstream s ;
                s << i ;
                return s.str() ;
                }
                 
                int gcd(int v,int b){
                	if(v>b) return gcd(b,v);
                	if(v==b) return b;
                	if(b%v==0) return v;
                	return gcd(v,b%v);
                }
                            double dot(Vector a,Vector b){
                            	return a.x*b.x+a.y*b.y;
                            }
                            double cross(Vector a,Vector b){
                            	return a.x*b.y-a.y*b.x;
                            }
                        
                double distans(double x1,double y1,double x2,double y2){
                	double rr=(x1-x2)*(x1-x2)+(y1-y2)*(y1-y2);
                	return sqrt(rr);
                	
                }


                         //----------------kokomade tenpure------------
int b[600][600];
int a[600][600]={0}; 


int aa(int z,int zz){
	if(z>0 && zz>0) return a[z][zz];
	else return 0;
	
}

signed main(){
int n,k;
	int H,W,h,w;
cin>>H>>W>>h>>w;

	if(H%h==0 && W%w==0) {
		cout<<"No"<<endl;
		return 0;
	}
	
	if(H%h==0){
		
		for(int i=1;i<=h;i++)for(int j=1;j<=w;j++){
		a[i][j]=-1;
		
		for(int k=i;k<=H;k+=h)for(int l=j;l<=W;l+=w){
			if(k==i && l==j) continue;
			a[k][l]=aa(k-h,l)+aa(k,l-w)-aa(k-h,l-w)-1;
			
		}
		
	}
		
		
		a[H][W]=1;
		for(int i=W-w;i>0;i-=w)a[H][i]=1000000*(W-i)/w;
		for(int i=H-h;i>0;i-=h)for(int j=W-w;j>0;j-=w){
			
		a[i][j]=-a[i+h][j+w]-1+a[i+h][j]+a[i][j+w];
		
		}
		
		
		
		
		cout<<"Yes"<<endl;
	for(int i=1;i<=H;i++){
		for(int j=1;j<=W;j++){
		b[i][j]=a[i][j]-a[i-1][j]-a[i][j-1]+a[i-1][j-1];
			cout<<b[i][j]<<" ";
		}
		cout<<endl;
	}
	return 0;
		
	}
	
	
	if(W%w==0){
		
		for(int i=1;i<=h;i++)for(int j=1;j<=w;j++){
		a[i][j]=-1;
		
		for(int k=i;k<=H;k+=h)for(int l=j;l<=W;l+=w){
			if(k==i && l==j) continue;
			a[k][l]=aa(k-h,l)+aa(k,l-w)-aa(k-h,l-w)-1;
			
		}
		
	}
		
		
		a[H][W]=1;
		for(int i=H-h;i>0;i-=h)a[i][W]=1000000*(H-i)/h;
		for(int i=H-h;i>0;i-=h)for(int j=W-w;j>0;j-=w){
			
		a[i][j]=-a[i+h][j+w]-1+a[i+h][j]+a[i][j+w];
		
		}
		
		
		
		
		cout<<"Yes"<<endl;
	for(int i=1;i<=H;i++){
		for(int j=1;j<=W;j++){
		b[i][j]=a[i][j]-a[i-1][j]-a[i][j-1]+a[i-1][j-1];
			cout<<b[i][j]<<" ";
		}
		cout<<endl;
	}
	return 0;
		
	}
	
	a[H][W]=1;
	for(int i=1;i<=h;i++)for(int j=1;j<=w;j++){
		a[i][j]=-1;
		
		for(int k=i;k<=H;k+=h)for(int l=j;l<=W;l+=w){
			if(k==i && l==j) continue;
			a[k][l]=aa(k-h,l)+aa(k,l-w)-aa(k-h,l-w)-1;
			
		}
		
	}
	a[H][W]=1;
	
	
	
	int HH=H,WW=W;
	

	
//	for(HH=H;HH>H-h;HH--)for(WW=W;WW>W-w;WW--){
	for(int i=HH-h;i>0;i-=h)for(int j=WW-w;j>0;j-=w){
		a[i][j]=-a[i+h][j+w]-1+a[i+h][j]+a[i][j+w];
		
	}
//	}
	/*
	for(int i=1;i<=H;i++){
		for(int j=1;j<=W;j++){
		
			cout<<a[i][j]<<" ";
		}
		cout<<endl;
	}
	cout<<endl;
	*/
	cout<<"Yes"<<endl;
	for(int i=1;i<=H;i++){
		for(int j=1;j<=W;j++){
		b[i][j]=a[i][j]-a[i-1][j]-a[i][j-1]+a[i-1][j-1];
			cout<<b[i][j]<<" ";
		}
		cout<<endl;
	}
	
	
	
	
	return 0;
}

Measured memory_bytes: 6912