[BZOJ4200][NOI2015]小园丁与老司机

bzoj
luogu
loj上可以交部分分程序

description

小园丁 Mr. S 负责看管一片田野,田野可以看作一个二维平面。田野上有 \(n\) 棵许愿树,编号 \(1,2,3, \ldots ,n\),每棵树可以看作平面上的一个点,其中第 \(i\) 棵树 (\(1 \leq i \leq n\)) 位于坐标 \((x_i,y_i)\)。任意两棵树的坐标均不相同。

老司机 Mr. P 从原点 \((0,0)\) 驾车出发,进行若干轮行动。每一轮,Mr. P 首先选择任意一个满足以下条件的方向:

  1. 为左、右、上、左上 \(45^{\circ}\) 、右上 \(45^{\circ}\) 五个方向之一。
  2. 沿此方向前进可以到达一棵他尚未许愿过的树。

完成选择后,Mr. P 沿该方向直线前进,必须到达该方向上距离最近的尚未许愿的树,在树下许愿并继续下一轮行动。如果没有满足条件的方向可供选择,则停止行动。他会采取最优策略,在尽可能多的树下许愿。若最优策略不唯一,可以选择任意一种。

不幸的是,小园丁 Mr. S 发现由于田野土质松软,老司机 Mr. P 的小汽车在每轮行进过程中,都会在田野上留下一条车辙印,一条车辙印可看作以两棵树(或原点和一棵树)为端点的一条线段。

在 Mr. P 之后,还有很多许愿者计划驾车来田野许愿,这些许愿者都会像 Mr. P 一样任选一种最优策略行动。Mr. S 认为非左右方向(即上、左上\(45^{\circ}\) 、右上 \(45^{\circ}\) 三个方向)的车辙印很不美观,为了维护田野的形象,他打算租用一些轧路机,在这群许愿者到来之前夯实所有「可能留下非左右方向车辙印」的地面。

「可能留下非左右方向车辙印」的地面应当是田野上的若干条线段,其中每条线段都包含在某一种最优策略的行进路线中。每台轧路机都采取满足以下三个条件的工作模式:

  1. 从原点或任意一棵树出发。
  2. 只能向上、左上 \(45^{\circ}\) 、右上 \(45^{\circ}\) 三个方向之一移动,并且只能在树下改变方向或停止。
  3. 只能经过「可能留下非左右方向车辙印」的地面,但是同一块地面可以被多台轧路机经过。

现在 Mr. P 和 Mr. S 分别向你提出了一个问题:

  1. 请给 Mr .P 指出任意一条最优路线。
  2. 请告诉 Mr. S 最少需要租用多少台轧路机。

sol

这道题又是个二合一啊。
对于第一问,为老司机求出一条最优路径,这个我们可以\(dp\)得到。
因为行走路线不能向下,我们令\(y_i\)相同的树为同一层,可以发现最优决策一定是:先在同一层中把所有可以走到的树全部走到,再向上走到达下一层。这么一来每一层就会有一个确定的入点和出点。(这里注意一下并不是所有层都会被经过)
先求出每个点向下、左下、右下三个方向到达的最近的树,这是这个点作为这一层的进入点的三个前驱。设\(f_i\)表示到达第\(i\)棵树作为当前这一层的出点时最多经过多少棵树,那么就可以枚举入点,分三种情况讨论:
1、入点就是出点,那么这一层就只能经过这一棵树。
2、入点在出点的左侧,那么可以经过出点左侧的所有树(在这一层中先一直往左走再一直往右走直到到达出点)
3、入点在出点的右侧,那么可以经过出点右侧的所有树。
然后这个\(dp\)是可以做到\(O(n)\)的。输出方案的时候,其实可以不记录\(dp\)转移的前驱,只要照着转移的方式再扫一遍所有前驱,看看是从哪个位置转移来的即可。

接下来就需要处理小园丁的问题。我们先可以求出那些边是必须经过的边,然后问题就转化为了要求这些边的流量都至少是\(1\)
这是一个上下界网络流的问题。我们先假设这些边的流量是\(1\),这样有的点就会不满足流量守恒。新建源点汇点对这些流量不平衡的点进行补流或者是分流。
因为求的是最小流,按照最小流的一般思路是先从新建源点向新建汇点跑最大流,在不建立从原汇点到原源点的\(inf\)的情况下,这样跑出的最大流就是原图中最多的可以承担的流量,所以答案就是\(\sum_{d_i>0}d_i\)-最大流。

code

码农题啊qwq

#include<cstdio>
#include<algorithm>
#include<cstring>
#include<vector>
#include<queue>
#include<map>
using namespace std;
int gi(){
    int x=0,w=1;char ch=getchar();
    while ((ch<'0'||ch>'9')&&ch!='-') ch=getchar();
    if (ch=='-') w=0,ch=getchar();
    while (ch>='0'&&ch<='9') x=(x<<3)+(x<<1)+ch-'0',ch=getchar();
    return w?x:-x;
}
#define cmax(a,b) (a<b?a=b:a)
const int N = 5e4+5;
const int inf = 1e9;
struct node{
    int x,y,id;
    bool operator < (const node &b) const{
        if (y==b.y) return x<b.x;return y<b.y;
    }
}p[N];
int n,o[N],len,d[N],ld[N],rd[N],f[N],g[N],tag[N],ans;
vector<int>v[N];map<int,int>M;queue<int>Q;
struct network{
    struct edge{int to,nxt,w;}a[N*100];
    int S,T,head[N],cnt,tot,du[N],dep[N],cur[N];
    void link(int u,int v){
        for (int e=head[u];e;e=a[e].nxt) if (a[e].to==v) return;
        a[++cnt]=(edge){v,head[u],inf};head[u]=cnt;
        a[++cnt]=(edge){u,head[v],0};head[v]=cnt;
        ++du[v];--du[u];tag[u]=1;
    }
    void link(int u,int v,int w){
        a[++cnt]=(edge){v,head[u],w};head[u]=cnt;
        a[++cnt]=(edge){u,head[v],0};head[v]=cnt;
    }
    bool bfs(){
        memset(dep,0,sizeof(dep));
        dep[S]=1;Q.push(S);
        while (!Q.empty()){
            int u=Q.front();Q.pop();
            for (int e=head[u];e;e=a[e].nxt)
                if (a[e].w&&!dep[a[e].to])
                    dep[a[e].to]=dep[u]+1,Q.push(a[e].to);
        }
        return dep[T];
    }
    int dfs(int u,int f){
        if (u==T) return f;
        for (int &e=cur[u];e;e=a[e].nxt)
            if (a[e].w&&dep[a[e].to]==dep[u]+1){
                int tmp=dfs(a[e].to,min(a[e].w,f));
                if (tmp) {a[e].w-=tmp;a[e^1].w+=tmp;return tmp;}
            }
        return 0;
    }
    int dinic(){
        int res=0;
        while (bfs()){
            for (int i=0;i<=T;++i) cur[i]=head[i];
            while (int tmp=dfs(S,inf)) res+=tmp;
        }
        return res;
    }
    int work(){
        S=n+1;T=n+2;
        for (int i=0;i<=n;++i){
            if (du[i]>0) link(S,i,du[i]),tot+=du[i];
            if (du[i]<0) link(i,T,-du[i]);
        }
        return tot-dinic();
    }
}G;
bool cmp(node i,node j){return i.id<j.id;}
int pre(int x){
    if (~d[x]&&g[x]==f[d[x]]) return d[x];
    if (~ld[x]&&g[x]==f[ld[x]]) return ld[x];
    if (~rd[x]&&g[x]==f[rd[x]]) return rd[x];
}
void dfs(int u){
    if (!u) return;int sz=v[p[u].y].size(),pos=0;
    while (v[p[u].y][pos]!=u) ++pos;
    if (f[u]==g[u]+1) {dfs(pre(u));printf("%d ",u);return;}
    for (int i=0;i<pos;++i){
        int x=v[p[u].y][i];
        if (f[u]==g[x]+pos+1){
            dfs(pre(x));
            for (int j=i;~j;--j) printf("%d ",v[p[u].y][j]);
            for (int j=i+1;j<=pos;++j) printf("%d ",v[p[u].y][j]);
            return;
        }
    }
    for (int i=pos+1;i<sz;++i){
        int x=v[p[u].y][i];
        if (f[u]==g[x]+sz-pos){
            dfs(pre(x));
            for (int j=i;j<sz;++j) printf("%d ",v[p[u].y][j]);
            for (int j=i-1;j>=pos;--j) printf("%d ",v[p[u].y][j]);
            return;
        }
    }
}
void build(int u){
    int sz=v[p[u].y].size(),pos=0;
    while (v[p[u].y][pos]!=u) ++pos;
    if (f[u]==g[u]+1){
        if (~d[u]&&g[u]==f[d[u]]) G.link(d[u],u);
        if (~ld[u]&&g[u]==f[ld[u]]) G.link(ld[u],u);
        if (~rd[u]&&g[u]==f[rd[u]]) G.link(rd[u],u);
    }
    for (int i=0;i<pos;++i){
        int x=v[p[u].y][i];
        if (f[u]==g[x]+pos+1){
            if (~d[x]&&g[x]==f[d[x]]) G.link(d[x],x);
            if (~ld[x]&&g[x]==f[ld[x]]) G.link(ld[x],x);
            if (~rd[x]&&g[x]==f[rd[x]]) G.link(rd[x],x);
        }
    }
    for (int i=pos+1;i<sz;++i){
        int x=v[p[u].y][i];
        if (f[u]==g[x]+sz-pos){
            if (~d[x]&&g[x]==f[d[x]]) G.link(d[x],x);
            if (~ld[x]&&g[x]==f[ld[x]]) G.link(ld[x],x);
            if (~rd[x]&&g[x]==f[rd[x]]) G.link(rd[x],x);
        }
    }
}
int main(){
    n=gi();G.cnt=1;
    for (int i=1;i<=n;++i) p[i]=(node){gi(),gi(),i},o[i]=p[i].y;
    sort(o+1,o+n+1);len=unique(o+1,o+n+1)-o-1;sort(p+1,p+n+1);
    memset(d,-1,sizeof(d));M.clear();M[0]=0;
    for (int i=1;i<=n;++i){
        if (M.find(p[i].x)!=M.end()) d[p[i].id]=M[p[i].x];
        M[p[i].x]=p[i].id;
    }
    memset(ld,-1,sizeof(ld));M.clear();M[0]=0;
    for (int i=1;i<=n;++i){
        if (M.find(p[i].x+p[i].y)!=M.end()) ld[p[i].id]=M[p[i].x+p[i].y];
        M[p[i].x+p[i].y]=p[i].id;
    }
    memset(rd,-1,sizeof(rd));M.clear();M[0]=0;
    for (int i=1;i<=n;++i){
        if (M.find(p[i].y-p[i].x)!=M.end()) rd[p[i].id]=M[p[i].y-p[i].x];
        M[p[i].y-p[i].x]=p[i].id;
    }
    for (int i=1;i<=n;++i){
        p[i].y=lower_bound(o+1,o+len+1,p[i].y)-o;
        v[p[i].y].push_back(p[i].id);
    }
    memset(f,-63,sizeof(f));memset(g,-63,sizeof(g));f[0]=0;
    for (int i=1;i<=len;++i){
        int sz=v[i].size(),mx;
        for (int j=0;j<sz;++j){
            int x=v[i][j];
            if (~d[x]) cmax(g[x],f[d[x]]);
            if (~ld[x]) cmax(g[x],f[ld[x]]);
            if (~rd[x]) cmax(g[x],f[rd[x]]);
            cmax(f[x],g[x]+1);
        }
        mx=-inf;
        for (int j=0;j<sz;++j){
            int x=v[i][j];
            cmax(f[x],mx+j+1);cmax(mx,g[x]);
        }
        mx=-inf;
        for (int j=sz-1;~j;--j){
            int x=v[i][j];
            cmax(f[x],mx+sz-j);cmax(mx,g[x]);
        }
    }
    for (int i=1;i<=n;++i) cmax(ans,f[i]);
    printf("%d\n",ans);sort(p+1,p+n+1,cmp);
    for (int i=1;i<=n;++i) if (f[i]==ans) {dfs(i);puts("");break;}
    for (int i=1;i<=n;++i) if (f[i]==ans) tag[i]=1;
    for (int i=len;i;--i)
        for (int j=0,sz=v[i].size();j<sz;++j)
            if (tag[v[i][j]]) build(v[i][j]);
    printf("%d\n",G.work());
    return 0;
}

转载于:https://www.cnblogs.com/zhoushuyu/p/9281144.html

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值