Libre 6005 「网络流 24 题」最长递增子序列 / Luogu 2766 最长递增子序列问题(网络流,最大流)
Libre 6005 「网络流 24 题」最长递增子序列 / Luogu 2766 最长递增子序列问题(网络流,最大流)
Description
问题描述:
给定正整数序列x1,...,xn 。
(1)计算其最长递增子序列的长度s。
(2)计算从给定的序列中最多可取出多少个长度为s的递增子序列。
(3)如果允许在取出的序列中多次使用x1和xn,则从给定序列中最多可取出多少个长度为s的递增子序列。
编程任务:
设计有效算法完成(1)(2)(3)提出的计算任务。
Input
第1 行有1个正整数n,表示给定序列的长度。接下来的1 行有n个正整数n:x1, ..., xn。
Output
第1 行是最长递增子序列的长度s。第2行是可取出的长度为s 的递增子序列个数。第3行是允许在取出的序列中多次使用x1和xn时可取出的长度为s 的递增子序列个数。
Sample Input
4
3 6 2 5
Sample Output
2
2
3
Http
Libre:https://loj.ac/problem/6005
Luogu:https://www.luogu.org/problem/show?pid=2766
Source
网络流,最大流
解决思路
看清题目,不是最长递增子序列是最长不下降子序列。
这道题目首先运用动态规划的方式求出最长不下降子序列,这也是第一问的内容。注意,本题不能使用单调队列的方式,因为要求出到每一个数的最长不下降子序列长度(后面记为F),这在后两问中要用。
那么如何求解第二问呢?
我们把每一个数拆成两个点入点和出点,在每一个数的入点和出点之间连容量为1的边,同时设置一个源点一个汇点。从前往后扫描每一个数,若发现第i个数的F[i]最长不下降子序列长度,则在源点与i的出点之间连一条容量为1的边。若F[i]1,则在其出点与汇点之间连一条容量为1的边。并且,对于任何数i,扫描其前面的每一个数j,若F[i]==F[j]+1且第j位的数<=第i位的数,则在i的出点与j的入点之间连一条容量为1的边。
这样建图,有点类似于分层图的思想,从最高层的F为最长不下降子序列长度,往下每一层长度减一,直到最底下一层长度为1。这样我们跑一边最短路就可以了。
至于第三问,我们只要在重新建图的时候把1到汇点,1的入点到出点,n的入点到出点,源点到n(如果存在的话)这几条边设置为无穷大即可。
但要注意,第三问要特判一下递减的情况,因为这样最长不下降子序列长度为1,跑最大流会出现无穷大的流的情况。
另:这里用Dinic求解最大流,具体请移步我的这篇文章
代码
#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<algorithm>
using namespace std;
const int maxN=2000;
const int maxM=maxN*maxN*4;
const int inf=147483647;
class Edge
{
public:
int u,v,flow;
};
int n;
int cnt=-1;
int F[maxN];
int Arr[maxN];
int Head[maxN];
int Next[maxM];
Edge E[maxM];
int depth[maxN];
int cur[maxN];
int Q[maxM];
void Add_Edge(int u,int v,int flow);
bool bfs();
int dfs(int u,int flow);
int main()
{
scanf("%d",&n);
for (int i=1;i<=n;i++)
scanf("%d",&Arr[i]);
for (int i=1;i<=n;i++)//动态规划求出最长不下降子序列
{
F[i]=1;
for (int j=1;j<i;j++)
if (Arr[j]<=Arr[i])
F[i]=max(F[i],F[j]+1);
}
int maxlength=1;
for (int i=1;i<=n;i++)
maxlength=max(maxlength,F[i]);//得出第一问答案
printf("%d\n",maxlength);
memset(Head,-1,sizeof(Head));
for (int i=1;i<=n;i++)//构造第二问的图
{
Add_Edge(i,i+n,1)//连接入点和出点
if (F[i]==maxlength)//若与最长长度相同,则连接源点
Add_Edge(0,i,1);
if (F[i]==1)//若为最小长度1,则连接汇点
Add_Edge(i+n,n*2+1,1);
for (int j=1;j<i;j++)
if ((F[j]==F[i]-1)&&(Arr[j]<=Arr[i]))//向前统计能连的
Add_Edge(i+n,j,1);
}
int Ans=0;//求解最大流
while (bfs())
{
for (int i=0;i<=2*n+1;i++)
cur[i]=Head[i];
while (int di=dfs(0,inf))
Ans+=di;
}
cout<<Ans<<endl;
memset(Head,-1,sizeof(Head));
cnt=-1;
for (int i=1;i<=n;i++)//构造第三问的图
{
int nowflow=1;
if ((i==1)||(i==n))//1和n的流量为无穷大
nowflow=inf;
if (maxlength==1)//注意这里特判递减序列的情况
Add_Edge(i,i+n,1);
else
Add_Edge(i,i+n,inf);
if (F[i]==maxlength)
Add_Edge(0,i,nowflow);
if (F[i]==1)
Add_Edge(i+n,n*2+1,nowflow);
for (int j=1;j<i;j++)
if ((F[j]==F[i]-1)&&(Arr[j]<=Arr[i]))
Add_Edge(i+n,j,1);
}
Ans=0;//求解最大流
while (bfs())
{
for (int i=0;i<=2*n+1;i++)
cur[i]=Head[i];
while (int di=dfs(0,inf))
Ans+=di;
}
cout<<Ans<<endl;
return 0;
}
void Add_Edge(int u,int v,int flow)
{
cnt++;
Next[cnt]=Head[u];
Head[u]=cnt;
E[cnt].u=u;
E[cnt].v=v;
E[cnt].flow=flow;
cnt++;
Next[cnt]=Head[v];
Head[v]=cnt;
E[cnt].v=u;
E[cnt].u=v;
E[cnt].flow=0;
}
bool bfs()
{
memset(depth,-1,sizeof(depth));
int h=1,t=0;
Q[1]=0;
depth[0]=1;
do
{
t++;
int u=Q[t];
//cout<<u<<endl;
for (int i=Head[u];i!=-1;i=Next[i])
{
int v=E[i].v;
if ((depth[v]==-1)&&(E[i].flow>0))
{
depth[v]=depth[u]+1;
h++;
Q[h]=v;
}
}
}
while (t!=h);
if (depth[n*2+1]==-1)
return 0;
return 1;
}
int dfs(int u,int flow)
{
if (u==n*2+1)
return flow;
for (int &i=cur[u];i!=-1;i=Next[i])
{
int v=E[i].v;
if ((depth[v]==depth[u]+1)&&(E[i].flow>0))
{
int di=dfs(v,min(flow,E[i].flow));
if (di>0)
{
E[i].flow-=di;
E[i^1].flow+=di;
return di;
}
}
}
return 0;
}
Libre 6005 「网络流 24 题」最长递增子序列 / Luogu 2766 最长递增子序列问题(网络流,最大流)的更多相关文章
- Libre 6013 「网络流 24 题」负载平衡 (网络流,最小费用最大流)
Libre 6013 「网络流 24 题」负载平衡 (网络流,最小费用最大流) Description G 公司有n 个沿铁路运输线环形排列的仓库,每个仓库存储的货物数量不等.如何用最少搬运量可以使n ...
- Libre 6012 「网络流 24 题」分配问题 (网络流,费用流)
Libre 6012 「网络流 24 题」分配问题 (网络流,费用流) Description 有n件工作要分配给n个人做.第i个人做第j件工作产生的效益为\(c_{ij}\).试设计一个将n件工作分 ...
- Libre 6011 「网络流 24 题」运输问题 (网络流,最小费用最大流)
Libre 6011 「网络流 24 题」运输问题 (网络流,最小费用最大流) Description W 公司有m个仓库和n个零售商店.第i个仓库有\(a_i\)个单位的货物:第j个零售商店需要\( ...
- Libre 6010「网络流 24 题」数字梯形 (网络流,最大费用最大流)
Libre 6010「网络流 24 题」数字梯形 (网络流,最大费用最大流) Description 给定一个由n 行数字组成的数字梯形如下图所示.梯形的第一行有m 个数字.从梯形的顶部的m 个数字开 ...
- Libre 6009 「网络流 24 题」软件补丁 / Luogu 2761 软件安装问题 (最短路径,位运算)
Libre 6009 「网络流 24 题」软件补丁 / Luogu 2761 软件安装问题 (最短路径,位运算) Description T 公司发现其研制的一个软件中有 n 个错误,随即为该软件发放 ...
- Libre 6008 「网络流 24 题」餐巾计划 (网络流,最小费用最大流)
Libre 6008 「网络流 24 题」餐巾计划 (网络流,最小费用最大流) Description 一个餐厅在相继的N天里,第i天需要Ri块餐巾(i=l,2,-,N).餐厅可以从三种途径获得餐巾. ...
- Libre 6007 「网络流 24 题」方格取数 / Luogu 2774 方格取数问题 (网络流,最大流)
Libre 6007 「网络流 24 题」方格取数 / Luogu 2774 方格取数问题 (网络流,最大流) Description 在一个有 m*n 个方格的棋盘中,每个方格中有一个正整数.现要从 ...
- Libre 6006 「网络流 24 题」试题库 / Luogu 2763 试题库问题 (网络流,最大流)
Libre 6006 「网络流 24 题」试题库 / Luogu 2763 试题库问题 (网络流,最大流) Description 问题描述: 假设一个试题库中有n道试题.每道试题都标明了所属类别.同 ...
- Libre 6004 「网络流 24 题」圆桌聚餐(网络流,最大流)
Libre 6004 「网络流 24 题」圆桌聚餐(网络流,最大流) Description 假设有来自n个不同单位的代表参加一次国际会议.每个单位的代表数分别为 ri.会议餐厅共有m张餐桌,每张餐桌 ...
随机推荐
- 20155229《网络对抗技术》Exp4:恶意代码分析
实验内容 使用schtasks指令监控系统运行 schtasks指令:允许管理员在本地或远程系统上创建计划任务. SCHTASKS /Create [/S system [/U username [/ ...
- 20155229《网络对抗技术》Exp3:免杀原理与实践
实验预习 免杀: 看为一种能使病毒木马避免被杀毒软件查杀的技术. 免杀的分类: 开源免杀:指在有病毒.木马源代码的前提下,通过修改源代码进行免杀.. 手工免杀:指在仅有病毒.木马的可执行文件(.exe ...
- QTableWidget 详细使用
QTableWidget 详细使用
- [CF1019D]Large Triangle[极角排序+二分]
题意 给出平面上 \(n\) 个点 \((x_i, y_i)\),问是否存在三个点构成的三角形的面积恰好为 \(S\) ,有的话,输出任意一组解即可. \(n\leq 2000\) 分析 BZOJ37 ...
- 蓝牙学习笔记三(Android Debug)
android 端可以通过两种方式去Debug: 一.在手机的设置功能里,开发者模式 Enable,如下图: http://blog.bluetooth.com/debugging-bluetoo ...
- 《Effective Java》学习笔记 —— 枚举、注解与方法
Java的枚举.注解与方法... 第30条 用枚举代替int常量 第31条 用实例域代替序数 可以考虑定义一个final int 代替枚举中的 ordinal() 方法. 第32条 用EnumSet代 ...
- App推荐 | Google Tasks
前不久,Google推出了一款移动任务管理应用Google Task,在使用2天后,写一下使用感受,并与Google同类产品Keep进行一个对比. 首先欣赏几张官方的App截图 然后来看一下官方的介绍 ...
- java 中的内部类总结
内部类不是很好理解,但说白了其实也就是一个类中还包含着另外一个类. 如同一个人是由大脑.肢体.器官等身体结果组成,而内部类相当于其中的某个器官之一,例如心脏:它也有自己的属性和行为(血液.跳动). 显 ...
- [沈航软工教学] 学生项目Coding地址汇总
同学们把自己的coding主页链接贴在评论里,要求格式"班号+学号+coding主页链接",如: "1301+13061193 + https://coding.net/ ...
- thinkphp5报错
thinkPHP5配置nginx环境无法打开(require():open_basedirrestrictionineffect.File(/mnt/hgfs/ro (2018-07-19 22:05 ...