loj2134/uoj132/洛谷P2304/bzoj4200 小园丁与老司机 DP+有上下界网络流

最新推荐文章于 2021-01-21 23:40:46 发布
最新推荐文章于 2021-01-21 23:40:46 发布
阅读量1.5k 收藏 2
点赞数 3
分类专栏： 网络流 动态规划 文章标签： DP 有上下界网络流 最小流
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/litble/article/details/80463466
版权

题目分析

老司机

关于老司机的这一部分，显而易见是个dp。将纵坐标相等的点组成的集合称为层，则我们要逐层dp。

将(0,0)这个点看作一棵树，算出答案后再减去它的贡献。

首先预处理出每棵树向左上，上，右上会到哪棵树，这个可以离散+桶完成，正上的处理就是将x坐标放桶里，左上右上的就把经过该点的直线 y=x+b1 y = x + b 1 和 y=x−b2 y = x − b 2 对应的 b1 b 1 ， b2 b 2 搞进去。

然后就逐层dp撒，如图，会发现假如我们从i点到达该层，从j点离开该层。如果 i<j i < j ，则该层在 j j 左边的树都会被走到。如果i>j，则该层在 j j 右边的树都会被走到。当然啦，如果i=j，就只有这棵树会被走到。

我才不会告诉你我一开始没想清楚必须是没许愿的树才能在此转向而GG了一个小时呢

灵魂画手litble

于是我们就可以设 f(i) f ( i ) 表示从i树进入i所在层，可以走的树的最大值。从y最大的层往y最小的处理，对每一层按x排序后，先从左到右扫一遍，再从右往左扫一遍，同时记录一种方案（记录以i进入该层，要以哪个点出去最优和以j出去，是走左上右上还是正上）。

利用我们记录的信息递归输出一种方案，即可获得40分。

小园丁

由于有每层树的个数不超过1000或只有一条路径的限制，所以我们可以从y最小的层往y最大的处理，对于每一个“可以作为这层起点”的树i，枚举这层所有树j，看是否可以从j出去。如果可以，那么j出去达到的树也要标记为“可以作为这层起点”，如此如此，可以弄出所有需要压路机压的路径，但是千万别弄出重边。

这是一个经典的有上下界网络流的模型，即每条边的流量上下界为 [1,inf] [ 1 , i n f ] ，然后建立源点 S S 和汇点T， S S 往每个点连一条[0,inf]的边，每个点往 T T 连一条[0,inf]的边。

什么？你不会有上下界网络流？还不去学!就是这篇文章讲的“有源汇最小流”的模型。

可是如果你打完，就TLE了。

以下没看懂文字请直接看代码。

这张图由于比较特殊，所以可以考虑，我们首先让每条边流量等于其下界，那么有些点进出不平衡了嘛，由于 S S 和T与所有点都有无限流量的边，所以在此时我们可以直接让“积蓄”了流量的点将流量排到 T T ，“缺少”流量的点从S拿一点来补充，那么此时的流为所有“积蓄”的流量之和。

然后我希望最小流，也就是想让这张图自己“消化调整”，那么求一次超级源汇之间的最大流 kans k a n s ，这就是这张图自己不经过S和T做调节的最大限度，“积蓄”流量之和减去 kans k a n s 就是答案。

void work2() {
    build();
    S=n+1,T=n+2,ss=n+3,tt=n+4;//呃，其实这个S和T可以删掉
    for(RI i=1;i<=n;++i) {
        if(du[i]>0) adde(ss,i,du[i]),ans+=du[i];//du:如果有[1,inf]边(x,y)，那么++du[y],--du[x]
        else adde(i,tt,-du[i]);
    }
    while(bfs(ss,tt)) ans-=dfs(ss,inf,tt);
    printf("%d\n",ans);
}

代码

哈哈哈哈我没疯放开我我怎么可能敲码农题敲了一下午就疯了呢哈哈哈哈

如果是考场上我肯定选择弃疗。

#include<bits/stdc++.h>
using namespace std;
int read() {
    int q=0,w=1;char ch=' ';
    while(ch!='-'&&(ch<'0'||ch>'9')) ch=getchar();
    if(ch=='-') w=-1,ch=getchar();
    while(ch>='0'&&ch<='9') q=q*10+ch-'0',ch=getchar();
    return q*w;
}
#define RI register int
const int N=50010,inf=0x3f3f3f;
int n,jsx,jsy,js1,js2;
struct node{int x,y,b1,b2,id;}p[N];
int bx[N],by[N],bb1[N],bb2[N];

int up[N],lup[N],rup[N],f[N],Tx[N],Tb1[N],Tb2[N],bj[N],tx[N],ty[N],gup[N],bup[N];
//gup:往上走最多走几棵树，bup:往上走走最多的树，下一步走到哪颗树,bj:记录方案,T开头的都是桶
bool cmpy(node a,node b) {return a.y>b.y;}
bool cmpx(node a,node b) {return a.x<b.x;}
vector<int> iny[N];
void prework() {//预处理每个点左上，右上，正上是那些点
    sort(bx+1,bx+1+n),sort(by+1,by+1+n);
    sort(bb1+1,bb1+1+n),sort(bb2+1,bb2+1+n);
    jsx=1;for(RI i=2;i<=n;++i) if(bx[i]!=bx[jsx]) bx[++jsx]=bx[i];
    jsy=1;for(RI i=2;i<=n;++i) if(by[i]!=by[jsy]) by[++jsy]=by[i];
    js1=1;for(RI i=2;i<=n;++i) if(bb1[i]!=bb1[js1]) bb1[++js1]=bb1[i];
    js2=1;for(RI i=2;i<=n;++i) if(bb2[i]!=bb2[js2]) bb2[++js2]=bb2[i];
    for(RI i=1;i<=n;++i) {
        p[i].x=lower_bound(bx+1,bx+1+jsx,p[i].x)-bx,tx[i]=p[i].x;
        p[i].y=lower_bound(by+1,by+1+jsy,p[i].y)-by,ty[i]=p[i].y;
        p[i].b1=lower_bound(bb1+1,bb1+1+js1,p[i].b1)-bb1;
        p[i].b2=lower_bound(bb2+1,bb2+1+js2,p[i].b2)-bb2;
    }
    sort(p+1,p+1+n,cmpy);
    for(RI i=1;i<=n;++i) {
        int k=p[i].id;
        up[k]=Tx[p[i].x],lup[k]=Tb1[p[i].b1],rup[k]=Tb2[p[i].b2];
        Tx[p[i].x]=Tb1[p[i].b1]=Tb2[p[i].b2]=k;
    }
    sort(p+1,p+1+n,cmpx);
    for(RI i=1;i<=n;++i) iny[p[i].y].push_back(p[i].id);
}
void print(int num) {//输出方案
    int y=ty[num],sz=iny[y].size(),nxt=bj[num];
    if(num!=n) printf("%d ",num);
    if(tx[nxt]<tx[num]) {
        for(RI i=0;i<sz;++i)
            if(tx[iny[y][i]]>tx[num]) printf("%d ",iny[y][i]);
        for(RI i=sz-1;i>=0;--i)
            if(tx[iny[y][i]]<tx[num]&&tx[iny[y][i]]>=tx[nxt]) printf("%d ",iny[y][i]);
    }
    else if(tx[nxt]>tx[num]) {
        for(RI i=sz-1;i>=0;--i)
            if(tx[iny[y][i]]<tx[num]) printf("%d ",iny[y][i]);
        for(RI i=0;i<sz;++i)
            if(tx[iny[y][i]]>tx[num]&&tx[iny[y][i]]<=tx[nxt]) printf("%d ",iny[y][i]);
    }
    if(bup[nxt]) print(bup[nxt]);
}
void work1() {//DP主体
    prework();
    for(RI y=jsy;y>=1;--y) {
        int kmx=0,kbj=0,sz=iny[y].size();
        for(RI i=0;i<sz;++i) {
            int k=iny[y][i];
            if(up[k]&&f[up[k]]>gup[k]) gup[k]=f[up[k]],bup[k]=up[k];
            if(lup[k]&&f[lup[k]]>gup[k]) gup[k]=f[lup[k]],bup[k]=lup[k];
            if(rup[k]&&f[rup[k]]>gup[k]) gup[k]=f[rup[k]],bup[k]=rup[k];
            f[k]=kmx,bj[k]=kbj;
            if(sz-i+gup[k]>kmx) kmx=sz-i+gup[k],kbj=k;
        }
        kmx=kbj=0;
        for(RI i=sz-1;i>=0;--i) {
            int k=iny[y][i];
            if(kmx>f[k]) f[k]=kmx,bj[k]=kbj;
            if(gup[k]+1>f[k]) f[k]=gup[k]+1,bj[k]=k;
            if(i+1+gup[k]>kmx) kmx=i+1+gup[k],kbj=k;
        }
    }
    printf("%d\n",f[n]-1),print(n),puts("");
}

int tot=1,S,T,ss,tt,ans;
int ok[N],du[N],h[N],ne[N*10],to[N*10],flow[N*10],lev[N],q[N];
void adde(int x,int y,int z) {
    to[++tot]=y,ne[tot]=h[x],h[x]=tot,flow[tot]=z;
    to[++tot]=x,ne[tot]=h[y],h[y]=tot,flow[tot]=0;
}
void add(int x,int y) {
    for(RI i=h[x];i;i=ne[i]) if(to[i]==y) return;//注意不要重复（由于数据特殊性不会很慢）
    ok[y]=1,++du[y],--du[x],adde(x,y,inf);
}
void calc(int y,int i) {//暴力查看符合条件的点
    int sz=iny[y].size(),num=iny[y][i];
    for(RI j=0;j<i;++j) {
        int k=iny[y][j];
        if(up[k]&&f[up[k]]+sz-j==f[num]) add(k,up[k]);
        if(lup[k]&&f[lup[k]]+sz-j==f[num]) add(k,lup[k]);
        if(rup[k]&&f[rup[k]]+sz-j==f[num]) add(k,rup[k]);
    }
    for(RI j=i+1;j<sz;++j) {
        int k=iny[y][j];
        if(up[k]&&f[up[k]]+j+1==f[num]) add(k,up[k]);
        if(lup[k]&&f[lup[k]]+j+1==f[num]) add(k,lup[k]);
        if(rup[k]&&f[rup[k]]+j+1==f[num]) add(k,rup[k]);
    }
    if(up[num]&&f[up[num]]+1==f[num]) add(num,up[num]);
    if(lup[num]&&f[lup[num]]+1==f[num]) add(num,lup[num]);
    if(rup[num]&&f[rup[num]]+1==f[num]) add(num,rup[num]);
}
void build() {//网络流建图
    ok[n]=1;
    for(RI y=1;y<=jsy;++y) {
        int sz=iny[y].size();
        for(RI i=0;i<sz;++i) if(ok[iny[y][i]]) calc(y,i);
    }
}
int dfs(int x,int liu,int t) {//dinic的dfs
    if(x==t) return liu;
    int kl,sum=0;
    for(RI i=h[x];i;i=ne[i])
        if(flow[i]>0&&lev[to[i]]==lev[x]+1) {
            kl=dfs(to[i],min(flow[i],liu-sum),t);
            sum+=kl,flow[i]-=kl,flow[i^1]+=kl;
            if(sum==liu) return sum;
        }
    if(!sum) lev[x]=-1;
    return sum;
}
int bfs(int s,int t) {//dinic的bfs
    for(RI i=1;i<=n+4;++i) lev[i]=0;
    int he=1,ta=1; lev[s]=1,q[1]=s;
    while(he<=ta) {
        int x=q[he];++he;
        if(x==t) return 1;
        for(RI i=h[x];i;i=ne[i])
            if(flow[i]>0&&!lev[to[i]])
                lev[to[i]]=lev[x]+1,q[++ta]=to[i];
    }
    return 0;
}
void work2() {
    build();
    S=n+1,T=n+2,ss=n+3,tt=n+4;
    for(RI i=1;i<=n;++i) {
        if(du[i]>0) adde(ss,i,du[i]),ans+=du[i];
        else adde(i,tt,-du[i]);
    }
    while(bfs(ss,tt)) ans-=dfs(ss,inf,tt);
    printf("%d\n",ans);
}

int main()
{
    n=read();
    for(RI i=1;i<=n;++i) {
        bx[i]=p[i].x=read(),by[i]=p[i].y=read();
        bb1[i]=p[i].b1=p[i].y-p[i].x,bb2[i]=p[i].b2=p[i].x+p[i].y;
        p[i].id=i;
    }
    ++n,p[n].id=n;
    work1(),work2();
    return 0;
}

题外话

今天litble的血泪史：
血泪史1
血泪史2

litble
关注
专栏目录
acwing-2189. 有源汇上下界最大流
m0_46656833的博客
05-12 172
给定一个包含 n 个点 m 条边的有向图，每条边都有一个流量下界和流量上界。 给定源点 S 和汇点 T，求源点到汇点的最大流。 输入格式 第一行包含四个整数 n,m,S,T。 接下来 m 行，每行包含四个整数 a,b,c,d 表示点 a 和 b 之间存在一条有向边，该边的流量下界为 c，流量上界为 d。 点编号从 1 到 n。 输出格式 输出一个整数表示最大流。 如果无解，则输出 No Solution。 数据范围 1≤n≤202, 1≤m≤9999, 1≤a,b≤n, 0≤c≤d≤105 输入样例： 10
题解 [SDOI2016] 征途（LOJ #2035 / 洛谷 P4072）【斜率优化DP】
破壁人五号的博客
06-01 229
题目链接：洛谷 P4072 / LOJ #2035 题意 将 nnn 个数的序列分为 mmm 段，求各段内数字之和的方差的最小值 ×m2\times m^2×m2。n,m≤3000n,m\leq 3000n,m≤3000。 题解 方差不好看，先推一下式子： m2s2=m2∑i=1m(xi−x‾)2m=m(∑i=1m(xi2−2xix‾+x‾2))=m((∑i=1mxi2)−2x‾n+m⋅n2m2)=m(∑i=1mxi2)−n2\begin{aligned}m^2s^2=&m^2\dfrac{\sum
[Noi2015]小园丁和老司机
weixin_30686845的博客
05-04 88
来自FallDream的博客，未经允许，请勿转载，谢谢。 小园丁 Mr. S 负责看管一片田野，田野可以看作一个二维平面。田野上有n棵许愿树，编号1,2,3,…,n，每棵树可以看作平面上的一个点，其中第 i棵树 （1≤i≤n） 位于坐标(xi,yi)。任意两棵树的坐标均不相同。老司机 Mr. P从原点(0,0)驾车出发，进行若干轮行动。每一轮，Mr.P首先选择任意一个满足以下条件的方向： ...
[NOI2015]小园丁与老司机(DP+上下界最小流)
weixin_30552811的博客
08-07 566
由于每行点的个数不超过1000，所以行内DP可以使用$O(n^2)$算法。 先找到每个点所能直接到达的所有点（x,y,x+y或x-y相同），用排序实现。 第一问：以行为阶段，对于每行，暴力枚举最有路径在这行上的起点和终点，g[x]记录当这行的最优路径以x为起点时，终点应在什么位置，f[x]记录走到x且这一行以x为起点，之后最多还能走到多少个点。 第二问：由于当一行的起点和终点都确定后，决策也...
【BZOJ4200】【NOI2015】小园丁与老司机（动态规划，网络流）
小蒟蒻yyb的博客
07-10 351
题面 BZOJ权限题，洛谷链接 题解 一道二合一的题目 考虑第一问。 先考虑如何计算六个方向上的第一个点。 左右上很好考虑，只需要按照xxx或者yyy轴排序就行了。 对于454545度的斜角，两点一定在同一条直线上。 这条直线是x+y=bx+y=bx+y=b或x−y=bx−y=bx-y=b 所以按照x+yx+yx+y和x−yx−yx-y的值分类考虑，再按照顺序在xxx轴扫一遍...
P2304 [NOI2015] 小园丁与老司机（网络流/上下界网络流）
xajd1906的博客
01-21 157
P2304 [NOI2015] 小园丁与老司机 平面上有n个点，每次可以向左、右、上、左上45度、右上45度移动，然后直线移动到达第一个没有到过的点，如果没有这样的点就不能移动，求解一条最长路，然后求解将所有可能不是左右移动的道路都覆盖需要多少路径。 那么首先最长路可以利用dp求解，然后需要覆盖可以先将这样的道路找到，然后连接一条边，利用上下界网络流求解最小流即可。 ...
[BZOJ5212][Zjoi2018]历史（树形dp+LCT）
女装的你，如此好看！
09-15 244
Address https://www.luogu.org/problemnew/show/P4338 https://www.lydsy.com/JudgeOnline/problem.php?id=5212 http://uoj.ac/problem/374 https://loj.ac/problem/2434 Solution 简版题意：给定树上每个点 Access 的次数...
bzoj5340/loj2552「CTSC2018」假面（期望与概率+背包dp）
Icefox的博客
05-08 694
这签到题送我见祖宗了啊qaq 全世界都A了这题，像我这样只拿了10分的傻子怕是不多了吧。我们考虑怎么求最后的期望生命值。 因为减到0血就不减了，所以我们不能直接每次期望算，只好考虑一个背包转移。 dp[i][j]表示i减了j滴血的概率。减了K[i]滴血就代表i死了。 每次攻击时对x做一个O(mi)的背包转移即可。 最后对每个人枚举减j滴血的情况计算即可。我们考虑“结界”询问。 即就是问活
题解 [十二省联考 2019] 春节十二响（LOJ #3052 / 洛谷 P5290）【启发式合并】
破壁人五号的博客
06-08 249
题目链接：洛谷 P5290 / LOJ #3052 题意 给定一棵有根树，节点有权值 MiM_iMi​，将节点分为多个集合（记为 SiS_iSi​），最小化 ∑imax⁡j∈siMj\sum\limits_{i}\max\limits_{j\in s_i}M_ji∑​j∈si​max​Mj​，要求 ∀i,j,k(i,j∈sk)\forall i,j,k(i,j\in s_k)∀i,j,k(i,j∈sk​)，i,ji,ji,j 无祖孙关系。 题解 假如有两个子树，要将它们合并起来，显然取它们的最大值放一起、次
有源汇有上下界的最小流（LOJ117）
最新发布
06-11
这是一个IT类问题。该问题可以使用网络流算法解决。首先用 Dinic 或 ISAP 算法求出最大流，...在新图上跑一遍费用流，如果总流量等于原图的最大流，则新图中的最小费用即为所求最小流，否则原图不存在满足上下界的流。
小园丁与老司机
某蒟蒻的博客
05-22 1183
小园丁与老司机题目描述 照例贴传送门http://www.lydsy.com/JudgeOnline/problem.php?id=4200 题解 前40分直接dp。 每个点预处理出下，左下，右下三个方向最近的点，直接转移。 对于同层的点，单独转移。由于看到这题同层点<=1000这个条件所以无脑写了n^2，后来发现有一个点没这个条件，于是在UOJ上被卡成了97分。然后就懒得改了。 对于第
[NOI2015]小园丁与老司机,洛谷P2304,恶心的建图+有源汇上下界最小流
Deep_Kevin的娱乐空间
09-17 183
正题 这题建图太恶心了,思路十分的简单,看到相同的显然可以把分层图建出来,一个流量代表一辆压路车,对于必须要走过的边就设下限流量为1,然后跑有源汇的上下界最小流就可以,但是怎么将那些边建出来呢? 每一层考虑dp,看看同层转移到这个点的答案,如果转移点在左边,那么它可以先走到最左边,然后再走到这个点,所以对于左边的任何转移点,转移的时候加上的都是左边的点数,对于右边的转移点同理. 把同层可能的转移点用vector存起来,把原来转移到该点值的下层转移点用vector存起...
[NOI2015]小园丁与老司机
yzyyylx的博客
08-29 471
题面 题意 给出n棵树，从原点开始，每次可以在左、右、上、左上 45∘、右上 45∘五个方向中选择一个，然后一直向这个方向走，直至走到一棵未经过的树，然后继续选择方向，直至5个方向都不存在未经过的树。 现在要经过尽可能多的树，输出最多可能经过的树的数量以及任意一条路径。 现在定义向上、左上 45∘、右上 45∘三个方向走形成的线段不优美，现在要覆盖所有最优路径中的不优美线段，每次从原点...
[UOJ#132][BZOJ4200][luogu_P2304][NOI2015]小园丁与老司机
weixin_30477293的博客
12-02 132
[UOJ#132][BZOJ4200][luogu_P2304][NOI2015]小园丁与老司机 试题描述 小园丁 Mr. S 负责看管一片田野，田野可以看作一个二维平面。田野上有 \(n\) 棵许愿树，编号 \(1,2,3, \cdots , n\)，每棵树可以看作平面上的一个点，其中第 \(i\) 棵树 \((1 \le i \le n)\) 位于坐标 \((x_i, y_i)\)。任意两棵树...
4200: [Noi2015]小园丁与老司机
CRZbulabula的博客
12-29 1991
4200: [Noi2015]小园丁与老司机 Time Limit: 20 Sec  Memory Limit: 512 MBSec  Special Judge Submit: 139  Solved: 80 [Submit][Status][Discuss] Description 小园丁 Mr. S 负责看管一片田野，田野可以看作一个二维平面。田野上有 nn 棵许愿树，编号 1
小园丁与老司机_疲倦的园丁
weixin_26742939的博客
08-19 307
小园丁与老司机An old man was walking along a road in the countryside, and came upon a stately mansion. He wandered in through the open gate and began admiring the garden: beautiful flowers in vibrant colours...
【BZOJ4200】【UOJ132】【NOI2015】小园丁与老司机
cz_xuyixuan的博客
05-28 545
【题目链接】BZOJUOJ【思路要点】将所有点按照纵坐标排序，分别处理同一纵坐标的点。显然，每个点在各个方向上的后继点若存在，是唯一的，先预处理。记\(f_i\)表示从节点\(i\)出发，能够经过的最多的点数。若不考虑左右的移动，\(f_i\)就是\(i\)在三个方向上后继结点的\(f\)值的最大值加1，记这个值为\(tmp_i\)。不妨令同一纵坐标的点横坐标递增，那么有\(f_i=max\{tm...
uoj132/BZOJ4200/洛谷P2304 [Noi2015]小园丁与老司机 【dp + 带上下界网络流】
weixin_30525825的博客
06-27 129
题目链接 uoj132 题解 真是一道大码题，，，肝了一个上午 老司机的部分是一个\(dp\)，观察点是按\(y\)分层的，而且按每层点的上限来看可以使用\(O(nd)\)的\(dp\)，其中\(d\)是每层的点数 我们设\(f[i]\)表示从\(i\)点进入该层，直到走完为止所经过的最多点的数量，我们把原点也看做一棵树，计算答案时减去即可 转移只需枚举出点\(j\)，假如\(i\)在\(j\)的...
BZOJ4200: [Noi2015]小园丁与老司机 最小流
リーインカーネイション
04-10 1385
题意：平面上有N（N<=50000）个点，从原点出发，每次可以走上、左、右、左上45度、右上45度五种方向，只有碰到一个未走过的点才能停止，求： 1.最多走多少点 2.输出一种方案 3.对于所有最优方案，只保留上、左上、右上的边，求最少条数的路径将所有边都覆盖 同一y坐标的点不超过1000个 前两问DP一下即可，将所有点按y为一关键字x为二关键字排序，然后对于每一行的每个点，枚举是从这行之
评论 1
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值