01分数规划
题目大意:给一张图,有点权 F[i] F [ i ] 和边权 T[i] T [ i ] ,求使得 ∑F[i]∑T[i] ∑ F [ i ] ∑ T [ i ] 最大的环的这个比值。
01分数规划的最大比例环问题。套路二分后SPFA跑最长路求正权环。
代码:
#include<cstdio>
#include<cstring>
#include<algorithm>
#define N 1005
#define eps 1e-8
using namespace std;
typedef double DB;
struct edge{ int next,to,d; }ed[N*5];
int n,m,k,h[N],F[N],q[N*N],sum[N];
DB d[N]; bool f[N];
#define addedge(x,y,z) ed[++k]=(edge){h[x],y,z},h[x]=k
bool pd(DB mid){
memset(sum,0,sizeof(sum));
memset(f,false,sizeof(f));
for (int i=2;i<=n;i++) d[i]=-1e100;
int l=0,r=1; d[1]=0,q[1]=1,f[1]=true;
while (l<=r){
int x=q[++l]; f[x]=false;
for (int i=h[x],v;i;i=ed[i].next)
if (d[x]+F[x]-mid*ed[i].d>d[v=ed[i].to]){
d[v]=d[x]+F[x]-mid*ed[i].d;
if (!f[v]){
q[++r]=v,f[v]=true;
if (++sum[v]>n) return true;
}
}
}
return false;
}
int main(){
scanf("%d%d",&n,&m);
for (int i=1;i<=n;i++) scanf("%d",&F[i]);
for (int i=1,x,y,z;i<=m;i++)
scanf("%d%d%d",&x,&y,&z),addedge(x,y,z);
DB l=0,r=1e6,mid;
while (r-l>-eps)
if (pd(mid=(l+r)/2)) l=mid+1e-6;
else r=mid-1e-6;
return printf("%.2f\n",r),0;
}