[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);
}