Advertisement
erfanul007

UVa 10245

Jun 27th, 2021
1,093
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
C++ 5.16 KB | None | 0 0
  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 <stack>
  12. // #include <ctime>
  13. // #include <cassert>
  14. // #include <complex>
  15. // #include <string>
  16. // #include <cstring>
  17. // #include <bitset>
  18. using namespace std;
  19.  
  20. // #pragma GCC optimize("Ofast,no-stack-protector")
  21. // #pragma GCC target("sse,sse2,sse3,ssse3,sse4,popcnt,abm,mmx,avx,tune=native")
  22. // #pragma GCC optimize("unroll-loops")
  23.  
  24. #define ll              long long int
  25. #define vi              vector< int >
  26. #define vll             vector< ll >
  27.  
  28. #define sc              scanf
  29. #define pf              printf
  30. #define cspf(i)         pf("Case %d:\n", i)
  31. #define spc             pf(" ")
  32. #define line            pf("\n")
  33.  
  34. #define ff              first
  35. #define ss              second
  36. #define mp              make_pair
  37. #define pb              push_back
  38. #define ppb             pop_back
  39. #define tp(v,j)         get<j>(v)
  40. #define Log(b,x)        (log(x)/log(b))
  41.  
  42. #define FOR(i,x,y)      for(int i = int(x); i < int(y); i++)
  43. #define ROF(i,x,y)      for(int i = int(x)-1; i >= int(y); i--)
  44. #define clr(arr,x)      memset(arr, x, sizeof arr)
  45. #define vout(v,sz)      for(int w=0;w<sz;w++){if(w) spc; cout<<v[w];}
  46. #define all(v)          v.begin(), v.end()
  47. #define rall(v)         v.rbegin(), v.rend()
  48. #define unq(v)          sort(all(v)),(v).resize(unique(all(v))-v.begin())
  49. #define fastIO          ios::sync_with_stdio(false); cin.tie(nullptr); cout.tie(nullptr)
  50.  
  51. #define sc1(x)          sc("%d",&x)
  52. #define sc2(x,y)        sc("%d %d", &x, &y)
  53. #define sc3(x,y,z)      sc("%d %d %d", &x, &y, &z)
  54. #define scl1(x)         sc("%lld",&x)
  55. #define scl2(x,y)       sc("%lld %lld", &x, &y)
  56. #define scf1(x)         sc("%lf",&x);
  57. #define scf2(x,y)       sc("%lf %lf", &x, &y)
  58.  
  59. #define pf1(x)          pf("%d",x)
  60. #define pf2(x,y)        pf("%d %d", x, y)
  61. #define pfl1(x)         pf("%lld\n",x)
  62. #define pfl2(x,y)       pf("%lld %lld", x, y)
  63.  
  64. #define MAX             1e9
  65. #define MaxN            10000
  66. #define inf             0x3f3f3f3f
  67. #define PI              acos(-1.0)  // 3.1415926535897932
  68. #define eps             1e-9
  69.  
  70. #ifdef ERFANUL007
  71.     #define debug(...) __f(#__VA_ARGS__, __VA_ARGS__)
  72.     template < typename Arg1 >
  73.     void __f(const char* name, Arg1&& arg1){
  74.         cout << name << " = " << arg1 << std::endl;
  75.     }
  76.     template < typename Arg1, typename... Args>
  77.     void __f(const char* names, Arg1&& arg1, Args&&... args){
  78.         const char* comma = strchr(names, ',');
  79.         cout.write(names, comma - names) << " = " << arg1 <<" | ";
  80.         __f(comma+1, args...);
  81.     }
  82. #else
  83.     #define debug(...)
  84. #endif
  85.  
  86. 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;}
  87. template <class T> inline T modInverse(T a,T M){return bigMod(a,M-2,M);}
  88. template <class T> inline T gcd(T a,T b){if(b==0)return a;return gcd(b,a%b);}
  89. template <class T> inline T lcm(T a,T b) {a=abs(a);b=abs(b); return (a/gcd(a,b))*b;}
  90.  
  91. int dx[] = { 1,-1, 0, 0};                //graph moves
  92. int dy[] = { 0, 0, 1,-1};               //graph moves
  93.  
  94. struct PT{
  95.     double x, y;
  96.     PT(){}
  97.     PT(double x, double y) : x(x), y(y) {}
  98.     void scan(){ scanf("%lf %lf",&x,&y);}
  99.     inline bool operator == (const PT &p) const {
  100.         return (abs(x - p.x) < eps && abs(y - p.y) < eps);
  101.     }
  102.     inline bool operator < (const PT &p) const {
  103.         return (((x + eps) < p.x) || (abs(x - p.x) < eps && (y + eps) < p.y));
  104.     }
  105. };
  106.  
  107. double SQ(double x){ return x*x;}
  108. /* return the euclidean distance of two point */
  109. double Dis(PT a, PT b){ return SQ(a.x-b.x) + SQ(a.y-b.y);}
  110.  
  111. double absDis(PT a, PT b){ return sqrt(Dis(a, b));}
  112.  
  113. /* receive a list of sorted by x point
  114. find the closest distance of any pair of points recursively
  115. O(n(log n)^2) */
  116. vector< PT > v;
  117. double ClosestDistance(int lo, int hi){
  118.     if(lo == hi) return MAX;
  119.     if(lo + 1 == hi) return absDis(v[lo], v[hi]);
  120.  
  121.     int mid = (lo + hi)/2;
  122.     double val = min( ClosestDistance(lo, mid), ClosestDistance(mid+1, hi));
  123.  
  124.     vector< PT > strip;
  125.     for(int i=lo; i<=hi; i++){
  126.         if((abs(v[mid].x - v[i].x) + eps) < val){
  127.             strip.pb(v[i]);
  128.         }
  129.     }
  130.  
  131.     sort(strip.begin(), strip.end(), [](const PT &a, const PT &b) -> bool {
  132.         return (a.y + eps) < b.y;
  133.     }); // sort by y coordinates
  134.  
  135.     for(int i=0; i<strip.size(); i++){
  136.         for(int j=i+1; j<strip.size() && (abs(strip[i].y - strip[j].y) + eps) < val; j++){
  137.             val = min(val, absDis(strip[i], strip[j]));
  138.         }
  139.     }
  140.     return val;
  141. }
  142.  
  143.  
  144. int main()
  145. {
  146.     #ifdef ERFANUL007
  147.         clock_t tStart = clock();
  148.         freopen("input.txt", "r", stdin);
  149.         freopen("output.txt", "w", stdout);
  150.     #endif
  151.  
  152.     int n;
  153.     while(~sc1(n) && n){
  154.         v.clear();
  155.         for(int i=0; i<n; i++){
  156.             PT a; a.scan();
  157.             v.pb(a);
  158.         }
  159.         sort(all(v));
  160.         double ans = ClosestDistance(0, n-1);
  161.         if(abs(ans - MaxN) < 1e-4 || (MaxN + eps) < ans) pf("INFINITY\n");
  162.         else pf("%.4f\n", ans);
  163.     }
  164.  
  165.     #ifdef ERFANUL007
  166.         fprintf(stderr, "\n>> Runtime: %.10fs\n", (double) (clock() - tStart) / CLOCKS_PER_SEC);
  167.     #endif
  168.  
  169.     return 0;
  170. }
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement