結果

問題 No.96 圏外です。
ユーザー mamekinmamekin
提出日時 2015-01-17 20:29:52
言語 C++11
(gcc 11.4.0)
結果
AC  
実行時間 243 ms / 5,000 ms
コード長 3,332 bytes
コンパイル時間 1,184 ms
コンパイル使用メモリ 112,512 KB
実行使用メモリ 103,116 KB
最終ジャッジ日時 2024-06-06 15:55:05
合計ジャッジ時間 6,656 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge3
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 81 ms
97,536 KB
testcase_01 AC 80 ms
97,664 KB
testcase_02 AC 80 ms
97,536 KB
testcase_03 AC 79 ms
97,536 KB
testcase_04 AC 82 ms
97,792 KB
testcase_05 AC 85 ms
97,920 KB
testcase_06 AC 86 ms
97,792 KB
testcase_07 AC 89 ms
98,176 KB
testcase_08 AC 95 ms
98,304 KB
testcase_09 AC 100 ms
98,432 KB
testcase_10 AC 105 ms
98,636 KB
testcase_11 AC 117 ms
98,932 KB
testcase_12 AC 134 ms
99,200 KB
testcase_13 AC 144 ms
99,328 KB
testcase_14 AC 160 ms
99,968 KB
testcase_15 AC 158 ms
100,096 KB
testcase_16 AC 218 ms
100,992 KB
testcase_17 AC 216 ms
102,016 KB
testcase_18 AC 210 ms
101,760 KB
testcase_19 AC 205 ms
101,760 KB
testcase_20 AC 181 ms
100,232 KB
testcase_21 AC 171 ms
103,116 KB
testcase_22 AC 243 ms
101,436 KB
testcase_23 AC 228 ms
101,436 KB
testcase_24 AC 77 ms
97,536 KB
testcase_25 AC 170 ms
100,608 KB
testcase_26 AC 198 ms
101,504 KB
testcase_27 AC 185 ms
101,120 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <cstdio>
#include <iostream>
#include <sstream>
#include <fstream>
#include <iomanip>
#include <algorithm>
#include <cmath>
#include <string>
#include <vector>
#include <list>
#include <queue>
#include <stack>
#include <set>
#include <map>
#include <bitset>
#include <numeric>
#include <limits>
#include <climits>
#include <cfloat>
#include <functional>
using namespace std;

class Point
{
public:
    int y, x;
    Point(){
        y = x = 0;
    }
    Point(int y0, int x0){
        y = y0;
        x = x0;
    }
    Point operator+(const Point& p) const{
        return Point(y + p.y, x + p.x);
    }
    Point operator-(const Point& p) const{
        return Point(y - p.y, x - p.x);
    }
    Point operator*(int a) const{
        return Point(y * a, x * a);
    }
    int length2() const{
        return y * y + x * x;
    }
    int dist2(const Point& p) const{
        return (y - p.y) * (y - p.y) + (x - p.x) * (x - p.x);
    }
    int dot(const Point& p) const{
        return y * p.y + x * p.x; // |a|*|b|*cosθ
    }
    int cross(const Point& p) const{
        return x * p.y - y * p.x; // |a|*|b|*sinθ
    }
};

void convexHull(const vector<Point>& p0, vector<Point>& cv)
{
    class Compare{
    public:
        bool operator()(const Point& a, const Point& b){
            return a.x < b.x || (a.x == b.x && a.y < b.y);
        }
    };

    vector<Point> p = p0;
    sort(p.begin(), p.end(), Compare());
    int n = p.size();
    if(n == 1){
        cv = p0;
        return;
    }

    int j = 0;
    cv.clear();
    for(int i=0; i<2*n-1; ++i){
        Point tmp = (i<n? p[i]:p[2*(n-1)-i]);
        while(j >= 2 && (tmp-cv[j-2]).cross(cv[j-1]-cv[j-2]) < 0){
            cv.pop_back();
            -- j;
        }
        cv.push_back(tmp);
        ++ j;
    }
    cv.pop_back();
}

int main()
{
    int n;
    cin >> n;
    vector<Point> station(n);
    vector<vector<vector<int> > > grid(2003, vector<vector<int> >(2003));
    for(int i=0; i<n; ++i){
        int y, x;
        cin >> x >> y;
        y += 10010;
        x += 10010;
        station[i].y = y;
        station[i].x = x;
        grid[y/10][x/10].push_back(i);
    }

    vector<bool> check(n, false);
    double ret = 1.0;
    for(int i=0; i<n; ++i){
        if(check[i])
            continue;

        check[i] = true;
        queue<int> q;
        q.push(i);
        vector<Point> select(1, station[i]);
        while(!q.empty()){
            Point p = station[q.front()];
            q.pop();

            int y = p.y / 10;
            int x = p.x / 10;
            for(int i=0; i<9; ++i){
                for(unsigned j=0; j<grid[y+i/3-1][x+i%3-1].size(); ++j){
                    int k = grid[y+i/3-1][x+i%3-1][j];
                    if(!check[k] && p.dist2(station[k]) <= 100){
                        check[k] = true;
                        q.push(k);
                        select.push_back(station[k]);
                    }
                }
            }
        }

        vector<Point> cv;
        convexHull(select, cv);

        int m = cv.size();
        int j = 0;
        for(int i=0; i<m; ++i){
            while(cv[i].dist2(cv[j]) < cv[i].dist2(cv[(j+1)%m]))
                j = (j + 1) % m;
            ret = max(ret, sqrt((double)cv[i].dist2(cv[j])) + 2.0);
        }
    }
    printf("%.10f\n", ret);

    return 0;
}
0