[知识点]Cantor展开
// 此博文为迁移而来，写于2015年3月14日，不代表本人现在的观点与看法。原始地址：http://blog.sina.com.cn/s/blog_6022c4720102vtyo.html
1、含义
一个很简单的概念哈，其实它的本质就是将你当前状态压缩成一个数，且状态与数一一对应，故一般用在哈希判重，因为有时哈希判重会存不下，或者根本不可能。这是一项辅助的知识点，故不详解。
2、公式
X=a[n]*(n-1)!+a[n-1]*(n-2)!+...+a[i]*(i-1)!+...+a[1]*0，其中,a[i]为整数,并且0<=a[i]<=i<=n。a[i]的意义参见举例中的解释部分。
举例
例如，3 5 7 4 1 2 9 6 8 展开为 98884。因为X=2*8!+3*7!+4*6!+2*5!+0*4!+0*3!+2*2!+0*1!+0*0!=98884.
排列的第一位是3，比3小的数有两个，以这样的数开始的排列有8!个，因此第一项为2*8! 排列的第二位是5，比5小的数有1、2、3、4，由于3已经出现，因此共有3个比5小的数，这样的排列有7!个，因此第二项为3*7!
以此类推，直至0*0!。
3、例题
很好理解的东西我们就直接通过例题来看如何实现。
-------------------------------------------------------------------------------------------------------
数字矩阵（number）
问题描述
还别说ZZD平常毛毛躁躁的，考场上的他可毫不含糊，孤身闯阵，奋勇杀敌，多般神勇，一鼓作气地攻下了第一题，第二题，第三题，真可谓一夫当关，万夫莫开啊！当下正浩浩荡荡地朝第四题前进呢。但是，敌人也不是吃干饭的，经受了一二三题的失败后，他们摸清了一个事实，那就是ZZD率领的银河站队所摆的战阵------银河九宫阵（指横、竖和对角线上的战士战斗力之和都等于15的九宫格，每个格子中的数都在1～9之间，并且不重复）非常厉害，所以在第三题之后他们就对这种战阵进行了超强的光子干扰，也就是使阵中每个战士的战斗资料相互间打乱，这样就让战阵的战斗力大打折扣。难怪刚才刷第四题的时候特别困难，ZZD恍然大悟。而他只会领兵作战，至于修复资料的工作自然就叫给了变成天才的你来解决。所以你只需告诉ZZD资料最少移动的步数（指战斗值上下或左右相交换），如果怎么移动都无法恢复战阵，则输出-1。
输入文件（number.in）
一共有50个战阵给你修复；
每个战阵用一个3*3的数字矩阵表示，每一个数表示一个战士的战斗值；
相邻两个数字距阵不用空行隔开。
输出文件（number.out）
应该有50行，每行对应一个距阵修复的最少步数，如无法修复，则为-1。
输入样例：
1 5 2
6 9 3
8 4 7
1 8 6
4 2 3
9 5 7
5 1 9
4 8 3
2 7 6
输出样例
5
6
6
样例说明
样例就只给出3个供大家参考啦！
-----------------------------------------------------------------------------------------------------
题解：这道题是搜索题，但是不是直接暴搜就搜得到的。不难看出，这道题如果不在BFS的过程中用哈希判重，那么将会爆栈。但是，我们总不可能用九维数组来存hash吧？正如今天的主题一样，在哈希判重时，它需要用到康拓展开。康拓展开为什么能够在这里使用？因为对于栈内每一个矩阵，都有且仅有1-9这9个数字且不重复，所以康拓展开式和矩阵一一对应，可以很方便的判重了。
代码：
#include<cstdio>
#include<cstring>
#define MAXN 100005
const int JC[9]={0,1,2,6,24,120,720,5040,40320};
struct Queue
{
int map[4][4],dep;
};
Queue q[MAXN];
int hash[3628810];
void Copy(int head,int tail)
{
for (int i=1;i<=3;i++)
for (int j=1;j<=3;j++)
q[tail].map[i][j]=q[head].map[i][j];
}
void Swap(int x1,int y1,int x2,int y2,int tail)
{
int temp;
temp=q[tail].map[x1][y1];
q[tail].map[x1][y1]=q[tail].map[x2][y2];
q[tail].map[x2][y2]=temp;
}
int Cantor(int tail)
{
int count[10],ans=0,now=8,tot=0,temp[10];
memset(count,0,sizeof(count));
for (int i=1;i<=3;i++)
for (int j=1;j<=3;j++)
{
tot++;
temp[tot]=q[tail].map[i][j];
}
for (int i=1;i<=8;i++)
{
for (int j=i+1;j<=9;j++)
if (temp[i]>temp[j]) count[now]++;
ans+=JC[now]*count[now];
now--;
}
return ans;
}
int Check(int tail)
{
if (q[tail].map[1][1]+q[tail].map[1][2]+q[tail].map[1][3]!=15) return 0;
if (q[tail].map[2][1]+q[tail].map[2][2]+q[tail].map[2][3]!=15) return 0;
if (q[tail].map[3][1]+q[tail].map[3][2]+q[tail].map[3][3]!=15) return 0;
if (q[tail].map[1][1]+q[tail].map[2][1]+q[tail].map[3][1]!=15) return 0;
if (q[tail].map[1][2]+q[tail].map[2][2]+q[tail].map[3][2]!=15) return 0;
if (q[tail].map[1][3]+q[tail].map[2][3]+q[tail].map[3][3]!=15) return 0;
if (q[tail].map[1][1]+q[tail].map[2][2]+q[tail].map[3][3]!=15) return 0;
if (q[tail].map[3][1]+q[tail].map[2][2]+q[tail].map[1][3]!=15) return 0;
return 1;
}
int BFS()
{
int head=1,tail=2;
while (head!=tail)
{
for (int i=1;i<=2;i++)
for (int j=1;j<=3;j++)
{
Copy(head,tail);
Swap(i,j,i+1,j,tail);
int NowCantor=Cantor(tail);
if (hash[NowCantor]!=1)
{
hash[NowCantor]=1;
q[tail].dep=q[head].dep+1;
if (Check(tail)==1) return q[tail].dep;
tail++;
}
}
for (int i=1;i<=3;i++)
for (int j=1;j<=2;j++)
{
Copy(head,tail);
Swap(i,j,i,j+1,tail);
int NowCantor=Cantor(tail);
if (hash[NowCantor]!=1)
{
hash[NowCantor]=1;
q[tail].dep=q[head].dep+1;
if (Check(tail)==1) return q[tail].dep;
tail++;
}
}
head++;
}
return -1;
}
int main()
{
freopen("number.in","r",stdin);
freopen("number.out","w",stdout);
for (int t=1;t<=50;t++)
{
for (int i=1;i<=3;i++)
for (int j=1;j<=3;j++)
scanf("%d",&q[1].map[i][j]);
printf("%d\n",BFS());
memset(hash,0,sizeof(hash));
}
return 0;
}
