【BubbleCup X】D. Exploration plan

这个题首先一眼能看出二分答案……

毕竟连可爱的边界都给你了。

下面就是怎么check

首先预处理跑一遍floyed，预处理出最短路。

用网络流判断能否达到即可。

#include<bits/stdc++.h>
typedef long long ll;
const int N=;
const int M=;
const int inf=1e9+;
using namespace std;
inline int read(){
    int f=,x=;char ch;
    do{ch=getchar();if(ch=='-')f=-;}while(ch<''||ch>'');
    do{x=x*+ch-'';ch=getchar();}while(ch>=''&&ch<='');
    return f*x;
}
namespace MaxFlow{

struct Edge{int u,v,next,f;}G[M];
int head[N],level[N],tot;
inline void init(){tot=;memset(head,,sizeof(head));}
inline void addedge(int u,int v,int f){
    G[++tot].u=u;G[tot].v=v;G[tot].f=f;G[tot].next=head[u];head[u]=tot;
    G[++tot].u=v;G[tot].v=u;G[tot].f=;G[tot].next=head[v];head[v]=tot;
}
inline bool bfs(int s,int t){
    memset(level,,sizeof(level));queue<int>q;
    q.push(s);level[s]=;
    while(!q.empty()){
        int u=q.front();q.pop();
        if(u==t)return ;
        for(int i=head[u];i;i=G[i].next){
            int v=G[i].v,f=G[i].f;
            if(f&&!level[v])level[v]=level[u]+,q.push(v);
        }
    }
    return ;
}
inline int dfs(int u,int maxf,int t){
    if(u==t)return maxf;int rat=;
    for(int i=head[u];i;i=G[i].next){
        int v=G[i].v,f=G[i].f;
        if(f&&level[v]==level[u]+){
            f=dfs(v,min(maxf-rat,f),t);
            rat+=f;G[i].f-=f;G[i^].f+=f;
        }
    }
    if(!rat)level[u]=inf;
    return rat;
}
inline int dinic(int s,int t){
    int ans=;
    while(bfs(s,t))ans+=dfs(s,inf,t);
    return ans;
}

}
int n,m,p,q,f[N],dis[][];
inline bool check(int x){
    MaxFlow::init();int s=,t=n*+;
    for(int i=;i<=n;i++)MaxFlow::addedge(s,i,f[i]);
    for(int i=;i<=n;i++)MaxFlow::addedge(i+n,t,);
    for(int i=;i<=n;i++)if(f[i])for(int j=;j<=n;j++)
    if(dis[i][j]<=x)MaxFlow::addedge(i,j+n,);
    int ans=MaxFlow::dinic(s,t);
    return ans>=q;
}
int main(){
    n=read();m=read();p=read();q=read();
    for(int i=;i<=p;i++){int x=read();f[x]++;}
    memset(dis,0x3f,sizeof(dis));
    for(int i=;i<=n;i++)dis[i][i]=;
    for(int i=;i<=m;i++){
        int u=read(),v=read(),w=read();dis[u][v]=dis[v][u]=min(dis[u][v],w);
    }
    for(int k=;k<=n;k++)for(int i=;i<=n;i++)for(int j=;j<=n;j++)
    dis[i][j]=min(dis[i][j],dis[i][k]+dis[k][j]);
    int l=,r=,ret=-;
    while(l<=r){
        int mid=(l+r)>>;
        if(check(mid))r=mid-,ret=mid;
        else l=mid+;
    }
    printf("%d\n",ret);
}