FZU2165 v11(带权的重复覆盖)

题意:有n个boss,m种武器,每种武器选用的时候需要有一定的花费ci,然后这个武器可以消灭掉其中一些BOSS,问你消灭完所有的BOSS,需要的最少花费是多少。

当时比赛的时候,看到这题以为是什么网络流的题,一种熟悉的感觉,后来才发现,购买一次武器可以消灭掉那么多怪物才不是什么费用流呢。赛后得知这个叫重复覆盖,然后只能用搜的办法,其中搜的比较机智的办法就是用DLX,然后前两天认真的学习了一下DLX,昨天看着别人的代码A了一道题练手。

今天做的时候首先就是将昨天的东西模板化了一下。DLX首先需要初始化,所以有init的函数,然后link,remove,resume就抄昨天的就好了。

今天覆盖的题目不一样就在于它是带权的,其实平时做的也是带权的,权值就是等于深度,这题的权值就不一定等于深度了,所以搜索的时候带个费用就好(其实可以直接将dep改成费用的)。然后搜到比当前最优解不优的情况就return,对于100的数据量就应该不会超时了。跑出来400多ms。

#pragma warning(disable:4996)
#include<iostream>
#include<string>
#include<cstring>
#include<cstdio>
#include<cmath>
#include<vector>
using namespace std;

#define maxnode 15000 // 最多的结点数
#define maxn 150  // 最多的行列数
struct DLX
{
	int L[maxnode], R[maxnode], U[maxnode], D[maxnode];
	int row[maxnode], col[maxnode];
	int S[maxn], H[maxn];
	int size;
	int n;
	int cost[maxn];
	int best;

	void init(int tot)
	{
		n = tot;
		for (int i = 0; i <= n; i++){
			S[i] = 0; U[i] = D[i] = i;
			L[i + 1] = i; R[i] = i+1;
		}R[n] = 0;
		memset(H, -1, sizeof(H)); size = n + 1;
		best = 1000000000;
	}
	void link(int r, int c)
	{
		S[c]++; row[size] = r; col[size] = c;
		U[size] = U[c]; D[U[c]] = size;
		D[size] = c; U[c] = size;

		if (H[r] == -1) H[r] = L[size] = R[size] = size;
		else{
			L[size] = L[H[r]]; R[L[H[r]]] = size;
			R[size] = H[r]; L[H[r]] = size;
		}
		size++;
	}

	void remove(int c)
	{
		for (int i = D[c]; i != c; i = D[i]){
			L[R[i]] = L[i]; R[L[i]] = R[i];
		}
	}

	void resume(int c)
	{
		for (int i = U[c]; i != c; i = U[i]){
			L[R[i]] = R[L[i]] = i;
		}
	}

	void dance(int dep,int curc)
	{
		if (curc> best) return;
		if (R[0] == 0){
			if (curc < best) best = curc; return;
		}
		int minv = maxn; int c;
		for (int i = R[0]; i; i = R[i]){
			if (S[i] < minv) minv = S[i], c = i;
		}
		for (int i = D[c]; i != c; i = D[i]){
			remove(i);
			for (int j = R[i]; j != i; j = R[j]){
				remove(j);
			}
			dance(dep + 1,curc+cost[row[i]]);
			for (int j = L[i]; j != i; j = L[j]){
				resume(j);
			}
			resume(i);
		}
		return;
	}
}dlx;

int n, m;
int vis[maxn];

int main()
{
	while (cin >> n >> m)
	{
		memset(vis, 0, sizeof(vis));
		dlx.init(n); int ci, ki;
		int tcost = 0;
		for (int i = 1; i <= m; i++){
			scanf("%d%d", &ci, &ki); tcost += ci;
			dlx.cost[i] = ci; int ti;
			for (int j = 1; j <= ki; j++){
				scanf("%d", &ti); vis[ti] = true;
				dlx.link(i, ti);
			}
		}
		bool flag = true;
		for (int i = 1; i <= n; i++){
			if (!vis[i]) {
				flag = false; break;
			}
		}
		if (!flag) {
			puts("-1"); continue;
		}
		dlx.best = tcost;
		dlx.dance(0, 0);
		printf("%d
", dlx.best);
	}
	return 0;
}
原文地址:https://www.cnblogs.com/chanme/p/3669576.html