結果
問題 | No.1265 Balloon Survival |
ユーザー |
|
提出日時 | 2020-10-23 22:29:48 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 311 ms / 2,000 ms |
コード長 | 3,846 bytes |
コンパイル時間 | 1,841 ms |
コンパイル使用メモリ | 174,184 KB |
実行使用メモリ | 19,916 KB |
最終ジャッジ日時 | 2024-07-21 11:18:00 |
合計ジャッジ時間 | 6,654 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 4 |
other | AC * 32 |
ソースコード
#include<bits/stdc++.h> using namespace std; #define rep(i,n) for(ll i=0;i<n;i++) #define repl(i,l,r) for(ll i=(l);i<(r);i++) #define per(i,n) for(ll i=n-1;i>=0;i--) #define perl(i,r,l) for(ll i=r-1;i>=l;i--) #define fi first #define se second #define pb push_back #define ins insert #define pqueue(x) priority_queue<x,vector<x>,greater<x>> #define all(x) (x).begin(),(x).end() #define CST(x) cout<<fixed<<setprecision(x) #define rev(x) reverse(x); using ll=long long; using vl=vector<ll>; using vvl=vector<vector<ll>>; using pl=pair<ll,ll>; using vpl=vector<pl>; using vvpl=vector<vpl>; const ll MOD=1000000007; const ll MOD9=998244353; const int inf=1e9+10; const ll INF=4e18; const ll dy[8]={1,0,-1,0,1,1,-1,-1}; const ll dx[8]={0,-1,0,1,1,-1,1,-1}; template<class T> inline bool chmin(T& a, T b) { if (a > b) { a = b; return true; } return false; } template<class T> inline bool chmax(T& a, T b) { if (a < b) { a = b; return true; } return false; } template <class T> class SegTree { int n;// 葉の数 vector<T> data;// データを格納するvector T def; // 初期値かつ単位元 function<T(T, T)> operation; // 区間クエリで使う処理 function<T(T, T)> update;// 点更新で使う処理 T find(int a, int b) { T val_left = def, val_right = def; for (a += (n - 1), b += (n - 1); a < b; a >>= 1, b >>= 1) { if ((a & 1) == 0){ val_left = operation(val_left, data[a]); } if ((b & 1) == 0){ val_right = operation(data[--b],val_right); } } return operation(val_left, val_right); } public: // _n:必要サイズ, _def:初期値かつ単位元, _operation:クエリ関数, // _update:更新関数 SegTree(size_t _n, T _def, function<T(T, T)> _operation, function<T(T, T)> _update=[](T a,T b){return b;}) : def(_def), operation(_operation), update(_update) { n = 1; while (n < _n) { n *= 2; } data = vector<T>(2 * n - 1, def); } void set(int i, T x) { data[i + n - 1] = x; } void build() { for (int k=n-2;k>=0;k--) data[k] = operation(data[2*k+1],data[2*k+2]); } // 場所i(0-indexed)の値をxで更新 void change(int i, T x) { i += n - 1; data[i] = update(data[i], x); while (i > 0) { i = (i - 1) / 2; data[i] = operation(data[i * 2 + 1], data[i * 2 + 2]); } } // [a, b)の区間クエリを実行 T query(int a, int b) { //return _query(a, b, 0, 0, n); return find(a,b); } // 添字でアクセス T operator[](int i) { return data[i + n - 1]; } }; int main(){ ll n;cin >> n; vector<ll> x(n),y(n); rep(i,n)cin >> x[i] >> y[i]; SegTree<ll> st(n*n,INF,[](ll a,ll b){return min(a,b);}); //vector<vector<double>> g(n,vector<double>(n)); rep(i,n){ rep(j,n){ ll len=(y[i]-y[j])*(y[i]-y[j])+(x[i]-x[j])*(x[i]-x[j]); if(i==j)len=INF; st.set(i*n+j,len); } } st.build(); ll ans=0; //rep(i,n*n)cout << st[i] <<endl; while(true){ ll p=st.query(0,n*n); if(p==INF)cout << ans <<endl,exit(0); ll ng=0,ok=n*n; while(ok-ng>1){ ll mid=(ok+ng)/2; if(st.query(0,mid)==p)ok=mid; else ng=mid; } ok--; //cout << ok <<endl; //break; ll l=ok/n,r=ok%n; if(l>r)swap(l,r); //cout << l <<" " << r <<endl; if(l==0){ ans++;l=r; } rep(i,n)st.change(l*n+i,INF); rep(i,n)st.change(n*i+l,INF); rep(i,n)st.change(r*n+i,INF); rep(i,n)st.change(n*i+r,INF); } }