首页 > 代码库 > 【poj3621】最优比率环
【poj3621】最优比率环
题意:
给定n个点,每个点有一个开心度F[i],每个点有m条单向边,每条边有一个长度d,要求一个环,使得它的 开心度的和/长度和 这个比值最大。
n<=1000,m<=5000
题解:
最优比率环,很像以前做过的一题最优比率生成树。
首先二分一个答案r=sigma(xi*fi)/sigma(xi*di),设z=r*sigma(xi*di)-sigma(xi*fi),若能找到一个环满足z<=0,则代表sigma(xi*fi)>=r*sigma(xi*di),则r可以比当前的r要大,则l=mid;否则r=mid;
判断能否有一个环满足让z<=0,我们就让所有的边权=r*di-f[a[i].x],然后用spfa判负环即可。
1 #include<cstdio>
2 #include<cstdlib>
3 #include<cstring>
4 #include<iostream>
5 #include<algorithm>
6 #include<queue>
7 #include<vector>
8 #include<cmath>
9 using namespace std;
10
11 const int N=1100,M=5100;
12 const double INF=(double)1e9;
13 int n,m,len,first[N],cnt[N];
14 double dis[N],F[N];
15 bool in[N];
16 struct node{
17 int x,y,next;
18 double d,w;
19 }a[2*M];
20 queue<int> q;
21
22 double mid;
23
24 double myabs(double x){return x>0 ? x:-x;}
25
26 void ins(int x,int y,double d)
27 {
28 a[++len].x=x;a[len].y=y;a[len].d=d;
29 a[len].next=first[x];first[x]=len;
30 }
31
32 double spfa(int st)
33 {
34 while(!q.empty()) q.pop();
35 for(int i=1;i<=n;i++) dis[i]=INF;
36 memset(cnt,0,sizeof(cnt));
37 memset(in,0,sizeof(in));
38 q.push(st);dis[st]=0;in[st]=1;
39 while(!q.empty())
40 {
41 int x=q.front();in[x]=0;cnt[x]++;q.pop();
42 if(cnt[x]==n) return 1;
43 for(int i=first[x];i;i=a[i].next)
44 {
45 int y=a[i].y;
46 if(dis[y]>=dis[x]+a[i].w)
47 {
48 dis[y]=dis[x]+a[i].w;
49 if(!in[y]) in[y]=1,q.push(y);
50 }
51 }
52 }
53 return 0;
54 }
55
56 bool check(double r)
57 {
58 for(int i=1;i<=len;i++) a[i].w=r*a[i].d-F[a[i].x];
59 return spfa(1);
60 }
61
62 int main()
63 {
64 freopen("a.in","r",stdin);
65 scanf("%d%d",&n,&m);
66 len=0;
67 memset(first,0,sizeof(first));
68 double l=0,r=0;
69 for(int i=1;i<=n;i++)
70 {
71 scanf("%lf",&F[i]);
72 r+=F[i];
73 }
74 for(int i=1;i<=m;i++)
75 {
76 int x,y;double d;
77 scanf("%d%d%lf",&x,&y,&d);
78 ins(x,y,d);
79 }
80 while(l<r)
81 {
82 mid=(l+r)/2;
83 if(check(mid)) l=mid;
84 else r=mid;
85 if(myabs(l-r)<=0.0000001) break;
86 }
87 mid=6.0;
88 if(myabs(l)<=0.0000001) printf("0\n");
89 else printf("%.2lf\n",l);
90 return 0;
91 }
【poj3621】最优比率环
声明:以上内容来自用户投稿及互联网公开渠道收集整理发布,本网站不拥有所有权,未作人工编辑处理,也不承担相关法律责任,若内容有误或涉及侵权可进行投诉: 投诉/举报 工作人员会在5个工作日内联系你,一经查实,本站将立刻删除涉嫌侵权内容。