BZOJ:4219: 跑得比谁都快 3007: 拯救小云公主
Posted Swm_sxt
tags:
篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了BZOJ:4219: 跑得比谁都快 3007: 拯救小云公主相关的知识,希望对你有一定的参考价值。
4219: 跑得比谁都快
3007: 拯救小云公主
后一道是前一道的弱化版。
如果考虑一下二分答案再判定,判定的过程无非就是把每个点变成一个以答案为半价的圆,然后如果左上和右下被切开了肯定就不合法,那就变成1001里面那种对偶图,跑一下spfa就好。
再看其实就不需要二分答案了,把各个boss连接起来,二分的答案达到边权一半时这两个圆就相交。整条路径会在达到最大边权的时候完全连通,那么最大边权最小的路径就是答案。
所以就三角剖分+最小生成树,至于三角剖分我选择拉模版。
貌似得用long double
#include<cmath> #include<cstdio> #include<algorithm> #define eps 1e-12 #define MN 110001 #define ld long double using namespace std; int read_p,read_ca; inline int read(){ read_p=0;read_ca=getchar(); while(read_ca<‘0‘||read_ca>‘9‘) read_ca=getchar(); while(read_ca>=‘0‘&&read_ca<=‘9‘) read_p=read_p*10+read_ca-48,read_ca=getchar(); return read_p; } struct po{ld x,y;}; struct na{po p;int n;}p[MN]; struct e{int x,y,ne,pre;}b[MN<<5]; struct bi{int x,y;ld z;}B[MN<<5]; bool operator < (bi a,bi b){return a.z<b.z;} bool operator < (na a,na b){return a.p.x==b.p.x?a.p.y<b.p.y:a.p.x<b.p.x;} int n,m,num=1,l[MN],st[MN],top,NUM=0,x,y,_l[MN],_num=0,fa[MN];ld _A,_B; inline void _add(int x,int y){b[++num].x=x;b[num].y=y;b[num].ne=l[x];b[l[x]].pre=num;l[x]=num;} inline void add(int x,int y){_add(x,y);_add(y,x);} inline void _del(int i){if(l[b[i].x]==i)l[b[i].x]=b[i].ne;else{b[b[i].pre].ne=b[i].ne;if(b[i].ne)b[b[i].ne].pre=b[i].pre;}} inline void del(int i){_del(i);_del(i^1);} ld cha(po a,po b,po c){return (b.x-a.x)*(c.y-b.y)-(c.x-b.x)*(b.y-a.y);} bool _cha(po a,po b,po c){return cha(a,b,c)<=0;} inline bool ju(po a,po b,po c,po d){return (_cha(a,b,c)^_cha(a,b,d))&&(_cha(c,d,a)^_cha(c,d,b));} ld sqr(ld x){return x*x;} inline ld min(ld a,ld b){return a<b?a:b;} bool incir(po a,po b,po c,po d){ ld m=(a.x*(b.y-c.y)+b.x*(c.y-a.y)+c.x*(a.y-b.y))*2.; ld x=(sqr(a.x)*(b.y-c.y)+sqr(b.x)*(c.y-a.y)+sqr(c.x)*(a.y-b.y)-(a.y-b.y)*(b.y-c.y)*(c.y-a.y))/m; ld y=((a.x-b.x)*(b.x-c.x)*(c.x-a.x)-sqr(a.y)*(b.x-c.x)-sqr(b.y)*(c.x-a.x)-sqr(c.y)*(a.x-b.x))/m; return sqrt(sqr(x-a.x)+sqr(y-a.y))+eps>sqrt(sqr(x-d.x)+sqr(y-d.y)); } void work(int L,int R){ if (R-L<=2){ for (int i=L;i<R;i++) for (int j=i+1;j<=R;j++) add(i,j); return; } int mid=L+R>>1; work(L,mid);work(mid+1,R); top=0; for (int i=L;i<=R;st[++top]=i++) while (top>1&&cha(p[st[top-1]].p,p[st[top]].p,p[i].p)<0) top--; int Ld,Rd,id; for (int i=1;i<top;i++) if (st[i]<=mid&&st[i+1]>mid) Ld=st[i],Rd=st[i+1]; add(Ld,Rd); for(;;){ id=-1; for (int i=l[Ld];i;i=b[i].ne) if (b[i].y!=Rd) if (cha(p[b[i].y].p,p[Ld].p,p[Rd].p)>0&&(id==-1||incir(p[Ld].p,p[Rd].p,p[id].p,p[b[i].y].p))) id=b[i].y; for (int i=l[Rd];i;i=b[i].ne) if (b[i].y!=Ld) if (cha(p[Ld].p,p[Rd].p,p[b[i].y].p)>0&&(id==-1||incir(p[Ld].p,p[Rd].p,p[id].p,p[b[i].y].p))) id=b[i].y; if (id==-1) break; if (id<=mid){ for (int i=l[Ld];i;i=b[i].ne) if (b[i].y!=Rd) if (ju(p[id].p,p[Rd].p,p[b[i].y].p,p[Ld].p)) del(i); Ld=id; }else{ for (int i=l[Rd];i;i=b[i].ne) if (b[i].y!=Ld) if (ju(p[id].p,p[Ld].p,p[b[i].y].p,p[Rd].p)) del(i); Rd=id; } add(Ld,Rd); } } int gf(int x){return fa[x]==x?x:fa[x]=gf(fa[x]);} int main(){ //freopen("a.in","r",stdin); //freopen("a.out","w",stdout); register int i,j; n=read();scanf("%Lf%Lf",&_A,&_B); for (i=0;i<n;i++) scanf("%Lf%Lf",&p[i].p.x,&p[i].p.y),p[i].n=i+1,fa[i]=i; sort(p,p+n); work(0,n-1); for (i=0;i<n;i++) for (j=l[i];j;j=b[j].ne) if (b[j].y>i) B[NUM].x=i,B[NUM].y=b[j].y,B[NUM++].z=sqrt(sqr(p[i].p.x-p[b[j].y].p.x)+sqr(p[i].p.y-p[b[j].y].p.y))/2; for (i=0;i<n;i++) B[NUM].x=i,B[NUM].y=n,B[NUM++].z=min(p[i].p.x,_B-p[i].p.y); for (i=0;i<n;i++) B[NUM].x=i,B[NUM].y=n+1,B[NUM++].z=min(_A-p[i].p.x,p[i].p.y); n+=2;fa[n-1]=n-1;fa[n-2]=n-2; sort(B,B+NUM); for (i=0;i<NUM;i++){ int x=gf(B[i].x),y=gf(B[i].y); if (x==y) continue; fa[x]=y; if (gf(n-2)==gf(n-1)) return printf("%.6Lf\n",B[i].z),0; } }
以上是关于BZOJ:4219: 跑得比谁都快 3007: 拯救小云公主的主要内容,如果未能解决你的问题,请参考以下文章