erfanul007

UVa 104

Apr 15th, 2021
607
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
  1. #include <bits/stdc++.h>
  2. // #include <iostream>
  3. // #include <cstdio>
  4. // #include <cstdlib>
  5. // #include <algorithm>
  6. // #include <cmath>
  7. // #include <vector>
  8. // #include <set>
  9. // #include <map>
  10. // #include <queue>
  11. // #include <ctime>
  12. // #include <cassert>
  13. // #include <complex>
  14. // #include <string>
  15. // #include <cstring>
  16. // #include <queue>
  17. // #include <bitset>
  18.  
  19. using namespace std;
  20.  
  21. // #pragma GCC optimize("Ofast,no-stack-protector")
  22. // #pragma GCC target("sse,sse2,sse3,ssse3,sse4,popcnt,abm,mmx,avx,tune=native")
  23. // #pragma GCC optimize("unroll-loops")
  24.  
  25.  
  26. #define ll              long long int
  27. #define vi              vector< int >
  28. #define vll             vector< ll >
  29.  
  30. #define sc              scanf
  31. #define pf              printf
  32. #define cspf(i)         pf("Case #%d\n", i)
  33. #define spc             pf(" ")
  34. #define line            pf("\n")
  35.  
  36. #define ff              first
  37. #define ss              second
  38. #define mp              make_pair
  39. #define pb              push_back
  40. #define ppb             pop_back
  41. #define tp(v,j)         get<j>(v)
  42. #define Log(b,x)        (log(x)/log(b))
  43.  
  44. #define FOR(i,x,y)      for(int i = int(x); i < int(y); i++)
  45. #define ROF(i,x,y)      for(int i = int(x)-1; i >= int(y); i--)
  46. #define clr(arr,x)      memset(arr, x, sizeof arr)
  47. #define vout(v)         for(int w=0;w<(int)v.size();w++){if(w) spc; cout<<v[w];}
  48. #define all(v)          v.begin(), v.end()
  49. #define rall(v)         v.rbegin(), v.rend()
  50. #define unq(v)          sort(all(v)),(v).resize(unique(all(v))-v.begin())
  51. #define fastIO          ios::sync_with_stdio(false); cin.tie(nullptr); cout.tie(nullptr)
  52.  
  53. #define sc1(x)          sc("%d",&x)
  54. #define sc2(x,y)        sc("%d %d", &x, &y)
  55. #define sc3(x,y,z)      sc("%d %d %d", &x, &y, &z)
  56. #define scl1(x)         sc("%lld",&x);
  57. #define scl2(x,y)       sc("%lld %lld", &x, &y)
  58. #define scf1(x)         sc("%lf",&x)
  59. #define scf2(x,y)       sc("%lf %lf", &x, &y)
  60.  
  61. #define pf1(x)          pf("%d",x)
  62. #define pf2(x,y)        pf("%d %d", x, y)
  63. #define pf3(x,y,z)      pf("%d %d %d", x, y, z)
  64. #define pfl1(x)         pf("%lld",x)
  65. #define pfl2(x,y)       pf("%lld %lld",x,y)
  66.  
  67. #define MOD             10056
  68. #define MaxN            1000001
  69. #define MAX             (int)(1e9+7)
  70. #define inf             0x3f3f3f3f
  71. #define PI              acos(-1.0)  // 3.1415926535897932
  72. #define eps             1e-6
  73.  
  74. template <class T> inline T bigMod(T p,T e,T M){T ret=1; for(;e>0;e>>=1){ if(e&1) ret=(ret*p)%M; p=(p*p)%M;} return (T)ret;}
  75. template <class T> inline T modInverse(T a,T M){return bigMod(a,M-2,M);}
  76. template <class T> inline T gcd(T a,T b){if(b==0)return a;return gcd(b,a%b);}
  77. template <class T> inline T lcm(T a,T b) {a=abs(a);b=abs(b); return (a/gcd(a,b))*b;}
  78. template <class T> inline T SQR(T a){return a*a;}
  79.  
  80. int dx[] = { 1,-1, 0, 0};
  81. int dy[] = { 0, 0, 1,-1};
  82.  
  83. double dp[55][55][55], g[55][55];
  84. int path[55][55][55], n;
  85.  
  86. void init(){
  87.     for(int k=1; k<=n; k++){
  88.         for(int i=1; i<=n; i++){
  89.             for(int j=1; j<=n; j++){
  90.                 path[i][j][k]=0;
  91.                 dp[i][j][k]=0;
  92.             }
  93.         }
  94.     }
  95. }
  96.  
  97. void pathp(int u, int v, int k){
  98.     if(k==0){
  99.         pf("%d",u);
  100.         return;
  101.     }
  102.     pathp(u, path[u][v][k], k-1);
  103.     pf(" %d",v);
  104. }
  105.  
  106. double dpsolve(int u, int v, int k){
  107.     if(k==0){
  108.         if(u==v) return 1;
  109.         return 0;
  110.     }
  111.     if(path[u][v][k]) return dp[u][v][k];
  112.     double mx = dp[u][v][k];
  113.     for(int i=1; i<=n; i++){
  114.         double val = dpsolve(u, i, k-1) * g[i][v];
  115.         if(val>mx){
  116.             mx = val;
  117.             path[u][v][k]=i;
  118.         }
  119.     }
  120.     return dp[u][v][k]=mx;
  121. }
  122.  
  123. bool Floyd(){
  124.     for(int k=2; k<=n; k++){
  125.         for(int i=1; i<=n; i++){
  126.             double val = dpsolve(i,i,k);
  127.             if(val>=1.01){
  128.                 pathp(i, i, k);
  129.                 return false;
  130.             }
  131.         }
  132.     }
  133.     return true;
  134. }
  135.  
  136. int main(){
  137.  
  138.     #ifndef ONLINE_JUDGE
  139.         clock_t tStart = clock();
  140.         freopen("input.txt", "r", stdin);
  141.         freopen("output.txt", "w", stdout);
  142.     #endif
  143.  
  144.     while(cin>>n){
  145.         init();
  146.  
  147.         for(int i=1; i<=n; i++){
  148.             for(int j=1; j<=n; j++){
  149.                 if(i==j) continue;
  150.                 cin>>g[i][j];
  151.             }
  152.         }
  153.        
  154.         if(Floyd()) pf("no arbitrage sequence exists");
  155.         line;
  156.     }
  157.  
  158.     #ifndef ONLINE_JUDGE
  159.         fprintf(stderr, "\n>> Runtime: %.10fs\n", (double) (clock() - tStart) / CLOCKS_PER_SEC);
  160.     #endif
  161.  
  162.     return 0;
  163. }
RAW Paste Data

Adblocker detected! Please consider disabling it...

We've detected AdBlock Plus or some other adblocking software preventing Pastebin.com from fully loading.

We don't have any obnoxious sound, or popup ads, we actively block these annoying types of ads!

Please add Pastebin.com to your ad blocker whitelist or disable your adblocking software.

×