hdu 4021(16数码)

         题意：给出一个board，上面有24个位置，其中23个位置上放置了标有数字1~23的方块，一个为空位（用数字0表示），现在可以把空位与它旁边的方块交换，给出board的起始状态，问是否可以达到指定的状态。


       思路：看起来很像著名的“八数码”问题，首先，针对八个特殊位置（死角），如果这里有空位就把它和相邻的位置交换，这样之后如果两个状态的对应死角上的数字不同，那么显然是不能达到指定状态的，因为无法把死角处的数字换出去。


      搞定了死角后就只剩下4×4的board了，这就变成了八数码问题的拓展——15数码。首先想想八数码是如何判断有解的：首先把所有数字（不包括空位的0）写成一行，就得到了一个1~8的排列，考虑空位的交换情况：1.左右交换，2.上下交换。对于左右交换而言，是不会改变写出的排列的逆序数的；而对上下交换，相当于在排列中向前或向后跳了两个位置，那么要么两个数都比它大或小，这样逆序数加2或减2，要么两个数一个比它大一个比它小，这样逆序数不变，综上，对于八数码问题，操作不会改变逆序数的奇偶性，所以只有初始状态和指定状态的逆序数奇偶性相同就有解。


      弄清楚了八数码，扩展起来就容易了，现在我们将其扩展到N维（即N*N的board，N*N-1数码问题）。


      首先无论N的奇偶，左右交换不改变逆序数，N为奇数时：上下交换逆序数增加N-1或减少N-1或不变，因为N为奇数，所以逆序数奇偶性不变；而N为偶数时：上下交换一次奇偶性改变一次。


     结论：N为奇数时，初始状态与指定状态逆序数奇偶性相同即有解；N为偶数时，先计算出从初始状态到指定状态，空位要移动的行数m，如果初始状态的逆序数加上m与指定状态的逆序数奇偶性相同，则有解。

好了，现在这道题就简单了，计算逆序数和空格要移动的行数即可。

代码如下：

#include <cstdio>
#include <stack>
#include <set>
#include <iostream>
#include <string>
#include <vector>
#include <queue>
#include <functional>
#include <cstring>
#include <algorithm>
#include <cctype>
#include <string>
#include <map>
#include <iomanip>
#include <cmath>
#define LL long long
#define ULL unsigned long long
#define SZ(x) (int)x.size()
#define Lowbit(x) ((x) & (-x))
#define MP(a, b) make_pair(a, b)
#define MS(arr, num) memset(arr, num, sizeof(arr))
#define PB push_back
#define F first
#define S second
#define ROP freopen("input.txt", "r", stdin);
#define MID(a, b) (a + ((b - a) >> 1))
#define LC rt << 1, l, mid
#define RC rt << 1|1, mid + 1, r
#define LRT rt << 1
#define RRT rt << 1|1
#define BitCount(x) __builtin_popcount(x)
#define BitCountll(x) __builtin_popcountll(x)
#define LeftPos(x) 32 - __builtin_clz(x) - 1
#define LeftPosll(x) 64 - __builtin_clzll(x) - 1
const double PI = acos(-1.0);
const int INF = 0x3f3f3f3f;
using namespace std;
const double eps = 1e-8;
const int MAXN = 300 + 10;
const int MOD = 1000007;
const int M=20010;
const int N=25;
typedef pair<int, int> pii;
int a[N],b[N],c[N],d[N];
void slove()
{
    int i,j;
    bool flag=true;
    if(a[0]==0)swap(a[0],a[3]);
    if(a[1]==0)swap(a[1],a[6]);
    if(a[2]==0)swap(a[2],a[3]);
    if(a[7]==0)swap(a[7],a[6]);
    if(a[16]==0)swap(a[16],a[17]);
    if(a[21]==0)swap(a[20],a[21]);
    if(a[22]==0)swap(a[17],a[22]);
    if(a[23]==0)swap(a[20],a[23]);

    if(b[0]==0)swap(b[0],b[3]);
    if(b[1]==0)swap(b[1],b[6]);
    if(b[2]==0)swap(b[2],b[3]);
    if(b[7]==0)swap(b[7],b[6]);
    if(b[16]==0)swap(b[16],b[17]);
    if(b[21]==0)swap(b[20],b[21]);
    if(b[22]==0)swap(b[17],b[22]);
    if(b[23]==0)swap(b[20],b[23]);

    if(a[0]!=b[0])flag=false;
    if(a[1]!=b[1])flag=false;
    if(a[2]!=b[2])flag=false;
    if(a[7]!=b[7])flag=false;
    if(a[16]!=b[16])flag=false;
    if(a[21]!=b[21])flag=false;
    if(a[22]!=b[22])flag=false;
    if(a[23]!=b[23])flag=false;

    if (!flag) {
        puts("Y");
        return;
    }
    for (i=0;i<4;i++){
        c[i]=a[i+3];
        d[i]=b[i+3];
    }
    for (i=4;i<12;i++){
        c[i]=a[i+4];
        d[i]=b[i+4];
    }
    for (i=12;i<16;i++){
        c[i]=a[i+5];
        d[i]=b[i+5];
    }
    int m1,m2,m,cnt1,cnt2;
    m1=m2=m=cnt1=cnt2=0;
    for (i=0;i<16;i++){
        if (c[i]==0) {
            m1=i;
            continue;
        }
        for (j=0;j<i;j++) {
            if (c[j]>c[i]) cnt1++;
        }
    }
    for (i=0;i<16;i++){
        if (d[i]==0) {
            m2=i;
            continue;
        }
        for (j=0;j<i;j++) {
            if (d[j]>d[i]) cnt2++;
        }
    }
    m1/=4;
    m2/=4;
    m=abs(m1-m2);
    if(((cnt1+m)&1)==((cnt2&1))) puts("N");
    else puts("Y");
}
int main()
{
    int t;
    cin>>t;
    while(t--)
    {
        for (int i=0;i<24;i++) scanf("%d",a+i);
        for (int i=0;i<24;i++) scanf("%d",b+i);
        slove();
    }
}