2019ICPC南昌邀请赛A （斯坦纳树）

做完HDU4085 来补这道 两道题目差不多 就是在合理状态的判断上有所区别。
题目链接：https://nanti.jisuanke.com/t/40253

大致题意

给定n个点，m条边，给定4对点，要求将这4对点对应联通，重复用边只要算一次。

思路

还是先跑一个斯坦纳树森林，然后合理状态有所不同，需要一一对应，也就是对应位置必须同时出现。但是有个坑点，4对点中可能有的点会重复出现，所以对于状态点用vector存一下相应的图上点。具体避免方式看代码。

这种模板题，大概就是先求出斯坦纳森林，然后用对应状态去更新dp。

代码

#include<bits/stdc++.h>
using namespace std;
#define maxn 35
#define maxm 1006
#define ll long long int
#define INF 0x3f3f3f3f
#define inc(i,l,r) for(int i=l;i<=r;i++)
#define dec(i,r,l) for(int i=r;i>=l;i--)
#define mem(a) memset(a,0,sizeof(a))
#define sqr(x) (x*x)
#define inf (ll)2e18+1
#define PI acos(-1)
#define ls x<<1
#define rs x<<1|1
ll read(){
    ll x=0,f=1ll;char ch=getchar();
    while(!isdigit(ch)){if(ch=='-')f=-1;ch=getchar();}
     while(isdigit(ch))x=x*10+ch-'0',ch=getchar();
     return f*x;
}
int n,m,cnt;
string s1,s2;
unordered_map<string ,int> mp;
struct edge{int v,w,nxt;}ed[maxm<<1];
int head[maxn],tot;
void add(int x,int y,int z){ed[++tot]={y,z,head[x]};head[x]=tot; }
int dp[1<<10],f[maxn][1<<10];
vector<int>id[maxn];
bool inq[maxn];
queue<int>q;
void spfa(int sta){
    while(!q.empty()){
        int u=q.front();q.pop();inq[u]=0;
        for(int i=head[u];i;i=ed[i].nxt){
            int v=ed[i].v;
            if(f[v][sta]>f[u][sta]+ed[i].w){
                f[v][sta]=f[u][sta]+ed[i].w;
                if(!inq[v]){
                    inq[v]=1;
                    q.push(v);
                }
            }
        }
    }
}
bool check(int sta){
    inc(i,0,3){
        int x=(((1<<(2*i))&sta)>0);
        int y=(((1<<(2*i+1))&sta)>0);
        if(x^y)return false;
    }
    return true;
}
int main()
{
    //n=read();m=read();
    scanf("%d%d",&n,&m);
    int s=(1<<8)-1;
    mem(inq);
    inc(i,1,n){
        //scanf("%s",&s1);
        cin>>s1;
        mp[s1]=i;
    }
    int x,y,z;
    inc(i,1,m){
        cin>>s1>>s2>>z;
        x=mp[s1];y=mp[s2];
        add(x,y,z);add(y,x,z);
    }
    inc(i,1,4){
        cin>>s1>>s2;
        x=mp[s1];y=mp[s2];
        id[x].push_back(++cnt);
        id[y].push_back(++cnt);
    }
    inc(i,1,n)inc(j,0,s)f[i][j]=INF;
    inc(i,1,n)for(int j=0;j<id[i].size();j++)f[i][(1<<(id[i][j]-1))]=0;
    inc(sta,0,s){
        inc(i,1,n){
            for(int j=sta;j;j=(j-1)&sta){
                f[i][sta]=min(f[i][sta],f[i][j]+f[i][sta^j]);
            }
            if(f[i][sta]<INF){
                inq[i]=1;
                q.push(i);
            }
        }
        spfa(sta);
    }
    inc(i,0,s){
        dp[i]=INF;
        inc(j,1,n)dp[i]=min(dp[i],f[j][i]);
    }
    inc(sta,0,s){
        if(!check(sta))continue;
        for(int j=sta;j;j=(j-1)&sta){
            if(!check(j))continue;
            dp[sta]=min(dp[sta],dp[j]+dp[sta^j]);
        }
    }
    printf("%d\n",dp[s]);
    return 0;
}