#include <stdio.h>
#include <string>
#include <cstring>
#include <stdlib.h>
#include <math.h>
#include <algorithm>
#include <vector>
#include <set>
#include <map>
#include <queue>
#include <stack>
#include <list>
#include <iterator>
#include <assert.h>
#pragma warning(disable:4996) 
 
typedef long long ll;
#define MIN(a, b) ((a)>(b)? (b): (a))
#define MAX(a, b) ((a)<(b)? (b): (a))
#define LINF 9223300000000000000
#define INF 2140000000
const long long MOD = 1000000007;
//const long long MOD = 998244353;
using namespace std;


int main(int argc, char* argv[])
{
    int n,K;
    scanf("%d%d", &n, &K);

    if(K==1) {
        printf("Yes\n");
        int i;
        for(i=0; i<n; i++) {
            printf("%d ", i+1);
        }
        printf("\n");
        return 0;
    }

    int m=n/K;
    if(K%2==0 && m%2) {
        printf("No\n"); return 0;
    }
    else if(K%2==0 || m%2==0) {
        printf("Yes\n");
        vector<vector<int> > z(K);
        int i,j;
        int cnt=0;
        for(i=0; i<m; i++) {
            for(j=0; j<K; j++) {
                int j2=(i%2==0? j: K-1-j);
                z[j2].push_back(cnt+1);
                cnt++;
            }
        }
        for(i=0; i<K; i++) {
            for(j=0; j<m; j++) {
                printf("%d ", z[i][j]);
            }
            printf("\n");
        }
    }
    else {
        if(m<=2) {
            printf("No\n"); return 0;
        }
        else {
            printf("Yes\n");
            vector<vector<int> > z(K);
            int i,j;
            int cnt=0;
            for(i=0; i<m; i++) {
                if(i<m-2) {
                    for(j=0; j<K; j++) {
                        int j2=(i%2==0? j: K-1-j);
                        z[j2].push_back(cnt+1);
                        cnt++;
                    }
                }
                else if(i==m-2) {
                    int K0=(K-1)/2;
                    for(j=0; j<K; j++) {
                        int j1;
                        if(j<K0) j1=K0-1-j;
                        else j1=K0*3-j;
                        int j2=(i%2==0? j1: K-1-j1);
                        z[j2].push_back(cnt+1);
                        cnt++;
                    }
                }
                else {
                    vector<pair<int,int> > v;
                    int p;
                    for(p=0; p<K; p++) {
                        int tmp=z[p][m-3]+z[p][m-2];
                        v.push_back(make_pair(tmp,p));
                    }
                    sort(v.rbegin(), v.rend());
                    for(j=0; j<K; j++) {
                        z[v[j].second].push_back(cnt+1);
                        cnt++;
                    }
                }
            }
            for(i=0; i<K; i++) {
                for(j=0; j<m; j++) {
                    printf("%d ", z[i][j]);
                }
                printf("\n");
            }
        }
    }


    return 0;
}