Arbitrage
Time Limit: 1000MS Memory Limit: 65536K
Total Submissions: 35390 Accepted: 14545
Description
Arbitrage is the use of discrepancies in currency exchange rates to transform one unit of a currency into more than one unit of the same currency. For example, suppose that 1 US Dollar buys 0.5 British pound, 1 British pound buys 10.0 French francs, and 1 French franc buys 0.21 US dollar. Then, by converting currencies, a clever trader can start with 1 US dollar and buy 0.5 * 10.0 * 0.21 = 1.05 US dollars, making a profit of 5 percent.
Your job is to write a program that takes a list of currency exchange rates as input and then determines whether arbitrage is possible or not.
Input
The input will contain one or more test cases. Om the first line of each test case there is an integer n (1<=n<=30), representing the number of different currencies. The next n lines each contain the name of one currency. Within a name no spaces will appear. The next line contains one integer m, representing the length of the table to follow. The last m lines each contain the name ci of a source currency, a real number rij which represents the exchange rate from ci to cj and a name cj of the destination currency. Exchanges which do not appear in the table are impossible.
Test cases are separated from each other by a blank line. Input is terminated by a value of zero (0) for n.
Output
For each test case, print one line telling whether arbitrage is possible or not in the format “Case case: Yes” respectively “Case case: No”.
Sample Input
3
USDollar
BritishPound
FrenchFranc
3
USDollar 0.5 BritishPound
BritishPound 10.0 FrenchFranc
FrenchFranc 0.21 USDollar
3
USDollar
BritishPound
FrenchFranc
6
USDollar 0.5 BritishPound
USDollar 4.9 FrenchFranc
BritishPound 10.0 FrenchFranc
BritishPound 1.99 USDollar
FrenchFranc 0.09 BritishPound
FrenchFranc 0.19 USDollar
0
Sample Output
Case 1: Yes
Case 2: No
spfa算法:(判断是否存在正环)
#include <iostream>
#include <cstdio>
#include <map>
#include <cstring>
#include <queue>
using namespace std;
const int N = 50, M = 1e5;
int n,m;
int h[N],e[M],ne[M],idx;
double w[M];
bool st[N];
int cnt[N];
double dist[N];
map<string,int> heap;
void add(int a,int b,double c){
e[idx] = b; w[idx] = c; ne[idx] = h[a]; h[a] = idx++;
}
bool spfa(){
queue<int> q;
for(int i = 1; i <= n; i++){
dist[i] = 1;
cnt[i] = 0;
q.push(i);
st[i] = true;
}
while(q.size()){
int t = q.front(); q.pop();
st[t] = false;
for(int i = h[t]; ~i; i = ne[i]){
int j = e[i];
if(dist[j] < dist[t] * w[i]){
dist[j] = dist[t] * w[i];
cnt[j] = cnt[t] + 1;
if(cnt[j] >= n) return true;
if(!st[j]){
q.push(j);
st[j] = true;
}
}
}
}
return false;
}
int main(int argc, char** argv) {
int T = 1;
while(cin>>n,n){
memset(h,-1,sizeof h);
for(int i = 1; i <= n; i++){
string s;
cin>>s;
heap[s] = i;
}
cin>>m;
for(int i = 0; i < m; i++){
string s1,s2;
double c;
cin>>s1>>c>>s2;
add(heap[s1],heap[s2],c);
}
if(spfa()) printf("Case %d: Yes\n",T++);
else printf("Case %d: No\n",T++);
}
return 0;
}
floyd算法:
#include <iostream>
#include <cstdio>
#include <map>
using namespace std;
const int N = 110;
int n,m;
double dist[N][N];
map<string,int> heap;
bool floyd(){
for(int k = 1; k <= n; k++)
for(int i = 1; i <= n; i++)
for(int j = 1; j <= n; j++)
if(dist[i][j] < dist[i][k] * dist[k][j]){
dist[i][j] = dist[i][k] * dist[k][j];
}
for(int i = 1; i <= n; i++)
if(dist[i][i] > 1) return true;
return false;
}
int main(int argc, char** argv) {
int T = 1;
while(cin>>n,n){
for(int i = 1; i <= n; i++){
string s;
cin>>s;
heap[s] = i;
dist[i][i] = 1; //初始化自身汇率为1
}
cin>>m;
for(int i = 0; i < m; i++){
string s1,s2;
double c;
cin>>s1>>c>>s2;
dist[heap[s1]][heap[s2]] = c;
}
if(floyd()) printf("Case %d: Yes\n",T++);
else printf("Case %d: No\n",T++);
}
return 0;
}
Bellman-ford算法:
当最后一次遍历m条边时,如果还能更新成功,那么说明存在最长路。
#include <iostream>
#include <cstdio>
#include <map>
#include <cstring>
using namespace std;
const int N = 10010;
struct Edge{
int a,b;
double c;
}edges[N];
int n,m;
double dist[N],backup[N];
map<string,int> heap;
bool Bellman_ford(){
for(int i = 1; i <= n; i++) dist[i] = 1;
for(int i = 0; i < n; i++){
memcpy(backup,dist,sizeof dist);
for(int j = 0; j < m; j++){
Edge e = edges[j];
if(dist[e.b] < backup[e.a] * e.c){
dist[e.b] = backup[e.a] * e.c;
if(i == n - 1) return true;
}
}
}
return false;
}
int main(int argc, char** argv) {
int T = 1;
while(cin>>n,n){
for(int i = 1; i <= n; i++){
string s;
cin>>s;
heap[s] = i;
}
cin>>m;
for(int i = 0; i < m; i++){
string s1,s2;
double c;
cin>>s1>>c>>s2;
edges[i] = {heap[s1],heap[s2],c};
}
if(Bellman_ford()) printf("Case %d: Yes\n",T++);
else printf("Case %d: No\n",T++);
}
return 0;
}