//#define _GLIBCXX_DEBUG #include "bits/stdc++.h" using namespace std; //------------------------------- Type Names -------------------------------// using i64 = int_fast64_t; using seika = string; //{akari : 1D, yukari : 2D, maki : 3D} vector template using akari = vector; template using yukari = akari>; template using maki = akari>; //{akane : ascending order, aoi : decending order} priority queue template using akane = priority_queue, greater>; template using aoi = priority_queue; //------------------------------- Libraries ---------------------------------// //------------------------------- Dubug Functions ---------------------------// inline void print() { cout << endl; } template void print(const First &first, const Rest &... rest) { cout << first << ' '; print(rest...); } //------------------------------- Solver ------------------------------------// void solve() { i64 n, k; cin >> n >> k; i64 l = n / k, x = n * (n + 1) / 2 / k; if (k * x * 2 != n * (n + 1)) { cout << "No" << endl; return; } yukari ans(k, akari(l)); if (l & 1) { if (l % k != 0) { cout << "No" << endl; return; } int cur = 1, s = 0; for (int i = 0; i < l; i++) { for (int j = 0; j < k; j++) { ans[(s + j) % k][i] = cur++; } s = (s + k - 1) % k; } } else { int cur = 1; for (int i = 0; i < l; i += 2) { for (int j = 0; j < k; j++) { ans[j][i] = cur++; } for (int j = k - 1; j >= 0; j--) { ans[j][i + 1] = cur++; } } } cout << "Yes\n"; for (int i = 0; i < k; i++) { for (int j = 0; j < l; j++) { cout << ans[i][j] << ' '; } cout << '\n'; } } int main() { solve(); return 0; }