P4542 [ZJOI2011]营救皮卡丘（Floyd+网络流）

P4542 [ZJOI2011]营救皮卡丘

乍一看似乎没啥题相似的

仔细一看，$N<=150$

边又是双向边，似乎可以用Floyd搞

 

先跑一遍Floyd处理出$dis[i][j]$

注意到走据点要先走小的才能走大的

也就是说，$i<j<k$时，$dis[i][j]$不能从$k$转移过来

并且实际走路径时，编号也必须从小到大

 

于是题目转化成了：

给定序列$0,1,2,3,.....,n-1,n$，给出每两个数字之间的转移代价$dis[i][j](i<j)$，

用$k$条从0开始的子序列覆盖整条序列的最小代价，且每个数（除0外）恰好被覆盖一次。（你闲着没事走两次干啥）

有没有可能$i$到$j$的最优路径中间经过$k$，使$k$被覆盖2次？不会，因为在Floyd中已经处理掉了。

 

卧槽这不是P2469 [SDOI2010]星际竞速吗！

于是就转化为一个最小路径覆盖问题写辣

#include<iostream>
#include<cstdio>
#include<cstring>
#include<queue>
using namespace std;
#define N 505
#define M 200005
int n,m,k,S,T0,T,tC,dis[N][N],d[N],a[N],p[N];
queue <int> h; bool inh[N];
int cnt=,hd[N],nxt[M],ed[N],poi[M],val[M],cst[M];
inline void adde(int x,int y,int v1,int v2){
    nxt[ed[x]]=++cnt, hd[x]=hd[x]?hd[x]:cnt,
    ed[x]=cnt, poi[cnt]=y, val[cnt]=v1, cst[cnt]=v2;
}
inline void link(int x,int y,int v1,int v2){adde(x,y,v1,v2),adde(y,x,,-v2);}
bool bfs(){
    memset(d,,sizeof(d)); int Inf=d[];
    h.push(S); inh[S]=; a[S]=Inf; d[S]=;
    while(!h.empty()){
        int x=h.front(); h.pop(); inh[x]=;
        for(int i=hd[x];i;i=nxt[i]){
            int to=poi[i];
            if(val[i]>&&d[to]>d[x]+cst[i]){
                d[to]=d[x]+cst[i]; p[to]=i;
                a[to]=min(a[x],val[i]);
                if(!inh[to]) h.push(to),inh[to]=;
            }
        }
    }if(d[T]==Inf) return ;
    tC+=a[T]*d[T];
    for(int i=T;i!=S;i=poi[p[i]^])
        val[p[i]]-=a[T],val[p[i]^]+=a[T];
    return ;
}
void Floyd(){
    memset(dis,,sizeof(dis));
    for(int i=;i<=n;++i) dis[i][i]=;
    for(int i=,u,v,w;i<=m;++i){
        scanf("%d%d%d",&u,&v,&w);
        dis[u][v]=dis[v][u]=min(dis[v][u],w);
    }
    for(int k=;k<=n;++k)
        for(int i=;i<=n;++i)
            for(int j=;j<=n;++j)
                if(k<=i||k<=j)//不能从大的转移回来
                    dis[i][j]=min(dis[i][j],dis[i][k]+dis[k][j]);
}
int main(){
    scanf("%d%d%d",&n,&m,&k); Floyd();
    S=n*+; T=S+;
    for(int i=;i<=n;++i)
        link(S,i,,),link(i+n+,T,,);
    link(S,,k,);//最大流量限制为k
    for(int i=;i<=n;++i)
        for(int j=i+;j<=n;++j)
                link(i,j+n+,,dis[i][j]);
    while(bfs());
    printf("%d",tC);
    return ;
}