洛谷P2868 [USACO07DEC]Sightseeing Cows G 题解
题目链接:P2868 [USACO07DEC]Sightseeing Cows G
题意:
给定有向图,有点权 F i F_i Fi 和边权 T i T_i Ti ,找一条回路(不一定是简单回路),使得
∑ F i ∑ T i \dfrac{\sum F_i}{\sum T_i} ∑Ti∑Fi
尽可能的大,求出这个值若回路的结点序列中含相同结点,只计算一次点权
2 ≤ n ≤ 1 0 3 , 1 ≤ m ≤ 5 × 1 0 3 , 1 ≤ F i , T i ≤ 1 0 3 2 \le n \le 10^3,1\le m \le 5 \times 10^3,1 \le F_i,T_i \le 10^3 2≤n≤103,1≤m≤5×103,1≤Fi,Ti≤103
比较经典的01分数规划问题
设
∑
F
i
×
1
∑
T
i
>
x
\sum F_i \times \dfrac{1}{\sum T_i} > x
∑Fi×∑Ti1>x
移项可得
∑
(
F
i
−
x
T
i
)
>
0
\sum (F_i - x T_i) > 0
∑(Fi−xTi)>0
取个反可得
∑
(
x
T
i
−
F
i
)
<
0
\sum (xT_i - F_i) < 0
∑(xTi−Fi)<0
然后就变成了判负环的问题。
二分一个 x x x ,然后每次跑个spfa检验一下
时间复杂度 O ( log 2 1 0 7 × n m ) O(\log_2 10^7 \times nm) O(log2107×nm)
为什么大于也可以求出答案呢?因为当 x x x 足够精确时,也就是答案了
这里有个问题,题目说的是回路,不是环(简单回路),
不过事实上答案中都是环(简单回路)
这里有个证明,比较繁琐,因此我来提供一种简单证明
证:考虑最简单的情况,
设两个环 A , B A,B A,B 的交集为结点 u u u (此时图可以是一个8字形)
设 A A A 的答案是 x x x , B B B 的答案是 y y y
-
若 x = y x=y x=y ,
则若合并 A , B A,B A,B ,多出来的几条连接 u u u 的边会使分母增加
则答案小于 x x x
只有这几条边的 T i T_i Ti 都为 0 0 0 时,
合并后的答案才能达到 x x x ,但是数据范围里没有这种情况。
故合并 A , B A,B A,B 一定不是最优的。
-
若 x > y x>y x>y 或 x < y x<y x<y
那显然合并了就更烂了哇。肯定不是最优的。
证毕。
代码:
#include <iostream>
#include <string>
#include <vector>
#include <algorithm>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <cmath>
#include <iomanip>
#include <queue>
using namespace std;
#define int long long
#define INF 0x3f3f3f3f3f3f3f3f
#define N (int)(1e3+15)
#define M (int)(5e3+15)
struct Edge
{
int u,v,T,next;
}e[M];
int n,m,F[N];
double d[N];
int head[N],pos=1,vis[N],cnt[N];
void addEdge(int u,int v,int T)
{
e[++pos]={u,v,T,head[u]};
head[u]=pos;
}
queue<int> q;
bool spfa(int st,double x)
{
while(!q.empty())q.pop();
q.push(st);vis[st]=1;d[st]=0;
while(!q.empty())
{
int u=q.front();q.pop();
vis[u]=0;
for(int i=head[u]; i; i=e[i].next)
{
int v=e[i].v;
double w=x*e[i].T-F[u];
if(d[v]>d[u]+w)
{
d[v]=d[u]+w;
if(!vis[v])
{
q.push(v);++cnt[v];vis[v]=1;
if(cnt[v]>=n)return 1;
}
}
}
}
return 0;
}
signed main()
{
ios::sync_with_stdio(0);
cin.tie(0);cout.tie(0);
// freopen("check.in","r",stdin);
// freopen("check.out","w",stdout);
cin >> n >> m;
for(int i=1; i<=n; i++) cin >> F[i];
for(int i=1,u,v,T; i<=m; i++)
{
cin >> u >> v >> T;
addEdge(u,v,T);
}
double l=0,r=1005;
while(fabs(r-l)>1e-4)
{
double mid=(l+r)/2;
bool ok=0;
for(int i=1; i<=n; i++)
{
d[i]=1e66;
cnt[i]=vis[i]=0;
}
for(int i=1; i<=n; i++)
{
if(!cnt[i])
{
if(spfa(i,mid))
{ok=1;break;}
}
}
ok?l=mid:r=mid;
}
cout << fixed << setprecision(2);
cout << l << '\n';
return 0;
}
转载请说明出处