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

23 篇文章 0 订阅
17 篇文章 0 订阅

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

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

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

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

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

#include <iostream>
#include <cstdio>
#include <cstring>
using namespace std;
#define INF 2000000000
int 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;
}



评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值