算法导论--学习笔记019

1.套汇是指利用货币汇兑率的差异将一个单位的某种货币转换为大于一个单位的同种货币的方法。

例如，假定1美元可以买0.7英镑，1英镑可以买9.5法郎，1法郎可以买到0.16美元。通过货币兑换，一个商人可以从1美元开始买入，得到0.7*9.5*0.16=1.064美元，从而获得6.4%的利润。

假定已知n种货币1c,2c,„,nc和有关兑换率的n*n的表R。其中，R[i,j]是一个单位货币ic可以买到的货币jc的单位数。

(1)试设计一个有效算法，用以确定是否存在一个货币序列1ic,2ic,„,ikc使得R[i1,i2]R[i2,i3]„R[ik,i1] >1 并分析算法的计算时间。

(2)试设计一个算法打印出满足本题(1)中条件的所有序列，并分析算法的计算时间。

二、算法问题

形式化表示此问题可转化为图问题。构造一个图G=(V,E)，其中顶点为ci，任意两点间都有两条有向边(vi, vj), (vj, vi)。则有|V| = n, |E| = n(n-1)/2。从题目条件可得：因此可设边(ci1,ci2) 权值为-lgR[i1, i2]，这样问题转化为在图中寻找负权回路的问题。可以现在图中添加一个新的顶点s, s指向所有入度为0的顶点，权值为0；然后以s为起点运行Bellman-Ford算法，即可确定是否有负权回路。

-----------------------------------------------我是分割线----------------------------------------------------------

为了防止非法套汇：

汇率兑换矩阵R=(rij)n*n是正反矩阵：

rii=1;(i=0,1,2,3,4.......);

rij=1/rji;

rij=rik/rjk

具有唯一非零的最大特征值。Xmax=n;

------------------------------------我是分割线------------------------------------------------------

[cpp] view plaincopy

1. #include <iostream>  
2. #include <cstdio>  
3. #include <algorithm>  
4. #include <cstring>  
5. using namespace std;  
6.   
7. const int N = 35;  
8. const int E = 905;  
9. const double MAX = 999999999;  
10.   
11. struct Edge  
12. {  
13.     int beg;  
14.     int end;  
15.     double rate;  
16. }edge[E];  
17. double mindis[N];  
18. int n, e;  
19. struct Name  
20. {  
21.     char str[55];  
22. }cname[N];  
23.   
24. bool cmp(const Name& n1, const Name& n2)  
25. {  
26.     return strcmp(n1.str, n2.str) < 0;  
27. }  
28.   
29. int bisearch(char* str)  
30. {  
31.     int low = 0, high = n - 1, mid;  
32.     int t;  
33.     while (low <= high)  
34.     {  
35.         mid = (low + high) >> 1;  
36.         t = strcmp(str, cname[mid].str);  
37.         if (t == 0) return mid;  
38.         else if (t < 0) high = mid - 1;  
39.         else low = mid + 1;  
40.     }  
41.     return -1;  
42. }  
43.   
44. bool relax(int pe)  
45. {  
46.     if (mindis[edge[pe].beg] * edge[pe].rate > mindis[edge[pe].end])   
47.     {  
48.         mindis[edge[pe].end] = mindis[edge[pe].beg] * edge[pe].rate;  
49.         return true;  
50.     }  
51.     return false;  
52. }  
53.   
54. bool bellman_ford()  
55. {  
56.     bool flag;  
57.     for (int i = 0; i < n; ++i) mindis[i] = 1.0;  
58.     for (int i = 0; i < n - 1; ++i)  
59.     {  
60.         flag = false;  
61.         for (int j = 0; j < e; ++j)  
62.         {  
63.             if (relax(j)) flag = true;  
64.         }  
65.         if (!flag) return false;  
66.     }  
67.     for (int i = 0; i < e; ++i)  
68.     {  
69.         if (relax(i)) return true;  
70.     }  
71.     return false;  
72. }  
73.   
74. int main()  
75. {  
76.     char sbeg[55], send[55];  
77.     double rate;  
78.     int T = 1;  
79.     while (scanf("%d", &n) != EOF && n)  
80.     {  
81.         for (int i = 0; i < n; ++i) scanf("%s", cname[i].str);  
82.         sort(cname, cname + n, cmp);  
83.         scanf("%d", &e);  
84.         for (int i = 0; i < e; ++i)  
85.         {  
86.             scanf("%s%lf%s", sbeg, &rate, send);  
87.             edge[i].beg = bisearch(sbeg);  
88.             edge[i].end = bisearch(send);  
89.             edge[i].rate = rate;  
90.         }  
91.         printf("Case %d: ", T++);  
92.         if (bellman_ford()) printf("Yes\n");  
93.         else printf("No\n");  
94.     }  
95.     return 0;  
96. }