这个题dfs的套路跟板子可以说是几乎无相同之处,但他的快是无法想象的。
1.迭代加深 我们寻找最小值,就直接从0开始找,也就是选0包能行吗,选1包能行吗,选2包能行吗 … 选m包能行吗
2.优先选择选择性最低的,我们这里对数据这样处理,一包糖就是一个二进制数,那么我们想凑够m种糖,我们就从含有第一种糖的包里挑一个,但是剩下的几种没选的糖,大概率是包数不一样的,例如 我们现在只有第一种糖和第二种糖没有,第一种糖,50包里有,第二种糖只有1(这里我们说这一包糖叫sniper)包里有,那么其实我们选sniper,就一定是最优解,并且这包里还可能含有第一种糖,可能直接优化掉下一步对50包糖的搜索。(这步优化是根据题目来的)
3.估值函数 估值估值,没必要太准确,大致推推就行,因为如果你弄的太细,会加高时间复杂度,么得必要
这里估值函数的意思就是估计一下当前的这个状态,最少选多少步就能达到最优解,那么也就是说你剩余的步数少于这个的话一定不对了。这里采用的估值逻辑是遍历剩下的糖的种类,让他们含有这些糖的糖包都算一个糖包,然后把这些糖包含的糖的种类都算到原状态里。
(其中第一点与第三点是IDA*的常规操作)
4.除此之外,下面的代码用到了大量二进制运算,以及lowbit优化,lowbit运算得到原二进制数第一位不为0的数,即lowbit(1010)= 0010,我们如果设四位二进制,代表四种糖,为1表示没选,为0表示选过了,那么我们lowbit(i) 就是第一种没选的糖,i - lowbit(i) 就是第二种没选的糖。
感想:怎么说呢,题挺难的,这个算法是专门针对重复覆盖问题的,所以记住也可,然后就是多做题,熟了就行了。这个算法其实第二个优化并不是主优化,主优化应该是第三步的估值函数,这是关键,那么针对不同问题,我们的估值函数到底写成什么样合适,那么其实就是一个对题目的感觉了,所以多做题!
1243. 糖果
糖果店的老板一共有 M 种口味的糖果出售。
为了方便描述,我们将 M 种口味编号 1∼M。
小明希望能品尝到所有口味的糖果。
遗憾的是老板并不单独出售糖果,而是 K 颗一包整包出售。
幸好糖果包装上注明了其中 K 颗糖果的口味,所以小明可以在买之前就知道每包内的糖果口味。
给定 N 包糖果,请你计算小明最少买几包,就可以品尝到所有口味的糖果。
输入格式
第一行包含三个整数 N,M,K。
接下来 N 行每行 K 这整数 T1,T2,⋅⋅⋅,TK,代表一包糖果的口味。
输出格式
一个整数表示答案。
如果小明无法品尝所有口味,输出 −1。
数据范围
1≤N≤100,
1≤M,K≤20,
1≤Ti≤M
输入样例:
6 5 3
1 1 2
1 2 3
1 1 3
2 3 5
5 4 2
5 1 2
输出样例:
2
#include <cstring>
#include <iostream>
#include <algorithm>
#include <vector>
using namespace std;
const int N = 110, M = 1 << 20;
int n, m, k;
vector<int> col[N];
int logg[M];//二进制向糖果种类的转化
int lowbit(int x)
{
return x & -x;
}
int h(int state) // 最少需要再选几行
{
int res = 0;
for (int i = (1 << m) - 1 - state; i; i -= lowbit(i))//logn
{
int c = logg[lowbit(i)];
res ++ ;
for (int row = 0; row < col[c].size(); row++) i &= ~col[c][row];//最优解作为估值
}
return res;
}
bool dfs(int depth, int state)
{
//这里dfs中的depth是倒着来的,depth为3的时候代表选3行达到最优解,但其实depth为3是三行中的第一行,而depth为1是最后一行
if (!depth || h(state) > depth) return state == (1 << m) - 1;//depth=0代表选完了,选完了比较比较返回就完事了
//h(state)估值函数,先预估一下当前这个状态,我再选最少选几步可以
// 找到选择性最少的一列
int t = -1;//特判
//当前可能是缺很多种糖果,但是我们选择,选择性最少的那一种
//意思就是:第一种糖果在50包里有,那第二种糖果只有一包有,那我们选择只有一包的
for (int i = (1 << m) - 1 - state; i; i -= lowbit(i))
{
int c = logg[lowbit(i)];
if (t == -1 || col[t].size() > col[c].size())
t = c;
}
// 枚举选哪行
for(int row = 0; row < col[t].size(); row ++)
if (dfs(depth - 1, state | col[t][row]))
return true;
return false;
}
int main()
{
cin >> n >> m >> k;
for (int i = 0; i < m; i ++ ) logg[1 << i] = i;
for (int i = 0; i < n; i ++ )
{
int state = 0;
for (int j = 0; j < k; j ++ )
{
int c;
cin >> c;
state |= 1 << c - 1;
}
for (int j = 0; j < m; j ++ )
if (state >> j & 1)
col[j].push_back(state);
}
int depth = 0;
while (depth <= m && !dfs(depth, 0)) depth ++ ;
if (depth > m) depth = -1;
cout << depth << endl;
return 0;
}