CSUSTOJ 搬砖(状压dp)

题目链接

emmm显然的状压dp
我们定义dp[状态][j] 表示该状态下以j为终点的最短时间
先把存图floyd一下，然后开始状压dp处理。
对于K个任务，只需要记录每个经过的点的二进制位就OK，把任务合成一个带有1号起点的状态，查询输出该状态下以目标终点结尾就OK

 #include<bits/stdc++.h>
#include<stdlib.h>
#include<algorithm>
#include<stdio.h>
#include<string.h>
#include<queue>
#include<time.h>
#include <cstdio>
#include <iostream>
#include <vector>
#define ll long long
#define int long long
#define inf 0x3f3f3f3f
#define mods 1000000007
#define modd 998244353
#define PI acos(-1)
#define fi first
#define se second
#define lowbit(x) (x&(-x))
#define mp make_pair
#define pb push_back
#define si size()
#define E exp(1.0)
#define fixed cout.setf(ios::fixed)
#define fixeds(x) setprecision(x)
#define IOS ios::sync_with_stdio(false);cin.tie(0)
 using namespace std;
 ll gcd(ll a,ll b){if(a<0)a=-a;if(b<0)b=-b;return b==0?a:gcd(b,a%b);}
template<typename T>void read(T &res){bool flag=false;char ch;while(!isdigit(ch=getchar()))(ch=='-')&&(flag=true);
for(res=ch-48;isdigit(ch=getchar());res=(res<<1)+(res<<3)+ch - 48);flag&&(res=-res);}
ll lcm(ll a,ll b){return a*b/gcd(a,b);}
ll qp(ll a,ll b,ll mod){ll ans=1;if(b==0){return ans%mod;}while(b){if(b%2==1){b--;ans=ans*a%mod;}a=a*a%mod;b=b/2;}return ans%mod;}//快速幂%
ll qpn(ll a,ll b, ll p){ll ans = 1;a%=p;while(b){if(b&1){ans = (ans*a)%p;--b;}a =(a*a)%p;b >>= 1;}return ans%p;}//逆元   (分子*qp(分母,mod-2,mod))%mod;

ll G[50][50];
const int N=(1<<19);
ll dp[N+20][22];
signed main()
{


    ll n,m,kk;
    read(n);
    read(m);
    read(kk);
    memset(G,inf,sizeof(G));
    for(int i=1; i<=m; i++)
    {
        ll u,v,w;
        read(u);
        read(v);
        read(w);
        u--;
        v--;
        G[u][v]=G[v][u]=min(G[u][v],w);
    }
    for(int i=0; i<n; i++)
    {
        G[i][i]=0;
    }

    for(int i=0; i<n; i++)
    {
        for(int j=0; j<n; j++)
        {
            for(int k=0; k<n; k++)
            {
                G[j][k]=min(G[j][k],G[j][i]+G[i][k]);

            }
        }
    }

    memset(dp,inf,sizeof(dp));
    dp[0][0]=0;
    dp[1][0]=0;
    for(int i=1; i<(1<<n); i++)
    {

        for(int j=0; j<n; j++)
        {

            if((1<<j)&i)
            {

                for(int k=0; k<n; k++)
                {

                    if(k!=j&&(!(i&(1<<k))))
                    {

                        dp[i+(1<<k)][k]=min(dp[i+(1<<k)][k],dp[i][j]+G[j][k]);



                    }


                }


            }


        }


    }



    for(int i=1; i<=kk; i++)
    {
        ll cnt=1;
        ll op;
        read(op);
        op--;
        if(!(cnt&(1<<op)))
            cnt=cnt+(1<<op);
        ll num;
        read(num);
        for(int j=1; j<=num; j++)
        {
            ll x;
            read(x);
            x--;
            if(!(cnt&(1<<x)))
                cnt=cnt+(1<<x);


        }



        printf("%lld\n",dp[cnt][op]);



    }



}