#include <bits/stdc++.h>
using namespace std;
#define int long long
#define REP(i,n) for(int i = 0;i < (int)(n);i++)
#define RREP(i,n) for(int i = (int)n-1;i >= 0;i--)
#define FOR(i,s,n) for(int i = s;i < (int)n;i++)
#define RFOR(i,s,n) for(int i = (int)n-1;i >= s;i--)
#define ALL(a) a.begin(),a.end()
#define IN(a, x, b) (a<=x && x<b)
template<class T>inline bool CHMAX(T&a,T b){if(a<b){a = b;return true;}return false;}
template<class T>inline bool CHMIN(T&a,T b){if(a>b){a = b;return true;}return false;}
constexpr long long INF = 1e18;

signed main(){
	int N,K;
	cin>>N>>K;
	int sum = N * (N + 1) / 2;
	if(sum % K){
		cout << "No" << endl;
		return 0;
	}
	vector<vector<int>>ans(K);
	if(K % 2){
		if(N == K){
			cout << "No" << endl;
			return 0;
		}
		REP(i,2*K){
			if(i % 2)ans[i % K].push_back(K + i / 2);
			else ans[i % K].push_back(i / 2);
		}
		REP(i,K){
			ans[i].push_back(3 * K - i - 1);
		}
		int turn = N / K - 3;
		REP(i,turn / 2){
			REP(j,K){
				ans[j].push_back(i * 2 * K + 3 * K + j);
				ans[j].push_back(i * 2 * K + 3 * K + 2 * K - j - 1);
			}
		}
	}else{
		int turn = N / K;
		REP(i,turn / 2){
			REP(j,K){
				ans[j].push_back(i * 2 * K + j);
				ans[j].push_back(i * 2 * K + 2 * K - j - 1);
			}
		}
	}
	cout << "Yes" << endl;
	REP(i,K){
		REP(j,ans[i].size()){
			if(j)cout << " ";
			cout << ans[i][j] + 1;
		}
		cout << endl;
	}
}