[USACO 2008 Jan. Silver]架设电话线 —— 最短路+二分
一道图论的最短路题。一开始连最短路都没想到，可能是做的题太少了吧，完全没有思路。
题目大意：
FJ的农场周围分布着N根电话线杆，任意两根电话线杆间都没有电话线相连。一共P对电话线杆间可以拉电话线，第i对电话线杆的两个端点分别为A_i、B_i，它们间的距离为L_i 。数据中保证每对{A_i，B_i}最多只出现1次。FJ的任务是找一条将1号和N号电话线杆连起来的路径。电信公司最终同意免费为FJ连结任意K对电话线杆。此外的电话线，FJ需要为它们付的费用，等于其中最长的电话线的长度（每根电话线仅连结一对电话线杆）。如果需要连结的电话线杆不超过K对，那么FJ的总支出为0。求FJ最少需要在电话线上花多少钱，如果任务不可能完成，输出-1。1 <= N <= 1,000
1 <= P <= 10,000
1 <= L_i <= 1,000,000
0 <= K < N
样例输入
5 7 1
1 2 5
3 1 4
2 4 8
3 2 3
5 2 9
3 4 7
4 5 6
样例输出
4
注意一句话题意 ： 在加权无向图上求出一条从 1号结点到n号结点的路径，使路径上第k+1大的边权尽量小。”
我们不妨将这条k+1大的边权设为\(P\),将图中大于等于\(P\)的边权类似离散的思想的将它们赋为1，反之赋为0。这里就不给出图了。我们再对离散后的图跑一遍最短路，设这时1-n的最短路长度为\(X\)，我们便可以发现\(P\)在这条最短录径上第\(X+1\)的边（因为比\(P\)大的边只有\(X\)条）
试问，如果这个\(X\)比\(K+1\)还大，请问这个\(P\)还可以作为我们的答案吗？肯定是不行的，因为我们跑的是最短路，也就是\(X\)最小的情况，如果这个\(X\)比\(K+1\)还大，那么这条边\(P\)就不可能是\(K+1\)大的边，自然就会舍去。
现在要解决的是\(P\)的问题，直接枚举吗？其实根本不需要，我们发现如果\(P\)跑出来的最短路小于等于\(K\)，那么\(P+1\)跑出来的最短路也一定小于等于\(K\),答案明显满足单调性，因此我们可以使用二分来做这道题。
（有的人可能要为为什么是小于等于\(K\),按照一句话题意不是要我们求\(X\)等于\(K\)的情况。我的理解是有的路径的条数可能比\(K\)还少，但我们依然可以选择这条路径上的边。另外，这样也适应了二分的情况。如果是等于，二分则跑不出来。）
关于对图的离散化，不要告诉我们你们会直接把所有边都赋为1和0，在代码中我们只需要这么做（我用的是dijkstra的堆优化）
v=G[now.u][i].v,w=G[now.u][i].w;
if(w>=P) w1=1;
else w1=0;
if(dis[now.u]+w1<dis[v]) {
dis[v]=dis[now.u]+w1;
q.push( node (dis[v],v));
}
代码
#include <iostream>
#include <cstdio>
#include <vector>
#include <algorithm>
#include <stack>
#include <cstring>
#include <queue>
using namespace std;
#define N 50010
#define inf 0x7f7f7f7f
struct node {
int u,dis;
node () {};
node (int D,int U) { u=U; dis=D; }
bool operator < (const node& h) const {
return dis>h.dis;
}
};
struct Gragh {
int v,w;
Gragh () {};
Gragh (int V,int W) {v=V;w=W;}
};
int P[10],n,m,k,ans=inf,vis[N],sum=1,dis[N];
vector<Gragh> G[N];
priority_queue<node> q;
int dijkstra(int st,int P) {
q.push( node (0,st) );
memset(dis,0x3f,sizeof(dis));
dis[st]=0;
memset(vis,0,sizeof(vis));
while(!q.empty() ) {
node now=q.top();
q.pop();
if(vis[now.u]) continue;
vis[now.u]=1;
for(int i=0,v,w,w1;i<G[now.u].size();i++) {
v=G[now.u][i].v,w=G[now.u][i].w;
if(w>=P) w1=1;
else w1=0;
if(dis[now.u]+w1<dis[v]) {
dis[v]=dis[now.u]+w1;
q.push( node (dis[v],v));
}
}
}
return dis[n];
}
bool check(int mid) {
if(dijkstra(1,mid)<=k) return 1;
return 0;
}
/*
4 3 3
1 2 1
2 3 2
3 4 3
*/
int main() {
int maxw=0;
cin>>n>>m>>k;
for(int i=1,u,v,w;i<=m;i++) {
cin>>u>>v>>w;
maxw=max(maxw,w);
G[u].push_back( Gragh(v,w) );
G[v].push_back( Gragh(u,w) );
}
int mid,l=0,r=maxw,flag=0;
while(l<r) {
mid=(l+r+1)/2;
if(check(mid)) r=mid-1,flag=1;
else l=mid;
}
if(!flag) cout<<"-1";
else cout<<l;
}
(话说二分每次都打不对orz)
