[ACM] POJ 2133 Cow Imposters, easy DP + BFS

简单DP + BFS
牌子之间两两组合出新牌子, 新产生的牌子只需与原始的号码两两组合即可,即,不用考虑两个新牌子组合的情况,因为两个新牌子组合的结果必然等价与一个新牌子和原始牌子组合的结果, 否则会超时。

很自然的就可以用类似BFS里的队列来做, 我由此引发的思考是, 这种新牌子和原始牌子组合到达新point, 和BFS的思想是十分一致的, 我可以把每一个原始牌子看作一个方向, 牌子号码的每一位是一维空间。 比如 01001 这个原始牌子, 对应了一个五维空间, 这个原始牌子的方向就是 ( 0, 1, 0, 0, 1 ) 。 题目就能抽象为, 给出一个B维的空间和可以走的E个方向, 问你能到达离某个点最近的维度

比如三维空间, 问你到达 (0,0,0) 最低的维度, 如果能到这个点本身那就是0维(0维就是一个点啦 — 我自己想的:),如果能到 (0,0,1) 就是最少1维, 能到 (0,1,1) 是最少二维。 有人问那 (0,0,0) 到 (0,1,1) 不是也能直接连线么, 这样想的话我们要考虑我们只能走给出的原始牌子对应的方向,我们甚至不是顺着xyz轴平移的,而是XOR运算的移动 ^^,所以从 (0,0,0) 到 (0,1,1) 必须顺着两个方向移动, 拉出一个二维空间。 相应的 (0,0,0) 到 (1,1,1) 就要穿越三维空间了。

在B>3的情况下就需要像我一样具有多维的想象力才能理解:)

细节:
变量名写明朗不会死斯基, 把B和E搞混了搞了好久
我输出的时候试图把2进制的每位乘上10来直接放在int里输出,但是没考虑到虽然它实际上是2^16的数,这样转换一下就变成10^16, 不止int的范围了

代码:

1#include <stdio.h>
2#include <stdlib.h>
3#include <algorithm>
4#include <string.h>
5#include <math.h>
6 
7using namespace std;
8 
9int inputBin(int B) {
10    while (getchar() != '\n')
11        ;
12    int temp = 0;
13    for (int i = B - 1; i >= 0; i--)
14        temp += (getchar() - '0') << i;
15    return temp;
16}
17 
18int main() {
19    int B, E;
20    int map[1 << 16];
21    int que[1 << 16];
22    int GOAL;
23 
24    scanf("%d%d", &B, &E);
25    GOAL = inputBin(B);
26 
27    //Initialize
28    for (int i = 0; i < 1 << B; i++)
29        map[i] = -1;
30 
31    for (int i = 0; i < E; i++) {
32        int temp = inputBin(B);
33        map[temp] = 0;
34        que[i] = temp;
35    }
36 
37    //BFS
38    int p0 = 0;
39    int p1 = E - 1;
40 
41    while (p0 <= p1) {
42        if (que[p0] == GOAL)
43            break;
44 
45        int temp = que[p0];
46        for (int i = 0; i < E; i++) {
47            int created = temp ^ que[i];
48            if (map[created] != -1)
49                continue;
50            //if(created==26)
51            //printf("**\n");
52            map[created] = map[temp] + 1;
53            que[++p1] = created;
54        }
55        p0++;
56    }
57 
58    //DetailFucker: Strip needs to be "created" as the text described.
59    int bla = map[0] == 0 ? 1 : 2;
60    for (int i = 0; i < E; i++)
61        map[que[i]] = bla;
62    map[0] = 1;
63 
64    int result[17];
65    int resultA[17];
66    memset(result, 0xff, sizeof(result));
67    for (int i = 0; i <= p1; i++) {
68        int diff = 0;
69        int temp = que[i];
70        for (int k = 0; k < B; k++)
71            diff += ((temp >> k) & 1) ^ ((GOAL >> k) & 1);
72        if (result[diff] == -1 || result[diff] > map[temp] || (result[diff]
73                == map[temp] && resultA[diff] > temp)) {
74            result[diff] = map[temp];
75            resultA[diff] = temp;
76        }
77    }
78 
79    int resultN;
80    for (resultN = 0;; resultN++) {
81        if (resultN > B)
82            exit(1);
83        if (result[resultN] != -1)
84            break;
85    }
86 
87    printf("%d\n", result[resultN]);
88    int notYet = resultA[resultN];
89    for (int i = B - 1; i >= 0; i--)
90        printf("%d", (notYet >> i) & 1);
91    printf("\n");
92 
93    return 0;
94}

Leave a Reply

Your email address will not be published. Required fields are marked *