[agc019c]Fountain Walk

前言

这是一道比较简单的题。
但是我没想好答案增大的情况，所以我一开始觉得细节讨论有点麻烦。。

题目大意

不想说，看网站

做法

先通过一系列翻转交换等满足x1<=x2,y1<=y2。
我们注意经过温泉有减少答案（走1/4圆）
和增加答案（走1/2圆）两种。
我们不希望增加答案，希望尽量减少答案。
由于减少答案很少，矩形边长很大，因此我们不能为了走温泉绕弯路，只能向右向上走。
然后对这个矩形内的点按x坐标递增，y坐标做个最长上升子序列，当然就是可以最多减少的答案。
手玩容易知道，只有我们每一行每一列都走了温泉，才会有不得不增加答案的时候。

#include<cstdio>#include<algorithm>#include<cmath>#define y1 yyy1#define y2 yyy2#define fo(i,a,b) for(i=a;i<=b;i++)using namespace std;typedef double db;const int maxn=200000+10,mx=100000000;const db pi=acos(-1);struct dong{    int x,y;    friend bool operator<(dong a,dong b){        return a.x<b.x||a.x==b.x&&a.y<b.y;    }} a[maxn];int b[maxn],p[maxn];int lis[maxn];int i,j,k,l,t,n,m,x1,y1,x2,y2,z1,z2,top;db ans;int main(){    scanf("%d%d%d%d",&x1,&y1,&x2,&y2);    scanf("%d",&n);    fo(i,1,n) scanf("%d%d",&a[i].x,&a[i].y);    if (x1>x2){        swap(x1,x2);        swap(y1,y2);    }    if (y1>y2){        y1=mx-y1;        y2=mx-y2;        fo(i,1,n) a[i].y=mx-a[i].y;    }    a[++n].x=x1;    a[n].y=y1;    a[++n].x=x2;    a[n].y=y2;    fo(i,1,n) b[i]=a[i].y;    sort(b+1,b+n+1);    top=unique(b+1,b+n+1)-b-1;    fo(i,1,n) a[i].y=lower_bound(b+1,b+top+1,a[i].y)-b;    z1=y1;z2=y2;    y1=lower_bound(b+1,b+top+1,y1)-b;    y2=lower_bound(b+1,b+top+1,y2)-b;    p[0]=y1;    fo(i,1,n) p[i]=n+1;    sort(a+1,a+n+1);    fo(i,1,n){        if (a[i].x<x1||a[i].y<y1||a[i].y>y2) continue;        if (a[i].x>x2) break;        if (a[i].x==x1&&a[i].y==y1){            lis[i]=0;            continue;        }        lis[i]=upper_bound(p,p+n+1,a[i].y)-p;        p[lis[i]]=a[i].y;        if (a[i].x==x2&&a[i].y==y2) break;    }    k=lis[i]-1;    ans=(x2-x1)+(z2-z1);    ans*=100;    ans+=(pi-4)*5*k;    if (k==min(x2-x1+1,z2-z1+1)&&k) ans+=pi*5;    printf("%.12lf\n",ans);}