#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
#define all(x) (x).begin(),(x).end()
const int mod=100000007,MAX=200005;
const ll INF=1LL<<60;

int main(){
    
    std::ifstream in("text.txt");
    std::cin.rdbuf(in.rdbuf());
    cin.tie(0);
    ios::sync_with_stdio(false);
    
    ll N,K;cin>>N>>K;
    if((1+N)*N/2%K==0){
        if(N/K%2==0){
            cout<<"YES"<<endl;
            for(int i=1;i<=K;i++){
                for(int j=1;j<=N/K;j++){
                    if(j%2==0){
                        cout<<j*K-(i-1)<<" ";
                    }else{
                        cout<<(j-1)*K+i<<" ";
                    }
                }
                cout<<endl;
            }
        }else if(N/K>=K&&(N/K-K)%2==0){
            cout<<"YES"<<endl;
            vector<vector<int>> S(K,vector<int>(N/K,0));
            
            int now=1;
            for(int j=0;j<K;j++){
                for(int i=0;i<K;i++){
                    S[(i+j)%K][j]=now;
                    now++;
                }
            }
            
            for(int j=K;j<N/K;j++){
                if((j-K)%2==0){
                    for(int i=0;i<K;i++){
                        S[i][j]=now;
                        now++;
                    }
                }else{
                    for(int i=0;i<K;i++){
                        S[K-1-i][j]=now;
                        now++;
                    }
                }
            }
            
            for(int i=0;i<K;i++){
                for(int j=0;j<N/K;j++){
                    cout<<S[i][j]<<" ";
                }
                cout<<endl;
            }
        }else if(N/K>=3&&abs(N/K-K)%2==0){
            cout<<"YES"<<endl;
            vector<vector<int>> S(K,vector<int>(N/K,0));
            
            int now=1;
            for(int i=0;i<K;i++){
                S[i][0]=now;
                now++;
            }
            
            for(int i=0;i<K;i++){
                S[(i+K/2)%K][1]=now;
                now++;
            }
            
            for(int i=0;i<K;i++){
                S[i][2]=(3*K+1)/2*3-S[i][0]-S[i][1];
                now++;
            }
            for(int j=3;j<N/K;j++){
                if((j-3)%2==0){
                    for(int i=0;i<K;i++){
                        S[i][j]=now;
                        now++;
                    }
                }else{
                    for(int i=0;i<K;i++){
                        S[K-1-i][j]=now;
                        now++;
                    }
                }
            }
            
            for(int i=0;i<K;i++){
                for(int j=0;j<N/K;j++){
                    cout<<S[i][j]<<" ";
                }
                cout<<endl;
            }
        }else cout<<"NO"<<endl;
    }else cout<<"NO"<<endl;
    
}