結果

問題 No.96 圏外です。
ユーザー kmjp
提出日時 2014-12-07 21:10:51
言語 C++11
(gcc 13.3.0)
結果
AC  
実行時間 161 ms / 5,000 ms
コード長 3,623 bytes
コンパイル時間 1,923 ms
コンパイル使用メモリ 182,620 KB
実行使用メモリ 18,908 KB
最終ジャッジ日時 2024-12-24 07:46:43
合計ジャッジ時間 5,234 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 28
権限があれば一括ダウンロードができます

ソースコード

diff #
プレゼンテーションモードにする

#include <bits/stdc++.h>
using namespace std;
typedef signed long long ll;
#undef _P
#define _P(...) (void)printf(__VA_ARGS__)
#define FOR(x,to) for(x=0;x<to;x++)
#define ITR(x,c) for(__typeof(c.begin()) x=c.begin();x!=c.end();x++)
#define ALL(a) (a.begin()),(a.end())
#define ZERO(a) memset(a,0,sizeof(a))
#define MINUS(a) memset(a,0xff,sizeof(a))
//-------------------------------------------------------
class UF {
    public:
    static const int ufmax=120200;
    int ufpar[ufmax],ufrank[ufmax],ufcnt[ufmax];
    UF() { init();}
    void init(){int i; FOR(i,ufmax) { ufpar[i]=i; ufrank[i]=0; ufcnt[i]=1; } }
    int find(int x) {   return (ufpar[x]==x)?(x):(ufpar[x] = find(ufpar[x]));}
    int operator[](int x) {return find(x);}
    int count(int x) {return ufcnt[find(x)];}
    void unite(int x,int y) {
        x = find(x); y = find(y);
        if(x==y) return;
        if(ufrank[x]<ufrank[y]) ufpar[x]=y, ufcnt[y]+=ufcnt[x];
        else {ufpar[y]=x; ufcnt[x]+=ufcnt[y]; if(ufrank[x]==ufrank[y]) ufrank[x]++;}
    }
};
ll veccross(pair<int,int> p1,pair<int,int> p2,pair<int,int> p3) {
    p3.first-=p1.first;p2.first-=p1.first;
    p3.second-=p1.second;p2.second-=p1.second;
    return p3.first*(ll)p2.second-p2.first*(ll)p3.second;
}
vector<int> convex_hull(vector< pair<int, int> >& vp) {
    vector<pair<pair<int, int>, int> > sorted;
    vector<int> res;
    int i,k=0,rb;
    
    if(vp.size()<=2) {
        if(vp.size()>=1) res.push_back(0);
        if(vp.size()>=2) res.push_back(1);
        return res;
    }
    
    FOR(i,vp.size()) sorted.push_back(make_pair(vp[i],i));
    sort(sorted.begin(),sorted.end());
    res.resize(vp.size()*2);
    /* bottom */
    FOR(i,vp.size()) {
        while(k>1 && veccross(vp[res[k-2]],vp[res[k-1]],sorted[i].first)<=0) k--;
        res[k++]=sorted[i].second;
    }
    /* top */
    for(rb=k, i=vp.size()-2;i>=0;i--) {
        while(k>rb && veccross(vp[res[k-2]],vp[res[k-1]],sorted[i].first)<=0) k--;
        res[k++]=sorted[i].second;
    }
    res.resize(k-1);
    return res;
}
int N;
int X[120001],Y[120001];
UF uf;
map<int,int> YS[20015];
vector<pair<int,int> > V[120001];
void solve() {
    int i,j,k,l,r,x,y,y2,x2;
    
    cin>>N;
    FOR(i,N) {
        cin>>X[i]>>Y[i], X[i]+=10001,Y[i]+=10001;
        YS[X[i]][Y[i]]=i;
    }
    
    FOR(x,20002) {
        ITR(it,YS[x]) {
            int cur=it->first;
            map<int,int>::iterator it2=it;
            for(it2++;it2!=YS[x].end();it2++) {
                if(it2->first-it->first>10) break;
                uf.unite(it->second,it2->second);
            }
            for(x2=x+1;x2<=x+10;x2++) {
                for(map<int,int>::iterator it2=YS[x2].lower_bound(it->first-10); it2!=YS[x2].end();it2++) {
                    if(it2->first<it->first && (x2-x)*(x2-x)+(it->first-it2->first)*(it->first-it2->first)>100) continue;
                    if((x2-x)*(x2-x)+(it->first-it2->first)*(it->first-it2->first)>100) break;
                    uf.unite(it->second,it2->second);
                }
            }
        }
    }
    
    FOR(i,N) V[uf[i]].push_back(make_pair(X[i],Y[i]));
    
    double ma=1;
    FOR(i,N) if(V[i].size()) {
        if(V[i].size()==1) ma=max(ma,2.0);
        else if(V[i].size()==2) ma=max(ma,2+sqrt((V[i][0].first-V[i][1].first)*(V[i][0].first-V[i][1].first)+(V[i][0].second-V[i][1].second)*(V[i][0]
            .second-V[i][1].second)));
        else {
            vector<int> ch=convex_hull(V[i]);
            FOR(x,ch.size()) for(y=x+1;y<ch.size();y++)
                ma=max(ma,2+sqrt((V[i][ch[x]].first-V[i][ch[y]].first)*(V[i][ch[x]].first-V[i][ch[y]].first)+(V[i][ch[x]].second-V[i][ch[y]].second
                    )*(V[i][ch[x]].second-V[i][ch[y]].second)));
        }
    }
    _P("%.12lf\n",ma);
}
int main(int argc,char** argv){
    string s;int i;
    if(argc==1) ios::sync_with_stdio(false);
    FOR(i,argc-1) s+=argv[i+1],s+='\n';
    FOR(i,s.size()) ungetc(s[s.size()-1-i],stdin);
    solve(); return 0;
}
הההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההה
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
0