#include <algorithm>
#include <climits>
#include <cmath>
#include <cstdio>
#include <cstdlib>
#include <ctime>
#include <iostream>
#include <sstream>
#include <functional>
#include <map>
#include <string>
#include <cstring>
#include <vector>
#include <queue>
#include <stack>
#include <deque>
#include <set>
#include <list>
#include <numeric>
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
typedef pair<ll,ll> P;
const double PI = 3.14159265358979323846;
const double EPS = 1e-12;
const ll INF = 1LL<<29;
const ll mod = 1e9+7;
#define rep(i,n) for(int (i)=0;(i)<(ll)(n);++(i))
#define repd(i,n,d) for(ll (i)=0;(i)<(ll)(n);(i)+=(d))
#define all(v) (v).begin(), (v).end()
#define pb(x) push_back(x)
#define mp(x,y) make_pair((x),(y))
#define mset(m,v) memset((m),(v),sizeof(m))
#define chmin(x,y) (x=min(x,y))
#define chmax(x,y) (x=max(x,y))
#define fst first
#define snd second
#define UNIQUE(x) (x).erase(unique(all(x)),(x).end())
template<class T> ostream &operator<<(ostream &os, const vector<T> &v){int n=v.size();rep(i,n)os<<v[i]<<(i==n-1?"":" ");return os;}

int main(){
	ll n, m;
	cin>>n>>m;
	ll n2 = n/2, n3 = n-n2;
	if(m<n-1 || n2*n3<m){
		cout<<"NO"<<endl;
		return 0;
	}
	cout<<"YES"<<endl;
	vector<ll> a(n);
	rep(i, n){
		a[i] = i%2?1e5-i:i;
	}
	cout<<a<<endl;
	rep(i, n-1) cout<<i+1<<" "<<i+2<<endl;
	m -= n-1;
	rep(i, n) rep(j, i-2){
		if(i%2==j%2) continue;
		if(m--){
			cout<<i+1<<" "<<j+1<<endl;
		} else return 0;
	}

	return 0;
}