It‘s not a Bug, it‘s a Feature! UVA - 658(dij + 隐式搜索)

原创于 2020-11-10 11:36:27 发布 · 212 阅读

0 ·

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

最短路专栏收录该内容

17 篇文章

订阅专栏

本文介绍了一个关于软件修复中寻找最优补丁应用顺序的问题。通过隐式图搜索算法，结合二进制操作技巧，在限定的补丁数量下找到修复所有bug所需的最短时间。算法采用Dijkstra实现，并详细展示了如何通过节点状态转移来枚举所有可能的补丁组合。

It’s not a Bug, it’s a Feature!

题意：

给你一个软件有n个bug。（n<=20)
现在提供补丁：（m<=100）

每个补丁有一个初态要求。(‘-’，表示该bug必须不存在，‘+’表示必须存在， 0 表示无所谓)。
打完补丁后，有一个终态。（‘-’表示不存在，’+'表示存在，0表示不变）
每个补丁都有个执行时间。

思路：

可以知道，一个软件的可能状态数是 $2^n$ 种，假如直接暴力。可能会TLE
这里引进一种方法（紫薯介绍的）

`隐式图搜索`

当需得到某个结点u出发的所有边时，不是去读G【u】，而是直接枚举所有m个补丁，看看是否能打得上。不管是 $D i j$ 还是 $B e l l m a n - F o r d$ .

反思：

本题也考察了一些二进制的操作。

对于要获得u里的1位数，先移位，再&1.
对于要判断u里是否有1位数， &（1<<i）.

AC

#include <iostream>
#include <cstring>
#include <vector>
#include <algorithm>
#include <queue>
#define mst(x,a) memset(x,a,sizeof(x))
#define For(i,x,y) for(int i=(x); i<=(y); i++)
#define fori(i,x,y) for(int i=(x); i<(y); i++)
#define fzhead EDGE(int _to, int _val, int _next)
#define fzbody to(_to), val(_val), next(_next)
#define mp make_pair
#define fi first
#define se second
#define pb push_back
using namespace std;
const int maxn=(1<<20)+5;
const int maxm=1e6+10;
typedef long long ll;
typedef pair<ll,ll>pai;
typedef pair<int,int>pa;
const int inf=0x3f3f3f3f;
struct EDGE{
    int to, val, next;
    EDGE(){}
    fzhead:fzbody{}
}e[maxm];
int head[maxn],tot,vis[maxn],dis[maxn],wei[200];
void add(int bg, int to, int val){
    e[tot]=EDGE(to, val, head[bg]);
    head[bg]=tot++;
}
int m,n;
string patch[210];
bool is_ok(int num, int u, int i){
    if(patch[num][i]=='0')return true;
    else if(patch[num][i]=='-'&& (u>>i)&1)return false;
    else if(patch[num][i]=='+'&& !((u>>i)&1))return false;
    else return true;
}
int modfy(int num, int u, int i){
    int res=(u>>i)&1;
    if(patch[num+m][i]=='0')return (res<<i);
    else if(patch[num+m][i]=='+')return (1<<i);
    else return 0;
}
void print(int u){
    fori(i,0,n){
        if(u&(1<<i))cout<<1;
        else cout<<0;
    }
}
inline void dij(int s){
    priority_queue<pa>q;
    dis[s]=0;q.push(mp(0,s));
    while(!q.empty()){
        int u=q.top().se;
        q.pop();
        if(vis[u])continue;
        //cout<<"u-----";print(u);cout<<endl;
        vis[u]=1;
        fori(num,0,m){
            int flag=1;
            fori(i,0,n){
                if(!is_ok(num,u,i)){
                    flag=0;
                    break;
                }
            }
            if(!flag)continue;
            int v=0;
           // cout<<"pre---v------";print(v);
           // cout<<endl;
            fori(i,0,n)v|=modfy(num,u,i);
            //cout<<"now---v------";print(v);
           // cout<<endl;
            if(dis[v]>dis[u]+wei[num]){
                dis[v]=dis[u]+wei[num];
                q.push(mp(-dis[v],v));
            }
        }
    }
}
void init(){
    mst(vis,0);mst(dis,0x3f);
    mst(head,-1);tot=0;
    fori(i,0,m){
        cin>>wei[i];///
        cin>>patch[i]>>patch[i+m];
    }
}
int main()
{
    int kase=0;
    ios::sync_with_stdio(0);cin.tie(0);cout.tie(0);
    while(cin>>n>>m,n||m){
        //if(kase)cout<<endl;
        init();
        int s=0;
        fori(i,0,n)s|=1<<i;
        dij(s);
        cout<<"Product "<<++kase<<endl;
        if(dis[0]==inf)cout<<"Bugs cannot be fixed."<<endl;
        else cout<<"Fastest sequence takes "<<dis[0]<<" seconds."<<endl;
        cout<<endl;
    }
    return 0;
}