結果
問題 | No.2376 障害物競プロ |
ユーザー |
![]() |
提出日時 | 2023-07-10 05:00:20 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 518 ms / 4,000 ms |
コード長 | 1,550 bytes |
コンパイル時間 | 1,801 ms |
コンパイル使用メモリ | 197,988 KB |
最終ジャッジ日時 | 2025-02-15 09:26:59 |
ジャッジサーバーID (参考情報) |
judge3 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 4 |
other | AC * 40 |
ソースコード
#define rep(i,n) for(int i=0;i<(int)(n);i++)#define ALL(v) v.begin(),v.end()typedef long long ll;#include <bits/stdc++.h>using namespace std;const double EPS=1e-9;const int MAX=310;const int INF=1e9+7;double dis[MAX][MAX];int n,m;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(){ios::sync_with_stdio(false);std::cin.tie(nullptr);cin>>n>>m;rep(i,2*n){rep(j,2*n){if(i==j) dis[i][j]=0;else dis[i][j]=INF;}}vector<double> X(2*n),Y(2*n);rep(i,2*n) cin>>X[i]>>Y[i];vector<int> A(m),B(m),C(m),D(m);rep(i,m){cin>>A[i]>>B[i]>>C[i]>>D[i];A[i]--,B[i]--,C[i]--,D[i]--;}auto cross=[&](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(cross(a,b,c)*cross(a,b,d)<0 && cross(c,d,a)*cross(c,d,b)<0) return false;return true;};auto g=[&](int i,int j)->double{return sqrt((X[i]-X[j])*(X[i]-X[j])+(Y[i]-Y[j])*(Y[i]-Y[j]));};rep(i,2*n){for(int j=i+1;j<2*n;j++){bool e=true;for(int k=0;k<n;k++){if(k==i/2 || k==j/2) continue;if(f(i,j,2*k,2*k+1)==false) e=false;}if(e){dis[i][j]=g(i,j);dis[j][i]=dis[i][j];}}}floyd();rep(i,m){printf("%.30f\n",dis[2*A[i]+B[i]][2*C[i]+D[i]]);}return 0;}