zoj 1156 Unscrambling Images 解题报告

---

1.题目：

zoj 1156 Unscrambling Images
题目大意：题意可以理解为将一个四分树的各个子女打乱顺序，再从0号开始依次为节点编号，然后表示为n x n的矩阵形式，最后加密。

---

2.解题思路：

这是一道模拟，重点在于理清题意。加密的方法是不变的，也就是说，不管四分树是哪种形式，最终会得到一个从0开始，从左到右，从上到下按递增顺序排列的n x n的intensity矩阵，如下图：

任务是还原出四分树，于是我们首先需要知道如何逆推。

在输入的模块中，第一部分是诸如“k intensity”的形式，利用上图的矩阵，我们可以将intensity矩阵转换为k矩阵，这样节点的编号，节点的顺序就知道了。需要记住的是，在k矩阵中四分树的节点是全部展开的（即使有四个子女的值相同）。
于是输入时，可以用sort对blk(pair)类型按intensity降序排序，然后填到矩阵：

for(int i=0;i<num;i++)
	scanf("%d%d",&b[i].k,&b[i].its);
sort(b,b+num);
//排布在矩阵上
int ind=0;
for(int i=0;i<n;i++)
	for(int j=0;j<n;j++)
		tree[i][j]=b[ind++].k;

第二部分是“k init_intensity”的形式，即给出的intensity是最初的像素强度。由此，我们可以由上述数据将k矩阵转换为init_intensity矩阵，即是最终答案。在这里，并非所有叶子节点的值都被给出，因为这里折叠了带有四个相同子女的节点，而给出父亲节点编号。所以当一个非叶子节点给出了一个值v，则以该节点为根的子树所有节点值都为v。此外，节点的编号有这样的规律：若父亲节点编号为k，则四个子节点编号依次为4k+1,4k+2,4k+3,4k+4。

综上所述，我们先利用简单哈希建立k和intensity的映射，对于未给出的强度，初始化为-1；然后使用dfs递归子树根节点编号root，若root不为-1，对[r,c]到(r+len,c+len)范围的矩阵赋值。若为-1，说明四个儿子不全相同，利用getroot函数根据叶子节点编号倒推子树根，得出递归的子树根的编号，然后进行四分。

void dfs(int root,int r,int c,int len)
{
	if(h[root]!=-1)
	{
		for(int i=r;i<r+len;i++)
			for(int j=c;j<c+len;j++)
				tree[i][j]=h[root];
		return;
	}
	dfs(getroot(root,r,c),r,c,len/2);
	dfs(getroot(root,r,c+len/2),r,c+len/2,len/2);
	dfs(getroot(root,r+len/2,c),r+len/2,c,len/2);
	dfs(getroot(root,r+len/2,c+len/2),r+len/2,c+len/2,len/2);
}

---

3.程序清单：

#include<bits/stdc++.h>
using namespace std;
const int MAX_N=20;
int N,n,cs,num;
int tree[MAX_N][MAX_N];
struct blk
{
	int k,its;
	friend bool operator<(const blk &b1,const blk &b2){return b1.its<b2.its;}
};
blk b[280];
int h[280];
void print()
{
	for(int i=0;i<n;i++)
	{
		for(int j=0;j<n;j++)
			printf("%4d",tree[i][j]);
		printf("\n");
	}
}
int getroot(int root,int r,int c)//获取子树根编号
{
	int t=tree[r][c];
	while(!((t<=4*root+4)&&(t>=4*root+1)))
		t=(t-1)/4;
	return t;
}
void dfs(int root,int r,int c,int len)
{
	if(h[root]!=-1)
	{
		for(int i=r;i<r+len;i++)
			for(int j=c;j<c+len;j++)
				tree[i][j]=h[root];
		return;
	}
	dfs(getroot(root,r,c),r,c,len/2);
	dfs(getroot(root,r,c+len/2),r,c+len/2,len/2);
	dfs(getroot(root,r+len/2,c),r+len/2,c,len/2);
	dfs(getroot(root,r+len/2,c+len/2),r+len/2,c+len/2,len/2);
}
int main()
{
	scanf("%d",&N);
	while(N--)
	{
		scanf("%d",&cs);//输入块数
		for(int csn=1;csn<=cs;csn++)
		{
			memset(h,-1,sizeof(h));
			scanf("%d",&n);//矩阵规模
			scanf("%d",&num);
			for(int i=0;i<num;i++)
				scanf("%d%d",&b[i].k,&b[i].its);
			sort(b,b+num);
			//排布在矩阵上
			int ind=0;
			for(int i=0;i<n;i++)
				for(int j=0;j<n;j++)
					tree[i][j]=b[ind++].k;
			scanf("%d",&num);
			int x,y;
			while(num--)
			{
				scanf("%d%d",&x,&y);
				h[x]=y;
			}
			printf("Case %d\n\n",csn);
			dfs(0,0,0,n);
			print();
			if(csn!=cs)
				printf("\n");
		}
		if(N!=0)
			printf("\n");
	}
	return 0;
}

评测结果：