hdu 4640 多校第四场 状态压缩+最短路+dp

来源:互联网 发布:淘宝我的仓库在哪里 编辑:程序博客网 时间:2024/05/21 04:02

应该注意到n很小,说明我们可以位压缩处理每种状态

预处理出一个人从1出发,已经走到哪些点(用状态记录),最后到达哪个点的最小花费,用广搜spfa

之后i个人走,是一个人的效果的i次叠加,用类似背包的过程,注意要从大到小枚举状态,避免转移重复,dp

标程错了,它处理成不能回到1,我wa了很久,主要手写队列搞挫了。。。

开始时不会做是因为看不到问题的本质,本质就是搞出一个人时候的最短路,然后dp叠加出来,中间状态压缩只要平时熟悉就不难的

#include <iostream>#include <cstdio>#include <cstring>using namespace std;#define INF 2000000000int q[4*(1<<17)][2];int maxQ=4*(1<<17)-3;int front,rear,n,m;int mat[20][20];int dp[1<<17],f[1<<17],dist[1<<17][18],vis[1<<17][18];int main (){    //freopen("input.txt","r",stdin);    //freopen("output.txt","w",stdout);    int T;scanf("%d",&T);    for(int kk=1;kk<=T;++kk)    {        scanf("%d%d",&n,&m);        memset(dist,127,sizeof(dist));        memset(dp,127,sizeof(dp));        memset(vis,0,sizeof(vis));        memset(mat,127,sizeof(mat));        front=rear=0;        int u,v,s;        for(int i=1;i<=m;++i)        {            scanf("%d%d%d",&u,&v,&s);            u--;v--;            if(mat[u][v]>s)                mat[u][v]=mat[v][u]=s;        }        dist[0][0]=0;        q[++rear][0]=0;q[rear][1]=0;        vis[0][0]=1;        while(front!=rear)        {            if(++front==maxQ) front=0;            int u_sta=q[front][0],u=q[front][1];            vis[u_sta][u]=0;            for(int i=0;i<n;++i)            {                int v_sta=u_sta|(1<<i);                if(dist[v_sta][i]>dist[u_sta][u]+mat[u][i])                {                    dist[v_sta][i]=dist[u_sta][u]+mat[u][i];                    if(!vis[v_sta][i])                    {                        if(++rear==maxQ) rear=0;                        vis[v_sta][i]=1;                        q[rear][0]=v_sta;q[rear][1]=i;                    }                }            }        }        for(int i=0;i<(1<<n);++i)            for(int j=0;j<n;++j)            dp[i>>1]=min(dp[i>>1],dist[i][j]);        int tot_sta=1<<n-1;        for(int i=0;i<tot_sta;++i)            f[i]=dp[i];        for(int k=1;k<=2;++k)        {            for(int i=tot_sta-1;i;--i)                for(u=i;u;u=(--u)&i)            {                v=i^u;                if(dp[i]>max(dp[u],f[v]))                {                    if(i==224 && max(dp[u],f[v])==5)                        cout<<"dd";                    dp[i]=max(dp[u],f[v]);                }            }        }        int query;scanf("%d",&query);        int sta=0;        while(query--)        {            scanf("%d",&u);            sta|=1<<u-2;        }        int ans=INF;        for(int i=0;i<tot_sta;++i)            if((i&sta)==sta)        {            ans=min(ans,dp[i]);        }        if(ans>=INF)            printf("Case %d: -1\n",kk);        else printf("Case %d: %d\n",kk,ans);    }    return 0;}


原创粉丝点击