結果

問題 No.2376 障害物競プロ
ユーザー umezoumezo
提出日時 2023-07-10 05:23:33
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 394 ms / 4,000 ms
コード長 1,182 bytes
コンパイル時間 3,218 ms
コンパイル使用メモリ 200,108 KB
実行使用メモリ 7,480 KB
最終ジャッジ日時 2023-10-09 22:31:48
合計ジャッジ時間 52,061 ms
ジャッジサーバーID
（参考情報） 		judge15 / judge11
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 1 ms
4,376 KB
testcase_01 AC 1 ms
4,372 KB
testcase_02 AC 1 ms
4,372 KB
testcase_03 AC 2 ms
4,372 KB
testcase_04 AC 202 ms
5,872 KB
testcase_05 AC 300 ms
6,884 KB
testcase_06 AC 117 ms
5,108 KB
testcase_07 AC 379 ms
7,456 KB
testcase_08 AC 366 ms
7,408 KB
testcase_09 AC 349 ms
7,340 KB
testcase_10 AC 362 ms
7,468 KB
testcase_11 AC 328 ms
7,468 KB
testcase_12 AC 307 ms
7,480 KB
testcase_13 AC 373 ms
7,424 KB
testcase_14 AC 384 ms
7,380 KB
testcase_15 AC 343 ms
7,184 KB
testcase_16 AC 371 ms
7,328 KB
testcase_17 AC 314 ms
7,448 KB
testcase_18 AC 308 ms
7,292 KB
testcase_19 AC 315 ms
7,476 KB
testcase_20 AC 335 ms
7,332 KB
testcase_21 AC 336 ms
7,344 KB
testcase_22 AC 293 ms
6,708 KB
testcase_23 AC 200 ms
5,668 KB
testcase_24 AC 210 ms
5,900 KB
testcase_25 AC 107 ms
4,692 KB
testcase_26 AC 244 ms
6,228 KB
testcase_27 AC 206 ms
5,892 KB
testcase_28 AC 132 ms
5,124 KB
testcase_29 AC 106 ms
4,792 KB
testcase_30 AC 130 ms
4,816 KB
testcase_31 AC 131 ms
5,160 KB
testcase_32 AC 18 ms
4,384 KB
testcase_33 AC 56 ms
4,400 KB
testcase_34 AC 67 ms
4,536 KB
testcase_35 AC 43 ms
4,384 KB
testcase_36 AC 234 ms
5,768 KB
testcase_37 AC 279 ms
6,700 KB
testcase_38 AC 102 ms
4,884 KB
testcase_39 AC 301 ms
6,524 KB
testcase_40 AC 105 ms
4,496 KB
testcase_41 AC 99 ms
4,696 KB
testcase_42 AC 394 ms
7,324 KB
testcase_43 AC 394 ms
7,344 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#define rep(i,n) for(int i=0;i<(int)(n);i++)
#include <bits/stdc++.h>
using namespace std;

const int MA=310,MB=2e5+7,INF=1e9+7;
double EPS=1e-9,dis[MA][MA],X[MA],Y[MA];
int n,m,A[MB],B[MB],C[MB],D[MB];

void floyd(){
  rep(k,2*n){
    rep(i,2*n){
      if(dis[i][k]>INF-EPS) continue;
      rep(j,2*n){
        if(dis[k][j]>INF-EPS) continue;
        dis[i][j]=min(dis[i][j],dis[i][k]+dis[k][j]);
      }
    }
  }
}

int main(){
  cin>>n>>m;
  rep(i,2*n){
    rep(j,2*n) if(i!=j) dis[i][j]=INF;
    cin>>X[i]>>Y[i];
  }
  rep(i,m) cin>>A[i]>>B[i]>>C[i]>>D[i];
  
  auto cr=[&](int a,int b,int c)->double{
    return (X[b]-X[a])*(Y[c]-Y[a])-(Y[b]-Y[a])*(X[c]-X[a]);
  };
  auto f=[&](int a,int b,int c,int d)->bool{
    if(cr(a,b,c)*cr(a,b,d)<0 && cr(c,d,a)*cr(c,d,b)<0) return false;
    return 1;
  }; 
  
  rep(i,2*n){
    for(int j=i+1;j<2*n;j++){
      bool e=1;
      rep(k,n){
        if(k==i/2 || k==j/2) continue;
        if(f(i,j,2*k,2*k+1)==0) e=0;
      }
      if(e){
        double p=X[i]-X[j],q=Y[i]-Y[j];
        dis[j][i]=dis[i][j]=sqrt(p*p+q*q);
      }
    }
  }
  floyd();
  rep(i,m) printf("%.10f\n",dis[2*(A[i]-1)+B[i]-1][2*(C[i]-1)+D[i]-1]);
  
  return 0;
}
0