#include <bits/stdc++.h>
using namespace std;
#define pb push_back
typedef long long ll;
#define rep(i,n) for(int i=0;i<(int)(n);i++)
#define rep1(i,n) for(int i=1;i<=(int)(n);i++)
const int n_max = 1e5+10;
vector<int> g[n_max];
// vector<int> color(n_max,-1);
// bool ok = true;
// void dfs(int cv, int cc){
//   color[cv] = cc;
//   int nc = (cc+1)%2;
//   for(auto i:g[cv]){
//     if(color[i]>=0){
//       // cout << i <<" " << color[i] << endl;
//       if(color[i] == nc) continue;
//       else ok = false;
//     }
//     else{
//       color[i] = nc;
//       dfs(i,nc);
//     }
//   }
// }

int main(){
  ios::sync_with_stdio(false);
  cin.tie(0);

  int n,m;
  cin >> n >> m;
  int poss_m = n/2 * ((n+1)/2);
  if(poss_m < m || m < n-1){
    cout << "NO" << endl;
    return 0;
  }
  cout << "YES" << endl;
  rep(i,(n+1)/2){
    if(i+(n+1)/2 < n){
      g[i].pb(i+(n+1)/2);
      g[i+(n+1)/2].pb(i);
    }
    if(i+(n+1)/2-1>=(n+1)/2){
      g[i].pb(i+(n+1)/2-1);
      g[i+(n+1)/2-1].pb(i);
    }
  }
  m-=n-1;
  rep(i,(n+1)/2){
    for(int j =(n+1)/2;j<n;j++){
      if(m==0) break;      
      if(j==i+(n+1)/2 || j==i+(n+1)/2-1) continue;
      else{
        m--;
        g[i].pb(j);
        g[j].pb(i);
      }      
    }
    if(m==0) break;
  }
//  dfs(0,1);
//  assert(ok);
  rep(i,n){
    cout << i+1 << " \n"[i==n-1];
  }
  rep(i,n){
    for(auto j:g[i]){
      if(i<j) cout << i+1 << " " << j+1 << endl;
    }
  }
  
  return 0;
    

}