【分数规划-最大比例环】POJ3621[Sightseeing Cows]题解

题目概述

给出 $n$ 个点， $m$ 条单向边。每个点有点权 $w_i$ ，每条边有边权 $t_i$ 。求点权比边权最大的环。

解题报告

在POJ做分数题简直就是作死，可以WA到飞起来，还是分数规划裸题（搞事啊，为什裸题我现在才做过啊），二分枚举答案 $mid$ ，然后用新边权 $w_x-mid\times t(x,y)$ 建边。只要存在一个环的边权和大于等于 $0$ ，就说明验证成功。

等等……这怎么找啊……其实边权和大于等于 $0$ 意味着这是一个正权回路，所以用spfa就行了。

解题报告

#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
using namespace std;
typedef double DB;
const int maxn=1000,maxm=5000;

int n,m,v[maxn+5];
int E,lnk[maxn+5],nxt[maxm+5],son[maxm+5],w[maxm+5];
int que[maxn+5],ti,vis[maxn+5],tot[maxn+5];DB dis[maxn+5];

inline int fcmp(DB a,DB b) {if (fabs(a-b)<1e-10) return 0;if (a<b) return -1;return 1;}
#define AM(x) (((x)+1)%maxn)
inline bool check(DB mid)
{
    for (int i=1;i<=n;i++) dis[i]=-1e100,tot[i]=0;ti++;
    int Head=0,Tail=0;que[++Tail]=1;dis[1]=0;vis[1]=ti;
    while (Head!=Tail)
    {
        int x=que[Head=AM(Head)];vis[x]=0;
        for (int j=lnk[x];j;j=nxt[j])
            if (dis[x]+v[x]-mid*w[j]>dis[son[j]])
            {
                dis[son[j]]=dis[x]+v[x]-mid*w[j];
                if (vis[son[j]]!=ti)
                {
                    que[Tail=AM(Tail)]=son[j];vis[son[j]]=ti;
                    if (++tot[son[j]]>n) return true;
                }
            }
    }
    return false;
}
#define Add(x,y,z) son[++E]=(y),w[E]=(z),nxt[E]=lnk[x],lnk[x]=E
int main()
{
    freopen("program.in","r",stdin);
    freopen("program.out","w",stdout);
    while (~scanf("%d%d",&n,&m))
    {
        for (int i=1;i<=n;i++) scanf("%d",&v[i]);E=0;memset(lnk,0,sizeof(lnk));
        for (int i=1,x,y,z;i<=m;i++) scanf("%d%d%d",&x,&y,&z),Add(x,y,z);
        DB L=0,R=1e6;
        for (DB mid=(L+R)/2;fcmp(L,R)<=0;mid=(L+R)/2)
            if (check(mid)) L=mid+1e-4; else R=mid-1e-4;
        printf("%.2f\n",R);
    }
    return 0;
}