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\n", dlx.best);

	}

	return 0;

}

 

你可能感兴趣的:(覆盖)