poj 3621 Sightseeing Cows

http://poj.org/problem?id=3621

分数规划+二分   最优比率环   

不是很难 本题中没有说明从哪个点开始 不过好像默认为1就可以过  poj数据又水了 

里面的用spfa判定部分还是不太懂 

代码及其注释：

#include<iostream>

#include<stdio.h>

#include<algorithm>

#include<math.h>

#include<queue>

#include<string.h>



using namespace std;



const double eps=1e-5;

const int N=1001;

const int M=5005;

int head[N];//邻接表头

struct node

{

    int j,k;

    int next;

}side[M];边

int value[N];//点的价值

bool in[N];//是否在队列中

int innum[N];//如队列次数

double dist[N];//记录

queue<int>str;

void build(int x,int t)//建树

{

    side[t].next=head[x];

    head[x]=t;

}

bool spfa(double L,int n)//看是否有正环  如果有返回true

{

    while(!str.empty())

    str.pop();

    for(int i=1;i<=n;++i)

    {

        dist[i]=0.0;

        str.push(i);

    }

    memset(in,true,sizeof(in));

    memset(innum,0,sizeof(innum));

    while(!str.empty())

    {

        int x=str.front();

        str.pop();

        in[x]=false;

        int t=head[x];

        while(t!=-1)

        {//cout<<t<<endl;

            int j=side[t].j;

            if(dist[j]<dist[x]+value[j]-L*side[t].k)

            {

                dist[j]=dist[x]+value[j]-L*side[t].k;

                if(!in[j])

                {

                    str.push(j);

                    in[j]=true;

                    ++innum[j];

                    if(innum[j]>=n)//有环

                    return true;

                }

            }

            t=side[t].next;

        }

    }

    return false;

}

int main()

{

   // freopen("data.txt","r",stdin);

    int n,m;

    while(scanf("%d %d",&n,&m)!=EOF)

    {

        memset(head,-1,sizeof(head));

        for(int i=1;i<=n;++i)

        scanf("%d",&value[i]);

        for(int i=1;i<=m;++i)

        {

            int x;

            scanf("%d %d %d",&x,&side[i].j,&side[i].k);

            build(x,i);

        }

        double l=1.0,r=200.0;

        while(r-l>eps)//用二分 不断向答案逼近

        {

            double mid=(l+r)/2.0;

            if(spfa(mid,n))

            l=mid;

            else

            r=mid;

        }

        printf("%.2f\n",l);

    }

    return 0;

}