CSU 2190 Fractal 模拟 规律

http://acm.csu.edu.cn:20080/csuoj/problemset/problem?pid=2190

Description

人民艺术家Wells的主要工作就是创作分形图，他有一个N*N的模板图，他制作分形的步骤如下：

1：将当前画纸分成N*N小块，在按照模板图染色(白或黑色)；

2：对于每个白色的小块，递归操作#1，即分成更小的N*N块，继续染色，而黑色的小块则不必操作。

你的任务是画出k次递归后的分形图。

Input

第一行两个数N,k，如题意。 接下来是一个NN的模板图，'.'表示白色，''表示黑色。 2≤n≤3,1≤k≤5

Output

一个N^k*N^k的矩阵表示答案。

Sample Input

2 3
.*
..

Sample Output

.*******
..******
.*.*****
....****
.***.***
..**..**
.*.*.*.*
........

Hint

Source

Author

Wells

思路：一开始一直在想递归的做法，反而限制了自己的思路，过了几周发现完全没必要啊，正着模拟就行了。用s[i][j][k]表示第i次分形后第j行第k列的符号，根据题意不难得出：s[i-1][j][k]的符号与第i次的一个矩形映射，我们很容易得到这个矩形表示的范围，行范围为：[n*(j-1)+1，nj]；列范围为：[n*(k-1)+1，nk]；那么我们遍历s[i-1]情况下的所有符号，就可以得到s[i]对应的图形。如果是'*'，映射到第i次的全部是'*'，否则映射到第i次的就是题目给定的模板图。

#include <bits/stdc++.h>
#define INF 0x3f3f3f3f
using namespace std;

char s[10][505][505];
int f[10];
int n,k;

int main()
{
	scanf("%d %d",&n,&k);
	for(int i=1;i<=n;i++)
		scanf("%s",s[1][i]+1);
	f[2]=n;
	for(int i=3;i<=6;i++)
		f[i]=f[i-1]*n;
	for(int i=2;i<=k;i++)
	{
		for(int j=1;j<=f[i];j++)
		{
			for(int k=1;k<=f[i];k++)
			{
				if(s[i-1][j][k]=='*')
				{
					for(int x=n*(j-1)+1;x<=n*j;x++)
						for(int y=n*(k-1)+1;y<=n*k;y++)
							s[i][x][y]='*';
				}
				else
				{
					for(int x=n*(j-1)+1;x<=n*j;x++)
						for(int y=n*(k-1)+1;y<=n*k;y++)
							s[i][x][y]=s[1][x-n*(j-1)][y-n*(k-1)];
				}
			}
		}
	}
	for(int i=1;i<=f[k+1];i++)
	{
		for(int j=1;j<=f[k+1];j++)
			putchar(s[k][i][j]);
		putchar('\n');
	}
	return 0;
}