题意:给定一个方格,然后告诉你一些规则,这些规则下方格中的有些位置是相冲突的。问在方格中最多能够放置多少个士兵。

分析:比赛时一开始就想着用状态压缩DP来搞,不过忘了相邻三行产生关系同样可以通过添加状态的维数来解,于是想着旋转45度之后再DP,这样就只有两行发生关系,不过写起来应该不太好写。后面就用二分匹配写了该题,结果一直WA,知道比赛之后才明白原来这题不能够进行拆点构图,因为从一个点引出去的边的端点相互之间是可以连边的,这就相当于一个成环了,而有环的图又怎么划分出两个内部没有边的二分图呢?其实拆点是用求树的最小路径覆盖的问题上的,这题显然不满足。

正确的解法是状态压缩DP或者是求反建图的最大团问题。DP解法与POJ1185相似。

代码如下:

状态压缩DP

#include <cstdlib>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <iostream>
using namespace std;int N, M, idx;
int sta[175], num[175];
int cnt[15], mp[105];
int f[105][175][175];void print(int x) {int rec[50] = {0}, p = 0;while (x) {rec[p++] = x % 2;x /= 2;}for (int i = 9; i >= 0; --i) {printf(i == 0 ? "%d\n" : "%d", rec[i]);}
}inline bool legal(int x) {if (x & (x << 2)) return false;return true;
}void init() {int lim;idx = 1;for (int k = 1; k <= 10; ++k) {lim = 1 << k;for (int i = 1<<(k-1); i < lim; ++i) {if (legal(i)) {sta[idx] = i;for (int j = 0; j < k; ++j) {if (i & (1 << j)) ++num[idx];}idx++;}}cnt[k] = idx;}
}inline bool place(int x, int y) {return (x & y) == x;
}inline bool judge_1(int x, int y) { // 判断是否斜线冲突 if ((x << 1) & y || (x >> 1) & y) return false;return true;
}inline bool judge_2(int x, int y) { // 判定是否直线冲突 if (x & y) return false;return true;
}int solve() {int ret = 0;memset(f, 0, sizeof (f));if (N == 1) {for (int i = 0; i < cnt[M]; ++i) {if (place(sta[i], mp[1])) {ret = max(ret, num[i]);}}return ret;}for (int i = 0; i < cnt[M]; ++i) {if (!place(sta[i], mp[2])) continue;for (int j = 0; j < cnt[M]; ++j) {if (!place(sta[j], mp[1])) continue;if (judge_1(sta[i], sta[j])) {f[2][i][j] = num[i] + num[j];}}}for (int i = 3; i <= N; ++i) {for (int j = 0; j < cnt[M]; ++j) {if (!place(sta[j], mp[i])) continue;for (int k = 0; k < cnt[M]; ++k) {if (!place(sta[k], mp[i-1]) || !judge_1(sta[j], sta[k])) continue;for (int h = 0; h < cnt[M]; ++h) {if (!place(sta[h], mp[i-2]) || !judge_1(sta[k], sta[h])) continue;if (judge_2(sta[j], sta[h])) {f[i][j][k] = max(f[i][j][k], f[i-1][k][h] + num[j]);    }}}}    }for (int i = 0; i < cnt[M]; ++i) {for (int j = 0; j < cnt[M]; ++j) {ret = max(ret, f[N][i][j]);}}return ret;
}int main() {init();int c;while (scanf("%d %d", &N, &M) != EOF) {memset(mp, 0, sizeof (mp));for (int i = 1; i <= N; ++i) {for (int j = 1; j <= M; ++j) {scanf("%d", &c);mp[i] <<= 1;mp[i] |= c;}}printf("%d\n", solve());}return 0;
}

最大团

#include <iostream>
#include <cstdlib>
#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;int N, M, idx, nd[1005];
char G[1005][1005];void check(int a, int b) {int x1 = nd[a]/M, y1 = nd[a]%M, x2 = nd[b]/M, y2 = nd[b]%M;if (abs(x1-x2) + abs(y1-y2) == 2) {G[a][b] = G[b][a] = 1;}
}void build() { // 对有冲突的点集构一个反图,0视为有边,1视为没有边 memset(G, 0, sizeof (G));for (int i = 0; i < idx; ++i) {for (int j = i+1; j < idx; ++j) {check(i, j);}    }
}int cnt[1005], ret;
int opt[1005], st[1005];void dfs(int x, int num) {int flag;for (int i = x+1; i < idx; ++i) {if (G[x][i]) continue;if (num + cnt[i] <= ret) return; // 可行性剪枝,就算后面完全子图全部加入进来也不会超过已知极大团 flag = true;for (int j = 0; j < num; ++j) {if (G[i][st[j]]) { // 如果没有边flag = false;break;}}if (flag) { // 说明该点能够并入正在搜查的完全子图st[num] = i;dfs(i, num + 1);}}if (num > ret) { // 只会在最里面一层进行更新ret = num;for (int i = 0; i < num; ++i) {opt[i] = st[i];}}
}void maxclique() {ret = 0;for (int i = idx-1; i >= 0; --i) {st[0] = i;dfs(i, 1);cnt[i] = ret;}
}int main() {int c;while (scanf("%d %d", &N, &M) != EOF) {idx = 0;for (int i = 0; i < N; ++i) {for (int j = 0; j < M; ++j) {scanf("%d", &c);if (c) {nd[idx++] = i*M+j;}}}build();maxclique();printf("%d\n", cnt[0]);}return 0;
}

四进制TLE代码留念:

View Code

#include <cstdlib>
#include <cstring>
#include <cstdio>
#include <iostream>
#include <algorithm>
using namespace std;/*
每一个空格使用两个二进制位进行表示:
00 代表当前行和上一行都没有士兵
01 代表当前行有士兵,而上一行没有
10 代表当前行没有士兵,而上一行有
11 代表当前行和上一行都有士兵
*/int N, M, LIM, idx, mp[105];
int sta[4000], num[4000], cnt[15];
const int hmask = 0x55555555; // 该掩码能够将2个二进制位的高位1赋值为0
const int lmask = 0xaaaaaaaa;int f[105][1050];void print(int x) {int rec[50] = {0}, idx = 0;while (x) {rec[idx++] = x % 2;x /= 2;}for (int i = 19; i >= 0; i-=2) {printf(i == 1 ? "%d\n" : "%d", rec[i]);}for (int i = 18; i >= 0; i-=2) {printf(i == 0 ? "%d\n" : "%d", rec[i]);}
}void ppt(int x) {int rec[50] = {0}, idx = 0;while (x) {rec[idx++] = x % 2;x /= 2;}for (int i = 18; i >= 0; i-=2) {printf(i == 0 ? "%d\n" : "%d", rec[i]);}puts("");
}bool judge(int cur, int pre) { // cur代表当前状态,pre代表上一行的综合状态// 首先判定低位的1是否合法int lcur = cur & hmask, hcur = cur & lmask;int lpre = pre & hmask, hpre = pre & lmask;if ((hcur >> 1) != lpre) { // 如果当前状态表示上一层有士兵,那么询问第i-1层的格子是否真有士兵 return false;    }if ((lcur << 1) & hpre) { // 如果当前图不能够满足放置的条件或者第i-2层放置了士兵与之冲突 return false;}if ((lcur << 2) & lpre || (lcur >> 2) & lpre) { // 如果当前行状态和第i-1行放置情况冲突return false;}return true;
}// legal用来判定一行是否合法,其中包括用来表示上一行的高位和低位本行以及互相不能冲突
bool legal(int x) {int lx = x & hmask, hx = x & lmask;if ((x << 4) & x || (lx << 3) & hx || (lx >> 1) & hx) {return false;}return true;
}void solve() {memset(f, 0, sizeof (f));
//    printf("idx = %d\n", idx);for (int i = 1; i <= N; ++i) {for (int j = 0; j < cnt[M]; ++j) { // 求解当前行的状态if (((sta[j]&hmask) & mp[i]) != (sta[j]&hmask)) continue;//    print(sta[j]);//    printf("i = %d, j = %d\n", i, j);//    getchar();for (int k = 0; k < cnt[M]; ++k) { // 遍历上一行的所有状态来生成这一行的状态if (judge(sta[j], sta[k])) {f[i][j] = max(f[i][j], f[i-1][k] + num[j]);/*    if (j == 9) {printf("f[%d][%d] = %d\n", i, j, f[i][j]);getchar();}*/}}}}
}void init() {idx = 1;for (int k = 1; k <= 10; ++k) {LIM = (1 << 2*k) - 1; // 每一个位需要2个二进制位来存储状态for (int i = 1<<2*(k-1); i <= LIM; ++i) {if (legal(i)) {num[idx] = 0;sta[idx] = i;for (int j = 0; j < M; ++j) {if (i & (1 << j*2)) {++num[idx];}}++idx;//    print(i);//    printf("%d\n", num[i]);
            }}cnt[k] = idx;//    printf("idx = %d\n", idx);
    }
}int main() {/*    for (int i = 0; i <= 1024; ++i) {if (!legal(i)) continue;for (int j = 0; j <= 1024; ++j) {    if (!legal(j)) continue;if (judge(i, j)) {print(j), ppt(i);getchar();puts("__________________");}}}*/while (scanf("%d %d", &N, &M) != EOF) {init();int c, Max = 0;memset(mp, 0, sizeof (mp));for (int i = 1; i <= N; ++i) {for (int j = 1; j <= M; ++j) {scanf("%d", &c);mp[i] <<= 2;if (c) mp[i] |= 1;}}//    printf("__%d\n", num[9]);
        solve();for (int i = 0; i < cnt[M]; ++i) {Max = max(Max, f[N][i]);    // 允许第N行以任意合法状态结束,没有更新的状态贡献为0
        }printf("%d\n", Max);}return 0;
}

HDU-4539 郑厂长系列故事——排兵布阵 状态压缩DP Or 最大团相关推荐

  1. HDU 4539 郑厂长系列故事――排兵布阵(状态压缩)

    郑厂长系列故事--排兵布阵 Time Limit: 10000/5000 MS (Java/Others)    Memory Limit: 65535/32768 K (Java/Others) T ...

  2. 郑厂长系列故事——排兵布阵 状态压缩DP

    郑厂长系列故事--排兵布阵 Time Limit: 10000/5000 MS (Java/Others)    Memory Limit: 65535/32768 K (Java/Others) T ...

  3. HDU 4539郑厂长系列故事――排兵布阵(状压DP)

    HDU 4539  郑厂长系列故事――排兵布阵 基础的状压DP,首先记录先每一行可取的所哟状态(一行里互不冲突的大概160个状态), 直接套了一个4重循环居然没超时我就呵呵了 1 //#pragma ...

  4. hdu 4539 郑厂长系列故事——排兵布阵

    郑厂长系列故事--排兵布阵 Time Limit: 10000/5000 MS (Java/Others)    Memory Limit: 65535/32768 K (Java/Others) T ...

  5. HDU 4539 郑厂长系列故事――排兵布阵

    /* 曼哈顿距离的定义是:两个点的坐标为(x1,y1),(x2,y2),两点的曼哈顿距离为|x1-x2|+|y1-y2| 题意:题上要求是两个士兵的距离不能是曼哈顿距离为2,意思就是这个点在同一行同一 ...

  6. 排兵布阵问题java语言_hdu 4539 郑厂长系列故事——排兵布阵

    郑厂长系列故事--排兵布阵 Time Limit: 10000/5000 MS (Java/Others)    Memory Limit: 65535/32768 K (Java/Others) T ...

  7. hdu4539 郑厂长系列故事——排兵布阵 + POJ1158 炮兵阵地

    题意: 郑厂长系列故事--排兵布阵 Time Limit: 10000/5000 MS (Java/Others)    Memory Limit: 65535/32768 K (Java/Other ...

  8. hdu-4539 郑厂长系列故事――排兵布阵(状态压缩)

    类似于poj-1185 //曼哈顿距离 d = |x1 - x2| + |y1 - y2| //解决曼哈顿距离冲突 s[i]&(s[j] << 1)||s[i]&(s[j] ...

  9. HDU4539-郑厂长系列故事——排兵布阵

    郑厂长不是正厂长 也不是副厂长 他根本就不是厂长 事实上 他是带兵打仗的团长 一天,郑厂长带着他的军队来到了一个nm的平原准备布阵. 根据以往的战斗经验,每个士兵可以攻击到并且只能攻击到与之曼哈顿距离 ...

最新文章

  1. GlideApp 引入不了问题
  2. 解决虚拟机 Virtualbox 中 linux 共享文件夹只读问题
  3. 【转】C++ Vector(向量容器)
  4. python 累加器_Python编程第5课:累加器,变量与赋值进阶练习
  5. ASP.NET-后台cookie与前台JQUERY解析cookie
  6. Informix IDS 11系统经管(918考试)认证指南,第 5 部分: 数据库做事器行使(5)
  7. 理解矩阵,矩阵背后的现实意义
  8. .Net Core 商城微服务项目系列(二):使用Ocelot + Consul构建具备服务注册和发现功能的网关...
  9. 前端学习(2733):重读vue电商网站43之使用 lodash 中 cloneDeep(obj) 来实现深拷贝
  10. Error Message ‘Column ows_ does not exist’ when update workflow association on a list
  11. 正点原子STM32F103学习笔记(二)
  12. 时间序列之平稳时间序列预测、趋势型序列预测、复合型序列预测
  13. 从学生宿舍到全球舞台: 米哈游的成长故事
  14. 计算机网络课程思政教学目标,计算机学院举行课程思政示范教学活动
  15. ae渲染文件服务器,在 After Effects 中自动执行渲染和网络渲染
  16. monthCalendar (日历)控件常用操作
  17. Java POI 导入Excell
  18. DeepFlow高效的光流匹配算法(下)
  19. 区块链技术研究热点有哪些
  20. 操作系统-进程调度实验报告

热门文章

  1. acdream 1201 SuSu's Power dp
  2. STM32与中显串口屏的通信
  3. Vue2基础-el与data的两种写法(HTML版)
  4. excel查找出不来了_Excel技巧:明明看到1了,为什么查找不到?
  5. LVM 的自定义报告(pvs vgs lvs各参数详解)
  6. mysql查询当前时间的前后几天时间
  7. 一生不可错过的世界50大美景[组图]
  8. Apache http设置反向代理和负载均衡
  9. 濡沫江湖一直显示获取服务器地址,濡沫江湖霜儿剧情任务怎么完成?霜儿专属剧情任务接取以及完成攻略...
  10. 大数据助力运营商创新转型 中国信息通信大数据大会圆满召开