QOJ.ac
QOJ
ID | Problem | Submitter | Result | Time | Memory | Language | File size | Submit time | Judge time |
---|---|---|---|---|---|---|---|---|---|
#735020 | #5104. Guardians of the Gallery | BananaMac | AC ✓ | 96ms | 4288kb | C++20 | 5.5kb | 2024-11-11 16:42:12 | 2024-11-11 16:42:14 |
Judging History
answer
#include <bits/stdc++.h>
//#include <ranges>
using namespace std;
const double EPS= 1e-9;
struct Point {
long double x, y;
Point(long double x = 0.0, long double y = 0.0) : x(x), y(y) {}
Point operator+(const Point& p) const { return {x+p.x, y+p.y}; }
Point operator-(const Point& p) const { return {x-p.x, y-p.y}; }
Point operator*(long double c) const { return {x*c, y*c}; }
Point operator/(long double c) const { return {x/c, y/c}; }
long double Len() const { return hypot(x, y); }
};
inline long double CrossProd(const Point& a, const Point& b) {
return a.x*b.y - a.y*b.x;
}
inline long double DotProd(const Point& a, const Point& b) {
return a.x*b.x + a.y*b.y;
}
//如果两个向量的叉积结果符号相反(一个为正,一个为负),则说明线段AC和线段AD在线段AB的两侧,即线段AC和线段AD相交。
double RayIntersect(const Point& a, const Point& b, const Point& c, const Point& d, int* sides = NULL) {
double cp1 = CrossProd(c-a, b-a), cp2 = CrossProd(d-a, b-a);
double dp1 = DotProd(c-a, b-a), dp2 = DotProd(d-a, b-a);
if (sides) *sides = (cp1 < -EPS || cp2 < -EPS) + 2 * (cp1 > EPS || cp2 > EPS);//sides 0,1,2,3,如果两个线段相交,一定是3 1,2是ab直线穿过c或者d的情况
if (cp1 < -EPS && cp2 < -EPS || cp1 > EPS && cp2 > EPS) return -1.0;//叉积结果同号,线段不相交
return (abs(cp1) < EPS && abs(cp2) < EPS) ? 0 : (dp1*cp2-dp2*cp1)/(cp2-cp1);//a到交点的距离乘以ab长度
}
bool PointOnLine(const Point& a, const Point& b, const Point& p) {
double ln = (b-a).Len(), cp = CrossProd(b-a, p-a), dp = DotProd(b-a, p-a);
return abs(cp/ln) < EPS && dp/ln > -EPS && dp/ln < ln+EPS;//cp=0 && 0<dp<ln^2
}
int main(){
int n;
cin>>n;
vector<Point> polygon(n);
Point s,g;
for (int i = 0; i < n; ++i) {
cin>>polygon[i].x>>polygon[i].y;
}
cin>>g.x>>g.y;
cin>>s.x>>s.y;
vector<Point> endpoints;
//找指向每个顶点方向距离最近的intersecting边
polygon.push_back(g);//g点方向也需要找
for(auto p:polygon){
vector<tuple<double,int>> intersections;//交点距离,遮挡面
if((p-s).Len()<EPS)continue;
p=(p-s)/(p-s).Len()+s;
//找出sculpture向各个顶点出发, 不被遮挡能走的最远距离maxd,即视野不被遮挡的区域的终点(在多边形的边上)
for(int i=0;i<n;++i){
int side=0;
auto d= RayIntersect(s,p,polygon[i],polygon[(i+1)%n],&side);
if(d>-EPS)
intersections.emplace_back(d,side);
}
ranges::sort(intersections);
double maxd=0;
int sides=0;
for(auto [d,s]:intersections){
maxd=d;
sides|=s;
if(sides==3)break;
}
endpoints.push_back(s+(p-s)*maxd);//从sculpture出发,沿着射线,视野不被遮挡的区域的终点(在多边形的边上)
//计算顶点p2,到sculpture发出的穿过顶点i的射线,的距离,以及交点
for(auto p2:polygon){
auto ortho=Point(s.y-p.y,p.x-s.x)*1e5;
auto d= RayIntersect(s,p,p2,p2+ortho);
if(d<-EPS)
d= RayIntersect(s,p,p2,p2-ortho);
if(d>EPS && d<maxd-EPS)
endpoints.push_back(s+(p-s)*d);
}
}
polygon.pop_back();
//dijkstra算法求最短长度
vector<Point> points(polygon.size()+2+endpoints.size());
int i=0;
points[i++]=g;
for(auto p:polygon)points[i++]=p;
points[i++]=s;
for(auto p:endpoints)points[i++]=p;
priority_queue<pair<double,int>, vector<pair<double,int>>, greater<>> que;
vector<double> dist(points.size(),1e20);//从g到每个节点的最短路径长度
que.emplace(0.0,0);//从0号guard的起始位置出发
for(i=que.top().second;i<=n;i=que.top().second){
auto d=que.top().first;
assert(que.size());
que.pop();
if(d>=dist[i])continue;
dist[i]=d;
for(int j=0;j<points.size();++j){
if(i==j)continue;
auto p1=points[i], p2=points[j];
auto ln=(p1-p2).Len();
if(ln<EPS ||
//如果j在poly顶点i关联的两条poly边上,可以直接沿着poly边走到j,不需要后面的包含性判断
i>=1&&i<=n && PointOnLine(polygon[i-1], polygon[(i) % n], p2) ||
i>=1&&i<=n && PointOnLine(polygon[i-1], polygon[(i + n - 2) % n], p2)) {
que.emplace(d+ln,j);
continue;
}
//判断p1p2是否被poly完全包含
bool fullContained=true;
p2=p1+(p2-p1)/ln;
for(int k=0;k<polygon.size();++k){
auto rd= RayIntersect(p1,p2,polygon[k],polygon[(k+1)%n]);
if(rd>EPS && rd<ln-EPS) {
fullContained=false;
break;
}
}
if(!fullContained)continue;
//判断一下p1p2上任意一点是不是在poly内部,来保证p1p2全段都在poly内部。通过p1p2中点画一任意射线, 看穿过poly为奇数还是偶数次
int cnt=0;
p1 = p1*2/3 + p2/ 3;
p2 = p1 + Point(cos(10), sin(10));
for(int k=0;k<n;++k){
auto rd= RayIntersect(p1,p2,polygon[k],polygon[(k+1)%n]);
cnt+=(rd>EPS);
}
if(cnt%2==1)
que.emplace(d+ln,j);
}
}
printf("%.12f\n", que.top().first);
}
Details
Tip: Click on the bar to expand more detailed information
Test #1:
score: 100
Accepted
time: 1ms
memory: 3976kb
input:
15 13 7 20 20 39 20 49 7 73 13 100 5 117 38 98 20 80 20 66 40 68 20 51 20 41 39 22 48 2 39 10 20 104 20
output:
29.000000000000
result:
ok found '29.0000000', expected '29.0000000', error '0.0000000'
Test #2:
score: 0
Accepted
time: 1ms
memory: 4044kb
input:
16 39 2 48 22 39 41 20 51 20 68 40 66 20 80 20 98 38 117 5 100 13 73 7 49 19 39 20 23 20 20 7 13 20 10 20 104
output:
13.000000000000
result:
ok found '13.0000000', expected '13.0000000', error '0.0000000'
Test #3:
score: 0
Accepted
time: 1ms
memory: 4068kb
input:
16 13 33 20 60 23 66 39 97 49 105 73 166 100 205 117 272 98 216 80 180 66 172 68 156 51 122 41 121 22 92 2 44 10 40 104 228
output:
140.872282582487
result:
ok found '140.8722826', expected '140.8722826', error '0.0000000'
Test #4:
score: 0
Accepted
time: 1ms
memory: 3964kb
input:
16 64 17 50 28 67 23 65 18 77 4 88 20 78 10 70 29 61 28 47 32 54 17 43 13 35 20 41 30 27 20 42 6 81 12 33 23
output:
64.204537702523
result:
ok found '64.2045377', expected '64.2045377', error '0.0000000'
Test #5:
score: 0
Accepted
time: 1ms
memory: 4024kb
input:
16 64 17 50 28 67 23 65 18 77 4 88 20 78 10 70 29 61 28 47 32 54 17 43 13 35 20 41 30 27 20 42 6 33 23 81 12
output:
72.283498041177
result:
ok found '72.2834980', expected '72.2834980', error '0.0000000'
Test #6:
score: 0
Accepted
time: 0ms
memory: 3992kb
input:
7 76 8 389 215 691 19 407 331 489 397 300 403 363 334 126 60 393 370
output:
6.657917756511
result:
ok found '6.6579178', expected '6.6579178', error '0.0000000'
Test #7:
score: 0
Accepted
time: 0ms
memory: 4012kb
input:
3 0 1000 1000 0 1000 1000 567 578 589 601
output:
0.000000000000
result:
ok found '0.0000000', expected '0.0000000', error '-0.0000000'
Test #8:
score: 0
Accepted
time: 0ms
memory: 4040kb
input:
3 0 1000 0 0 1000 0 366 366 367 366
output:
0.000000000000
result:
ok found '0.0000000', expected '0.0000000', error '-0.0000000'
Test #9:
score: 0
Accepted
time: 0ms
memory: 4044kb
input:
5 50 93 278 178 199 300 596 362 208 519 421 388 142 153
output:
175.169759391735
result:
ok found '175.1697594', expected '175.1697594', error '0.0000000'
Test #10:
score: 0
Accepted
time: 0ms
memory: 4016kb
input:
7 50 93 278 178 199 300 401 312 483 162 596 362 208 519 488 252 142 153
output:
289.682139876890
result:
ok found '289.6821399', expected '289.6821399', error '0.0000000'
Test #11:
score: 0
Accepted
time: 0ms
memory: 4008kb
input:
8 10 10 40 25 20 25 20 35 12 23 30 23 10 20 5 40 15 15 19 26
output:
25.000000000000
result:
ok found '25.0000000', expected '25.0000000', error '0.0000000'
Test #12:
score: 0
Accepted
time: 0ms
memory: 3992kb
input:
9 5 1000 6 3 5 999 0 1000 0 0 500 2 500 0 1000 0 1000 1000 1 4 993 1
output:
5.101047906986
result:
ok found '5.1010479', expected '5.1010479', error '0.0000000'
Test #13:
score: 0
Accepted
time: 26ms
memory: 4060kb
input:
100 695 43 538 87 463 208 597 329 750 306 812 481 960 555 912 344 983 450 987 573 994 852 941 985 801 855 792 800 849 806 792 696 924 701 939 672 710 546 722 668 723 807 715 767 624 524 634 554 547 503 357 352 627 458 651 495 937 558 932 545 864 509 753 489 509 397 341 335 300 495 199 528 380 688 48...
output:
1695.186573023642
result:
ok found '1695.1865730', expected '1695.1865730', error '0.0000000'
Test #14:
score: 0
Accepted
time: 1ms
memory: 4036kb
input:
20 840 854 839 45 996 905 959 938 852 938 730 423 425 493 136 481 213 778 527 740 691 941 22 830 83 313 462 155 636 21 462 321 360 324 238 422 402 492 806 406 952 822 410 838
output:
1424.384201454849
result:
ok found '1424.3842015', expected '1424.3842015', error '0.0000000'
Test #15:
score: 0
Accepted
time: 24ms
memory: 3980kb
input:
74 89 395 52 622 124 832 115 698 95 598 199 491 190 356 191 398 132 315 94 371 34 221 91 0 153 139 220 465 260 283 312 30 409 15 338 50 343 52 437 69 359 89 332 213 377 505 375 396 405 199 657 90 658 50 360 50 618 23 642 7 824 191 688 417 795 227 709 286 662 321 646 175 485 210 381 357 420 329 441 3...
output:
2036.755709876637
result:
ok found '2036.7557099', expected '2036.7557099', error '0.0000000'
Test #16:
score: 0
Accepted
time: 14ms
memory: 4228kb
input:
100 380 626 511 639 548 551 651 476 706 462 636 604 652 617 776 577 794 566 821 433 765 410 778 276 735 345 700 329 448 550 283 482 537 332 706 213 741 204 833 152 657 182 626 173 568 225 602 213 673 203 537 286 459 317 609 261 493 344 334 430 468 338 331 400 350 326 512 197 553 155 424 120 446 179 ...
output:
307.850710857283
result:
ok found '307.8507109', expected '307.8507109', error '0.0000000'
Test #17:
score: 0
Accepted
time: 96ms
memory: 4108kb
input:
100 425 641 614 667 719 714 598 761 548 727 505 713 415 832 505 856 724 762 764 767 803 755 773 727 826 633 832 509 842 570 829 456 742 430 706 513 604 527 942 208 912 569 959 330 975 605 977 878 882 609 844 694 869 789 930 896 992 894 763 937 699 930 701 854 732 810 709 820 657 881 507 896 342 805 ...
output:
1941.568735726889
result:
ok found '1941.5687357', expected '1941.5687357', error '0.0000000'
Test #18:
score: 0
Accepted
time: 78ms
memory: 4184kb
input:
100 845 528 842 889 837 997 809 663 786 746 793 554 782 470 769 798 709 992 520 993 95 983 191 897 250 666 136 715 139 745 32 979 32 918 5 916 0 740 31 283 10 238 36 177 102 740 141 635 145 353 132 435 106 607 130 383 41 66 139 12 403 11 330 45 225 48 153 216 251 342 233 374 289 424 266 99 334 62 34...
output:
1863.571740263405
result:
ok found '1863.5717403', expected '1863.5717403', error '0.0000000'
Test #19:
score: 0
Accepted
time: 0ms
memory: 4044kb
input:
4 0 0 1000 0 1000 1000 0 1000 4 939 27 58
output:
0.000000000000
result:
ok found '0.0000000', expected '0.0000000', error '-0.0000000'
Test #20:
score: 0
Accepted
time: 4ms
memory: 4088kb
input:
94 5 5 995 5 995 995 5 995 990 990 5 990 970 970 5 970 950 950 5 950 930 930 5 930 910 910 5 910 890 890 5 890 870 870 5 870 850 850 5 850 830 830 5 830 810 810 5 810 790 790 5 790 770 770 5 770 750 750 5 750 730 730 5 730 710 710 5 710 690 690 5 690 670 670 5 670 650 650 5 650 630 630 5 630 610 610...
output:
620.291060712630
result:
ok found '620.2910607', expected '620.2910607', error '0.0000000'
Test #21:
score: 0
Accepted
time: 0ms
memory: 4040kb
input:
8 0 0 20 0 20 30 60 30 60 0 80 0 80 50 0 50 70 30 70 10
output:
0.000000000000
result:
ok found '0.0000000', expected '0.0000000', error '-0.0000000'
Test #22:
score: 0
Accepted
time: 0ms
memory: 4008kb
input:
5 2 0 10 0 0 10 0 3 5 3 1 8 5 2
output:
5.000000000000
result:
ok found '5.0000000', expected '5.0000000', error '0.0000000'
Test #23:
score: 0
Accepted
time: 0ms
memory: 4008kb
input:
5 2 0 10 0 0 10 0 3 5 3 1 4 5 2
output:
4.000000000000
result:
ok found '4.0000000', expected '4.0000000', error '0.0000000'
Test #24:
score: 0
Accepted
time: 0ms
memory: 4020kb
input:
12 0 0 2 0 2 4 1 4 1 5 2 5 2 7 0 7 0 3 1 3 1 2 0 2 1 6 1 1
output:
2.000000000000
result:
ok found '2.0000000', expected '2.0000000', error '0.0000000'
Test #25:
score: 0
Accepted
time: 0ms
memory: 3904kb
input:
10 0 0 2 0 2 4 1 4 2 5 2 7 0 7 0 3 1 2 0 2 1 6 1 1
output:
2.000000000000
result:
ok found '2.0000000', expected '2.0000000', error '0.0000000'
Test #26:
score: 0
Accepted
time: 0ms
memory: 4012kb
input:
12 0 0 2 0 3 3 5 0 6 3 7 0 8 2 7 6 5 2 4 6 2 2 0 2 1 1 7 1
output:
6.478708664619
result:
ok found '6.4787087', expected '6.4787087', error '0.0000000'
Test #27:
score: 0
Accepted
time: 0ms
memory: 4048kb
input:
8 10 0 11 3 12 0 14 0 15 3 16 0 18 4 8 4 10 1 16 1
output:
5.876122922140
result:
ok found '5.8761229', expected '5.8761229', error '0.0000000'
Test #28:
score: 0
Accepted
time: 0ms
memory: 4108kb
input:
8 10 0 11 3 12 0 16 3 14 0 17 0 17 4 8 4 10 1 15 1
output:
7.236067977500
result:
ok found '7.2360680', expected '7.2360680', error '0.0000000'
Test #29:
score: 0
Accepted
time: 0ms
memory: 3960kb
input:
8 0 0 20 0 20 30 60 30 60 0 80 0 80 50 0 50 70 10 10 10
output:
58.137767414995
result:
ok found '58.1377674', expected '58.1377674', error '0.0000000'
Test #30:
score: 0
Accepted
time: 0ms
memory: 4088kb
input:
11 0 0 4 0 4 1 5 1 5 0 7 0 7 2 3 2 3 1 2 2 0 2 6 1 1 1
output:
2.000000000000
result:
ok found '2.0000000', expected '2.0000000', error '0.0000000'
Test #31:
score: 0
Accepted
time: 0ms
memory: 3904kb
input:
3 31 41 59 26 53 58 36 41 56 31
output:
0.000000000000
result:
ok found '0.0000000', expected '0.0000000', error '-0.0000000'
Test #32:
score: 0
Accepted
time: 7ms
memory: 3996kb
input:
97 6 10 8 10 8 12 6 12 6 14 4 8 2 6 2 10 4 10 2 12 4 12 4 14 0 14 0 0 4 0 4 2 2 2 2 4 4 4 4 6 6 6 6 8 8 4 8 2 6 4 6 0 8 0 10 2 10 0 16 0 12 2 16 2 18 0 22 0 12 4 22 4 22 6 24 6 24 2 18 2 24 0 26 2 26 6 22 8 14 10 32 10 24 8 26 8 28 6 28 2 26 0 30 0 32 2 36 4 36 2 34 2 32 0 38 0 38 6 30 2 32 4 30 4 3...
output:
72.481035080187
result:
ok found '72.4810351', expected '72.4810351', error '0.0000000'
Test #33:
score: 0
Accepted
time: 28ms
memory: 4180kb
input:
97 180 20 190 20 170 30 160 30 160 40 170 40 180 30 190 30 190 40 180 40 170 50 160 50 160 60 170 60 180 50 190 50 180 60 190 60 190 70 150 60 100 50 110 50 110 30 100 30 100 40 90 40 90 30 80 40 80 30 70 20 70 30 50 30 60 40 70 40 50 60 80 60 70 50 90 50 90 60 140 60 180 70 60 70 30 60 50 70 20 70 ...
output:
240.478385698309
result:
ok found '240.4783857', expected '240.4783857', error '0.0000000'
Test #34:
score: 0
Accepted
time: 14ms
memory: 4004kb
input:
89 10 15 10 10 5 10 5 15 0 15 0 0 5 0 5 5 15 5 10 0 100 0 20 5 25 5 15 10 85 10 30 5 90 5 105 0 285 0 295 5 350 5 290 0 445 0 355 5 380 5 450 0 735 0 715 5 665 5 710 10 730 10 720 5 735 5 740 0 745 0 740 5 745 5 745 15 740 10 735 10 740 15 685 15 705 10 630 10 660 5 590 5 610 10 625 10 680 15 615 15...
output:
306.457898336281
result:
ok found '306.4578983', expected '306.4578983', error '0.0000000'
Test #35:
score: 0
Accepted
time: 10ms
memory: 4164kb
input:
94 5 5 995 5 995 995 5 995 990 990 5 990 970 970 5 970 950 950 5 950 930 930 5 930 910 910 5 910 890 890 5 890 870 870 5 870 850 850 5 850 830 830 5 830 810 810 5 810 790 790 5 790 770 770 5 770 750 750 5 750 730 730 5 730 710 710 5 710 690 690 5 690 670 670 5 670 650 650 5 650 630 630 5 630 610 610...
output:
1306.733331632740
result:
ok found '1306.7333316', expected '1306.7333316', error '0.0000000'
Test #36:
score: 0
Accepted
time: 12ms
memory: 3984kb
input:
100 1000 20 1 40 1000 60 1 80 1000 100 1 120 1000 140 1 160 1000 180 1 200 1000 220 1 240 1000 260 1 280 1000 300 1 320 1000 340 1 360 1000 380 1 400 1000 420 1 440 1000 460 1 480 1000 500 1 520 1000 540 1 560 1000 580 1 600 1000 620 1 640 1000 660 1 680 1000 700 1 720 1000 740 1 760 1000 780 1 800 ...
output:
47913.598737540771
result:
ok found '47913.5987375', expected '47913.5987375', error '0.0000000'
Test #37:
score: 0
Accepted
time: 77ms
memory: 4172kb
input:
100 0 0 1000 0 1000 1000 0 1000 0 3 998 2 998 998 2 998 2 5 996 4 996 996 4 996 4 7 994 6 994 994 6 994 6 9 992 8 992 992 8 992 8 11 990 10 990 990 10 990 10 13 988 12 988 988 12 988 12 15 986 14 986 986 14 986 14 17 984 16 984 984 16 984 16 19 982 18 982 982 18 982 18 21 980 20 980 980 20 980 20 23...
output:
46834.005636191228
result:
ok found '46834.0056362', expected '46834.0056362', error '0.0000000'
Test #38:
score: 0
Accepted
time: 8ms
memory: 4216kb
input:
100 7 409 21 321 39 254 62 198 74 177 79 169 96 146 105 135 118 120 154 90 181 72 188 68 199 62 210 57 219 53 239 45 311 24 320 22 356 15 376 12 408 8 483 2 499 1 500 1 516 2 570 6 581 7 630 13 746 40 749 41 768 48 780 53 823 75 834 82 845 90 854 97 903 146 913 159 925 177 964 267 972 295 973 299 98...
output:
0.000000000000
result:
ok found '0.0000000', expected '0.0000000', error '-0.0000000'
Test #39:
score: 0
Accepted
time: 1ms
memory: 4016kb
input:
13 100 100 110 110 100 120 110 130 100 140 110 150 100 160 90 155 100 145 90 135 100 125 90 115 100 105 100 106 100 159
output:
34.000000000000
result:
ok found '34.0000000', expected '34.0000000', error '0.0000000'
Test #40:
score: 0
Accepted
time: 1ms
memory: 3968kb
input:
13 100 100 110 90 120 100 130 90 140 100 150 90 160 100 155 110 145 100 135 110 125 100 115 110 105 100 106 100 159 100
output:
34.000000000000
result:
ok found '34.0000000', expected '34.0000000', error '0.0000000'
Test #41:
score: 0
Accepted
time: 1ms
memory: 3964kb
input:
13 100 100 120 100 120 120 140 120 140 140 160 140 160 160 145 165 145 145 125 145 125 125 105 125 105 105 106 106 159 159
output:
48.083261120685
result:
ok found '48.0832611', expected '48.0832611', error '0.0000000'
Test #42:
score: 0
Accepted
time: 1ms
memory: 4048kb
input:
13 100 100 100 120 80 120 80 140 60 140 60 160 40 160 35 145 55 145 55 125 75 125 75 105 95 105 94 106 41 159
output:
48.083261120685
result:
ok found '48.0832611', expected '48.0832611', error '0.0000000'
Test #43:
score: 0
Accepted
time: 16ms
memory: 4288kb
input:
100 7 409 21 321 39 254 62 198 74 177 79 169 96 146 105 135 118 120 154 90 181 72 188 68 199 62 210 57 219 53 239 45 311 24 320 22 356 15 376 12 408 8 483 2 499 1 500 1 516 2 570 6 581 7 630 13 746 40 749 41 768 48 780 53 823 75 834 82 845 90 854 97 903 146 913 159 925 177 964 267 972 295 981 335 99...
output:
467.001070662584
result:
ok found '467.0010707', expected '467.0010707', error '0.0000000'
Extra Test:
score: 0
Extra Test Passed