HDU 2467 Déjà vu IDA*

题意：一个机器内有M个开关，每个开关表示$0$~ $2^n - 1$之间的数字，并且每个每个数字的二进制表达中，有且只有三个1。一个数字Y通过M个开关得到的结果为：
$Y \oplus X_1 \oplus X_2\oplus ...\oplus X_M$
现在给出初始的数字S，和目标的数字T，想用最小的开关个数，使S变成T，求出最小的开关数。如果不可能输出，impossible.
思路：求出最小的步数的方法有很多，图上的最短路，bfs,IDA*,A*dp;
但是题意中提示了：WARNING: a naive algorithm might not be sufficient to solve this problem.这就基本否定了：图上最短路，BFS,DP，原因是这三者的复杂度很高。
这样，看来只能用IDA*和A*
从编程复杂度上，IDA*更简单，所以我们先考虑用IDA*。
在IDA*中，最关键的是对启发式函数h()的设计。
首先，注意到，每个开关至多用一次。因为用两次异或的话，相当于没用。这就明确了，在搜索中我，我们直接向后搜索，无需回溯。
然后，注意到，题目中告诉每个开关的二进制表达中有且只有3个1。
这样，我们可以求出当前数字Z和目标T有多少位不同，作为h函数的自变量。
最重要的内容是对h函数的设计，因为在IDA*中没有判重，必须要有高效的启发式函数h，使搜索更快的到达目标状态。
所以，我们要求出对于有i位不同的两个二进制数字，至少需要多少个开关才能从一个变成另外一个，这就是我们的启发式函数h，这个可以用暴力得到
这样，我们就完成了IDA*的设计。

还有一个比较重要的问题，什么情况下是impossible情况。如果M个开关都搜索，但是没有找到解，这就表明了这样情况下就是impossible。
代码如下：

#include <cstdio>
#include <algorithm>
#include <cstring>

using namespace std;

int N,M,S,T;

int h[]={0,3,2,1,2,3,2,3,4,3,4,5,4,5,6,5,6,7,6,7,8};
int p[1<<21];
int v[1<<21];

int count(int x){
    x = (x & 0x55555555) + ((x >> 1) & 0x55555555);
    x = (x & 0x33333333) + ((x >> 2) & 0x33333333);
    x = (x & 0x0F0F0F0F) + ((x >> 4) & 0x0F0F0F0F);
    x = (x & 0x00FF00FF) + ((x >> 8) & 0x00FF00FF);
    x = (x & 0x0000FFFF) + ((x >> 16) & 0x0000FFFF);
    return x;
}

bool dfs(int state,int id,int d, int maxd)
{
    if(d == maxd && state == T)
        return true;
    if(h[p[state ^ T]] + d > maxd || M - id + d < maxd)
        return false;
    for(int i = id; i < M; ++i)
        if(dfs(state ^ v[i],i+1,d+1,maxd)) return true;
    return false;
}

int main(void)
{
    //freopen("input.txt","r",stdin);
    int cas = 1;
    for(int i = 0; i <1100000; ++i)
        p[i] = count(i);
    while(scanf("%d%d",&N,&M),N||M){
        scanf("%d%d",&S,&T);
        for(int i = 0; i < M; ++i)
            scanf("%d",&v[i]);
        int dep = 0;
        bool sig = false;
        for(;dep <= M;++dep){
            if(dfs(S,0,0,dep)){
                sig = true;
                break;
            }
        }
        printf("Case #%d: ",cas++);
        if(sig)
            printf("%d\n",dep);
        else
            puts("Impossible");
    }
    return 0;
}