poj2240Arbitrage

题目链接：http://poj.org/problem?id=2240

floyd:

 1 #include<cstdio>
 2 #include<cstring>
 3 #include<string>
 4 #include<map>
 5 #include<iostream>
 6 using namespace std;
 7 const int maxn=32;
 8 int n,m,u,v;
 9 double r;
10 double pic[maxn][maxn];
11
12 void floyd()
13 {
14     for(int k=1;k<=n;k++)
15         for(int i=1;i<=n;i++)
16             for(int j=1;j<=n;j++)
17                 pic[i][j]=max(pic[i][j],pic[i][k]*pic[k][j]);
18 }
19
20 int main()
21 {
22     char s1[100],s2[100],s[100];
23     int cas=0;
24     while(scanf("%d",&n)!=EOF&&n)
25     {
26         int ok=0;
27         map<string,int> mapp;
28         mapp.clear();
29         for(int i=1;i<=n;i++){
30             scanf("%s",s);
31             mapp[s]=i;
32         }
33         scanf("%d",&m);
34
35          for(int i=1;i<=n;i++)
36             for(int j=1;j<=n;j++)
37             pic[i][j]=(i==j?1:0);
38
39         for(int i=0;i<m;i++)
40         {
41             scanf("%s%lf%s",s1,&r,s2);
42             u=mapp[s1];
43             v=mapp[s2];
44             pic[u][v]=r;
45         }
46
47         floyd();
48         for(int i=1;i<=n;i++)
49           if(pic[i][i]>1) {ok=1;break;}
50         printf("Case %d: ",++cas);
51         if(ok) puts("Yes");
52         else puts("No");
53     }
54     return 0;
55 }

也可用bellman-ford判"正环"。

 1 #include<cstdio>
 2 #include<cstring>
 3 #include<string>
 4 #include<map>
 5 using namespace std;
 6 const int maxn=32;
 7 double dis[maxn];
 8 double r;
 9 int n,m,u,v;
10 struct node
11 {
12     int u,v;
13     double r;
14 }e[1010];
15 bool bellman_ford()
16 {
17     int ok=0;
18     for(int i=1;i<=n;i++)
19         dis[i]=0;
20     dis[1]=1;
21     for(int i=0;i<n;i++)
22     {
23         ok=0;
24         for(int i=0;i<m;i++)
25             if(dis[e[i].v]<dis[e[i].u]*e[i].r){
26                     ok=1;
27                 dis[e[i].v]=dis[e[i].u]*e[i].r;
28             }
29         if(!ok) break;
30         if(i==n-1) return 1;
31     }
32     return 0;
33 }
34
35 int main()
36 {
37     int cas=0;
38     char s[100],s1[100],s2[100];
39     while(scanf("%d",&n)!=EOF&&n)
40     {
41         map<string,int> mapp;
42         mapp.clear();
43         for(int i=1;i<=n;i++)
44         {
45             scanf("%s",s);
46             mapp[s]=i;
47         }
48         scanf("%d",&m);
49         for(int i=0;i<m;i++)
50         {
51             scanf("%s%lf%s",s1,&r,s2);
52             u=mapp[s1];
53             v=mapp[s2];
54             e[i].u=u;
55             e[i].v=v;
56             e[i].r=r;
57         }
58         printf("Case %d: ",++cas);
59         if(bellman_ford()) puts("Yes");
60         else puts("No");
61     }
62 }

或者spfa。

 1 #include<cstdio>
 2 #include<cstring>
 3 #include<string>
 4 #include<map>
 5 #include<queue>
 6 using namespace std;
 7
 8 const int maxn=32;
 9 double r,pic[maxn][maxn];
10 int n,m,u,v;
11 double dis[maxn];
12
13 bool inq[maxn];
14
15 int spfa(int s)
16 {
17     for(int i=1;i<=n;i++){
18             inq[i]=0;
19         dis[i]=0;
20     }
21
22     queue<int> q;
23     while(!q.empty()) q.pop();
24     dis[s]=1;
25     q.push(s);
26     inq[s]=1;
27     while(!q.empty())
28     {
29         int x=q.front();
30         q.pop();
31         inq[x]=0;
32         for(int i=1;i<=n;i++)
33             if(dis[i]<dis[x]*pic[x][i])
34         {
35             dis[i]=dis[x]*pic[x][i];
36             if(dis[s]>1.0) return 1;
37             if(!inq[i])
38             {
39                 inq[i]=1;
40                 q.push(i);
41             }
42         }
43     }
44     return 0;
45
46 }
47 int main()
48 {
49     int cas=0;
50     char s[100],s1[100],s2[100];
51     while(scanf("%d",&n)!=EOF&&n)
52     {
53         map<string,int> mapp;
54         mapp.clear();
55         for(int i=1;i<=n;i++){
56             scanf("%s",s);
57             mapp[s]=i;
58         }
59         for(int i=1;i<=n;i++)
60             for(int j=1;j<=n;j++)
61             pic[i][j]=(i==j?1:0);
62         scanf("%d",&m);
63         for(int i=0;i<m;i++)
64         {
65             scanf("%s%lf%s",s1,&r,s2);
66             u=mapp[s1];
67             v=mapp[s2];
68             pic[u][v]=r;
69         }
70         int ok=0;
71         for(int i=1;i<=n;i++)
72         {
73             if(spfa(i)==1) {ok=1;break;}
74         }
75         printf("Case %d: ",++cas);
76         if(ok) puts("Yes");
77         else puts("No");
78     }
79 }