毒药？解药？[vijos]

原创已于 2022-08-23 22:34:32 修改 · 430 阅读

0 ·

本内容遵循CC 4.0 BY-SA版权协议

标签

#算法 #c++ #宽度优先 #广度优先

于 2022-08-12 14:55:03 首次发布

本文介绍了一道来自Vijos的算法题目，涉及到治疗病症的药剂组合问题。每种药剂对病症有不同的效果，目标是找出治愈所有病症所需的最少药剂数。题目提供病症总数、药剂种类及每种药剂的效果清单，要求通过宽度优先搜索（BFS）或Shortest Path Faster Algorithm（SPFA）来解决。文中还分享了解题者的代码实现。

我在洛谷载的题目
 原题目 ~~（双倍经验）~~

题目背景

羽毛笔和im是抽签到同一个考场的，她们突然闻到一阵刺鼻的化学试剂的气味。

机灵鼠：（头都不抬）你们是考生么？还在门口磨蹭什么？快进来帮我忙！！……怎么还不进来？你们拖赛，拖赛，把你们的青春都拖掉赛……

羽毛笔：哎呀~机灵鼠大人要我们帮什么忙？quq

题目描述

机灵鼠：你们看这里的这些药，都是我研制的对付各种症状的解药。可是我一个不小心，每种药都小小地配错了一点原料，所以这些药都有可能在治愈某些病症的同时又使人患上某些别的病症……（im：那…那是解药还是毒药啊？！）……经过我天才的努力，终于弄清了每种药的具体性能，我会把每种药能治的病症和能使人患上的病症列一张清单给你们，然后你们要根据这张清单找出能治愈所有病症的最少药剂组合……我的药是用不完的，就是说每种药剂都可以被重复使用。给你们的单子里第一行是病症的总数n，第二行是药剂的种类m，以下有m行，每行有n个数字用空格隔开，文件的第i+2行的n个数字中，如果第j个数为1，就表示第i种药可以治愈病症j（如果患有这种病的话则治愈，没有这种病则无影响），如果为0表示无影响，如果为-1表示反而能使人得上这种病（无病患上，有病无影响）。我制的药任何两种性能都不同。你们只要给我用的最少的药剂数就可以了。其实还有可能用尽了所有的药也不能将所有病治愈,那就输出：“The patient will be dead.”

输入格式

第一行是病症的总数n，病症的数目不超过10种;

第二行是药剂的种类m，0< m<=100.

以下有m行，每行有n个数字用空格隔开，文件的第i+2行的n个数字中，如果第j个数为1，就表示第i种药可以治愈病症j。

输出格式

最少的药剂数 (不可能将所有病治愈,输出：The patient will be dead.)

样例 #1

样例输入 #1

样例输出 #1

本蒟蒻的解法

理论上来说状压bfs比较好，但是可以spfa。
代码奉上：

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cstring>
#include <queue>
using namespace std;

const int MAXN=12;
const int MAXM=105;
const int MAX=(1<<10)+5;
const int INF=(1<<30)-1;
queue<int> q;
bool in[MAX];   int d[MAX];//SPFA相关
int a[MAXM];//记录治愈药性
int b[MAXM];//记录毒性
int n,m;
int s;//初始状态

void init()
{
    int x;
    cin>>n>>m;
    for(int i=1;i<=m;i++)
        for(int j=0;j<n;j++)
        {
            cin>>x;
            if(x==1)
                a[i]|=(1<<j);
            else    if(x==-1)
                b[i]|=(1<<j);
        }
    s=(1<<n)-1; //初始时所有患病
    for(int i=0;i<=(1<<n)-1;i++)
        d[i]=INF;
}

inline int update(int x,int k)
{
    x&=(~a[k]);//a[k]取反之后不能治愈的变为了1,只有两个都为1才会患病
    x|=b[k];//有一个为1则这种病会患上
    return x;
}

void SPFA()
{
    d[s]=0; q.push(s);  in[s]=1;
    while(!q.empty())
    {
        int u=q.front();    in[u]=0;    q.pop();
        for(int i=1;i<=m;i++)
        {
            int v=update(u,i);
            if(d[v]>d[u]+1)
            {
                d[v]=d[u]+1;
                if(!in[v])
                    q.push(v),in[v]=1;
            }
        }
    }
    if(d[0]==INF)
        cout<<"The patient will be dead."<<endl;
    else
        cout<<d[0]<<endl;
}

int main()
{
    init();
    SPFA();
    return 0;
}