km算法+slack优化 poj3565

km算法在百度上讲的很清楚

讲一些针对这道题的

有n只蚂蚁要吃n个苹果，求总路径和最小，且不交叉的方案

网上有很多人说是计算几何题

其实，研究一下题意可以得出很重要的性质————最优方案一定不交叉。

根据三角形不等式

ad+bc<ac+bd    dc+ab<ac+bd

因此，我们只需要构造二分图，求最小权和匹配

将所有蚂蚁看做一个集合，所有苹果看做一个集合，求出各蚂蚁到各苹果的距离并取反，在这个二分图上用km算法

容易发现这个模型是完全二分图，因此，加上slack优化可以快很多（在其他情况下slack优化优化为O（n^3），不过也有人说只优化常数）。

这道题的相等判断要注意精度，否则会死循环。

const eps=1e-6; var lx,ly,slack:array[1..100]of real; x1,x2,y1,y2:array[1..100]of longint; vx,vy:array[1..100]of boolean; cost:array[1..100,1..100]of real; b,c:array[1..100]of longint; n:longint; function km(x:longint):boolean; var i:longint; tmp:real; begin if vx[x] then exit(false); vx[x]:=true; for i:=1 to n do begin tmp:=lx[x]+ly[i]-cost[x,i]; if not(vy[i])and(tmp<eps) then begin vy[i]:=true; if (b[i]=0)or(km(b[i])) then begin b[i]:=x;c[x]:=i; exit(true) end end else if tmp<slack[i] then slack[i]:=tmp end; exit(false) end; procedure init; var i,j,k:longint; d:real; begin readln(n); for i:=1 to n do readln(x1[i],y1[i]); for i:=1 to n do readln(x2[i],y2[i]); fillchar(ly,sizeof(ly),0); for i:=1 to n do lx[i]:=-maxlongint; for i:=1 to n do for j:=1 to n do begin cost[i,j]:=-sqrt(sqr(x1[i]-x2[j])+sqr(y1[i]-y2[j])); if cost[i,j]>lx[i] then lx[i]:=cost[i,j] end; fillchar(b,sizeof(b),0); for i:=1 to n do begin fillchar(vx,sizeof(vx),false); fillchar(vy,sizeof(vy),false); fillchar(slack,sizeof(slack),127); while not(km(i)) do begin d:=maxlongint; for k:=1 to n do if not(vy[k]) then if slack[k]<d then d:=slack[k]; for j:=1 to n do begin if vx[j] then lx[j]:=lx[j]-d; if v