poj1860-Currency Exchange

第一个单源最短路径算法，因为题目可能存在负边，所以用Bellman Ford算法，原始Bellman Ford可以用来求负环，这题需要改进一下用来求正环。

前几次试验一直被图的数据结构所困扰：一开始采用链表，因为BF算法要遍历所有边，而链表对边的遍历占速度优势；后来发现链表操作太繁琐，总共100x100的矩阵，用数组也行，于是想用数组；后来参考了一下大牛们的代码，发现其实用一维数组就能解决问题了。

数据结构：用value[]储存每个点（exchange point）经过exchange操作之后可以转换得到的值，初值统设为0（求负环相应设为正无穷）。rate[][]和com[][]用来储存汇率和回扣。

算法：对原始的Bellman Ford算法进行小改进即可。

另外从discuss中可以知道需要注意精度。

以下是AC代码：Memory:308K，Time:0MS

 1 #include <stdio.h>
 2 #include <string.h>
 3 #define MAX 101
 4 #define eps 1e-8
 5 int N, M, S;
 6 double value[MAX];
 7 double rate[MAX][MAX];
 8 double com[MAX][MAX];
 9 
 int Bellman_Ford()
 {
     int i, j, k, relaxed;
     for (k=0;k<N-1;k++){//N-1 passes
         relaxed=0;
         for (i=1;i<=N;i++){
             for (j=1;j<=N;j++){
                 if (i!=j && rate[i][j]>eps && value[j]+eps<(value[i]-com[i][j])*rate[i][j]){
                     //RELAX
                     value[j]=(value[i]-com[i][j])*rate[i][j];
                     relaxed=1;
                 }
             }
         }
         //如果没有松弛操作,即不存在正环
         if (!relaxed){
             return 0;
         }
     }
     //查找正环
     for (i=1;i<=N;i++){
         for (j=1;j<=N;j++){
             if (i!=j && rate[i][j]>eps && value[j]+eps<(value[i]-com[i][j])*rate[i][j]){
                 return 1;        
             }
         }
     }
     return 0;
 }
 int main()
 {
     int a, b, i;
     double Rab, Cab, Rba, Cba, v;
     //freopen("input.txt", "r", stdin);
     while (scanf("%d %d %d %lf", &N, &M, &S, &v)!=EOF){
         value[S]=v;
         for (i=0;i<M;i++){
             scanf("%d %d %lf %lf %lf %lf", &a, &b, &Rab, &Cab, &Rba, &Cba);
             rate[a][b]=Rab;
             rate[b][a]=Rba;
             com[a][b]=Cab;
             com[b][a]=Cba;
         }
         if (Bellman_Ford()){
             printf("YES\n");
         }
         else{
             printf("NO\n");
         }
     }
     return 0;
 }

END