# include <bits/stdc++.h>
using namespace std;
int n,s,t,W,H,C,cnt,du[250010],bel[510][510],a[510][510],h[250010],dis[250010],cur[250010],head[250010],to[500010],edge[500010],flow[500010],cap[500010],nxt[500010],tot=1;
vector<int> G[250010];
char s1[510],s2[510];
bool vis[250010];
struct node
{
	int v,w;
	node(int _v=0,int _w=0):v(_v),w(_w){}
	bool operator<(node a)const{return w>a.w;}
};
struct node2
{
	int x,y,dir;// 0 : ·, 1 : | , 2 : —
	node2(int _x=0,int _y=0,int _d=0):x(_x),y(_y),dir(_d){}
}b[250010];
int id(int i,int j){return (i-1)*H+j;}
void addedge(int u,int v,int w,int c)
{
	to[++tot]=v;
	edge[tot]=w;
	cap[tot]=c;
	nxt[tot]=head[u];
	head[u]=tot;
}
void add(int u,int v,int w,int c){addedge(u,v,w,c);addedge(v,u,-w,0);}
bool Dij()
{
	priority_queue<node> que;
	for(int i=1;i<=t;i++) dis[i]=1e9,vis[i]=0;
	dis[s]=0;que.emplace(s,0);
	while(!que.empty())
	{
		node t=que.top();que.pop();
		if(vis[t.v]) continue;
		int u=t.v;vis[u]=1;
		for(int i=head[u];i;i=nxt[i])
		{
			int v=to[i],w=edge[i]+h[u]-h[v];
			
			if(cap[i]>flow[i] && dis[v]>dis[u]+w) dis[v]=dis[u]+w,que.emplace(v,dis[v]);
		}
	}
	return h[t]+dis[t]<0;
}
int dfs(int u,int minn)
{
	if(u==t || !minn) return minn;
	vis[u]=1;
	int f,ans=0;
	for(int &i=cur[u];i;i=nxt[i])
	{
		int v=to[i],w=edge[i]+h[u]-h[v];
		if(!vis[v] && cap[i]>flow[i] && dis[v]==dis[u]+w && (f=dfs(v,min(minn,cap[i]-flow[i])))>0)
		{
			flow[i]+=f;flow[i^1]-=f;
			ans+=f;minn-=f;
		}
		if(!minn) break;
	}
	vis[u]=0;
	return ans;
}
int main()
{
	cin>>n>>C;
	scanf("%s%s",s1+1,s2+1);
	for(int i=1;i<=n;i++) (s1[i]=='1'?W:H)++;
	int x=0,y=0;
	for(int i=1;i<=n;i++)
		if(s1[i]=='1') x++;
		else
		{
			y++;
			for(int i=x+1;i<=W;i++) a[i][y]++;
		}
	x=y=0;
	for(int i=1;i<=n;i++)
		if(s2[i]=='1') x++;
		else
		{
			y++;
			for(int i=x+1;i<=W;i++) a[i][y]--;
		}
	s=W*H+1;t=s+1;
	for(int i=1;i<=W;i++)
		for(int j=1;j<=H;j++)
			if(a[i][j])
			{
				if((i+j)&1) add(s,id(i,j),0,1);
				else add(id(i,j),t,0,1);
			}
	for(int i=1;i<=W;i++)
		for(int j=1;j<=H;j++)
			if(a[i][j] && (i+j)&1)
			{
				if(a[i-1][j]) add(id(i,j),id(i-1,j),-C-2,1);
				if(a[i+1][j]) add(id(i,j),id(i+1,j),-C-2,1);
				if(a[i][j-1]) add(id(i,j),id(i,j-1),-C-1,1);
				if(a[i][j+1]) add(id(i,j),id(i,j+1),-C-1,1);
			}
	while(Dij())
	{
		for(int i=1;i<=t;i++) cur[i]=head[i],vis[i]=0;
		dfs(s,1e9);
		for(int i=1;i<=t;i++)
			if(dis[i]<1e9) h[i]+=dis[i];
	}
	for(int i=1;i<=W;i++)
		for(int j=1;j<=H;j++)
			if(a[i][j] && (i+j)&1)
			{
				int u=id(i,j),v=0;
				for(int k=head[u];k;k=nxt[k])
					if(to[k]<s && flow[k]==cap[k])
					{
						v=to[k];
						break;
					}
				if(v)
				{
					int ti=(v-1)/H+1,tj=(v-1)%H+1,t=(ti==i?2:1);
					b[++cnt]=node2(min(i,ti),min(j,tj),t);
					bel[i][j]=bel[ti][tj]=cnt;
				}
				else b[++cnt]=node2(i,j,0),bel[i][j]=cnt;
			}
	for(int i=1;i<=W;i++)
		for(int j=1;j<=H;j++)
			if(a[i][j] && !((i+j)&1) && !bel[i][j])
				b[++cnt]=node2(i,j,0),bel[i][j]=cnt;
	auto add=[&](int u,int v)
	{
		if(!u || !v || u==v) return;
		G[u].push_back(v);du[v]++;
	};
	auto print=[&](int i)
	{
		int x=b[i].x,y=b[i].y;
		if(!b[i].dir) printf("%d %d\n",x+y-1,x+y);
		else printf("%d %d\n",x+y-1,x+y+1);
	};
	for(int i=1;i<=W;i++)
		for(int j=1;j<=H;j++)
			if(a[i][j]>0) add(bel[i-1][j],bel[i][j]),add(bel[i][j+1],bel[i][j]);
			else add(bel[i][j],bel[i-1][j]),add(bel[i][j],bel[i][j+1]);
	queue<int> que;
	for(int i=1;i<=cnt;i++)
		if(!du[i]) que.push(i);
	cout<<cnt<<endl;
	while(!que.empty())
	{
		int u=que.front();que.pop();
		print(u);
		for(int v:G[u])
			if(!(--du[v])) que.push(v);
	}
	return 0;
}