Problem: People do many crazy things to stand out in a crowd. Some of them dance, some learn by heart rules of Russian language, some try to become an outstanding competitive programmers, while others collect funny math objects.

Alis is among these collectors. Right now she wants to get one of k-special tables. In case you forget, the table n × n is called k-special if the following three conditions are satisfied:

  * every integer from 1 to n2 appears in the table exactly once; 
  * in each row numbers are situated in increasing order; 
  * the sum of numbers in the k-th column is maximum possible. 



Your goal is to help Alice and find at least one k-special table of size n × n. Both rows and columns are numbered from 1 to n, with rows numbered from top to bottom and columns numbered from left to right.

Input

The first line of the input contains two integers n and k (1 ≤ n ≤ 500, 1 ≤ k ≤ n) — the size of the table Alice is looking for and the column that should have maximum possible sum.

Output

First print the sum of the integers in the k-th column of the required table.

Next n lines should contain the description of the table itself: first line should contains n elements of the first row, second line should contain n elements of the second row and so on.

If there are multiple suitable table, you are allowed to print any.

Examples

Input

4 1


Output

28
1 2 3 4
5 6 7 8
9 10 11 12
13 14 15 16


Input

5 3


Output

85
5 6 17 18 19
9 10 23 24 25
7 8 20 21 22
3 4 14 15 16
1 2 11 12 13

Solution (in CPP):

Solution:
#include <bits/stdc++.h>
using namespace std;
template <class T>
inline T bigmod(T p, T e, T M) {
  if (e == 0) return 1;
  if (e % 2 == 0) {
    T t = bigmod(p, e / 2, M);
    return (t * t) % M;
  }
  return (bigmod(p, e - 1, M) * p) % M;
}
template <class T>
inline T gcd(T a, T b) {
  if (b == 0) return a;
  return gcd(b, a % b);
}
template <class T>
inline T modinverse(T a, T M) {
  return bigmod(a, M - 2, M);
}
template <class T>
inline T multimod(T a, T b, T m) {
  T ans = 0ll;
  a %= m, b %= m;
  while (b) {
    if (b & 1ll) ans = m - ans > a ? (ans + a) : (ans + a - m);
    b >>= 1ll;
    a = (m - a) > a ? a + a : a + a - m;
  }
  return (T)ans;
}
const int maxBufSize = (6000000);
struct Input {
  int bufSize, bufEnd, bufPos;
  char buffer[maxBufSize];
  void grabBuffer() {
    bufSize = (maxBufSize);
    bufPos = 0;
    bufEnd = fread(buffer, sizeof(char), bufSize, stdin);
    buffer[bufEnd] = '\0';
  }
  bool bufEof() { return bufPos == bufEnd; }
  int getChar() { return buffer[bufPos++]; }
  void skipWS() {
    while ((buffer[bufPos] == '\n' || buffer[bufPos] == ' ' ||
            buffer[bufPos] == '\t'))
      bufPos++;
  }
  int rInt() {
    int n = 0, x;
    skipWS();
    bool neg = false;
    if ((x = getChar()) == '-') {
      neg = true;
      x = getChar();
    }
    for (; x <= '9' && x >= '0'; x = getChar())
      n = (n << 3) + (n << 1) + (x - '0');
    if (neg)
      return -n;
    else
      return n;
  }
  inline bool isWhiteSpace(char x) {
    return x == ' ' || x == '\n' || x == '\t';
  }
  char rChar() {
    skipWS();
    return getChar();
  }
} input;
long long arr[509][509];
int main() {
  int a, b;
  int n, e, i, u, v;
  string str;
  int num;
  cin >> n >> num;
  int rest = (num - 1) * n;
  long long sum = 0ll;
  int row = 1;
  for (i = rest + 1; i <= n * n; i += (n - num + 1)) {
    sum += (long long)i;
    arr[row][num] = i;
    for (int j = num + 1; j <= n; j++) {
      arr[row][j] = arr[row][j - 1] + 1;
    }
    row++;
  }
  cout << sum << endl;
  int have = 0;
  for (i = 1; i <= n; i++) {
    for (int j = 1; j <= n; j++) {
      if (!arr[i][j])
        cout << (++have) << " ";
      else
        cout << arr[i][j] << " ";
    }
    cout << endl;
  }
  return 0;
}