uva658(最短路径+隐式图+状态压缩)

Posted cryingrain

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了uva658(最短路径+隐式图+状态压缩)相关的知识,希望对你有一定的参考价值。

题目连接(vj):https://vjudge.net/problem/UVA-658

题意:补丁在修正 bug 时,有时也会引入新的 bug。假定有 n(n≤20)个潜在 bug 和 m(m≤100) 个补丁,每个补丁用两个长度为 n 的字符串表示,其中字符串的每个位置表示一个 bug。第一 个串表示打补丁之前的状态(“-” 表示该 bug 必须不存在,“+” 表示必须存在,0 表示无所 谓),第二个串表示打补丁之后的状态(“-” 表示不存在,“+” 表示存在,0 表示不变)。每 个补丁都有一个执行时间,你的任务是用最少的时间把一个所有 bug 都存在的软件通过打补 丁的方式变得没有 bug。一个补丁可以打多次。 在任意时刻,每个 bug 可能存在也可能不存在,所以可以用一个 n 位二进制串表示当前软 件的 “状态”。打完补丁之后,bug 状态会发生改变,对应 “状态转移”。是不是很像动态规 划?可惜动态规划是行不通的,因为状态经过多次转移之后可能会回到以前的状态,即状态 图并不是 DAG。如果直接用记忆化搜索,会出现无限递归。 正确的方法是把状态看成结点,状态转移看成边,转化成图论中的最短路径问题,然后 使用 Dijkstra 或 Bellman-Ford 算法求解。不过这道题和普通的最短路径问题不一样:结点很 多,多达 2 n 个,而且很多状态根本遇不到(即不管怎么打补丁,也不可能打成那个状态), 所以没有必要像前面那样先把图储存好。 还记得第 7 章中介绍的 “隐式图搜索” 吗?这里也可以用相同的方法:当需要得到某个结 点 u 出发的所有边时,不是去读 G[u],而是直接枚举所有 m 个补丁,看看是否能打得上。不管 是 Dijsktra 算法还是 Bellman-Ford 算法,这个方法都适用。本题很经典,强烈建议读者编程实现。 

以上题意以及思路来自紫书

解题思路:抽象成最短路问题,因为状态(节点)太多,而且很多不会出现,所以不能单独建完图再跑最短路,所以采用隐式图,每次扩展遍历一遍已知补丁,匹配成功则看是否可以松弛,记录松弛后的状态。

当前状态唯一确定,所以压缩为一个n位二进制串,0表示无bug,1表示有,而补丁的状态不唯一确定,所以用两个二进制串表示,has串和no串,has串1表示该位置一定有,no串1表示该位置一定无。

匹配方法:cur表示当前状态。如果 cur&no==0并且cur|yes=cur 则匹配成功,匹配成功后松弛后的状态为 v=(u|p.to_has)&(~p.to_no)。然后跑dijkstra即可。

代码如下:

#include<bits/stdc++.h>
#define MAX (1<<20)+10
#define INF 0x3fffffff
using namespace std;

struct data
{
    int from_no,from_has,to_no,to_has,dist;
} patch[105];

struct heapnode
{
    int d,u;
    bool operator < (const heapnode& rhs) const
    {
        return d>rhs.d;
    }
};

int l,m;
bool vis[MAX];
long long d[MAX];

void init(void)
{
    for(int i=0; i<m; i++)
    {
        patch[i].from_has=0;
        patch[i].from_no=0;
        patch[i].to_has=0;
        patch[i].to_no=0;
    }
    memset(vis,0,sizeof(vis));
    for(int i=0; i<MAX; i++)
        d[i]=INF;
}

bool dijkstra(void)
{
    priority_queue<heapnode> Q;
    int s=(1<<l)-1;
    d[s]=0;;
    Q.push((heapnode){0,s});
    while(!Q.empty())
    {
        heapnode x=Q.top();
        Q.pop();
        int u=x.u;
        if(vis[u]) continue;
        vis[u]=1;
        for(int i=0; i<m; i++)
        {
            data &p=patch[i];
            if(!(u&p.from_no)&&((u|p.from_has)==u))//匹配成功
            {
                int v=(u|p.to_has)&(~p.to_no);
                if(d[v]>d[u]+p.dist)
                {
                    d[v]=d[u]+p.dist;
                    Q.push((heapnode){d[v],v});
                }
            }
        }
    }
    if(d[0]==INF)
        return 0;
    return 1;
}

int main()
{
    int T=0;
    while(~scanf("%d%d",&l,&m)&&l)
    {
        init();
        for(int i=0; i<m; i++)
        {
            char from[25],to[25];
            int dis;
            scanf("%d%s%s",&dis,from,to);
            patch[i].dist=dis;
            for(int j=0; j<l; j++)
            {
                if(from[j]==+) patch[i].from_has=patch[i].from_has|(1<<j);
                if(from[j]==-) patch[i].from_no=patch[i].from_no|(1<<j);
                if(to[j]==+) patch[i].to_has=patch[i].to_has|(1<<j);
                if(to[j]==-) patch[i].to_no=patch[i].to_no|(1<<j);
            }
        }
        printf("Product %d\n",++T);
        if(dijkstra())
            printf("Fastest sequence takes %d seconds.\n",d[0]);
        else
            printf("Bugs cannot be fixed.\n");
        printf("\n");
    }
}

 

以上是关于uva658(最短路径+隐式图+状态压缩)的主要内容,如果未能解决你的问题,请参考以下文章

UVa658 It's not a Bug, it's a Feature! (最短路,隐式图搜索)

uva658 It's not a Bug, it's a Feature!

倒水问题 (FillUVa 10603) 隐式图

UVa 658 - It's not a Bug, it's a Feature!(Dijkstra + 隐式图搜索)

路径寻找(隐式图遍历)

UVa10603 Fill (隐式图搜索+Dijkstra)