#include <bits/stdc++.h>
using namespace std;
using ll = long long;
#define rep(i, s, t) for (ll i = s; i < (ll)(t); i++)
#define rrep(i, s, t) for(ll i = (ll)(t) - 1; i >= (ll)(s); i--)
#define all(x) begin(x), end(x)
#define rall(x) rbegin(x), rend(x)

#define TT template<typename T>
TT using vec = vector<T>;
TT bool chmin(T &x, T y) { return x > y ? (x = y, true) : false; }
TT bool chmax(T &x, T y) { return x < y ? (x = y, true) : false; }

struct io_setup {
	io_setup() {
		ios::sync_with_stdio(false);
		std::cin.tie(nullptr);
		cout << fixed << setprecision(15);
	}
} io_setup;

int main(){
	int a,b;
	cin>>a>>b;
	int d=gcd(a,b);
	a/=d;
	b/=d;
	set<pair<int,int>> va,vb;
	int l=a*a+b*b;
	{
		int x=0,y=0;
		rep(lp,0,l){
			va.insert({x,y});
			x=(x+a)%l;
			y=(y+b)%l;
		}
		assert(va.size()==l);
	}
	{
		int x=a,y=0;
		rep(lp,0,l){
			vb.insert({x,y});
			x=(x+a)%l;
			y=(y+b)%l;
		}
		assert(vb.size()==l);
	}
	rep(dx,0,d) rep(dy,0,d){
		for(auto[x,y]:va) cout<<x*d+dx*l*d<<" "<<y*d+dy*l*d<<"\n";
	}
	rep(dx,0,d) rep(dy,0,d){
		for(auto[x,y]:vb) cout<<x*d+dx*l*d<<" "<<y*d+dy*l*d<<"\n";
	}
}