結果

問題 No.2353 Guardian Dogs in Spring
ユーザー Taiki0715Taiki0715
提出日時 2023-06-16 21:48:16
言語 C++17(clang)
(17.0.6 + boost 1.87.0)
結果
MLE  
実行時間 -
コード長 4,752 bytes
コンパイル時間 5,710 ms
コンパイル使用メモリ 166,132 KB
実行使用メモリ 535,152 KB
最終ジャッジ日時 2024-06-24 13:43:12
合計ジャッジ時間 8,007 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 2
other AC * 3 MLE * 1 -- * 36
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <iostream>
#include <vector>
#include <set>
#include <queue>
#include <stack>
#include <math.h>
#include <cassert>
#include <bitset>
#include <map>
#include <algorithm>
#include <iterator>
#include <string>
#include <utility>
#include <numeric>
#include <regex>
#include <complex>
#include <random>
#include <atcoder/all>
using namespace atcoder;
using namespace std;
using ll=long long;
using ull=unsigned long long;
using P=pair<ll,ll>;
template<typename T>using minque=priority_queue<T,vector<T>,greater<T>>;
template<typename T>bool chmax(T &a,const T &b){return (a<b?(a=b,true):false);}
template<typename T>bool chmin(T &a,const T &b){return (a>b?(a=b,true):false);}
template<typename T1,typename T2>ostream &operator<<(ostream &os,const pair<T1,T2>&p){os<<p.first<<' '<<p.second;return os;}
template<typename T1,typename T2>istream &operator>>(istream &is,pair<T1,T2>&p){is>>p.first>>p.second;return is;}
template<typename T> ostream &operator<<(ostream &os,const vector<T> &a){cout<<'{';for(int i=0;i<a.size();i++){cout<<a[i];if(i+1!=a.size())cout<<',';};cout<<'}';return os;}
template<typename T>istream &operator>>(istream &is,vector<T> &a){for(auto &i:a)is>>i;return is;}
#define reps(i,a,n) for(int i=(a);i<(n);i++)
#define rep(i,n) reps(i,0,n)
#define all(x) x.begin(),x.end()
#define pcnt(x) __builtin_popcount(x)
ll myceil(ll a,ll b){return (a+b-1)/b;}
template<typename T1,typename ...T2>
void debug(const T1& x,const T2& ...y){
  #ifdef LOCAL
  cout<<x;
  ((cout<<' '<<y),...);
  cout<<'\n';
  #endif
}
void SOLVE();
int main(){
  #ifdef LOCAL
  clock_t start=clock();
  #endif
  int testcase=1;
  //cin>>testcase;
  for(int i=0;i<testcase;i++){
    SOLVE();
  }
  #ifdef LOCAL
  cout<<"time:";
  cout<<(clock()-start)/1000;
  cout<<"ms\n";
  #endif
}
using R=double;
using Point=complex<R>;
constexpr R eps=1e-12;
constexpr R pi=3.141592653589793238;
inline bool equal(const R &a,const R &b){
  return fabs(a-b)<eps;
}
Point unit_vector(const Point &a){return a/abs(a);}
R dot(const Point &a,const Point &b){//内積
  return a.real()*b.real()+a.imag()*b.imag();
}
R cross(const Point &a,const Point &b){//外積
  return a.real()*b.imag()-a.imag()*b.real();
}
Point rotate(const Point &a,const R t){
  return Point(cos(t)*a.real()-sin(t)*a.imag(),sin(t)*a.real()+cos(t)*a.imag());
}
struct Line{
  Point a,b;
  Line():a(Point(0,0)),b(Point(0,0)){};
  Line(Point a,Point b):a(a),b(b){}
  Line(R A,R B,R C){
    assert(!(equal(A,0)&&equal(B,0)));
    if(equal(A,0)){
      a=Point(0,C/B),b=Point(1,C/B);
    }
    else if(equal(B,0)){
      a=Point(C/A,0),b=Point(C/A,1);
    }
    else{
      a=Point(C/A,0),b=Point(0,C/B);
    }
  }
};
Point projection(const Line &line,const Point &p){
  R t=dot(p-line.a,line.a-line.b)/norm(line.a-line.b);
  return line.a+(line.a-line.b)*t;
}
Point reflection(const Line &line,const Point &p){
  return p+(projection(line,p)-p)*(R)2;
}
bool ON(const Point &a,const Point &b,const Point &p){
  if(a==p)return true;
  if(b==p)return true;
  if(!equal(cross(a-p,b-p),0))return false;
  return dot(a-p,b-p)<0;
}
Point crosspoint(const Line &a,const Line &b){
  R c=cross(a.b-a.a,b.b-b.a);
  R d=cross(a.b-a.a,a.b-b.a);
  if(equal(abs(c),0)&&equal(abs(d),0))return b.a;
  return b.a+(b.b-b.a)*d/c;
}
Line suityokunitoubunsenn(const Point &a,const Point &b){
  Point c=(a+b)*(R)0.5;
  Point d=rotate(a-c,pi/2)+c;
  Line ret(c,d);
  return ret;
}
struct Circle{
  Point p;
  R r;
  Circle():p(Point(0,0)),r(0){};
  Circle(Point p,R r):p(p),r(r){}
  Circle(Point a,Point b,Point c){
    assert(!equal(cross(b-a,c-a),0));
    Line x=suityokunitoubunsenn(a,b),y=suityokunitoubunsenn(a,c);
    p=crosspoint(x,y);
    r=abs(a-p);
  }
  bool in(const Point &a){
    return abs(a-p)<=r+eps;
  }
};
struct E{
  int u,v;
  R dst;
};
bool comp(E a,E b){return a.dst<b.dst;}
void SOLVE(){
  int n;
  cin>>n;
  vector<Point>a(n);
  rep(i,n){
    R x,y;
    cin>>x>>y;
    a[i]={x,y};
  }
  dsu uf(n);
  vector<E>edge;
  vector<vector<int>>to(n);
  vector<int>cnt(n,0);
  vector<bool>seen(n,false);
  rep(i,n-1)reps(j,i+1,n){
    R cu=abs(a[i]-a[j]);
    edge.push_back({i,j,cu});
  }
  sort(all(edge),comp);
  for(auto [u,v,dst]:edge){
    if(uf.same(u,v))continue;
    uf.merge(u,v);
    to[u].push_back(v);
    to[v].push_back(u);
    cnt[u]++,cnt[v]++;
    debug(u,v);
  }
  queue<int>que;
  rep(i,n)if(cnt[i]==1)que.push(i),debug(i);
  vector<pair<int,int>>ans;
  while(!que.empty()){
    int x=que.front();
    que.pop();
    for(auto i:to[x]){
      if(!seen[i]){
        if(--cnt[i]<=1){
          if(!seen[x]){
            ans.push_back({i+1,x+1});
            seen[x]=true;
            seen[i]=true;
          }
        }
      }
    }
  }
  cout<<ans.size()<<endl;
  for(auto i:ans)cout<<i<<endl;
}
0