POJ-1860 Currency Exchange SPFA判断环
题目链接：http://poj.org/problem?id=1860
题目大意就是判断一个有向图中的是否存在正权环。利用Bellman-ford判断负环的思想，改为求最长路，如果在n-1次后再来一次bellman，如果还能更新，那么就是存在正权环。
SPFA：
1 //STATUS:G++_AC_0MS_728KB 2 #include<stdio.h> 3 #include<stdlib.h> 4 #include<string.h> 5 #include<math.h> 6 #include<iostream> 7 #include<string> 8 #include<algorithm> 9 #include<vector> 10 #include<queue> 11 #include<stack> 12 using namespace std; 13 #define LL __int64 14 #define pii pair<int,int> 15 #define Max(a,b) ((a)>(b)?(a):(b)) 16 #define Min(a,b) ((a)<(b)?(a):(b)) 17 #define mem(a,b) memset(a,b,sizeof(a)) 18 #define lson l,mid,rt<<1 19 #define rson mid+1,r,rt<<1|1 20 const int N=110,INF=0x3f3f3f3f,MOD=1999997; 21 const LL LLNF=0x3f3f3f3f3f3f3f3fLL; 22 23 struct Edge{ 24 int u,v; 25 double r,s; 26 }e[N*N]; 27 28 double d[N]; 29 double v; 30 int first[N],next[N*N],inq[N],cou[N]; 31 int n,m,s,mt; 32 33 void adde(int a,int b,double r,double s) 34 { 35 e[mt].u=a,e[mt].v=b; 36 e[mt].r=r,e[mt].s=s; 37 next[mt]=first[a],first[a]=mt++; 38 } 39 40 int spfa() 41 { 42 int x,i; 43 double t; 44 queue<int> q; 45 mem(inq,0);mem(cou,0); 46 for(i=1;i<=n;i++)d[i]=-1; 47 d[s]=v; 48 cou[s]=1; 49 q.push(s); 50 while(!q.empty()) 51 { 52 x=q.front();q.pop(); 53 inq[x]=0; 54 for(i=first[x];i!=-1;i=next[i]){ 55 t=(d[x]-e[i].s)*e[i].r; 56 if(t>=0 && t>d[e[i].v]){ 57 d[e[i].v]=t; 58 if(!inq[e[i].v]){ 59 if(++cou[e[i].v]>=n)return 1; 60 inq[e[i].v]=1; 61 q.push(e[i].v); 62 } 63 } 64 } 65 } 66 return 0; 67 } 68 69 int main() 70 { 71 // freopen("in.txt","r",stdin); 72 int i,j,a,b; 73 double r1,s1,r2,s2; 74 while(~scanf("%d%d%d%lf",&n,&m,&s,&v)) 75 { 76 mt=0; 77 mem(first,-1); 78 for(i=0;i<m;i++){ 79 scanf("%d%d%lf%lf%lf%lf",&a,&b,&r1,&s1,&r2,&s2); 80 adde(a,b,r1,s1); 81 adde(b,a,r2,s2); 82 } 83 84 printf("%s\n",spfa()?"YES":"NO"); 85 } 86 return 0; 87 }
