HDU 3498 whosyourdaddy DLX重复覆盖
2017-08-20 19:57
525 查看
题目:
http://acm.hdu.edu.cn/showproblem.php?pid=3498题意:
有n个敌人,其中有m对敌人互为邻居,当你攻击杀死一个敌人时,同时会杀死它所有的邻居,问最少需要多少次攻击才能杀死所有敌人思路:
DLX重复覆盖,题很裸,很容易建图,对于第i个敌人,如果它和j互为邻居,则有matrix[i][j]=matrix[j][i]=1,否则为0,然后跑一下重复覆盖求最少次数即可。对比精确覆盖和重复覆盖,二者有些不同,精确覆盖选定一列后,删除此列以及能覆盖此列的行,然后枚举选取某行作为解集一部分的时候,要把覆盖此行的所有列以及覆盖这些列的行删除(读起来有点拗口),因为这些行和当前选取的行已经冲突了,这样保证最终答案每列只被覆盖1次,对于重复覆盖,选定一列后,只删除当前列,然后枚举某行作为解集一部分的时候,只删除覆盖此行的所有列,因为不需要保证每列仅被覆盖1次。两者模板稍有不同。另外重复覆盖每次删除的行变少了,意味着速度变慢了,因此需要强剪枝,用一个估价函数计算在当前状况下至少还需要多少次才能完全覆盖,加上已用的次数,和答案对比,大于等于答案就可以剪枝了
#include <iostream> #include <cstdio> #include <cstring> #include <algorithm> #include <cmath> using namespace std; const int X = 10000 + 10, N = 100 + 10, M = 100 + 10, INF = 0x3f3f3f3f; bool g [M]; int ans; struct DLX { int U[X], D[X], L[X], R[X], row[X], col[X]; int H , S[M]; int head, sz, tot, n, m; bool vis[M]; void init(int _n, int _m) { n = _n, m = _m; for(int i = 0; i <= m; i++) L[i] = i-1, R[i] = i+1, U[i] = D[i] = i, S[i] = 0; head = 0, tot = 0, sz = m; L[head] = m, R[m] = head; for(int i = 1; i <= n; i++) H[i] = -1; } void link(int r, int c) { ++S[col[++sz]=c]; row[sz] = r; D[sz] = D[c], U[D[c]] = sz; U[sz] = c, D[c] = sz; if(H[r] < 0) H[r] = L[sz] = R[sz] = sz; else R[sz] = R[H[r]], L[R[H[r]]] = sz, L[sz] = H[r], R[H[r]] = sz; } void del(int x)//这里的x是点编号,精确覆盖中这里是列号 { for(int i = D[x]; i != x; i = D[i]) R[L[i]] = R[i], L[R[i]] = L[i]; } void recover(int x)//同上 { for(int i = U[x]; i != x; i = U[i]) R[L[i]] = L[R[i]] = i; } int fun_f() //估价函数,估计当前情况下最少还需要多少次才能完成覆盖,只是估计值,只会比真实次数少,从代码中可以明显看出 { memset(vis, 0, sizeof vis); int num = 0; for(int i = R[head]; i != head; i = R[i]) if(! vis[i]) { vis[i] = true; num++; for(int j = D[i]; j != i; j = D[j]) for(int k = R[j]; k != j; k = R[k]) vis[col[k]] = true; } return num; } void dance(int dep) { if(R[head] == head) { ans = min(ans, dep-1); return; } if(dep-1 + fun_f() >= ans) return ;//剪枝 int c = R[head]; for(int i = R[head]; i != head; i = R[i]) if(S[i] < S[c]) c = i; for(int i = D[c]; i != c; i = D[i]) { del(i); for(int j = R[i]; j != i; j = R[j]) del(j); dance(dep + 1); for(int j = L[i]; j != i; j = L[j]) recover(j); recover(i); } } }dlx; int main() { int n, m; while(~ scanf("%d%d", &n, &m)) { memset(g, 0, sizeof g); int a, b; for(int i = 1; i <= m; i++) { scanf("%d%d", &a, &b); g[a][b] = g[b][a] = true; } dlx.init(n, n); for(int i = 1; i <= n; i++) g[i][i] = true; for(int i = 1; i <= n; i++) for(int j = 1; j <= n; j++) if(g[i][j]) dlx.link(i, j); ans = INF; dlx.dance(1); printf("%d\n", ans); } return 0; }
相关文章推荐
- HDU 3498 whosyourdaddy 重复覆盖 DLX+A*
- [DLX重复覆盖] hdu 3498 whosyourdaddy
- hdu - 3498 - whosyourdaddy(重复覆盖DLX)
- hdu 3498 whosyourdaddy (重复覆盖,DLX+迭代加深A*)
- 【HDU】3498 whosyourdaddy 重复覆盖入门题
- hdu 3498 whosyourdaddy(重复覆盖+估价函数剪枝)
- HDU 3498 whosyourdaddy(DLX+A*||多重覆盖)
- HDU 3498 whosyourdaddy (可重复覆盖舞蹈链)
- HDU 3498 whosyourdaddy【Dancing Links重复覆盖】
- HDU 4398 whosyourdaddy 精确覆盖,允许重复覆盖
- hdu 3498(DLX 重复覆盖)
- hdu 3498 whosyourdaddy 重复覆盖+DLX 每次攻击i和i的所有邻居,最少攻击多少人才能杀死所有敌人
- hdu 3498 whosyourdaddy
- HDU 3498 whosyourdaddy
- hdu3498 whosyourdaddy--可重复覆盖舞蹈链
- 【HDU】3498 whosyourdaddy
- HDU 2295 Radar 重复覆盖 DLX
- HDU 5046 Airport(DLX重复覆盖)
- HDU 4735 Little Wish~ lyrical step~(DLX , 重复覆盖)
- hdu 2295 dlx重复覆盖+二分答案