Tromino谜题 ---- 分治法（C++版&OpenCV版）

一：问题描述

Tromino是一个由棋盘上的三个邻接方块组成的L型瓦片。我们的问题是，如何用Tromino覆盖少了一个方块（可以在棋盘上的任何位置）的2^k*2^k棋盘，除了这个缺失的方块，Tromino应该覆盖棋盘上的所有方块。而且不能有重叠。

例如：

如果是2*2的棋盘，（0,0）是缺失的位置，则输出下图：

0 1

1 1

其中1代表L型瓦片。

如果是4*4的棋盘，（0,0）是缺失的位置，则输出下图：

0   2   3   3
2   2   1   3
4   1   1   5
4   4   5   5

二：分析理解

K>0时，可将2^k*2^k的棋盘划分为4个2^(k-1)*2^(k-1)的子棋盘。这样划分后，由于原棋盘只有一个特殊方格，所以，这4个子棋盘中只有1个子棋盘中有特殊方格，其余3个子棋盘中没有特殊方格。为了将这3个没有特殊方格的子棋盘转化成为特殊棋盘，以便采用递归方法求解，可以用一个L型骨牌覆盖这3个较小的棋盘的会合处，从而将原问题转化为4个较小规模的棋盘覆盖问题。递归地使用这种划分策略，直至将棋盘分割为1×1的子棋盘。

三：C++代码与运行图

#define _CRT_SECURE_NO_DEPRECATE 
#define _CRT_SECURE_CPP_OVERLOAD_STANDARD_NAMES 1 

#include<iostream>
#include<cmath>

using namespace std;

int a[1000][1000];//棋盘d
int c;//骨牌的号码

/* tr和tc是棋盘左上角的下标,dr和dc是特殊方格的下标,size是棋盘的大小 */
void Tromino(int tr, int tc, int dr, int dc, int size)
{
	if (size == 1)
		return;

	int t = c++;
	int s = size / 2;

	//左上角
	if (dr < tr + s&&dc < tc + s)
		Tromino(tr, tc, dr, dc, s);
	else
	{
		a[tr + s - 1][tc + s - 1] = t;
		Tromino(tr, tc, tr + s - 1, tc + s - 1, s);
	}

	//右上角
	if (dr < tr + s&&dc >= tc + s)
		Tromino(tr, tc + s, dr, dc, s);
	else
	{
		a[tr + s - 1][tc + s] = t;
		Tromino(tr, tc + s, tr + s - 1, tc + s, s);
	}

	//左下角
	if (dr >= tr + s&&dc < tc + s)
		Tromino(tr + s, tc, dr, dc, s);
	else
	{
		a[tr + s][tc + s - 1] = t;
		Tromino(tr + s, tc, tr + s, tc + s - 1, s);
	}

	//右下角
	if (dr >= tr + s&&dc >= tc + s)
		Tromino(tr + s, tc + s, dr, dc, s);
	else
	{
		a[tr + s][tc + s] = t;
		Tromino(tr + s, tc + s, tr + s, tc + s, s);
	}
}

int main()
{
	int k;
	int x, y;
	printf("输入数据请自行保证k值大于等于1，下标不越界!!!\n");
	printf("请分别输入k值，特殊方格下标，输入3个0表示结束：\n");
	while (scanf("%d%d%d", &k,&x,&y) && (k||x||y))
	{
		//init
		c = 1;
		int num = pow(2, k);
		memset(a, 0, sizeof(a));

		Tromino(0, 0, x, y, num);

		//print
		for (int i = 0; i < num; i++)
		{
			for (int j = 0; j < num; j++)
				printf("%-3d ", a[i][j]);
			printf("\n");
		}

		printf("\n");
	}

	return 0;
}

运行如下：

四：OpenCV代码与运行图

总觉得有些东西应该多实践才能理解，虽然之前学了一点OpenCV，但是一直找不到机会实战，这次正好看到同学用OpenCV实现这题，便来试试

#define _CRT_SECURE_NO_DEPRECATE 
#include<iostream>
#include<cmath>
#include<opencv2/core.hpp>
#include<opencv2/imgproc.hpp>
#include<opencv2/highgui.hpp>

using namespace std;
using namespace cv;

RNG rng(1000);

/* tr和tc是棋盘左上角的下标,dr和dc是特殊方格的下标,size是棋盘的大小 */
void Tromino(int tr, int tc, int dr, int dc, int size, Mat mat)
{
	if (size == 1)
		return;

	int s = size / 2;
	int r = rng.uniform(0, 255);
	int g = rng.uniform(0, 255);
	int b = rng.uniform(0, 255);

	//左上角
	if (dr < tr + s&&dc < tc + s)
		Tromino(tr, tc, dr, dc, s, mat);
	else
	{
		//a[tr + s - 1][tc + s - 1] = t;
		mat.at<Vec3b>(tr + s - 1, tc + s - 1) = Vec3b(r, g, b);
		Tromino(tr, tc, tr + s - 1, tc + s - 1, s, mat);
	}

	//右上角
	if (dr < tr + s&&dc >= tc + s)
		Tromino(tr, tc + s, dr, dc, s, mat);
	else
	{
		//a[tr + s - 1][tc + s] = t;
		mat.at<Vec3b>(tr + s - 1, tc + s) = Vec3b(r, g, b);
		Tromino(tr, tc + s, tr + s - 1, tc + s, s, mat);
	}

	//左下角
	if (dr >= tr + s&&dc < tc + s)
		Tromino(tr + s, tc, dr, dc, s, mat);
	else
	{
		//a[tr + s][tc + s - 1] = t;
		mat.at<Vec3b>(tr + s, tc + s - 1) = Vec3b(r, g, b);
		Tromino(tr + s, tc, tr + s, tc + s - 1, s, mat);
	}

	//右下角
	if (dr >= tr + s&&dc >= tc + s)
		Tromino(tr + s, tc + s, dr, dc, s, mat);
	else
	{
		//a[tr + s][tc + s] = t;
		mat.at<Vec3b>(tr + s, tc + s) = Vec3b(r, g, b);
		Tromino(tr + s, tc + s, tr + s, tc + s, s, mat);
	}
}

int main()
{
	int k;
	int x, y;
	printf("输入数据请自行保证k值大于等于1，下标不越界!!!\n");
	printf("请分别输入k值，特殊方格下标，输入3个0表示结束：\n");

	while (scanf("%d%d%d", &k, &x, &y) && (k || x || y))
	{
		//init
		int num = pow(2, k);
		Mat myMat(num, num, CV_8UC3);
		myMat.at<Vec3b>(x, y) = Vec3b(0, 0, 0);

		Tromino(0, 0, x, y, num, myMat);

		//enlarge the image and show it
		resize(myMat, myMat, Size(512, 512), 0.0, 0.0, 3);
		imshow("OpenCV显示", myMat);
		waitKey(0);
	}
	return 0;
}

运行如下：