借鉴这里:http://www.cnblogs.com/wangfang20/archive/2013/07/25/3213628.html
有n个城市,m个边界线,p名士兵。现在士兵要按一定顺序攻占城市,不能穿过一些线段。士兵有一个容量为K的背包装粮食,士兵如果攻占城市,就能装满背包。从城市到城市消耗的粮食等于两城市的距离,如果距离大于士兵当前的背包的容量,士兵就不能走这条路。求p个士兵攻占完所有城市所需要的最小背包容量k。
二分答案 , 将所有的直线的两个端点和城市混在一起,将能直接到达的两个点连线,求一次floyd最短路径。二分枚举bag容量,只需要判断最小路径覆盖是否小于等于p,就是最少需要多少个士兵才能全部占领
#include<iostream>
#include<cstring>
#include<cstdio>
#include<cmath>
using namespace std;
#define clr(x) memset(x,0,sizeof(x));
#define INF 0x1f1f1f1f
#define maxn 510
#define maxm 510
#define eps 1e-6
struct Point{
double x;
double y;
}city[maxn];
struct Edge{
Point a,b;
}edge[maxm];
int ord[maxn];
double d[maxn+maxm*2][maxn+maxm*2];
int n,m,p;
double multi(Point p0, Point p1, Point p2)
{ return (p1.x-p0.x)*(p2.y-p0.y)-(p2.x-p0.x)*(p1.y-p0.y);}
double Dis(Point a,Point b)
{ return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));}
bool isIntersected(Point s1,Point e1,Point s2,Point e2){
//判断线段是否相交(非规范相交)
return multi(s1,e1,s2)*multi(s1,e1,e2)<-eps&&
multi(s2,e2,s1)*multi(s2,e2,e1)<-eps;
}
int g[maxn][maxn];
int v[maxn];
int link[maxn];
int find(int k){
int i;
for(i=1;i<=n;i++){
if(g[k][i]&&!v[i]){
v[i]=1;
if(link[i]==0||find(link[i])){
link[i]=k;
return 1;
}
}
}
return 0;
}
int maxmatch(){
int i;
clr(link);
int res=0;
for(i=1;i<=n;i++){
clr(v);
if(find(i))res++;
}
return n-res;
}
int main(){
int T,tot,i,j,k;
scanf("%d",&T);
while(T--){
scanf("%d%d%d",&n,&m,&p);
for(i=1;i<=n;i++)
scanf("%lf%lf",&city[i].x,&city[i].y);
tot=n;
for(i=1;i<=m;i++){
scanf("%lf%lf%lf%lf",&edge[i].a.x,&edge[i].a.y,&edge[i].b.x,&edge[i].b.y);
city[++tot].x=edge[i].a.x;city[tot].y=edge[i].a.y;
city[++tot].x=edge[i].b.x;city[tot].y=edge[i].b.y;
}
for(i=1;i<=n;i++)
scanf("%d",&ord[i]);
for(i=1;i<=tot;i++)
for(j=1;j<=tot;j++)
d[i][j]=INF;
for(i=1;i<tot;i++){
for(j=i+1;j<=tot;j++){
int flag=0;
for(k=1;k<=m;k++)
if(isIntersected(city[i],city[j],edge[k].a,edge[k].b)){
flag=1;break;
}
if(flag==0)d[i][j]=d[j][i]=Dis(city[i],city[j]);
}
}
for(k=1;k<=tot;k++)
for(i=1;i<=tot;i++)
for(j=1;j<=tot;j++){
d[i][j]=min(d[i][j],d[i][k]+d[k][j]);
}
double maxlen=0.0;
for(i=1;i<=n;i++)
for(j=1;j<=n;j++){
if(d[i][j]<INF)maxlen=max(d[i][j],maxlen);
}
double mid,l=0,r=maxlen;
while(r-l>eps){
mid=(l+r)/2;
clr(g);
for(i=1;i<n;i++)
for(j=i+1;j<=n;j++){
if(d[ord[i]][ord[j]]<=mid)
g[ord[i]][ord[j]]=1;
}
if(maxmatch()<=p)r=mid;
else l=mid;
}
printf("%.2lf\n",r);
}
return 0;
}