国产探花免费观看_亚洲丰满少妇自慰呻吟_97日韩有码在线_资源在线日韩欧美_一区二区精品毛片,辰东完美世界有声小说,欢乐颂第一季,yy玄幻小说排行榜完本

首頁 > 學院 > 開發設計 > 正文

【JZOJ3297】【SDOI2013】逃考(escape)

2019-11-06 06:05:21
字體:
來源:轉載
供稿:網友

Mission

高考又來了,對于不認真讀書的來講真不是個好消息。為了小楊能在家里認真讀書,他的親戚決定駐扎在他的家里監督他學習,有爺爺奶奶、外公外婆、大舅、大嫂、阿姨…… 小楊實在是忍無可忍了,這種生活跟監獄有什么區別!為了他親愛的小紅,為了他的dota,他決定越獄! 假設小楊的家是個n*m 的矩陣,左下角坐標為(0,0),右上角坐標為(x1,y1)。小楊有n 個親戚,駐扎在矩陣里(位置不同,且不在矩陣的邊上)。小楊家里的每個地方都被親戚監控著,而且只被距離最近的親戚監控: 也就是說假設小楊所在的位置是(3,3),親戚A 在(3,0),A 距離小楊距離是3;親戚B 在(6,7),則B 距離小楊距離是5。距離A < 距離B,所以(3,3)位置由A 監控。 如果“最近距離”出現同時有幾個親戚,那么那個位置同時被那幾個親戚監控。 給出小楊的坐標(x0,y0)。因為被發現的人數越少,越獄成功的機會越大,所以小楊需要你設計一條越獄路線到達矩形的邊上,且被發現的人數最少。 Ps:小楊做的方向是任意的,也就是說路線上的任意位置只需要是實數。 保證一開始小楊只被一個親戚監控著。 n<=600

Solution

顯然,每個點與其他點連線的中垂線形成的半平面求交,就是這個點的監視范圍。 然后最短路即可。

Code

#include<iostream>#include<algorithm>#include<stdio.h>#include<math.h>#include<string.h>#define ll long long#define db double using namespace std;const char* fin="ex3297.in";const char* fout="ex3297.out";const db eps=10e-10;const int inf=0x7fffffff;const int maxn=607,maxm=maxn*maxn;bool equ(db v,db x){return fabs(v-x)<=eps;}int sgn(db x){return equ(x,0.0)?0:(x<0?-1:1);}struct P{ db x,y; P(db _x=0,db _y=0){x=_x;y=_y;} P Operator +(P b){return P(x+b.x,y+b.y);} P operator -(P b){return P(x-b.x,y-b.y);} P operator *(db b){return P(x*b,y*b);} P per(){return P(y,-x);} db operator ^(P b){return x*b.y-y*b.x;} db arg(){return atan2(y,x);} void read(){scanf("%lf",&x);scanf("%lf",&y);}}a[maxn];struct L{ P p,v; int id; L(){} L(P _p,P _v,int _id=0){p=_p;v=_v;id=_id;} db arg(){return v.arg();} P operator &(L b){return b.p+b.v*((v^(p-b.p))/(v^b.v));}}b[maxn],c[maxn];int num;bool cmp(L a,L b){return a.arg()<b.arg();}bool in(P p,L l){return sgn((p-l.p)^l.v)>=0;}P mid(P a,P b){return P((a.x+b.x)/2.0,(a.y+b.y)/2.0);}int n,t,i,j,k,X,Y,sx,sy,st,tot,head,tail;int fi[maxn],ne[maxm],la[maxm],va[maxm];int B[maxn*10],dis[maxn];bool bz[maxn];void add_line(int a,int b,int c){ tot++; ne[tot]=fi[a]; la[tot]=b; va[tot]=c; fi[a]=tot;}void add(int v,int u){ if (dis[v]>u){ dis[v]=u; if (!bz[v]){ B[++tail]=v; bz[v]=true; } }}void spfa(int v){ int i,j,k; memset(dis,127,sizeof(dis)); head=tail=0; add(v,0); while (head++<tail){ for (k=fi[B[head]];k;k=ne[k]) add(la[k],dis[B[head]]+va[k]); bz[B[head]]=false; }}int main(){ scanf("%d",&t); while (t--){ scanf("%d",&n); scanf("%d%d%d%d",&X,&Y,&sx,&sy); tot=0; st=0; P stp(sx,sy); memset(fi,0,sizeof(fi)); for (i=1;i<=n;i++) a[i].read(); for (i=1;i<=n;i++){ num=0; for (j=1;j<=n;j++){ if (i==j) continue; b[++num]=L(mid(a[i],a[j]),(a[j]-a[i]).per(),j); } b[++num]=L(P(0,0),P(0,1),0); b[++num]=L(P(0,Y),P(1,0),0); b[++num]=L(P(X,Y),P(0,-1),0); b[++num]=L(P(X,0),P(-1,0),0); sort(b+1,b+num+1,cmp); int N=0; for (j=1;j<=num;j++){ if (!N) b[++N]=b[j]; else if (!equ(b[N].arg(),b[j].arg())) b[++N]=b[j]; else if (in(b[j].p,b[N])) b[N]=b[j]; } num=N; int head=1,tail=2; bool noans=false; c[1]=b[1],c[2]=b[2]; for (j=3;j<=num;j++){ while (head<tail && !in(c[tail]&c[tail-1],b[j])) tail--; while (head<tail && !in(c[head]&c[head+1],b[j])) head++; if (head==tail && sgn(c[head].v^b[j].v)<=0){noans=true;break;} c[++tail]=b[j]; } if (noans) continue; while (head<tail && !in(c[tail]&c[tail-1],c[head])) tail--; for (j=head;j<=tail;j++) add_line(i,c[j].id,1); if (!st){ bool ST=true; for (j=head;j<=tail;j++) if (!in(stp,c[j])){ ST=false; break; } if (ST) st=i; } } spfa(st); Warning

1.半平面交的流程 1)加入半平面以及四個特殊的邊界半平面; 2)按半平面的極角排序,極角相同的半平面取較內側的一個; 3)利用雙端隊列求出半平面交,頭尾都要檢查交點是否在新加入的半平面內; 注意:判斷無解(隊列中只剩一個半平面) 4)隊列尾交點是否在頭半平面中


發表評論 共有條評論
用戶名: 密碼:
驗證碼: 匿名發表
主站蜘蛛池模板: 靖西县| 潞城市| 灵川县| 确山县| 昭平县| 胶州市| 仁化县| 新津县| 赣榆县| 瓮安县| 江孜县| 金昌市| 合阳县| 兴业县| 长岭县| 盖州市| 西华县| 武强县| 穆棱市| 收藏| 广德县| 泗阳县| 福贡县| 晋城| 巍山| 五原县| 于都县| 双桥区| 鄱阳县| 卢湾区| 大余县| 多伦县| 咸丰县| 龙海市| 拉萨市| 来安县| 方正县| 青田县| 伊金霍洛旗| 正安县| 乐业县|