BZOJ_1486_[HNOI2009]最小圈_01分数规划
Description
Input
Output
Sample Input
4 5
1 2 5
2 3 5
3 1 5
2 4 3
4 1 3
1 2 5
2 3 5
3 1 5
2 4 3
4 1 3
Sample Output
3.66666667
二分答案，边权减去答案，判负环即可。
然而spfa判负环会T掉，于是我使用了dfs判负环的方法。
dfs判负环代码：
void dfs(int x)
{
vis[x]=1;
int i;
for(i=head[x];i&&!ok;i=nxt[i]){
if(dis[to[i]]>dis[x]+cost[i]){
dis[to[i]]=dis[x]+cost[i];
if(vis[to[i]]){
ok=1;return ;
}
dfs(to[i]);
}
}
vis[x]=0;
}
不断更新最小值，直到更新了一圈回来，则说明有负环存在。
总之是比spfa快到不知哪里去。
dis数组清零即可。
代码：
#include <stdio.h>
#include <string.h>
#include <algorithm>
using namespace std;
#define N 3050
#define M 10050
#define du double
#define eps 1e-9
int head[N],to[M],nxt[M],val[M],cnt,n,m;
int dep[N],inq[N],Q[N],l,r,ok,vis[N];
du dis[N],cost[M];
inline void add(int u,int v,int w)
{
to[++cnt]=v;nxt[cnt]=head[u];head[u]=cnt;val[cnt]=w;
}
void dfs(int x)
{
vis[x]=1;
int i;
for(i=head[x];i&&!ok;i=nxt[i]){
if(dis[to[i]]>dis[x]+cost[i]){
dis[to[i]]=dis[x]+cost[i];
if(vis[to[i]]){
ok=1;return ;
}
dfs(to[i]);
}
}
vis[x]=0;
}
bool check(du x)
{
int i;
for(i=1;i<=m;i++) cost[i]=val[i]-x;
for(i=1;i<=n;i++) dis[i]=vis[i]=0;
for(i=1,ok=0;i<=n;i++)
{
dfs(i);
if(ok) return 1;
}
return 0;
}
int main()
{
scanf("%d%d",&n,&m);
int i,x,y,z;
for(i=1;i<=m;i++)
{
scanf("%d%d%d",&x,&y,&z);
add(x,y,z);
}
du L=-100000,R=100000;
while(R-L>eps)
{
du mid=(L+R)*0.5;
if(check(mid)) R=mid;
else L=mid;
}
printf("%.8lf\n",L);
}