由于每行点的个数不超过1000,所以行内DP可以使用$O(n^2)$算法。

先找到每个点所能直接到达的所有点(x,y,x+y或x-y相同),用排序实现。

第一问:以行为阶段,对于每行,暴力枚举最有路径在这行上的起点和终点,g[x]记录当这行的最优路径以x为起点时,终点应在什么位置,f[x]记录走到x且这一行以x为起点,之后最多还能走到多少个点。

第二问:由于当一行的起点和终点都确定后,决策也是确定的,故只需要沿着DFS一遍即可得到最优路径。

第三问:先考虑怎么找到所有可能在最优路径上的边,同样暴力枚举每行的起点和终点即可,注意处理(0,0)的情况。找到符合要求的边后,问题就是用最少的路径数覆盖这些边,使用有上下界的最小流求解。

 #include<map>
#include<cstdio>
#include<algorithm>
#define rep(i,l,r) for (int i=(l); i<=(r); i++)
#define For(i,x) for (int i=h[x],k; i; i=nxt[i])
using namespace std; const int N=,M=,inf=;
bool va[N],vb[N];
int n,mx,cnt=,s0,ss,tt,S,T,d[N],f[N],h[N],dis[N],idx[N],q[N],cur[N],a[N],mp[N][];
int Lst[N],Nxt[N],pre[N],suf[N],premx[N],sufmx[N],to[M],fl[M],nxt[M],Do[N],Do2[N];
struct P{ int x,y; }p[N]; bool cmp1(int a,int b){ return (p[a].y==p[b].y) ? p[a].x<p[b].x : p[a].y<p[b].y; }
bool cmp2(int a,int b){ return (p[a].x==p[b].x) ? p[a].y<p[b].y : p[a].x<p[b].x; }
bool cmp3(int a,int b){ return (p[a].x-p[a].y==p[b].x-p[b].y) ? p[a].y<p[b].y : p[a].x-p[a].y<p[b].x-p[b].y; }
bool cmp4(int a,int b){ return (p[a].x+p[a].y==p[b].x+p[b].y) ? p[a].y<p[b].y : p[a].x+p[a].y<p[b].x+p[b].y; }
bool cmp5(int a,int b){ return (p[a].y==p[b].y) ? p[a].x<p[b].x : p[a].y>p[b].y; } inline void add(int u,int v,int w){
to[++cnt]=v; fl[cnt]=w; nxt[cnt]=h[u]; h[u]=cnt;
to[++cnt]=u; fl[cnt]=; nxt[cnt]=h[v]; h[v]=cnt;
} void work(int u,int v,int low,int high){ add(u,v,high-low); d[u]-=low; d[v]+=low; } bool bfs(){
rep(i,,n+) dis[i]=; dis[S]=; q[]=S;
for (int st=,ed=; st<ed; ){
int x=q[++st];
For(i,x) if (fl[i] && !dis[k=to[i]]) dis[k]=dis[x]+,q[++ed]=k;
}
return dis[T];
} int dfs(int x,int lim){
if (x==T) return lim;
int c=;
for (int i=cur[x],k; i; i=nxt[i])
if (fl[i] && dis[k=to[i]]==dis[x]+){
int t=dfs(k,min(lim-c,fl[i]));
c+=t; fl[i]-=t; fl[i^]+=t;
if (fl[i]) cur[x]=i;
if (c==lim) return lim;
}
if (!c) dis[x]=-;
return c;
} int dinic(){
int res=;
while (bfs()){
rep(i,,n+) cur[i]=h[i];
res+=dfs(S,inf);
}
return res;
} void Print(int x){
if (!x) return;
printf("%d ",x);
if (Do[x]==) { Print(mp[x][Do2[x]]); return; }
if (idx[x]>idx[Do[x]]){
int t=Do[x];
for (int i=Nxt[x]; i; i=Nxt[i]) printf("%d ",i);
for (x=Lst[x]; x!=t; x=Lst[x]) printf("%d ",x);
printf("%d ",t); Print(mp[t][Do2[t]]); return;
}
if (idx[x]<idx[Do[x]]){
int t=Do[x];
for (int i=Lst[x]; i; i=Lst[i]) printf("%d ",i);
for (x=Nxt[x]; x!=t; x=Nxt[x]) printf("%d ",x);
printf("%d ",t); Print(mp[t][Do2[t]]); return;
}
} void solve1(){
rep(i,,n) a[i]=i;
sort(a+,a+n+,cmp1);
rep(i,,n) if (p[a[i-]].y==p[a[i]].y) Lst[a[i]]=a[i-],Nxt[a[i-]]=a[i];
sort(a+,a+n+,cmp2);
rep(i,,n) if (p[a[i-]].x==p[a[i]].x) mp[a[i-]][]=a[i];
sort(a+,a+n+,cmp3);
rep(i,,n) if (p[a[i-]].x-p[a[i-]].y==p[a[i]].x-p[a[i]].y) mp[a[i-]][]=a[i];
sort(a+,a+n+,cmp4);
rep(i,,n) if (p[a[i-]].x+p[a[i-]].y==p[a[i]].x+p[a[i]].y) mp[a[i-]][]=a[i];
sort(a+,a+n+,cmp5);
rep(i,,n) idx[a[i]]=i;
} void W(int x){
rep(i,,) if (mp[x][i]){
int to=mp[x][i];
if (f[to]+>f[x]) Do2[x]=i,f[x]=f[to]+;
}
} int calc(int x,int y){
return (idx[x]<idx[y]) ? pre[y]+f[mp[y][Do2[y]]] : suf[y]+f[mp[y][Do2[y]]];
} bool jud(P &p){ return p.x== || p.x==p.y || p.x==-p.y; } void solve2(){
rep(k,,n) if (p[a[k]].y!=p[a[k-]].y) pre[a[k]]=; else pre[a[k]]=pre[a[k-]]+;
for (int k=n; k; k--) if (p[a[k]].y!=p[a[k+]].y) suf[a[k]]=; else suf[a[k]]=suf[a[k+]]+;
rep(k,,n) if (p[a[k]].y!=p[a[k-]].y){
for (int x=a[k]; x; x=Nxt[x]) f[x]=,W(x);
for (int x=a[k]; x; x=Nxt[x])
for (int y=a[k]; y; y=Nxt[y])
if (y!=x && calc(x,y)>f[x]) f[x]=calc(x,y),Do[x]=y;
}
rep(i,,n) if (jud(p[i])) mx=max(mx,f[i]);
printf("%d\n",mx);
rep(i,,n) if (f[i]==mx && jud(p[i])) { Print(i); break; }
puts("");
} int solve3(){
rep(i,,n) if (f[i]==mx && jud(p[i])) work(s0,i,,inf),va[i]=;
for (int i=n; i; i--) if (p[a[i]].y!=p[a[i+]].y){
for (int x=a[i]; x; x=Lst[x]) if (va[x]){
if (f[mp[x][Do2[x]]]+==f[x]) vb[x]=;
for (int y=a[i]; y; y=Lst[y]) if (y!=x && calc(x,y)==f[x]) vb[y]=;
}
for (int x=a[i]; x; x=Lst[x])
if (vb[x]) rep(j,,) if (mp[x][j] && f[mp[x][j]]==f[mp[x][Do2[x]]])
va[mp[x][j]]=,work(x,mp[x][j],,inf);
}
rep(i,,n+) work(ss,i,,inf),work(i,tt,,inf);
rep(i,,n+) if (d[i]>) add(S,i,d[i]); else add(i,T,-d[i]);
add(tt,ss,inf); dinic(); S=tt; T=ss; return inf-dinic();
} int main(){
freopen("driver.in","r",stdin);
freopen("driver.out","w",stdout);
scanf("%d",&n); s0=n+; ss=n+; tt=n+; S=n+; T=n+;
rep(i,,n) scanf("%d%d",&p[i].x,&p[i].y);
solve1(); solve2(); printf("%d\n",solve3());
return ;
}

[NOI2015]小园丁与老司机(DP+上下界最小流)的更多相关文章

  1. luogu P2304 [NOI2015]小园丁与老司机 dp 上下界网络流

    LINK:小园丁与老司机 苦心人 天不负 卧薪尝胆 三千越甲可吞吴 AC的刹那 真的是泪目啊 很久以前就写了 当时记得特别清楚 写到肚子疼.. 调到胳膊疼.. ex到根不不想看的程度. 当时wa了 一 ...

  2. 【BZOJ4200】[Noi2015]小园丁与老司机 DP+最小流

    [BZOJ2839][Noi2015]小园丁与老司机 Description 小园丁 Mr. S 负责看管一片田野,田野可以看作一个二维平面.田野上有 nn 棵许愿树,编号 1,2,3,…,n1,2, ...

  3. [BZOJ4200][Noi2015]小园丁与老司机

    4200: [Noi2015]小园丁与老司机 Time Limit: 20 Sec  Memory Limit: 512 MBSec  Special JudgeSubmit: 106  Solved ...

  4. [UOJ#132][BZOJ4200][luogu_P2304][NOI2015]小园丁与老司机

    [UOJ#132][BZOJ4200][luogu_P2304][NOI2015]小园丁与老司机 试题描述 小园丁 Mr. S 负责看管一片田野,田野可以看作一个二维平面.田野上有 \(n\) 棵许愿 ...

  5. bzoj 4200: [Noi2015]小园丁与老司机【dp+有上下界最小流】

    洛谷上有个点死活卡不过去,不知道是哪里写丑了orz 参考:https://www.cnblogs.com/ditoly/p/BZOJ4200.html 从上往下dp,设f为不向左右走直接上去的值,g为 ...

  6. BZOJ4200 NOI2015小园丁与老司机(动态规划+上下界网络流)

    一看上去就是一个二合一的题.那么先解决第一部分求最优路线(及所有可能在最优路线上的线段). 由于不能往下走,可以以y坐标作为阶段.对于y坐标不同的点,我们将可以直接到达的两点连边,显然这样的边的个数是 ...

  7. [Noi2015]小园丁和老司机

    来自FallDream的博客,未经允许,请勿转载,谢谢. 小园丁 Mr. S 负责看管一片田野,田野可以看作一个二维平面.田野上有n棵许愿树,编号1,2,3,…,n,每棵树可以看作平面上的一个点,其中 ...

  8. UOJ#132&bzoj4200[Noi2015]小园丁与老司机

    看,这是一个传送门 Part A 把坐标离散化,按照纵坐标为第一关键字,横坐标为第二关键字排序 以$f_i$记录来到$i$这个点最多经过点数,那么答案显而易见就是$f_i$加上该层点数 转移的话就是分 ...

  9. 并不对劲的loj2134:uoj132:p2304:[NOI2015]小园丁与老司机

    题目大意 给出平面直角坐标系中\(n\)(\(n\leq5*10^4\))个点,第\(i\)个点的坐标是\(x_i,y_i(|x_i|\leq10^9,1\leq y_i\leq10^9)\),只有朝 ...

随机推荐

  1. 【BZOJ】1697: [Usaco2007 Feb]Cow Sorting牛排序

    [算法]数学置换 [题意]给定n个数,要求通过若干次交换两个数的操作得到排序后的状态,每次交换代价为两数之和,求最小代价. [题解] 考虑置换的定义:置换就是把n个数做一个全排列. 从原数组到排序数组 ...

  2. bzoj 1503 平衡树

    我们可以用一颗平衡树维护每个人的工资,因为工资的变化会影响到后面所有的人,所以我们打一个标签,向平衡树里插入的时候减去这个标签的值,这样代表改变了之后的零点,,这样维护这个标签就好了,输出的时候要加上 ...

  3. Angular2.0 基础: Form

    对于Angular2.0 的Form表单中的隐藏和验证,个人觉得还是挺有意思的. 1.通过ngModel 跟踪修改状态与验证. 在表单中使用 ngModel 可以获得更多的控制权,包括一些常用的验证. ...

  4. ribbon使用eureka的meta进行动态路由

    序 使用eureka的元数据信息,再配上ribbon的路由功能,就可以在api-gateway实现很多功能,比如灰度测试.生产调试等等.下面介绍一下,怎么使用jmnarloch大神提供的ribbon- ...

  5. Lucene7.2.1系列(一)快速入门

    系列文章: Lucene系列(一)快速入门 Lucene系列(二)luke使用及索引文档的基本操作 Lucene系列(三)查询及高亮 Lucene是什么? Lucene在维基百科的定义 Lucene是 ...

  6. Python3安装Celery模块后执行Celery命令报错

    1 Python3安装Celery模块后执行Celery命令报错 pip3 install celery # 安装正常,但是执行celery 命令的时候提示没有_ssl模块什么的 手动在Python解 ...

  7. [How to] 使用Xib来创建view

    1.简介 代码库 正如之前博客介绍的,xib可定义页面的某个部分,特别当此部分区域的view集中并且还有一些相互关联性(如隐藏等)是i特别适合使用xib来进行封装. 本文为[How to]使用自定义c ...

  8. keil问题:Error: failed to execute 'C:\Keil\ARM\BIN\ArmCC'

    1.打开cmd,进入相应的路径下 cd C:\Keil\ARM\BIN\ 输入 armcc  若显示如下界面则说明keil已经注册 2.若注册成功还没解决问题,则说明软件是在网上下载的破解版的,建议卸 ...

  9. HDU 2829 Lawrence(四边形优化DP O(n^2))

    题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=2829 题目大意:有一段铁路有n个站,每个站可以往其他站运送粮草,现在要炸掉m条路使得粮草补给最小,粮草 ...

  10. 使用GitLab进行落地项目的管理,并且自动更新、重启、回滚

      Gitlab 清空项目历史commit,节省空间 http://blog.csdn.net/dounine/article/details/77840416?locationNum=6&f ...