BZOJ 1486 HNOI2009 最小圈 二分答案+DFS

题目大意：裸的最优比例环

直接二分答案+SPFA 这样会T 因为数据卡SPFA

SPFA在负环非常小的时候会退化成Bellman-Ford 时间复杂度是O(nm) (好像是O(n*m^2)?我忘了)的

换一种方法 枚举每个点 从每个点开始DFS 只沿着能将指向的点dis减小的边搜索 搜到栈中的点就返回true

期望复杂度O(n^2) 最坏复杂度O(2^n) 这种东西能过我也是醉了- -

#include <cstdio>#include <cstring>#include <iomanip>#include <iostream>#include <algorithm>#define M 3030#define EPS 1e-9#define INF 1e7using namespace std;struct abcd{int to,next;double f;}table[10100];int head[M],tot;int n,m,T;double a[10100],dis[M];void Add(int x,int y){table[++tot].to=y;table[tot].next=head[x];head[x]=tot;}void DFS(int x){static int v[M];int i;if(v[x]==T) throw(true);v[x]=T;for(i=head[x];i;i=table[i].next)if(dis[x]+table[i].f<dis[table[i].to]){dis[table[i].to]=dis[x]+table[i].f;DFS(table[i].to);}v[x]=0;}bool Judge(double x){int i;for(i=1;i<=m;i++)table[i].f=a[i]-x;for(i=1;i<=n;i++){try{++T;memset(dis,0,sizeof dis);DFS(i);}catch(bool){return true;}}return false;}double Bisection(){int i;double l=INF,r=-INF;for(i=1;i<=m;i++)l=min(l,a[i]),r=max(r,a[i]);while(r-l>EPS){double mid=(l+r)/2;if( Judge(mid) )r=mid;else l=mid;}return (l+r)/2;}int main(){int i,x,y;cin>>n>>m;for(i=1;i<=m;i++)scanf("%d%d%lf",&x,&y,&a[i]),Add(x,y);cout<<fixed<<setprecision(8)<<Bisection()<<endl;return 0;}