洛谷3001 [USACO10DEC]巨无霸（spfa）

题目

洛谷3001 [USACO10DEC]巨无霸

题解

spfa
d[y]=min{ d[x]*e[k].c }
这题的坑点是有“小于1环”，所以不但不能用dijkstra，还要加一个判“小于1环”。
还卡了精度，要用long double，long double的输入输出用%Lf。

代码

#include<cstdio>
#include<cstring>
#include<algorithm>
#define double long double
using namespace std;
typedef long long ll;
const int END=2005;
const double INF=1e15;
const int MAXN=2010,MAXM=25010;

int n,m,S,T;double W;

struct E{int y;double c;int next;}e[MAXM];int len=0,last[MAXN];
void ins(int x,int y,double c)
{
    e[++len]=(E){y,c,last[x]};last[x]=len;
}

double d[MAXN];int cnt[MAXN];
int head,tail,q[MAXN];bool inq[MAXN];
void spfa()
{
    for(int i=1;i<=n;i++) d[i]=INF;d[S]=W;
    head=0,tail=1;q[0]=S;
    while(head!=tail)
    {
        int x=q[head++];if(head==END) head=0;
        inq[x]=false;
        for(int k=last[x];k;k=e[k].next)
        {
            int y=e[k].y;
            if(d[y]>d[x]*e[k].c)
            {
                d[y]=d[x]*e[k].c;
                cnt[y]=cnt[x]+1;if(cnt[y]>=n){puts("0");exit(0);}
                if(!inq[y])
                {
                    q[tail++]=y;if(tail==END) tail=0;
                    inq[y]=true;
                }
            }
        }
    }
}

int main()
{
    scanf("%d%d%Lf%d%d",&n,&m,&W,&S,&T);
    for(int i=1;i<=m;i++)
    {
        int x,y;double c;
        scanf("%d%d%Lf",&x,&y,&c);
        ins(x,y,c);
    }
    spfa();
    if(d[T]==INF) puts("0");
    else printf("%.6Lf\n",d[T]);
    return 0;
}