結果

問題 No.96 圏外です。
ユーザー yaoshimaxyaoshimax
提出日時 2015-03-08 23:28:53
言語 C++11
(gcc 11.4.0)
結果
WA  
実行時間 -
コード長 4,285 bytes
コンパイル時間 1,557 ms
コンパイル使用メモリ 105,328 KB
実行使用メモリ 18,372 KB
最終ジャッジ日時 2024-06-24 16:36:34
合計ジャッジ時間 5,381 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge3
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 4 ms
8,516 KB
testcase_01 AC 5 ms
8,460 KB
testcase_02 AC 4 ms
8,488 KB
testcase_03 AC 4 ms
8,640 KB
testcase_04 AC 7 ms
8,896 KB
testcase_05 AC 8 ms
9,028 KB
testcase_06 WA -
testcase_07 AC 16 ms
9,484 KB
testcase_08 AC 20 ms
9,796 KB
testcase_09 WA -
testcase_10 AC 44 ms
10,648 KB
testcase_11 AC 52 ms
11,972 KB
testcase_12 WA -
testcase_13 AC 115 ms
12,708 KB
testcase_14 AC 126 ms
13,908 KB
testcase_15 AC 194 ms
14,532 KB
testcase_16 AC 204 ms
16,220 KB
testcase_17 AC 198 ms
18,372 KB
testcase_18 AC 274 ms
17,864 KB
testcase_19 AC 274 ms
17,860 KB
testcase_20 AC 150 ms
16,912 KB
testcase_21 AC 169 ms
16,460 KB
testcase_22 AC 241 ms
17,704 KB
testcase_23 AC 156 ms
17,840 KB
testcase_24 AC 5 ms
8,528 KB
testcase_25 WA -
testcase_26 WA -
testcase_27 AC 146 ms
16,360 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <cstdio>
#include <cstdlib>
#include <cmath>
#include <climits>
#include <cfloat>
#include <map>
#include <utility>
#include <set>
#include <iostream>
#include <memory>
#include <string>
#include <vector>
#include <algorithm>
#include <functional>
#include <sstream>
#include <complex>
#include <stack>
#include <queue>
#include <cstring>

using namespace std;

class UnionFindTree
{
   private:
      int N;
      int *rootnum;
      int *rank;
      int *unionnum;

   public:
      UnionFindTree(int n){
         /*0からn-1までの要素をunionとしてもつ*/
         N = n;
         rootnum = new int[N];
         rank = new int[N];
         unionnum = new int[N];
         for( int i = 0 ; i < N ; i++ ){
            rootnum[i] = i;
            rank[i] = 0;
            unionnum[i] = 1;
         }
      }
      
      int Find( int j ){
         if( rootnum[j] == j ){
            return j;
         }
         else{
            return rootnum[j] = Find(rootnum[j]);
         }
      }

      void Union( int i , int j ){
         int x = Find(i);
         int y = Find(j);
         if( x == y ) return;
         N--;
         if( rank[x] < rank[y] ){
            rootnum[y] = x;
            unionnum[x] += unionnum[y];
         }else{
            rootnum[x] = y;
            unionnum[y] += unionnum[x];
            if( rank[x]==rank[y] ) rank[x]++;
         }
      }

      int Num( int j ){
         /*与えられた要素の入っている集合の要素数を答える*/
         int x = Find(j);
         return unionnum[x];
      }

      int UnionNum(){
         /*集合の数を答える*/
         return N;
      }

      ~UnionFindTree(){
         delete rootnum;
         delete rank;
         delete unionnum;
      }

};

long long product( pair<int,int> p1, pair<int,int> p2, pair<int,int> p3 ){
   pair<int,int> v1 = make_pair( p2.first-p1.first, p2.second-p1.second);
   pair<int,int> v2 = make_pair( p3.first-p1.first, p3.second-p1.second);
   return v1.first*(long long)v2.second-v1.second*(long long)v2.first;
}

vector< pair<int,int> > convex_hull( vector<pair<int,int> > &V){
   if( V.size() <=2 ) return V;
   int size = V.size();
   sort(V.begin(),V.end());
   vector<pair<int,int> > ret(size*2);
   int retlen=0; 
   for( int i = 0 ; i <size; i++ ){
      while( retlen >= 2 && product( ret[retlen-2],ret[retlen-1],V[i] ) >= 0 ) retlen--;
      ret[retlen++]=V[i];
   }
   for( int i = size-1 ; i > 0 ; i-- ){
      while( retlen >= 2 && product( ret[retlen-2],ret[retlen-1],V[i] ) >= 0 ) retlen--;
      ret[retlen++]=V[i];
   }
   ret.resize(retlen);
   return ret;
}

pair<int,int> points[120001];
vector< pair<int,int > > V[120001];
map<int, int> m[20011];
UnionFindTree uft(120001);

int main(){
   int N;
   cin >> N;

   for( int i = 0 ; i < N ; i++ ){
      V[i].clear();
      cin >> points[i].first >> points[i].second;
      m[points[i].first+10000][points[i].second+10000]=i+1;
   }
   if( N == 0 ){
      printf("%.10f\n",1.0);
      return 0;
   }
   for( int i = 0 ; i < N ; i++ ){
      int x = points[i].first+10000;
      int y = points[i].second+10000;
      for( int x2 = x ; x2 <= x+10; x2++ ){
         map<int,int>::iterator it=m[x2].lower_bound(y-10);
         map<int,int>::iterator itEnd=m[x2].lower_bound(y+11);
         for( ;it != itEnd; it++ ){  
            int y2=it->first;
            int dx=x-x2;
            int dy=y-y2;
            if( dx*dx+dy*dy <= 100 ){
               //cout << i <<" "<< it->second-1<<endl;
               uft.Union(i,it->second-1);
            }
         }
      }
   }
   for( int i = 0 ; i < N ; i++ ){
      V[uft.Find(i)].push_back(points[i]);
   }
   
   long long maxDist = 0;
   for( int i = 0 ; i <N; i++ ){
      vector<pair<int,int> > convPoints = convex_hull(V[i]);  
      for( int j = 0 ; j < (int) convPoints.size(); j ++ ){
         //cout << convPoints[j].first << "-" << convPoints[j].second<<", ";
         for( int k = j+1 ; k < (int) convPoints.size(); k++ ){
            long long dx = convPoints[j].first-convPoints[k].first;
            long long dy = convPoints[j].second-convPoints[k].second;
            maxDist = max(maxDist,dx*dx+dy*dy);
         }
      }
      //cout << endl;
   }
   printf("%.10f\n",2.0+sqrt((double)maxDist));
   return 0;
}
0