Description
Alice是一个奇怪的画家。她想对一副有N*N个像素点组成的画进行着色,N是2的幂(1,2,4,8,16等等)。每个像素点可以着成黑色或白色。
Alice着色方案不是唯一的,她采用以下不确定的规则:
•如果画作只有一个像素点,那可以直接着白色或黑色;
•否则,把画平均分成四块,然后进行以下操作:
(1) 选择一块全部着白色;
(2) 选择一块全部着黑色;
(3) 把剩下的两块当作是独立的画作并采用同样的方法进行着色。
对于每一幅画作,Alice心目中已经有一个蓝图,接下来请你帮她采用上述方法着色,要求选择跟心目中的蓝图差异最小的着色方案,当然要遵循上述的着色规则,两幅图的差异是指对应位置颜色不相同的像素点的个数。
Input
输入第一行包含整数N(1<=N<=512),表示画作的尺寸为N*N,N保证是2的幂。
接下来N行每行包含N个0或1,描述Alice心目中的蓝图,0表示白色,1表示黑色。
Output
第一行输出最小的差异是多少。
Sample Input
输入1:
4
0001
0001
0011
1110
输入2:
4
1111
1111
1111
1111
输入3:
8
01010001
10100011
01010111
10101111
01010111
10100011
01010001
10100000
Sample Output
输出1:
1
输出2:
6
输出3:
16
Hint
【数据范围】
50%的数据N<=8
思路
DP
设f[i][j][k][0…2]表示以(i,j)为左下角的正方形,边长为2^k,染色状态为w(0…2),与脑补图的不同点数
0代表全白,1代表全黑,2代表按要求自由染色
0和1很容易处理(把四个矩阵加起来)
2根据题意有12种情况。。。
具体看代码
代码
#include<bits/stdc++.h>
const int N=577;
using namespace std;
int f[N][N][13][3],mi[13],n;
int main()
{
memset(f,127,sizeof(f));
scanf("%d\n",&n);
mi[0]=1;
for(int i=1; i<=12; i++)
mi[i]=mi[i-1]*2;
for(int i=1; i<=n; i++)
{
for(int j=1; j<=n; j++)
{
char ch;
scanf("%c",&ch);
int t=ch-48;
f[i][j][0][0]=t;
f[i][j][0][1]=1-t;
f[i][j][0][2]=0;
}
scanf("\n");
}
for(int k=1; k<=log2(n); k++)
{
for(int i=1; i<=n-mi[k]+1; i++) for(int j=1; j<=n-mi[k]+1; j++)
{
f[i][j][k][0]=f[i][j][k-1][0]+f[i+mi[k-1]][j][k-1][0]+f[i][j+mi[k-1]][k-1][0]+f[i+mi[k-1]][j+mi[k-1]][k-1][0];
f[i][j][k][1]=f[i][j][k-1][1]+f[i+mi[k-1]][j][k-1][1]+f[i][j+mi[k-1]][k-1][1]+f[i+mi[k-1]][j+mi[k-1]][k-1][1];
f[i][j][k][2]=min(f[i][j][k][2],f[i][j][k-1][0]+f[i+mi[k-1]][j][k-1][1]+f[i][j+mi[k-1]][k-1][2]+f[i+mi[k-1]][j+mi[k-1]][k-1][2]);
f[i][j][k][2]=min(f[i][j][k][2],f[i][j][k-1][0]+f[i+mi[k-1]][j][k-1][2]+f[i][j+mi[k-1]][k-1][1]+f[i+mi[k-1]][j+mi[k-1]][k-1][2]);
f[i][j][k][2]=min(f[i][j][k][2],f[i][j][k-1][0]+f[i+mi[k-1]][j][k-1][2]+f[i][j+mi[k-1]][k-1][2]+f[i+mi[k-1]][j+mi[k-1]][k-1][1]);
f[i][j][k][2]=min(f[i][j][k][2],f[i][j][k-1][1]+f[i+mi[k-1]][j][k-1][0]+f[i][j+mi[k-1]][k-1][2]+f[i+mi[k-1]][j+mi[k-1]][k-1][2]);
f[i][j][k][2]=min(f[i][j][k][2],f[i][j][k-1][1]+f[i+mi[k-1]][j][k-1][2]+f[i][j+mi[k-1]][k-1][0]+f[i+mi[k-1]][j+mi[k-1]][k-1][2]);
f[i][j][k][2]=min(f[i][j][k][2],f[i][j][k-1][1]+f[i+mi[k-1]][j][k-1][2]+f[i][j+mi[k-1]][k-1][2]+f[i+mi[k-1]][j+mi[k-1]][k-1][0]);
f[i][j][k][2]=min(f[i][j][k][2],f[i][j][k-1][2]+f[i+mi[k-1]][j][k-1][0]+f[i][j+mi[k-1]][k-1][1]+f[i+mi[k-1]][j+mi[k-1]][k-1][2]);
f[i][j][k][2]=min(f[i][j][k][2],f[i][j][k-1][2]+f[i+mi[k-1]][j][k-1][0]+f[i][j+mi[k-1]][k-1][2]+f[i+mi[k-1]][j+mi[k-1]][k-1][1]);
f[i][j][k][2]=min(f[i][j][k][2],f[i][j][k-1][2]+f[i+mi[k-1]][j][k-1][1]+f[i][j+mi[k-1]][k-1][0]+f[i+mi[k-1]][j+mi[k-1]][k-1][2]);
f[i][j][k][2]=min(f[i][j][k][2],f[i][j][k-1][2]+f[i+mi[k-1]][j][k-1][1]+f[i][j+mi[k-1]][k-1][2]+f[i+mi[k-1]][j+mi[k-1]][k-1][0]);
f[i][j][k][2]=min(f[i][j][k][2],f[i][j][k-1][2]+f[i+mi[k-1]][j][k-1][2]+f[i][j+mi[k-1]][k-1][0]+f[i+mi[k-1]][j+mi[k-1]][k-1][1]);
f[i][j][k][2]=min(f[i][j][k][2],f[i][j][k-1][2]+f[i+mi[k-1]][j][k-1][2]+f[i][j+mi[k-1]][k-1][1]+f[i+mi[k-1]][j+mi[k-1]][k-1][0]);
}
}
int pp=log2(n);
printf("%d",f[1][1][pp][2]);
}